BeeTheory – Galaxias SPARC – Ajuste – 2025

Teoría de la abeja aplicada a 20 galaxias externas:
Fórmula ajustada y metodología de prueba ciega

El catálogo SPARC proporciona 175 galaxias con perfiles bariónicos y curvas de rotación medidos. Aplicamos la ecuación de masa oscura de BeeTheory -ajustando su ley de escala para que se ajuste a la población de galaxias- e informamos del resultado: 18 de 20 galaxias predichas dentro del 20% de su velocidad de rotación plana observada, con χ²/dof = 0,93.

0. Resultados

Mejor ajuste – 20 galaxias SPARC, Q = 1, alta calidad

Con la fórmula BeeTheory modificada Kd = K0/Rd y ℓd = c – Rd, dos constantes universales se ajustan simultáneamente a las 20 galaxias.

La densidad de masa oscura en cada galaxia se predice únicamente a partir de los parámetros de su disco bariónico – sin ajuste por galaxia.

Parámetros de mejor ajuste: K0 = 0,3759, adimensional, y c = 6,40, adimensional. Resultado: 18/20 galaxias dentro del 20% del Vf observado, χ²/dof = 0,93. Dos valores atípicos, CamB y NGC 3741, son enanas dominadas por el gas en las que la modelización del disco estelar se rompe.

18/20

Dentro del 20% de Vf

6.8%

Error medio

0.93

χ²/dof

2

Constantes universales

0.91

Pearson r, Tully-Fisher

Fórmula BeeTheory modificada ajustada a la población de galaxias \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{R_{\mathrm{max}}}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_d D)e^{-\alpha_d D}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+R’^2},\qquad \alpha_d=\frac{1}{cR_d},\qquad K_0=0.3759,\qquad c=6.40\)

1. Las 20 galaxias SPARC – Datos y predicciones

La muestra SPARC abarca galaxias que abarcan cinco décadas en luminosidad, desde irregulares enanas hasta espirales masivas. Para cada galaxia, los parámetros de entrada se toman directamente de la tabla 1 de Lelli et al. 2016: radio de escala del disco Rd, brillo de la superficie central Σd, masa de gas HI MHI y velocidad de rotación plana Vf.

La masa estelar se calcula como M★ = Υ★ × L3,6, con Υ★ = 0,5 M⊙/L⊙. La masa de gas se calcula como Mgas = 1,33 × MHI.

Galaxia Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Error Estado
Cargando tabla de galaxias…

Todas las velocidades en km/s. Error = (VBTVf)/Vf. Parámetros: Kd = 0,3759/Rd, ℓd = 6,40 × Rd. Predicción de BeeTheory evaluada enReval = 5Rd.

VBT frente a Vf observado – 20 galaxias SPARC, BeeTheory ajustada
Dentro del 20%, 18 galaxias Valores atípicos: CamB, NGC 3741 Predicción perfecta, 1:1 Envolvente de ±20

2. La fórmula modificada – Por qué K ∝ 1/Rd

El ajuste original de BeeTheory para la Vía Láctea utilizaba una única constante de acoplamiento K = 0,02365 kpc-¹ con una longitud de coherencia ℓ = 3,17Rd. Cuando se aplicó a ciegas a 20 galaxias SPARC, esto subestimó sistemáticamente Vf en aproximadamente un 50%.

El análisis por galaxia reveló un patrón claro: la constante de acoplamiento necesaria varía como K ∝ 1/Rd.

2.1 De una constante a una ley de escala

La idea clave es dimensional. La densidad oscura BeeTheory en el radio r de un disco de escala Rd y densidad superficial Σ0 es, en el régimen asintótico de rotación plana r ≪ ℓ:

Densidad oscura asintótica, Rd ≪ r ≪ ℓ \(\rho_{\mathrm{dark}}(r)\approx \frac{K\Sigma_0 R_d^2}{r^2}f\left(\frac{r}{\ell}\right)\) \(M_{mathrm{dark}(<r)\propto K\Sigma_0R_d^2r\)

La velocidad de rotación plana escala entonces como:

Velocidad de rotación plana a partir de la masa oscura de BeeTheory \(V_f^2=\frac{GM_{\mathrm{dark}}(<R_{\mathrm{eval}})}{R_{\mathrm{eval}}}\propto GK\Sigma_0R_d^2=GK\frac{M_\star}{2\pi}\) \(V_f\propto (K M_\star)^{1/2}\)

La relación bariónica Tully-Fisher observada establece Vf4Mbar, lo que significa VfMbar1/4. Para que esto sea reproducido por la Teoría de la Abeja, necesitamos Vf2M★/Rd, la densidad superficial media del disco. Esto requiere

Escala necesaria para reproducir la pendiente de Tully-Fisher \(V_f^2\propto GK\frac{M_\star}{2\pi}\propto\frac{M_\star}{R_d}\implies K\propto\frac{1}{R_d}\) \(\boxed{K_d=\frac{K_0}{R_d},\qquad K_0=0.3759}\)

Este escalado no es un parche ad hoc, es lo que exige la relación Tully-Fisher. Un acoplamiento K ∝ 1/Rd significa que los discos más compactos generan un campo oscuro más intenso por unidad de masa.

2.2 La longitud de coherencia – Por qué c = 6,40 ≠ 3,17

El ajuste de la Vía Láctea dio cMW = ℓd/Rd = 3,17. La muestra SPARC da cSPARC = 6,40, aproximadamente el doble. Dos explicaciones son posibles:

  • Sesgo de selección: las 20 galaxias SPARC se eligieron por sus curvas de rotación extendidas de alta calidad, lo que supone un sesgo hacia las galaxias con discos HI más extendidos.
  • Contribución del disco de gas: en muchas galaxias SPARC, el disco de gas HI tiene un radio de escala RHI ≈ 1,7Rd. Incluir el gas como una fuente de disco separada aumentaría el tamaño efectivo de la fuente.

Ambos efectos son reales y mensurables. El valor definitivo de c requiere modelizar por separado el gas y los discos estelares.

Fórmula completa de la Teoría de la Abeja modificada – dos constantes universales \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{8R_d}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(\alpha_d=\frac{1}{cR_d},\qquad D=\sqrt{r^2+R’^2},\qquad K_0=0.3759,\qquad c=6.40\)

3. El cálculo – Paso a paso

Para cada galaxia SPARC, la predicción BeeTheory procede en cinco pasos. No se ajustan parámetros libres por galaxia.

1
Leer las entradas bariónicas de la tabla 1 del SPARC

Rd, Σd, MHI y Vf. Convierta Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc², y Mgas = 1,33 × MHI.

2
Calcule los parámetros de BeeTheory a partir de Rd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd, y αd = 1/ℓd. No hay ajuste.

3
Integrar la densidad oscura en r = 5Rd
\(\rho_{\mathrm{dark}}(5R_d)=K_d\sum_{i=1}^{60}\Sigma_0e^{-R_i’/R_d}\frac{(1+\alpha_dD_i)e^{-\alpha_dD_i}}{D_i^2}2\pi R_i’\Delta R’\)

Integración numérica con 60 anillos, R′ ∈ [0, 8Rd]. Luego integrar esféricamente para obtener la masa oscura encerrada Mdark(<5Rd).

4
Calcule la velocidad circular bariónica
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Predecir la velocidad circular total
\(V_{\mathrm{BT}}=\sqrt{V_{\mathrm{bar}}^2+V_{\mathrm{dark}}^2},\qquad V_{\mathrm{dark}}=\sqrt{\frac{GM_{\mathrm{dark}}(<5R_d)}{5R_d}}\)

Compárelo con el Vf observado. Error = (VBTVf)/Vf.

Relación Vdark/Vbar – ¿Hasta qué punto es dominante el componente oscuro?

4. Por qué la prueba ciega es la única prueba honesta

Un modelo que reproduce los datos con los que fue calibrado no prueba nada. Todo modelo, incluso uno erróneo, puede ajustarse a sus datos de entrenamiento. La única prueba científicamente significativa es una predicción a ciegas: aplicar el modelo a datos que nunca ha visto, con los parámetros congelados de la calibración, e informar del resultado, sea cual sea.

4.1 Qué significa «ciego» aquí

Los parámetros K0 = 0,3759 y c = 6,40 de BeeTheory se determinaron ajustando simultáneamente las 20 galaxias SPARC. Ahora están fijados.

La prueba a ciegas sería: aplicar estos parámetros a las 155 galaxias SPARC restantes, que no se utilizaron en el ajuste, e informar del resultado antes de observar sus curvas de rotación observadas. Esta prueba aún no se ha realizado: es el siguiente paso.

Los parámetros originales de la Vía Láctea, Kd = 0,02365 y ℓd = 3,17Rd, se determinaron en una sola galaxia. Aplicarlos al SPARC sin ajustes dio 0/20 galaxias correctas – un fallo honesto e importante. Ese fallo reveló el escalado K ∝ 1/Rd.

4.2 Significado estadístico de la calidad del ajuste

Con χ²/dof = 0,93 en 20 galaxias, el modelo se ajusta aproximadamente al nivel esperado de las incertidumbres de velocidad del 15% asumidas.

Interpretación Chi-cuadrado \(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{BT}}(i)-V_f(i)}{0.15V_f(i)}\right)^2=0.93\) \(N=19\ \text{(excluyendo CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

Un valor de 0,93 está muy cerca del ideal 1,0. El modelo da cuenta de la dispersión al nivel de la incertidumbre de medición.

4.3 Los dos valores atípicos

CamB – enana de gas puro, Vf = 2,0 km/s

CamB casi no tiene masa estelar, M★ ≈ 2×10⁷ M⊙. La fórmula de BeeTheory utiliza Σ0e-R/Rd como fuente – pero en CamB, los bariones son casi en su totalidad gas HI, no estrellas. El modelo de disco estelar es inaplicable.

NGC 3741 – sobreestimada en un 47

NGC 3741 es una pequeña enana de bajo brillo superficial con un disco HI muy extendido. La fuente de BeeTheory, el disco estelar, subestima la extensión bariónica real. La inclusión del disco de gas como componente separado de la fuente con un radio de escala mayor reduciría la masa oscura prevista y corregiría la sobreestimación.

Las otras 18 – predicciones auténticas

Para las 18 galaxias dentro del 20%, el error medio es del 6,8%, muy dentro de las incertidumbres observacionales. Éstas abarcan Rd de 1,3 a 5,8 kpc y Vf de 76 a 278 km/s. BeeTheory predice correctamente este intervalo de un factor de 3,7 en la velocidad -la pendiente Tully-Fisher- con dos constantes universales.

5. Significado físico – Lo que revela la escala

5.1 El acoplamiento universal adimensional

Con Kd = K0/Rd y ℓd = cRd, el acoplamiento adimensional BeeTheory es:

Acoplamiento efectivo – escala con el tamaño de la galaxia \(\lambda_{\mathrm{eff}}=K_d\ell_d^2=\frac{K_0}{R_d}(cR_d)^2=K_0c^2R_d\) \(\lambda_{\mathrm{eff}}=0.3759\times41.0\times R_d=15.4R_d\ \text{(kpc)}\)

λeff crece con Rd. Las galaxias más grandes generan proporcionalmente más masa oscura. Esta es la predicción de la Teoría de la Abeja de por qué las espirales masivas están más dominadas por la materia oscura que las enanas.

5.2 Conexión con la relación de aceleración radial

McGaugh et al. descubrieron que la aceleración centrípeta observadagobs = Vc2/R es una función universal de la contribución bariónica gbar = GMbar/R². En la Teoría de la Abeja, esta relación surge porque:

La predicción BeeTheory de la RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

El escalado gdarkgbar1/2 produce la forma RAR observada. Derivar la curva RAR exacta a partir de BeeTheory es la siguiente tarea teórica inmediata.

Datos: Lelli, F., McGaugh, S. S., Schombert, J. M., SPARC: Modelos de masa para 175 galaxias de disco con fotometría Spitzer y curvas de rotación precisas, AJ 152, 157, 2016.

RAR: McGaugh, S. S., Lelli, F., Schombert, J. M., PRL 117, 201101, 2016.

BTFR: Lelli, F. et al., ApJ 816, 2016.

BeeTheory: Dutertre, X., BeeTheory.com v2, 2023, ampliado 2025.

Masa a luz: Υ★ = 0,5 M⊙/L⊙ a 3,6 μm, McGaugh & Schombert 2014.