Matematisk grundlag for galaktisk manglende masse: Skive, kugle, tæthed, potentiale og radial skalering
TL;DR: Problemet med den manglende masse opstår, når den masse, der udledes af galaktiske rotationskurver, overstiger den masse, der observeres direkte i stjerner, gas og støv. Matematisk set kræver det, at man forbinder overfladetætheden på en skive, rumtætheden i tre dimensioner, gravitationspotentialet, den radiale acceleration og den indesluttede masse.
1. Radiale koordinater og geometri
Vi skelner mellem to geometrier:
- Skivegeometri: Synligt galaktisk stof er hovedsageligt fordelt i en tynd roterende skive.
- Sfærisk geometri: Mørk eller manglende masse modelleres ofte som en nogenlunde sfærisk halo.
Diskområde-element:
\[ dA = R\,dR\,d\phi \]
Sfærisk volumenelement:
\[ dV = 4\pi r^2\,dr \]
Det samme symbol \(r\) bruges ofte for galaktocentrisk radius, men betydningen afhænger af geometrien. I en skive er \(R\) en cylindrisk radius. I en halo er \(r\) normalt en sfærisk radius.
2. Synlig masse på en galaktisk skive
Den synlige skive er ofte tilnærmet med en eksponentiel overfladetæthed:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Massen af en ring mellem \(R\) og \(R+dR\) er:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Den samlede synlige diskmasse er derfor:
\[ 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] \]
Ved stor radius:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Den synlige skives masse nærmer sig en endelig værdi.
3. Sfærisk masse og volumendensitet
For en sfærisk massefordeling bestemmer volumendensiteten \(\rho(r)\) den indesluttede masse:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Den omvendte relation er:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dette forhold er centralt for problemet med den manglende masse. Hvis den udledte masse vokser lineært med radius, så falder den tilsvarende sfæriske tæthed som \(1/r^2\).
4. Dynamisk masse fra cirkulær bevægelse
For cirkulær bevægelse opfylder tyngdeaccelerationen:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Derfor:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
For en flad rotationskurve:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Dette giver standardskalering:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definition af manglende masse
Den manglende masse er forskellen mellem den dynamiske masse og den synlige masse:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
For en eksponentiel synlig disk:
\[ 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] \]
For \(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] \]
Ved stor radius mættes skiveudtrykket, mens det dynamiske udtryk fortsætter med at vokse omtrent som \(r\).
6. Tyngdekraftpotentiale i 3D
Det newtonske gravitationspotentiale, der genereres af en punktmasse, er:
\[ \Phi(r)=-\frac{GM}{r} \]
Det tilsvarende tyngdefelt er den radiale afledning af potentialet:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Det forklarer forholdet mellem \(1/r\) og \(1/r^2\): Potentialet i en lokaliseret masse falder som \(1/r\), mens kraften eller accelerationen falder som \(1/r^2\).
7. Poisson-ligningen
Massetæthed og tyngdekraftspotentiale er forbundet gennem Poissons ligning:
\[ \nabla^2\Phi=4\pi G\rho \]
I sfærisk symmetri bliver dette til:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \højre) = 4\pi G\rho(r) \]
Denne ligning forbinder tre størrelser:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Potentialet i en udvidet 3D-tæthed
For en generel 3D-tæthedsfordeling er potentialet:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Kernen \(1/|\mathbf{x}-\mathbf{x}’|\) er den matematiske oprindelse af \(1/r\)-potentialet i tre dimensioner.
9. Potentialet i en tynd disk
For en tynd skive med overfladetæthed \(\Sigma(R’)\) kan tyngdepotentialet i skiveplanet skrives 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}} \]
Afstanden mellem et feltpunkt ved radius \(R\) og et kildepunkt ved radius \(R’\) er:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Den radiale acceleration i skiven fås ved at differentiere potentialet:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Rotationshastigheden følger af:
\[ v^2(R)=R\,|g_R(R)|. \]
10. Projektion af en 3D-interaktion på disken
Hvis en interaktion udbreder sig i tre dimensioner, men evalueres i diskens plan, indfører den radiale projektion en geometrisk faktor. For to punkter i skiven, der er adskilt af afstanden \(d\), er den radiale projektionsfaktor:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
En generisk 3D-radial kerne \(K(d)\), projiceret på disken, ser således ud:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
For eksempel har en Newtonsk kraftlignende kerne:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
En eksponentiel 3D-kerne kan skrives som:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Eksponentielle kerner i tre dimensioner
En rent eksponentiel radial faktor har formen:
\[ e^{-r/\lambda} \]
I tredimensionel feltteori optræder et eksponentielt skærmet potentiale ofte i Yukawa-lignende form:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Den tilhørende radiale kraft indeholder både \(1/r^2\) og eksponentielle udtryk:
\[ 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) \]
Dette viser, hvorfor eksponentiel radial opførsel i 3D ikke er uafhængig af \(1/r\)-geometrien. Eksponentialet styrer dæmpningen, mens \(1/r\) og \(1/r^2\) stammer fra tredimensionel spredning.
12. Radiale skaleringslove
Problemet med den manglende masse er stærkt knyttet til radial skalering. Flere vigtige radiale love optræder gentagne gange:
| Mængde | Typisk skalering | Betydning |
|---|---|---|
| Potentiale for punktmasse | \(\Phi(r)\sim 1/r\) | 3D grøn funktion af tyngdekraften |
| Kraft af punktmasse | \(g(r)\sim 1/r^2\) | Afledt af \(1/r\) |
| Flad rotationshastighed | \(v(r)\sim konstant\) | Observeret i ydre galaktiske skiver |
| Dynamisk masse | \(M(r)\sim r\) | Påkrævet af flad rotation |
| Halo-tæthed | \(\rho(r)\sim 1/r^2\) | Giver \(M(r)\sim r\) |
| Eksponentiel disk | \(\Sigma(R)\sim e^{-R/R_d}\) | Synlig disk forsvinder hurtigt |
| Skærmet 3D-potentiale | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Eksponentiel dæmpning plus 3D-spredning |
13. Fra tæthed til rotationskurve
For en sfærisk halo med tæthed:
\[ \rho(r)\propto \frac{1}{r^2} \]
den indesluttede masse er:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Og så:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approx konstant \]
\[ v(r)\approx konstant \]
Dette er den matematiske bro mellem en \(1/r^2\) halotæthed og en flad galaktisk rotationskurve.
14. Fra diskmasse til manglende masse
Den synlige skivemasse vokser hurtigt i starten og mættes derefter:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Den dynamiske masse, der udledes af en flad rotationskurve, bliver ved med at vokse:
\[ M_{\rm dyn}(r)\propto r \]
Derfor opfører den manglende masse sig omtrent 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 \]
Ved tilstrækkelig stor radius:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Matematisk advarsel: Skive og kugle kan ikke udskiftes
Ligningen
\[ M(r)=\frac{v(r)^2r}{G} \]
er eksakt for sfærisk symmetri. For en fladtrykt skive skal man beregne potentialet ved at integrere over skiven og derefter udlede den radiale acceleration. Det sfæriske udtryk bruges ofte som en effektiv tilnærmelse, især når man diskuterer den masse, der er nødvendig for at understøtte en given rotationskurve.
16. Oversigt over vigtige ligninger
Diskens overfladetæthed:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Diskens masseelement:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Synlig diskmasse:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Sfærisk volumenmasse:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Densitet fra indesluttet masse:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dynamisk masse:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Manglende masse:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Newtonsk potentiale:
\[ \Phi(r)=-\frac{GM}{r} \]
Poisson-ligningen:
\[ \nabla^2\Phi=4\pi G\rho \]
3D-potentialeintegral:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potentiale for tynde skiver:
\[ \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}} \]
Projiceret radial kerne:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Eksponentielt 3D-skærmpotentiale:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Konklusion
Problemet med den manglende masse er et matematisk misforhold mellem to radiale adfærdsmønstre. Den synlige skive følger en eksponentiel overfladetæthed og når en endelig kumulativ masse. Den dynamiske masse, der udledes af omtrent flade rotationskurver, vokser nogenlunde lineært med radius. Hvis den fortolkes som en sfærisk halo, svarer det til en tæthed, der falder omtrent som \(1/r^2\). Ligningerne for diskintegration, sfæriske skaller, gravitationspotentiale og radial projektion giver det matematiske sprog, der er nødvendigt for at analysere denne uoverensstemmelse.