BeeTheory – Análisis residual – 2025

Leyendo los residuales:
Por qué algunas galaxias están por encima y otras por debajo de la predicción

En el ajuste SPARC de 20 galaxias, 15 galaxias están subestimadas y 3 sobreestimadas, excluyendo los dos grandes valores atípicos. Esta asimetría no es aleatoria. Dos parámetros físicos determinan el patrón, y apuntan directamente a la física que falta en el modelo BeeTheory actual.

0. Los dos grupos – Expuestos en primer lugar

Por encima de la línea – el modelo sobreestima VBT > Vf

BeeTheory predice más masa oscura de la que realmente tiene la galaxia. Las 3 galaxias de la muestra del núcleo son:

  • NGC 4051 – +15,8% – baja fracción de gas, AGN Seyfert
  • NGC 0100 – +6,7% – borde de bajo brillo superficial
  • NGC 0300 – +0,9% – espiral floculenta, aislada

Rasgo común: baja fracción de gas, fgas ≈ 0,37, sin protuberancia, disco estelar puro. El modelo atribuye demasiada masa oscura porque sólo ve el disco estelar y no el hecho de que el disco por sí solo no puede generar tanto campo como se supone.

Por debajo de la línea – el modelo subestima VBT < Vf

BeeTheory predice menos masa oscura de la que realmente tiene la galaxia. Esto afecta a 15 galaxias, incluidas las 7 con protuberancias:

  • Las 7 galaxias con bulbo – subestimadas sin excepción
  • NGC 3621 – -16,2% – muy rica en gas, fgas = 0,82
  • NGC 3521 – -17,1% – gran disco de gas extendido
  • NGC 0925, NGC 3198 – espirales tardías ricas en gas

Rasgos comunes: o bien tienen una protuberancia, no modelada completamente, o una fracción de gas elevada, fgas ≈ 0,50-0,82. El disco HI extendido no se incluye como fuente de BeeTheory.

1. Los dos impulsores físicos

-0.68

Pearson r: error frente a fracción de gas, galaxias no turbulentas

100%

Galaxias de bulbo subestimadas, 7 de 7

+0.04

Pearson r: error vs Σd, sin señal

Dos efectos físicos independientes determinan los residuos: uno para las galaxias con protuberancias y otro para los sistemas ricos en gas. La densidad superficial Σd no tiene esencialmente ningún poder predictivo para los residuales una vez que se controla la presencia de protuberancias.

Controlador 1 – Presencia de protuberancias

Todas las galaxias con protuberancia detectada están infravaloradas. El modelo actual asigna el 15% de la masa estelar al bulbo, una estimación aproximada. La fracción real de masa del bulbo varía desde aproximadamente el 5% en las espirales tardías hasta cerca del 40% en las galaxias Sa.

Y lo que es más importante, el bulbo genera un campo oscuro BeeTheory con una longitud de coherencia corta, ℓbℓd, intenso a r pequeño y que contribuye significativamente a Vf en el radio de rotación plano. Esto no se capta plenamente en el modelo actual.

En la Vía Láctea, el análisis de dos regímenes mostró que el bulbo contribuye con cerca del 35% de la masa oscura total a r = 8 kpc, a pesar de tener sólo el 18% de la masa bariónica. En estas galaxias se produce la misma amplificación, pero el modelo utiliza unaKb genérica y una fracción del bulbo en lugar de valores específicos de cada galaxia.

Conductor 2 – Fracción de gas

En las galaxias no bulbosas, el residuo se correlaciona con la fracción de gas a r = -0,68. La dirección es clara: más gas significa que el modelo subestima la velocidad.

El modelo BeeTheory actual utiliza el disco estelar como fuente, con escala Rd. Pero en las galaxias ricas en gas, el disco HI se extiende hasta RHI ≈ 1,7-3 × Rd. Este disco de gas extendido es una fuente BeeTheory que no se incluyó.

Su radio de escala mayor implica un ℓgas más grande y un perfil de campo oscuro diferente. Cuando el gas domina la masa bariónica, ignorar el disco de gas subestima significativamente el campo oscuro total.

La aparente paradoja: ¿por qué más gas significa menos masa oscura prevista?

El modelo utiliza Kd = K0/Rd, donde Rd es el radio de escala del disco estelar. Cuando fgas es alto, el disco de gas se extiende más allá del disco estelar, pero el modelo sólo ve el disco estelar.

El disco estelar genera por sí solo un campo oscuro calibrado a Vf. En realidad, el disco de gas también contribuye, y como el gas se extiende más, su ℓ efectiva es mayor, creando un perfil de campo oscuro diferente y más extendido. El modelo atribuye toda la masa oscura a una única fuente estelar y subestima el total.

2. Análisis cuantitativo – Las correlaciones

Error residual frente a la fracción de gas fgas – Sólo galaxias no abombadas
Sin protuberancia – utilizada en la correlación Tiene protuberancia – excluida de la correlación de gas Tendencia de mejor ajuste
Residuos ordenados – La señal del abombamiento
Tiene protuberancia, siempre subestimada Sin protuberancia Línea cero

Comparación: Por encima frente a por debajo – Propiedades medias

Propiedad Sobreestimado ↑ Subestimada ↓ Interpretación
N galaxias315Predomina la subestimación sistemática
Error medio+7.8%-7.1%Magnitud simétrica, recuento asimétrico
Media fgas0.370.50Subestimación de las medias ricas en gas
Media Σd, L⊙/pc²146247La media de los discos más densos está subestimada, principalmente un efecto de protuberancia
Tiene protuberancia0 / 37 / 15Todas las galaxias con protuberancia subestimadas
Media del tipo de Hubble T5,7, Sc5,1, ScSin señal. T no es un impulsor.

3. El mecanismo – Lo que le falta a la fórmula de la Teoría de la Abeja

3.1 El disco de gas como fuente ausente de la teoría de la abeja

La fórmula actual utiliza el disco estelar como única fuente de campo oscuro:

Fórmula actual – sólo fuente del disco estelar \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d^\star}\int \Sigma_0^\star e^{-R’/R_d^\star}\,\mathrm{kernel}\,dR’\)

Pero cada elemento de masa bariónica es una fuente de Teoría de la Abeja. En las galaxias ricas en gas, el disco HI contiene tanta masa como el disco estelar y se extiende hasta RHI ≈ 1,7Rd★. La fórmula correcta debería ser

Fórmula corregida – fuentes estelares + disco de gas \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d^\star}\int\Sigma_\star e^{-R’/R_d^\star}\mathrm{kernel}_d\,dR’+\frac{K_0}{R_{\mathrm{HI}}}\int\Sigma_{\mathrm{HI}}e^{-R’/R_{\mathrm{HI}}}\mathrm{kernel}_{\mathrm{gas}}\,dR’\) \(R_{\mathrm{HI}}\approx1.7R_d^\star,\qquad \mathrm{kernel}_{\mathrm{gas}}=\frac{(1+\alpha_{\mathrm{gas}}D)e^{-\alpha_{\mathrm{gas}}D}}{D^2},\qquad \alpha_{\mathrm{gas}}=\frac{1}{c_{\mathrm{gas}}R_{\mathrm{HI}}}\)

La fuente del disco de gas tiene una longitud de coherencia mayor que la del disco estelar. Su campo oscuro es más extenso y contribuye de forma diferente a gran r.

  • Debería aumentar la masa oscura prevista para las galaxias ricas en gas, reduciendo la subestimación.
  • Debería reducir la sobreestimación para las galaxias de disco estelar puro con poco gas.
  • Explica por qué NGC 0925, NGC 3198 y NGC 3621 están subestimadas.

3.2 El campo oscuro de la protuberancia: una fuente compacta e intensa

En la Teoría de la Abeja, una fuente compacta genera un campo oscuro más intenso por unidad de masa, porque el núcleo de Yukawa da más peso a las distancias pequeñas. La protuberancia, concentrada en un radio de aproximadamente 1-2 kpc, genera un campo oscuro con una longitud de coherencia corta ℓbℓd.

Campo oscuro del bulbo – fuente compacta de alta intensidad \(\rho_{\mathrm{dark,bulge}}(r)=K_b\int_0^{R_b}\rho_{0,b}e^{-r’/r_b}\frac{(1+\alpha_bD)e^{-\alpha_bD}}{D^2}4\pi r’^2\,dr’\) \(\alpha_b=\frac{1}{\ell_b}\gg \alpha_d=\frac{1}{\ell_d}\)

El modelo actual utiliza Kb = 1,055 kpc-¹, calibrado en la Vía Láctea, y asigna el 15% de la masa estelar al bulbo, ambas estimaciones aproximadas.

Por qué es difícil modelizar el bulbo en las galaxias SPARC

SPARC proporciona el radio Rd de la escala del disco estelar y la luminosidad total, pero no una descomposición fiable entre bulbo y disco para todas las galaxias. Separar el bulbo del disco requiere un ajuste fotométrico 2D, disponible para algunas galaxias SPARC pero no de manera uniforme.

4. Lo que esto predice – El modelo corregido

Si se incluyen los dos efectos físicos identificados -el disco de gas como fuente separada de BeeTheory, y el bulbo con una fracción de masa específica de la galaxia- el patrón residual debería desaparecer.

Corrección de las galaxias sobreestimadas
  • Para NGC 4051 y NGC 0100, la fórmula basada únicamente en el disco estelar es casi correcta porque el gas es escaso. La corrección es pequeña.
  • La sobreestimación puede proceder de una ligera sobreestimación de la relación masa estelar/luz Υ★.
  • NGC 0300 ya es esencialmente correcta en un +0,9%.
Corrección de las galaxias subestimadas
  • Las galaxias ricas en gas como NGC 3621, NGC 0925 y NGC 3198 deberían mejorar cuando se incluye una fuente de disco HI.
  • El uso de RHI = 1,7Rd y KHI = K0/RHI puede añadir alrededor de un 10-15% de campo oscuro.
  • Las galaxias de bulbo como NGC 3521 y NGC 0891 requieren fracciones de masa de bulbo específicas para cada galaxia.
Fórmulacompleta corregida de BeeTheory – tres fuentes \(\rho_{\mathrm{dark}}(r)=\underbrace{\frac{K_0}{R_d^\star}\int\Sigma_\star e^{-R’/R_d^\star}K_d\,dR’}_{\mathrm{stellar\ disk}}+\underbrace{\frac{K_0}{R_{\mathrm{HI}}}\int\Sigma_{\mathrm{HI}}e^{-R’/R_{\mathrm{HI}}}K_{\mathrm{gas}}\,dR’}_{\mathrm{HI\ gas\ disk}}+\underbrace{K_b\int\rho_{0,b}e^{-r’/r_b}K_{\mathrm{bulge}}\,dr’}_{\mathrm{bulge}}\)

Esta fórmula de tres fuentes sigue utilizando sólo dos constantes universales -K0 = 0,3759 y c = 6,40- porque RHI y rb son propiedades bariónicas medidas de cada galaxia, no parámetros libres.

Conclusión científica positiva

El patrón residual no es ruido aleatorio. Está estructurado, es explicable y apunta a una física ausente bien definida.

Un modelo que falla de forma estructurada es más valioso que uno que falla al azar: le indica exactamente qué debe mejorar. En este caso, el marco de BeeTheory tiene una estructura correcta; lo que falta es la inclusión de todas las fuentes bariónicas, no sólo el disco estelar dominante.

La predicción es clara: añada el disco de gas HI y la descomposición adecuada del bulbo a la fórmula, y el porcentaje de aciertos de 18/20 debería mejorar, incluidos los dos valores atípicos CamB y NGC 3741.

Datos: Lelli, McGaugh, Schombert, AJ 152, 157, 2016.

Modelo BeeTheory: Dutertre, 2023, ampliado 2025.

Escalado del disco HI: RHI/Rd ≈ 1,7, Broeils & Rhee 1997; Swaters et al. 2009; Lelli et al. 2014.