Fondamenti matematici della massa mancante galattica: Disco, sfera, densità, potenziale e scala radiale
TL;DR: Il problema della massa mancante si presenta quando la massa dedotta dalle curve di rotazione galattica supera la massa osservata direttamente nelle stelle, nel gas e nella polvere. Matematicamente, questo richiede di collegare la densità di superficie su un disco, la densità di volume in tre dimensioni, il potenziale gravitazionale, l’accelerazione radiale e la massa racchiusa.
1. Coordinate radiali e geometria
Distinguiamo due geometrie:
- Geometria del disco: la materia galattica visibile è distribuita principalmente in un sottile disco rotante.
- Geometria sferica: la massa oscura o mancante è spesso modellata come un alone approssimativamente sferico.
Elemento dell’area del disco:
\[ dA = R\,dR\,d\phi \]
Elemento di volume sferico:
\[ dV = 4\pi r^2\,dr \]
Lo stesso simbolo \(r\) è spesso utilizzato per il raggio galattocentrico, ma il significato dipende dalla geometria. In un disco, \(R\) è un raggio cilindrico. In un alone, \(r\) è solitamente un raggio sferico.
2. Massa visibile su un disco galattico
Il disco visibile è spesso approssimato da una densità superficiale esponenziale:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
La massa di un anello compreso tra \(R\) e \(R+dR\) è:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
La massa cumulativa del disco visibile è quindi:
\[ 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 \sinistra[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \destra) \right] \]
A grande raggio:
\[ M_{\\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
La massa del disco visibile si avvicina a un valore finito.
3. Massa sferica e densità di volume
Per una distribuzione di massa sferica, la densità di volume \(\rho(r)\) determina la massa racchiusa:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
La relazione inversa è:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Questa relazione è centrale per il problema della massa mancante. Se la massa dedotta cresce linearmente con il raggio, allora la densità sferica corrispondente diminuisce come \(1/r^2\).
4. Massa dinamica dal moto circolare
Per il moto circolare, l’accelerazione gravitazionale soddisfa:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Pertanto:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Per una curva di rotazione piatta:
\[ v(r)\i}, approssimazione di v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Questo fornisce la scala standard:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definizione di massa mancante
La massa mancante è la differenza tra la massa dinamica e la massa visibile:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Per un disco visibile esponenziale:
\[ M_{{{rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \sinistra[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \destra) \right] \]
Per \(v(r)\ approssimativamente v_0\):
\[ M_{{{rm miss}(r)\approssimativamente \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \sinistra[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \destra) \right] \]
A grandi raggi, il termine del disco si satura, mentre il termine dinamico continua a crescere approssimativamente come \(r\).
6. Potenziale gravitazionale in 3D
Il potenziale gravitazionale newtoniano generato da una massa puntiforme è:
\[ \Phi(r)=-\frac{GM}{r} \]
Il campo gravitazionale corrispondente è la derivata radiale del potenziale:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Questo spiega la relazione tra \(1/r\) e \(1/r^2\): il potenziale di una massa localizzata cade come \(1/r\), mentre la forza o l’accelerazione cade come \(1/r^2\).
7. Equazione di Poisson
La densità di massa e il potenziale gravitazionale sono collegati attraverso l’equazione di Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Nella simmetria sferica, questo diventa:
\[ \frac{1}{r^2}\frac{d}{dr} \sinistra( r^2\frac{d\Phi}{dr} \destra) = 4\pi G\rho(r) \]
Questa equazione collega tre quantità:
\[ \rho(r) \freccia dritta M(r) \freccia retta \Phi(r) \longrightarrow v(r) \]
8. Potenziale di una densità 3D estesa
Per una distribuzione generale della densità in 3D, il potenziale è:
\[ \Phi(\mathbf{x}) = -Gint \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Il kernel \(1/|mathbf{x}-\mathbf{x}’|\) è l’origine matematica del potenziale \(1/r\) in tre dimensioni.
9. Potenziale di un disco sottile
Per un disco sottile con densità superficiale \(\Sigma(R’)\), il potenziale gravitazionale nel piano del disco può essere scritto come:
\[ \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}} \]
La distanza tra un punto di campo di raggio \(R\) e un punto sorgente di raggio \(R’\) è:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
L’accelerazione radiale nel disco si ottiene differenziando il potenziale:
\[ g_R(R)=-\frac{\parziale \Phi}{\parziale R} \]
La velocità di rotazione deriva da:
\[ v^2(R)=R\,|g_R(R)| \]
10. Proiezione di un’interazione 3D sul disco
Se un’interazione si propaga in tre dimensioni ma viene valutata nel piano del disco, la proiezione radiale introduce un fattore geometrico. Per due punti nel disco separati da una distanza \(d\), il fattore di proiezione radiale è:
\[ \coseta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Quindi un generico kernel radiale 3D \(K(d)\), proiettato sul disco, appare come:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Per esempio, un kernel simile a una forza newtoniana ha:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{{rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Un kernel 3D esponenziale può essere scritto come:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{{rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Kernel esponenziali in tre dimensioni
Un fattore radiale puramente esponenziale ha la forma:
\[ e^{-r/\lambda} \]
Nella teoria di campo tridimensionale, un potenziale schermato esponenzialmente appare spesso nella forma di Yukawa:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
La forza radiale associata contiene sia termini \(1/r^2\) che esponenziali:
\[ g_Y(r) = -frac{d}{dr} \sinistra( \frac{e^{-r/\lambda}}{r} \destra) \]
\[ g_Y(r)\propto e^{-r/\lambda} \a sinistra( \frac{1}{r^2} + \frac{1}{\lambda r} \destra) \]
Questo dimostra perché il comportamento radiale esponenziale in 3D non è indipendente dalla geometria di \(1/r\). L’esponenziale controlla l’attenuazione, mentre \(1/r\) e \(1/r^2\) derivano dalla diffusione tridimensionale.
12. Leggi di scala radiale
Il problema della massa mancante è fortemente legato alla scalatura radiale. Diverse leggi radiali importanti appaiono ripetutamente:
| Quantità | Scala tipica | Significato |
|---|---|---|
| Potenziale della massa puntiforme | \(\Phi(r)\sim 1/r\) | Funzione verde 3D della gravità |
| Forza della massa puntiforme | \(g(r)\sim 1/r^2\) | Derivata di \(1/r\) |
| Velocità di rotazione piana | \(v(r)\sim costante\) | Osservata nei dischi galattici esterni |
| Massa dinamica | \(M(r)\sim r\) | Richiesto dalla rotazione piatta |
| Densità dell’alone | \(\rho(r)\sim 1/r^2\) | Dà \(M(r)\sim r\) |
| Disco esponenziale | \(\Sigma(R)\sim e^{-R/R_d}\) | Il disco visibile svanisce rapidamente |
| Potenziale 3D schermato | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Attenuazione esponenziale più diffusione 3D |
13. Dalla densità alla curva di rotazione
Per un alone sferico con densità:
\[ \rho(r)\propto \frac{1}{r^2} \]
la massa racchiusa è:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Poi:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approssimativamente costante \]
\[ v(r)\approssimativamente costante \]
Questo è il ponte matematico tra una densità dell’alone \(1/r^2\) e una curva di rotazione galattica piatta.
14. Dalla massa del disco alla massa mancante
La massa del disco visibile cresce rapidamente all’inizio e poi si satura:
\[ M_{\rm disk}(R)\rightarrow M_d \]
La massa dinamica desunta da una curva di rotazione piatta continua a crescere:
\[ M_{\rm dyn}(r)\propto r \]
Pertanto, la massa mancante si comporta approssimativamente come:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\approssimativamente \frac{v_0^2}{G}r-M_d \]
Con un raggio sufficientemente grande:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Avviso matematico: disco e sfera non sono intercambiabili.
L’equazione
\[ M(r)=\frac{v(r)^2r}{G} \]
è esatto per la simmetria sferica. Per un disco appiattito, si dovrebbe calcolare il potenziale integrando sul disco e poi ricavare l’accelerazione radiale. L’espressione sferica è spesso utilizzata come approssimazione efficace, soprattutto quando si parla della massa necessaria per sostenere una determinata curva di rotazione.
16. Riassunto delle equazioni chiave
Densità della superficie del disco:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Elemento di massa del disco:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Massa del disco visibile:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \Sinistra[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \destra) \right] \]
Massa volumetrica sferica:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Densità dalla massa racchiusa:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Massa dinamica:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Massa mancante:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Potenziale newtoniano:
\[ \Phi(r)=-\frac{GM}{r} \]
Equazione di Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Integrale di potenziale 3D:
\[ \Phi(\mathbf{x}) = -G´int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potenziale del disco sottile:
\[ \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}} \]
Kernel radiale proiettato:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Potenziale esponenziale 3D schermato:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Conclusione
Il problema della massa mancante è una mancata corrispondenza matematica tra due comportamenti radiali. Il disco visibile segue una densità superficiale esponenziale e raggiunge una massa cumulativa finita. La massa dinamica dedotta dalle curve di rotazione approssimativamente piatte cresce in modo approssimativamente lineare con il raggio. Se interpretato come un alone sferico, questo corrisponde a una densità che diminuisce approssimativamente come \(1/r^2\). Le equazioni di integrazione del disco, i gusci sferici, il potenziale gravitazionale e la proiezione radiale forniscono il linguaggio matematico necessario per analizzare questo disallineamento.