BeeTheory – Física de la gravedad y de las ondas
Las ecuaciones radiales de la masa oculta de la Vía Láctea
De los perfiles de densidad a las integrales de anillo y las curvas de rotación: un tratamiento matemático de la masa oculta en función del radio galáctico R.
Esta página presenta las ecuaciones radiales utilizadas para describir la masa oculta de la Vía Láctea. Compara los perfiles clásicos de densidad de la materia oscura, las integrales de anillo y cáscara, las ecuaciones de la masa oculta, las curvas de rotación y la interpretación BeeTheory de la masa oculta como un posible efecto de interferencia de ondas.
~10¹² M⊙
Estimación clásica de la masa total del halo de la Vía Láctea.
0,39 GeV/cm³
Densidad local típica de materia oscura cerca del Sol.
R⊙ ≈ 8 kpc
Distancia aproximada del Sol al centro galáctico.
~200 kpc
Escala exterior aproximada utilizada para las estimaciones del halo de la Vía Láctea.
Contenido
- ¿Por qué falta masa?
- Perfiles de densidad ρ(R)
- Masa anular y anillo dM/dR
- Masa de materia oscura encerrada M(<R)
- La curva de rotación V(R)
- Estimaciones observacionales actuales
- Hipótesis contrapuestas
- La perspectiva de BeeTheory
1. ¿Por qué falta masa?
En 1933, Fritz Zwicky observó que las galaxias del cúmulo de Coma se movían demasiado rápido para mantenerse unidas sólo por su masa visible. En la década de 1970, Vera Rubin y Kent Ford midieron las curvas de rotación de las galaxias espirales y descubrieron algo igualmente sorprendente: las estrellas situadas en los grandes radios orbitan casi tan rápido como las situadas cerca del centro, mientras que la gravedad newtoniana de la materia visible predice que deberían ralentizarse.
Para una órbita kepleriana simple alrededor de una masa central, esperamos:
\(V(R)\propto \frac{1}{\sqrt{R}}\)Lo que se observa en su lugar es una curva de rotación aproximadamente plana, o sólo lentamente decreciente:
\(V(R)\capprox V_{{infty}}={mathrm{const}{qquad}{mathrm{for}{ R\gtrsim 5\,\mathrm{kpc}\)Reconciliar estos hechos con la gravedad newtoniana requiere un componente de masa invisible adicional cuya densidad cae aproximadamente como:
\(\rho(r)\propto r^{-2}\)Esto produce una masa total encerrada proporcional al radio:
\(M(<R)\propto R\)y por lo tanto:
\(V\propto \sqrt{\frac{M}{R}}=\mathrm{const}\)Enigma cuantitativo clave
La masa bariónica luminosa de la Vía Láctea es de aproximadamente 5 × 10¹⁰ M⊙. La masa dinámica total inferida a partir de la cinemática hasta aproximadamente 200 kpc es de alrededor de 10¹² M⊙. Esto implica una relación de masa oscura a luminosa de aproximadamente 10 a 1.
2. Perfiles de densidad ρ(R)
Un perfil de densidad es una función matemática que describe cómo varía la densidad ρ de la materia oscura con el radio galactocéntrico r, o con el radio cilíndrico R en el plano galáctico.
2.1 Perfil de la NFW
El perfil NFW, introducido por Navarro, Frenk y White, se deriva de simulaciones cosmológicas de N-cuerpos. Produce una ley de potencia doble característica con una cúspide central.
\(\rho_{\mathrm{NFW}}(r)=\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}\)| Parámetro | Símbolo | Estimación de la Vía Láctea | Papel |
|---|---|---|---|
| Radio de escala | rs | 15-25 kpc | Transición entre las vertientes interior y exterior |
| Densidad característica | ρ0 | Calibrada con la densidad local de materia oscura | Normalización global |
| Pendiente interior | γ | -1 | Comportamiento cáustico |
| Pendiente exterior | – | -3 | Rápido declive en radios grandes |
2.2 Perfil de Einasto
El perfil de Einasto evita una divergencia central estricta y utiliza un parámetro de forma α que permite que la pendiente de la densidad cambie suavemente con el radio.
\(\rho_{\mathrm{Ein}}(r)=\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)| Parámetro | Símbolo | Estimación de la Vía Láctea | Papel |
|---|---|---|---|
| Índice de forma | α | Depende del modelo | Controla la rapidez con la que cambia la pendiente |
| Radio de escala | r-2 | ~18-22 kpc | Radio donde la pendiente logarítmica es igual a -2 |
| Densidad en r-2 | ρ-2 | Calibrado a la densidad local | Normalización |
Tensión observacional reciente
Estudios recientes basados en Gaia sugieren que la curva de rotación de la Vía Láctea puede declinar más rápidamente más allá del radio solar de lo que predeciría un halo NFW estándar. Esto hace que los perfiles con núcleo o de variación suave como el de Einasto sean especialmente importantes en los debates actuales.
2.3 Perfil pseudoisotérmico
El perfil pseudoisotérmico se utiliza a menudo como una simple aproximación analítica para un halo con núcleo.
\(\rho_{\mathrm{iso}}(r)=\frac{\rho_0}{1+\left(\frac{r}{r_s}\right)^2}\)En radios pequeños, la densidad se aproxima a un valor constante. A radio grande, cae como r-² y produce una curva de rotación plana.
\(V_{\infty}=\sqrt{4\pi G\rho_0 r_s^2}\)Problema de la cúspide frente al núcleo
Las simulaciones de cuerpo N a menudo predicen perfiles de NFW cuspidados, mientras que muchas galaxias enanas observadas parecen preferir perfiles de densidad con núcleo. Este problema cúspide-núcleo sigue siendo una de las principales cuestiones sin resolver en la física de la materia oscura.
3. Masa del anillo y del anillo – dM/dR
Para calcular cuánta materia oscura hay en cada corte radial de la galaxia, integramos la densidad sobre un fino caparazón o anillo. La geometría depende de si el halo se trata como esférico o aplanado.
3.1 Caparazón esférico delgado
Para un halo esféricamente simétrico, la masa en una envoltura de espesor dr en el radio r es:
\(\frac{dM_{\mathrm{shell}}}{dr}=4\pi r^2\rho(r)\)3.2 Anillo anular disco-plano
Para un anillo situado en el plano galáctico, con radio cilíndrico R y semiespesor efectivo H(R), la masa anular es:
\(dM_{\mathrm{ann}}=2\pi R\,\rho(R,0)\,2H(R)\,dR\)Para un halo esférico, esto puede escribirse como una integral sobre la altura z:
\(\frac{dM}{dR}=2\pi R\int_{-\infty}^{+\infty}\rho\left(\sqrt{R^2+z^2}\right)dz\)En la aproximación esférica, esto conecta de nuevo con:
\(\frac{dM}{dR}\approx4\pi R^2\rho(R)\)3.3 Masa de NFW por concha
\(\frac{dM_{\mathrm{NFW}}}{dr}=4\pi r^2\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}=\frac{4\pi\rho_0 r_s r}{\left(1+\frac{r}{r_s}\right)^2}\)Esta función alcanza su máximo en torno al radio de escala rs, lo que significa que gran parte de la masa de materia oscura por envoltura se deposita en el halo intermedio y no sólo en el centro o en la periferia.
3.4 Masa de einasto por concha
\(\frac{dM_{\mathrm{Ein}}}{dr}=4\pi r^2\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)La masa encerrada Einasto suele evaluarse numéricamente.
Significado físico
La función dM/dr nos indica qué radio galáctico contribuye más al presupuesto de masa oculta. Un perfil exterior más pronunciado reduce la masa total inferida del halo, mientras que un perfil menos pronunciado la aumenta.
4. Masa de materia oscura encerrada M(
Integrando el elemento de envoltura de 0 a R se obtiene la masa total de materia oscura encerrada en el radio R:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)
4.1 Masa cerrada NFW
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)
4.2 Einasto Masa Encerrada
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)
4.3 Descomposición de la masa total
La masa dinámica total encerrada puede descomponerse en componentes visibles y ocultos:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)
~3 × 10¹⁰ M⊙
Masa oculta aproximada dentro del círculo solar.
~1-2 × 10¹¹ M⊙
Masa oculta aproximada dentro de 20 kpc.
5-12 × 10¹¹ M⊙
Masa oculta aproximada dentro de la región virial, fuertemente dependiente del modelo.
El perfil de masa sigue dependiendo del modelo.
Las estimaciones de la masa del halo de la Vía Láctea dependen en gran medida de cómo se extrapole el halo exterior más allá de la región con fuertes limitaciones observacionales.
Integrando el elemento de envoltura de 0 a R se obtiene la masa total de materia oscura encerrada en el radio R:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)4.1 Masa cerrada NFW
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)4.2 Einasto Masa Encerrada
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)4.3 Descomposición de la masa total
La masa dinámica total encerrada puede descomponerse en componentes visibles y ocultos:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)~3 × 10¹⁰ M⊙
Masa oculta aproximada dentro del círculo solar.
~1-2 × 10¹¹ M⊙
Masa oculta aproximada dentro de 20 kpc.
5-12 × 10¹¹ M⊙
Masa oculta aproximada dentro de la región virial, fuertemente dependiente del modelo.
El perfil de masa sigue dependiendo del modelo.
Las estimaciones de la masa del halo de la Vía Láctea dependen en gran medida de cómo se extrapole el halo exterior más allá de la región con fuertes limitaciones observacionales.
5. La curva de rotación V(R)
La velocidad circular en el radio R es fijada por la masa total encerrada mediante el equilibrio de la atracción gravitatoria y la aceleración centrípeta:
\(V_c(R)=\sqrt{\frac{G\,M_{\mathrm{tot}}(<R)}{R}}\)Dado que los componentes de masa independientes contribuyen al potencial gravitatorio, sus contribuciones de velocidad se suman a menudo en cuadratura:
\(V_c^2(R)=V_{\mathrm{bulge}}^2(R)+V_{\mathrm{thin\,disk}}^2(R)+V_{\mathrm{thick\,disk}}^2(R)+V_{\mathrm{gas}}^2(R)+V_{\mathrm{DM}}^2(R)\)5.1 Contribución del disco bariónico
El disco fino estelar sigue un perfil de densidad superficial exponencial:
\(\Sigma(R)=\Sigma_0\exp\left(-\frac{R}{R_d}\right)\)La velocidad circular correspondiente para un disco exponencial es:
\(V_{\mathrm{disk}}^2(R)=\frac{2GM_d}{R_d}y^2\left[I_0(y)K_0(y)-I_1(y)K_1(y)\right],\qquad y=\frac{R}{2R_d}\)AquíIn yKn son funciones de Bessel modificadas. Los parámetros típicos del disco fino de la Vía Láctea son Rd ≈ 2,6 kpc y Md ≈ 3,5 × 10¹⁰ M⊙.
5.2 Contribución de la materia oscura
\(V_{\mathrm{DM,NFW}}(R)=\sqrt{\frac{4\pi G\rho_0r_s^3}{R}\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]}\)5.3 La relación bariónica Tully-Fisher
La relación bariónica Tully-Fisher relaciona la velocidad de rotación plana de una galaxia con su masa bariónica total:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0,\qquad a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)~230 km/s
Velocidad circular cerca del radio solar.
~170-180 km/s
Posible valor decreciente en el disco exterior, en función de los datos del trazador.
~150 km/s
Escala aproximada de velocidad del halo exterior a partir de trazadores del halo.
6. Estimaciones observacionales actuales
La tabla siguiente resume los valores representativos de la densidad y la masa de la materia oscura en radios galácticos clave. Los valores exactos varían según el conjunto de datos, la población de trazadores y el modelo de halo.
| Radio R | Densidad de materia oscura | Masa oscura cerrada | Método |
|---|---|---|---|
| Centro | Divergente en NFW, finito en modelos con núcleo | Depende del modelo | Simulaciones de cuerpo N y modelización del interior de la galaxia |
| R⊙ ≈ 8 kpc | ~0,39 GeV/cm³ | ~3 × 10¹⁰ M⊙ | Curva de rotación y cinemática vertical |
| 20 kpc | ~0,05 GeV/cm³ | ~1-2 × 10¹¹ M⊙ | Gaia y los trazadores espectroscópicos |
| 50 kpc | ~5 × 10-³ GeV/cm³ | ~3-5 × 10¹¹ M⊙ | Cúmulos globulares y halo estelar |
| 100-200 kpc | ≤10-³ GeV/cm³ | ~5-12 × 10¹¹ M⊙ | Galaxias satélites y métodos de velocidad de escape |
La combinación de la cinemática de los cúmulos globulares, las estrellas del halo, las galaxias satélite y la astrometría de la era Gaia sugiere que el perfil del halo exterior de la Vía Láctea sigue siendo incierto. Esta incertidumbre es fundamental para el problema de la masa oculta.
7. Hipótesis contrapuestas para la masa desaparecida
Varias familias principales de explicaciones siguen activas. Ninguna ha sido definitivamente confirmada o descartada en todas las escalas de observación.
7.1 Partículas frías de materia oscura
La materia oscura fría sigue siendo el paradigma principal. Entre las partículas candidatas se incluyen las WIMP, los neutrinos estériles y otras posibilidades más allá del modelo estándar. Estos candidatos forman halos extendidos a menudo modelados con perfiles NFW o Einasto.
\(m_{\chi}\sim10\text{–}1000\,\mathrm{GeV}\)La principal tensión es experimental: la detección directa aún no ha encontrado una partícula de materia oscura confirmada.
7.2 Materia oscura ultraligera o difusa
La materia oscura difusa utiliza partículas ultraligeras cuya longitud de onda de Broglie puede llegar a ser astrofísicamente grande, suprimiendo la estructura a pequeña escala.
\(m_a\sim10^{-22}\,\mathrm{eV}\) \(\lambda_{\mathrm{dB}}\sim\mathrm{kpc}\)Este marco produce de forma natural núcleos de densidad interna más suaves, pero está limitado por los datos del bosque Lyman-alfa y la estructura de las galaxias enanas.
7.3 Dinámica newtoniana modificada
MOND modifica la aceleración gravitatoria efectiva por debajo de una escala característica:
\(a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)En el régimen profundo de Mond, la aceleración efectiva pasa a ser:
\(g_{\mathrm{eff}}=\sqrt{g_Na_0}\)MOND predice la relación bariónica Tully-Fisher:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0\)Funciona bien para muchas curvas de rotación de galaxias, pero los cúmulos de galaxias y la cosmología siguen siendo difíciles.
7.4 Materia oscura autointeractiva
La autointeracción de la materia oscura propone que las partículas de materia oscura interactúen entre sí con la fuerza suficiente para remodelar los perfiles de densidad del halo interior.
\(\frac{\sigma}{m}\sim1\text{–}100\,\mathrm{cm^2/g}\)Esto puede ayudar a explicar la diversidad de los núcleos del halo, pero aún no se ha confirmado ninguna partícula candidata específica.
7.5 Agujeros negros primordiales
Los agujeros negros primordiales formados en el universo primitivo podrían constituir parte de la masa oculta. Muchas ventanas de masa están fuertemente constreñidas por observaciones de microlentes, fondo cósmico de microondas y ondas gravitacionales.
\(10^{-16}\text{–}10^{-11}\,M_\odot\)Siguen siendo especulativas como explicación completa de la masa oculta de la Vía Láctea.
8. La perspectiva de BeeTheory
La Teoría de la Abeja propone que la gravedad puede entenderse como un efecto emergente derivado del comportamiento ondulatorio y no como una fuerza fundamental transportada sólo por una partícula o producida únicamente por la curvatura del espaciotiempo.
En este marco, todo sistema masivo está asociado a una función de onda ψ(r,t). Un punto de partida cuántico básico es la ecuación tridimensional de Schrödinger:
\(i\hbar\frac{\partial\psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\psi+V(\mathbf{r})\psi\)Cuando dos distribuciones de masa se aproximan, sus funciones de onda se superponen. La convolución de estas funciones de onda puede escribirse como:
\(\Psi_{12}(\mathbf{r})=(\psi_1*\psi_2)(\mathbf{r})=\int\psi_1(\mathbf{r}’)\psi_2(\mathbf{r}-\mathbf{r}’)\,d^3r’\)La Teoría de la Abeja interpreta la atracción gravitatoria como una manifestación a gran escala de la superposición de ondas estructuradas, la resonancia y la coherencia de campo.
8.1 Reinterpretación de la teoría de la abeja sobre la masa oculta
En la Teoría de la Abeja, lo que suele denominarse materia oscura puede interpretarse como el efecto gravitatorio acumulativo de la interferencia de ondas de muchos sistemas oscilatorios distribuidos por el halo galáctico.
\(\rho_{\mathrm{eff}}(R)=\rho_{\mathrm{bar}}(R)+\Delta\rho_{\mathrm{wave}}(R)\)Aquí Δρwave(R) representa una densidad gravitatoria efectiva adicional que surge de la estructura coherente del campo de ondas y no de la materia bariónica directamente visible.
Este término tendría que reproducir el comportamiento radial atribuido normalmente a la materia oscura. En particular, tendría que generar curvas de rotación aproximadamente planas en el rango galáctico pertinente.
\(\rho_{\mathrm{wave}}(R)\propto R^{-2}\)Desafío cuantitativo abierto
BeeTheory debe demostrar si un modelo de interferencia basado en ondas puede reproducir el perfil de densidad radial preciso que exigen las curvas de rotación observadas. También debe explicar por qué la masa oculta efectiva es a menudo mucho mayor que la masa bariónica visible.
Referencias
- Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – El perfil de materia oscura de la Vía Láctea inferido a partir de su curva de velocidad circular, MNRAS 528, 693-710, 2024.
- Navarro, J. F., Frenk, C. S., White, S. D. M. – Un perfil de densidad universal a partir de la agrupación jerárquica, ApJ 490, 493, 1997.
- Einasto, J. – Sobre la construcción de un modelo compuesto para la Galaxia, Trudy 5, 87, 1965.
- Watkins, L. L., van der Marel, R. P. et al. – Pruebas de una anticorrelación entre las masas de la Vía Láctea y las galaxias de Andrómeda, ApJ 873, 111, 2019.
- Milgrom, M. – Una modificación de la dinámica newtoniana como posible alternativa a la hipótesis de la masa oculta, ApJ 270, 365, 1983.
- McGaugh, S. S. et al. – Radial Acceleration Relation in Rotationally Supported Galaxies, PRL 117, 201101, 2016.
Nota: las referencias recientes o con fecha futura deben verificarse antes de su publicación si la página se utiliza como fuente de citas científicas.
Perspectiva de BeeTheory
El problema de la masa oculta no es sólo una cuestión de cuánta materia falta. Es una cuestión de qué tipo de estructura física produce gravedad a escala galáctica.
Los modelos clásicos de materia oscura interpretan la masa desaparecida como materia invisible. La Teoría de la Abeja explora una posibilidad complementaria: parte del efecto gravitatorio oculto puede surgir de la coherencia estructurada de las ondas.
El siguiente paso es matemático: definir el término de densidad de onda radial, derivar su curva de rotación y compararlo directamente con los datos de la Vía Láctea de la era Gaia.