Fundamentos matemáticos de la masa faltante galáctica: Disco, esfera, densidad, potencial y escala radial
TL;DR: El problema de la masa ausente aparece cuando la masa inferida a partir de las curvas de rotación galáctica supera la masa observada directamente en las estrellas, el gas y el polvo. Matemáticamente, esto requiere conectar la densidad de superficie en un disco, la densidad de volumen en tres dimensiones, el potencial gravitatorio, la aceleración radial y la masa encerrada.
1. Coordenadas radiales y geometría
Distinguimos dos geometrías:
- Geometría del disco: la materia galáctica visible se distribuye principalmente en un fino disco giratorio.
- Geometría esférica: la masa oscura o ausente se modela a menudo como un halo aproximadamente esférico.
Elemento del área del disco:
\[ dA = R\,dR\,d\phi \]
Elemento de volumen esférico:
\[ dV = 4\pi r^2\,dr \]
El mismo símbolo \(r\) se utiliza a menudo para el radio galactocéntrico, pero su significado depende de la geometría. En un disco, \(R\) es un radio cilíndrico. En un halo, \(r\) suele ser un radio esférico.
2. Masa visible en un disco galáctico
El disco visible suele aproximarse mediante una densidad de superficie exponencial:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
La masa de un anillo entre \(R\) y \(R+dR\) es:
\[ dM_{\rm disco}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
La masa acumulada del disco visible es, por tanto:
\[ 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 \izquierda[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \derecha) \derecha] \]
A gran radio:
\[ M_{{rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
La masa del disco visible se aproxima a un valor finito.
3. Masa esférica y densidad volumétrica
Para una distribución de masa esférica, la densidad de volumen \(\rho(r)\) determina la masa encerrada:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
La relación inversa es:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Esta relación es fundamental para el problema de la masa desaparecida. Si la masa inferida crece linealmente con el radio, entonces la densidad esférica correspondiente disminuye como \(1/r^2\).
4. Masa dinámica a partir del movimiento circular
Para un movimiento circular, la aceleración gravitatoria satisface:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Por lo tanto:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Para una curva de rotación plana:
\[ v(r)\aprox v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Esto da la escala estándar:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{rm dyn}(r)\propto \frac{1}{r^2} \]
5. Definición de masa faltante
La masa que falta es la diferencia entre la masa dinámica y la masa visible:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Para un disco visible exponencial:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \izquierda[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \derecha) \derecha] \]
Para \(v(r)\aprox v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \izquierda[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \derecha) \derecha] \]
A radios grandes, el término disco se satura, mientras que el término dinámico sigue creciendo aproximadamente como \(r\).
6. Potencial gravitatorio en 3D
El potencial gravitatorio newtoniano generado por una masa puntual es:
\[ \Phi(r)=-\frac{GM}{r} \]
El campo gravitatorio correspondiente es la derivada radial del potencial:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Esto explica la relación entre \(1/r\) y \(1/r^2\): el potencial de una masa localizada cae como \(1/r\), mientras que la fuerza o aceleración cae como \(1/r^2\).
7. Ecuación de Poisson
La densidad de masa y el potencial gravitatorio están conectados a través de la ecuación de Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
En simetría esférica, esto se convierte en:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \derecha) = 4\pi G\rho(r) \]
Esta ecuación relaciona tres cantidades:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Potencial de una densidad 3D ampliada
Para una distribución de densidad tridimensional general, el potencial es:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
El núcleo \(1/||mathbf{x}-\mathbf{x}’|\) es el origen matemático del potencial \(1/r\) en tres dimensiones.
9. Potencial de un disco delgado
Para un disco delgado con densidad superficial \(\Sigma(R’)\), el potencial gravitatorio en el plano del disco puede escribirse como:
\[ \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 distancia entre un punto campo de radio \(R\) y un punto fuente de radio \(R’\) es:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
La aceleración radial en el disco se obtiene diferenciando el potencial:
\[ g_R(R)=-\frac{\parcial \Phi}{\parcial R} \]
La velocidad de rotación se deduce de:
\[ v^2(R)=R\,|g_R(R)|| \]
10. Proyección de una interacción 3D sobre el disco
Si una interacción se propaga en tres dimensiones pero se evalúa en el plano del disco, la proyección radial introduce un factor geométrico. Para dos puntos del disco separados por una distancia \(d\), el factor de proyección radial es:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Así, un núcleo radial tridimensional genérico \(K(d)\), proyectado sobre el disco, aparece como:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Por ejemplo, un núcleo similar a una fuerza newtoniana tiene:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disco}\propto \frac{R-R’\cos\phi}{d^3} \]
Un núcleo exponencial 3D puede escribirse como:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disco}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Núcleos exponenciales en tres dimensiones
Un factor radial puramente exponencial tiene la forma:
\[ e^{-r/\lambda} \]
En la teoría de campos tridimensionales, un potencial exponencialmente apantallado aparece a menudo en forma yukawa:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
La fuerza radial asociada contiene tanto términos \(1/r^2\) como exponenciales:
\[ g_Y(r) = -\frac{d}{dr} \( \frac{e^{-r/\lambda}}{r} \derecha) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \derecha) \]
Esto demuestra por qué el comportamiento radial exponencial en 3D no es independiente de la geometría \(1/r\). La exponencial controla la atenuación, mientras que \(1/r\) y \(1/r^2\) surgen de la propagación tridimensional.
12. Leyes de escala radial
El problema de la masa perdida está fuertemente ligado a la escala radial. Varias leyes radiales importantes aparecen repetidamente:
| Cantidad | Escala típica | Significado |
|---|---|---|
| Potencial de masa puntual | \(\Phi(r)\sim 1/r\) | Función verde 3D de la gravedad |
| Fuerza de la masa puntual | \(g(r)\sim 1/r^2\) | Derivada de \(1/r\) |
| Velocidad de rotación plana | \(v(r)\sim constante\) | Observada en los discos galácticos exteriores |
| Masa dinámica | \(M(r)\sim r\) | Requerida por la rotación plana |
| Densidad del halo | \(\rho(r)\sim 1/r^2\) | Da \(M(r)\sim r\) |
| Disco exponencial | \(\Sigma(R)\sim e^{-R/R_d}\) | El disco visible se desvanece rápidamente |
| Potencial 3D tamizado | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Atenuación exponencial más dispersión 3D |
13. De la densidad a la curva de rotación
Para un halo esférico con densidad
\[ \rho(r)\propto \frac{1}{r^2} \]
la masa encerrada es:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Entonces:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\aprox constante \]
\[ v(r)\aprox constante \]
Este es el puente matemático entre una densidad de halo \(1/r^2\) y una curva de rotación galáctica plana.
14. De la masa del disco a la masa desaparecida
La masa visible del disco crece rápidamente al principio y luego se satura:
\[ M_{{rm disk}(R)\rightarrow M_d \]
La masa dinámica deducida de una curva de rotación plana sigue creciendo:
\[ M_{\rm dyn}(r)\propto r \]
Por lo tanto la masa faltante se comporta aproximadamente como:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\frac{v_0^2}{G}r-M_d \]
Con un radio suficientemente grande:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{rm miss}(r)\sim \frac{1}{r^2} \]
15. Advertencia matemática: disco y esfera no son intercambiables
La ecuación
\[ M(r)=\frac{v(r)^2r}{G} \]
es exacta para la simetría esférica. Para un disco aplanado, hay que calcular el potencial integrando sobre el disco y luego derivar la aceleración radial. La expresión esférica se utiliza a menudo como una aproximación efectiva, especialmente cuando se habla de la masa necesaria para soportar una curva de rotación dada.
16. Resumen de ecuaciones clave
Densidad de la superficie del disco:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Elemento de masa del disco:
\[ dM_{\rm disco}=2\pi R\Sigma(R)dR \]
Masa visible del disco:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \izquierda[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \derecha) \derecha] \]
Masa volumétrica esférica:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Densidad de la masa encerrada:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Masa dinámica:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Masa perdida:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Potencial newtoniano:
\[ \Phi(r)=-\frac{GM}{r} \]
Ecuación de Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Integral de potencial 3D:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Potencial de disco delgado:
\[ \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}} \]
Núcleo radial proyectado:
\[ \cos\theta= \frac{R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Potencial apantallado exponencial 3D:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Conclusión
El problema de la masa faltante es un desajuste matemático entre dos comportamientos radiales. El disco visible sigue una densidad superficial exponencial y alcanza una masa acumulada finita. La masa dinámica inferida a partir de curvas de rotación aproximadamente planas crece de forma aproximadamente lineal con el radio. Si se interpreta como un halo esférico, esto corresponde a una densidad que disminuye aproximadamente como \(1/r^2\). Las ecuaciones de integración del disco, las envolturas esféricas, el potencial gravitatorio y la proyección radial proporcionan el lenguaje matemático necesario para analizar este desajuste.