Wiskundige grondslagen van ontbrekende massa in de galactica: Schijf, bol, dichtheid, potentieel en radiale schaling
TL;DR: Het probleem van de ontbrekende massa treedt op wanneer de massa die wordt afgeleid uit galactische rotatiecurves groter is dan de massa die direct is waargenomen in sterren, gas en stof. Wiskundig gezien vereist dit het verbinden van oppervlaktedichtheid op een schijf, volumedichtheid in drie dimensies, zwaartekrachtpotentiaal, radiale versnelling en ingesloten massa.
1. Radiale coördinaten en meetkunde
We onderscheiden twee geometrieën:
- Schijfgeometrie: zichtbare galactische materie is voornamelijk verdeeld in een dunne roterende schijf.
- Bolvormige geometrie: donkere of ontbrekende massa wordt vaak gemodelleerd als een ruwweg bolvormige halo.
Schijfruimte-element:
\[ dA = R,dR,dphi \]
Bolvormig volume-element:
\[ dV = 4\pi r^2,dr \]
Hetzelfde symbool wordt vaak gebruikt voor de galactocentrische straal, maar de betekenis hangt af van de geometrie. In een schijf is \(r) een cilindrische straal. In een halo is \(r) meestal een bolstraal.
2. Zichtbare massa op een galactische schijf
De zichtbare schijf wordt vaak benaderd door een exponentiële oppervlaktedichtheid:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
De massa van een annulus tussen \(R) en \(R+dR) is:
\[ dM_{schijf}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm schijf}=2\pi R\Sigma_0e^{-R/R_d}, dR \]
De cumulatieve zichtbare schijfmassa is dus:
\[ M_{schijf}(R)=2\pi\int_0^R R’\Sigma_0e^{-R’/R_d},dR’ \]
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \links[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Bij grote straal:
\[ M_{\rm disk}(R)\rrow 2\piSigma_0R_d^2 \]
De zichtbare schijfmassa nadert een eindige waarde.
3. Bolmassa en volumedichtheid
Voor een bolvormige massaverdeling bepaalt de volumedichtheid \rho(r)\ de ingesloten massa:
\[ M(r)=4\piint_0^r \rho(r’)r’^2,dr’ \]
De omgekeerde relatie is:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Deze relatie staat centraal in het probleem van de ontbrekende massa. Als de afgeleide massa lineair toeneemt met de straal, dan neemt de bijbehorende bolvormige dichtheid af als \(1/r^2).
4. Dynamische massa van cirkelvormige beweging
Voor een cirkelvormige beweging voldoet de gravitatieversnelling:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Daarom:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Voor een vlakke rotatiecurve:
\[ v(r)approx v_0 \]
\[ M_{\rm dyn}(r)approx frac{v_0^2}{G}r \]
Dit geeft de standaard schaalverdeling:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definitie van ontbrekende massa
De ontbrekende massa is het verschil tussen de dynamische massa en de zichtbare massa:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Voor een exponentiële zichtbare schijf:
\[ M_{miss}(r)= \frac{v(r)^2r}{G} – 2piSigma_0R_d^2 \links[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Voor (v(r)eter dan v_0):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2piSigma_0R_d^2 \links[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Bij een grote straal verzadigt de schijfterm, terwijl de dynamische term blijft groeien met ongeveer \(r).
6. Zwaartekrachtpotentiaal in 3D
De Newtoniaanse zwaartekrachtpotentiaal gegenereerd door een puntmassa is:
\[ \Phi(r)=-\frac{GM}{r} \]
Het bijbehorende gravitatieveld is de radiale afgeleide van de potentiaal:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Dit verklaart het verband tussen \(1/r) en \(1/r2): de potentiaal van een gelokaliseerde massa daalt als \(1/r), terwijl de kracht of versnelling daalt als \(1/r2).
7. Poissonvergelijking
Massadichtheid en zwaartekrachtpotentiaal zijn met elkaar verbonden door de vergelijking van Poisson:
\[ \nabla^2\Phi=4pi G\rho \]
In sferische symmetrie wordt dit:
\[ \frac{1}{r^2}\frac{d}{dr} \links( r^2{frac{d{Phi}{dr} \rechts) = 4pi Gho(r) \]
Deze vergelijking koppelt drie grootheden:
\[ \rho(r) \rechtopstaand M(r) \rechtopstaand \Phi(r) \longrightarrow v(r) \]
8. Potentieel van een uitgebreide 3D dichtheid
Voor een algemene 3D dichtheidsverdeling is de potentiaal:
\[ \Phi(\mathbf{x}) = -G \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \^3x’ \]
De kern \(1/|mathbf{x}-\mathbf{x}’|) is de wiskundige oorsprong van de \(1/r) potentiaal in drie dimensies.
9. Potentieel van een dunne schijf
Voor een dunne schijf met oppervlaktedichtheid \(\Sigma(R’)\) kan de zwaartekrachtpotentiaal in het schijfvlak geschreven worden als:
\[ \pi(R) = -G \int_0^^infty \int_0^{2\pi} \frac{Sigma(R’)R’^, dR’^, d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
De afstand tussen een veldpunt met straal \ en een bronpunt met straal \ is:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
De radiale versnelling in de schijf wordt verkregen door de potentiaal te differentiëren:
\[ g_R(R)=-\frac{deel \Phi}{deel R} \]
De rotatiesnelheid volgt uit:
\[ v^2(R)=R, |g_R(R)| \]
10. Projectie van een 3D interactie op de schijf
Als een interactie zich in drie dimensies voortplant, maar geëvalueerd wordt in het schijfvlak, introduceert de radiale projectie een geometrische factor. Voor twee punten in de schijf die gescheiden zijn door een afstand \(d) is de radiale projectiefactor:
\[ \cos. = \frac{R-R’\cosheta}{d(R,R’,\phi)} \]
Een generieke 3D radiale kern (K(d)\), geprojecteerd op de schijf, ziet er dus als volgt uit:
\[ K_{schijf}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Een Newtoniaanse krachtachtige kernel heeft bijvoorbeeld:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{schijf}{propto \frac{R-R’\cos\phi}{d^3} \]
Een exponentiële 3D kernel kan worden geschreven als:
\[ K(d)\propto e^{-d/ambda} \]
\[ K_{schijf}}propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Exponentiële kernels in drie dimensies
Een zuiver exponentiële straalfactor heeft de vorm:
\[ e^{-r/\lambda} \]
In driedimensionale veldtheorie verschijnt een exponentieel afgeschermde potentiaal vaak in een Yukawa-achtige vorm:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
De bijbehorende radiale kracht bevat zowel \(1/r^2) als exponentiële termen:
\[ g_Y(r) = -frac{d}{dr} \links( \frac{e^{-r/\lambda}}{r} \rechts) \]
\[ g_Y(r)\propto e^{-r/\lambda}} \links \frac{1}{r^2} + \frac{1}{lambda r} \rechts) \]
Dit laat zien waarom het exponentiële radiale gedrag in 3D niet onafhankelijk is van de geometrie van \(1/r). De exponentieel controleert de demping, terwijl \(1/r) en \(1/r^2) ontstaan door driedimensionale spreiding.
12. Radiale schalingswetten
Het probleem van de ontbrekende massa is sterk verbonden met radiale schaling. Verschillende belangrijke radiale wetten komen herhaaldelijk voor:
| Hoeveelheid | Typische schaling | Betekenis |
|---|---|---|
| Potentiaal van puntmassa | \Phi(r)\sim 1/r | 3D Groene functie van zwaartekracht |
| Kracht van puntmassa | \g(r)\1/r^2) | Afgeleide van \(1/r) |
| Vlakke rotatiesnelheid | \(v(r)\constante) | Waargenomen in buitenste galactische schijven |
| Dynamische massa | \M(r)ijkconstante) | Vereist door vlakke rotatie |
| Halodichtheid | \(\rho(r)\sim 1/r^2) | Geeft \(M(r)\sim r) |
| Exponentiële schijf | \Sigma(R)\sim e^{-R/R_d}}) | Zichtbare schijf vervaagt snel |
| Afgeschermde 3D potentiaal | \Phi(r)\sim e^{-r/\lambda}/r) | Exponentiële verzwakking plus 3D-spreiding |
13. Van dichtheid naar rotatiecurve
Voor een bolvormige halo met dichtheid:
\[ \rho(r)\frac{1}{r^2} \]
de ingesloten massa is:
\[ M(r)=4\piint_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Dan:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)ijk constant \]
\[ v(r)eter constant \]
Dit is de wiskundige brug tussen een \(1/r^2) halodichtheid en een vlakke galactische rotatiecurve.
14. Van schijfmassa naar ontbrekende massa
De zichtbare schijfmassa groeit eerst snel en verzadigt dan:
\[ M_{rm disk}(R)rightarrow M_d \]
De dynamische massa die afgeleid wordt uit een vlakke rotatiecurve blijft groeien:
\[ M_{\rm dyn}(r)\propto r \]
Daarom gedraagt de ontbrekende massa zich ongeveer als:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)çfrac{v_0^2}{G}r-M_d \]
Bij voldoende grote straal:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Wiskundige waarschuwing: schijf en bol zijn niet verwisselbaar
De vergelijking
\[ M(r)=\frac{v(r)^2r}{G} \]
is exact voor sferische symmetrie. Voor een afgeplatte schijf moet men de potentiaal berekenen door over de schijf te integreren en dan de radiale versnelling af te leiden. De sferische uitdrukking wordt vaak gebruikt als een effectieve benadering, vooral bij het bespreken van de massa die nodig is om een gegeven rotatiecurve te ondersteunen.
16. Samenvatting van sleutelvergelijkingen
Schijf oppervlaktedichtheid:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Schijfmassa-element:
\[ dM_{schijf}=2\pi R\Sigma(R)dR \]
Zichtbare schijfmassa:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \links[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Bolmassa:
\[ M(r)=4\piint_0^r\rho(r’)r’^2dr’ \]
Dichtheid van ingesloten massa:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dynamische massa:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Ontbrekende massa:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Newtoniaanse potentiaal:
\[ \Phi(r)=-\frac{GM}{r} \]
Poisson vergelijking:
\[ \nabla^2\Phi=4pi G\rho \]
3D potentiaalintegraal:
\[ \Phi(\mathbf{x}) = -G \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Dunne schijf potentiaal:
\[ \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}} \]
Geprojecteerde radiale kern:
\[ \cosheta= \frac{R-R’coshephi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Exponentiële 3D afgeschermde potentiaal:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Conclusie
Het probleem van de ontbrekende massa is een wiskundige mismatch tussen twee radiale gedragingen. De zichtbare schijf volgt een exponentiële oppervlaktedichtheid en bereikt een eindige cumulatieve massa. De dynamische massa die wordt afgeleid uit ongeveer vlakke rotatiekrommen groeit ruwweg lineair met de straal. Als dit geïnterpreteerd wordt als een bolvormige halo, komt dit overeen met een dichtheid die ongeveer afneemt als \(1/r^2). De vergelijkingen van schijfintegratie, sferische schillen, zwaartekrachtpotentiaal en radiale projectie bieden de wiskundige taal die nodig is om deze wanverhouding te analyseren.