Mathematische Grundlagen der fehlenden galaktischen Masse: Scheibe, Sphäre, Dichte, Potential und radiale Skalierung
TL;DR: Das Problem der fehlenden Masse tritt auf, wenn die aus galaktischen Rotationskurven abgeleitete Masse die in Sternen, Gas und Staub direkt beobachtete Masse übersteigt. Mathematisch gesehen erfordert dies eine Verbindung zwischen der Oberflächendichte einer Scheibe, der Volumendichte in drei Dimensionen, dem Gravitationspotential, der Radialbeschleunigung und der eingeschlossenen Masse.
1. Radiale Koordinaten und Geometrie
Wir unterscheiden zwei Geometrien:
- Scheibengeometrie: Die sichtbare galaktische Materie ist hauptsächlich in einer dünnen rotierenden Scheibe verteilt.
- Sphärische Geometrie: Dunkle oder fehlende Masse wird oft als ein grob kugelförmiger Halo modelliert.
Plattenbereich-Element:
\[ dA = R\,dR\,d\phi \]
Kugelförmiges Volumenelement:
\[ dV = 4\pi r^2\,dr \]
Das gleiche Symbol \(r\) wird oft für den galaktozentrischen Radius verwendet, aber die Bedeutung hängt von der Geometrie ab. In einer Scheibe ist \(R\) ein zylindrischer Radius. In einem Halo ist \(r\) normalerweise ein kugelförmiger Radius.
2. Sichtbare Masse auf einer galaktischen Scheibe
Die sichtbare Scheibe wird oft durch eine exponentielle Oberflächendichte approximiert:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Die Masse eines Ringes zwischen \(R\) und \(R+dR\) ist:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Die kumulative sichtbare Scheibenmasse ist also:
\[ 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] \]
Bei großem Radius:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Die sichtbare Scheibenmasse nähert sich einem endlichen Wert.
3. Sphärische Masse und Volumendichte
Bei einer kugelförmigen Massenverteilung bestimmt die Volumendichte \(\rho(r)\) die eingeschlossene Masse:
\[ M(r)=4\pi\int_0^r \rho(r‘)r’^2\,dr‘ \]
Die umgekehrte Beziehung lautet:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Diese Beziehung ist von zentraler Bedeutung für das Problem der fehlenden Masse. Wenn die abgeleitete Masse linear mit dem Radius wächst, dann nimmt die entsprechende kugelförmige Dichte mit \(1/r^2\) ab.
4. Dynamische Masse aus Kreisbewegung
Für eine Kreisbewegung erfüllt die Gravitationsbeschleunigung folgende Bedingungen:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Deshalb:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Für eine flache Rotationskurve:
\[ v(r)\ca. v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Dies ergibt die Standard-Skalierung:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definition der fehlenden Masse
Die fehlende Masse ist die Differenz zwischen der dynamischen Masse und der sichtbaren Masse:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Für eine exponentielle sichtbare Scheibe:
\[ 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] \]
Für \(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] \]
Bei großen Radien sättigt der Scheibenterm, während der dynamische Term weiterhin ungefähr mit \(r\) wächst.
6. Gravitationspotential in 3D
Das Newtonsche Gravitationspotenzial, das von einer Punktmasse erzeugt wird, ist:
\[ \Phi(r)=-\frac{GM}{r} \]
Das entsprechende Gravitationsfeld ist die radiale Ableitung des Potentials:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Dies erklärt die Beziehung zwischen \(1/r\) und \(1/r^2\): Das Potenzial einer lokalisierten Masse fällt mit \(1/r\), während die Kraft oder Beschleunigung mit \(1/r^2\) fällt.
7. Poisson-Gleichung
Massendichte und Gravitationspotential sind durch die Poisson-Gleichung miteinander verbunden:
\[ \nabla^2\Phi=4\pi G\rho \]
In sphärischer Symmetrie wird dies zu:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Diese Gleichung verbindet drei Größen:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Das Potenzial einer erweiterten 3D-Dichte
Für eine allgemeine 3D-Dichteverteilung ist das Potenzial:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}‘)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x‘ \]
Der Kernel \(1/|\mathbf{x}-\mathbf{x}’|\) ist der mathematische Ursprung des \(1/r\) Potentials in drei Dimensionen.
9. Potential einer dünnen Scheibe
Für eine dünne Scheibe mit einer Oberflächendichte \(\Sigma(R‘)\) kann das Gravitationspotential in der Scheibenebene wie folgt geschrieben werden:
\[ \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}} \]
Der Abstand zwischen einem Feldpunkt am Radius \(R\) und einem Quellpunkt am Radius \(R’\) ist:
\[ d(R,R‘,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Die radiale Beschleunigung in der Scheibe erhält man durch Differenzierung des Potentials:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Die Rotationsgeschwindigkeit ergibt sich aus:
\[ v^2(R)=R\,|g_R(R)| \]
10. Projektion einer 3D-Interaktion auf die Scheibe
Wenn sich eine Wechselwirkung in drei Dimensionen ausbreitet, aber in der Scheibenebene ausgewertet wird, führt die radiale Projektion einen geometrischen Faktor ein. Für zwei Punkte in der Scheibe, die durch den Abstand \(d\) getrennt sind, ist der radiale Projektionsfaktor:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R‘,\phi)} \]
Ein allgemeiner radialer 3D-Kernel \(K(d)\), der auf die Scheibe projiziert wird, sieht also so aus:
\[ K_{\rm disk}(R,R‘,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Ein Newtonscher kraftähnlicher Kern hat zum Beispiel:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Ein exponentieller 3D-Kernel kann wie folgt geschrieben werden:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Exponentiale Kerne in drei Dimensionen
Ein rein exponentieller Radialfaktor hat die Form:
\[ e^{-r/\lambda} \]
In der dreidimensionalen Feldtheorie erscheint ein exponentiell abgeschirmtes Potential oft in Yukawa-ähnlicher Form:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Die zugehörige Radialkraft enthält sowohl \(1/r^2\) als auch exponentielle Terme:
\[ g_Y(r) = -\frac{d}{dr} \left( \frac{e^{-r/\lambda}}{r} \rechts) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \rechts) \]
Dies zeigt, warum das exponentielle radiale Verhalten in 3D nicht unabhängig von der Geometrie von \(1/r\) ist. Der Exponentialwert steuert die Dämpfung, während \(1/r\) und \(1/r^2\) durch die dreidimensionale Ausbreitung entstehen.
12. Radiale Skalierungsgesetze
Das Problem der fehlenden Masse ist eng mit der radialen Skalierung verbunden. Mehrere wichtige radiale Gesetze tauchen immer wieder auf:
| Menge | Typische Skalierung | Bedeutung |
|---|---|---|
| Potential einer Punktmasse | \(\Phi(r)\sim 1/r\) | 3D Green-Funktion der Schwerkraft |
| Kraft der Punktmasse | \(g(r)\sim 1/r^2\) | Ableitung von \(1/r\) |
| Flache Rotationsgeschwindigkeit | \(v(r)\sim konstant\) | Beobachtet in äußeren galaktischen Scheiben |
| Dynamische Masse | \(M(r)\sim r\) | Erforderlich für die flache Rotation |
| Halo-Dichte | \(\rho(r)\sim 1/r^2\) | Ergibt \(M(r)\sim r\) |
| Exponentiale Scheibe | \(\Sigma(R)\sim e^{-R/R_d}\) | Sichtbare Scheibe verblasst schnell |
| Abgeschirmtes 3D-Potenzial | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Exponentielle Abschwächung plus 3D-Ausbreitung |
13. Von der Dichte zur Rotationskurve
Für einen kugelförmigen Halo mit Dichte:
\[ \rho(r)\propto \frac{1}{r^2} \]
die eingeschlossene Masse ist:
\[ M(r)=4\pi\int_0^r \rho(r‘)r’^2dr‘ \]
\[ M(r)\propto r \]
Dann:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\ca. konstant \]
\[ v(r)\ annähernd konstant \]
Dies ist die mathematische Brücke zwischen einer \(1/r^2\) Halo-Dichte und einer flachen galaktischen Rotationskurve.
14. Von der Scheibenmasse zur fehlenden Masse
Die sichtbare Scheibenmasse wächst zunächst schnell und sättigt dann:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Die dynamische Masse, die aus einer flachen Rotationskurve abgeleitet wird, wächst weiter:
\[ M_{\rm dyn}(r)\propto r \]
Daher verhält sich die fehlende Masse ungefähr so:
\[ 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 \]
Bei ausreichend großem Radius:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Mathematische Warnung: Scheibe und Kugel sind nicht austauschbar
Die Gleichung
\[ M(r)=\frac{v(r)^2r}{G} \]
ist für sphärische Symmetrie exakt. Für eine abgeflachte Scheibe sollte man das Potenzial durch Integration über die Scheibe berechnen und dann die radiale Beschleunigung ableiten. Der kugelförmige Ausdruck wird oft als effektive Annäherung verwendet, insbesondere wenn es um die Masse geht, die zur Unterstützung einer bestimmten Rotationskurve benötigt wird.
16. Zusammenfassung der wichtigsten Gleichungen
Dichte der Scheibenoberfläche:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Massenelement der Scheibe:
\[ dM_{\rm Scheibe}=2\pi R\Sigma(R)dR \]
Sichtbare Scheibenmasse:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Sphärische Volumenmasse:
\[ M(r)=4\pi\int_0^r\rho(r‘)r’^2dr‘ \]
Dichte aus eingeschlossener Masse:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dynamische Masse:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Fehlende Masse:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Newtonsches Potential:
\[ \Phi(r)=-\frac{GM}{r} \]
Poisson-Gleichung:
\[ \nabla^2\Phi=4\pi G\rho \]
3D-Potentialintegral:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}‘)}{|\mathbf{x}-\mathbf{x}’|} d^3x‘ \]
Dünnes Scheibenpotential:
\[ \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}} \]
Projizierter radialer Kernel:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Exponentielles 3D abgeschirmtes Potential:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Fazit
Das Problem der fehlenden Masse ist eine mathematische Unstimmigkeit zwischen zwei radialen Verhaltensweisen. Die sichtbare Scheibe folgt einer exponentiellen Oberflächendichte und erreicht eine endliche kumulative Masse. Die dynamische Masse, die aus annähernd flachen Rotationskurven abgeleitet wird, wächst ungefähr linear mit dem Radius. Wenn man sie als kugelförmigen Halo interpretiert, entspricht dies einer Dichte, die ungefähr mit \(1/r^2\) abnimmt. Die Gleichungen der Scheibenintegration, der Kugelschalen, des Gravitationspotenzials und der radialen Projektion liefern die mathematische Sprache, die für die Analyse dieses Missverhältnisses erforderlich ist.