BeeTheory – Galáxias SPARC – Ajuste ajustado – 2025

BeeTheory aplicada a 20 galáxias externas:
Fórmula ajustada e metodologia de teste cego

O catálogo SPARC fornece 175 galáxias com perfis bariônicos e curvas de rotação medidos. Aplicamos a equação de massa escura BeeTheory – ajustando sua lei de escala para corresponder à população de galáxias – e relatamos o resultado: 18 de 20 galáxias previstas dentro de 20% de sua velocidade de rotação plana observada, com χ²/dof = 0,93.

0. Resultados – Primeiro, o senhor disse

Melhor ajuste – 20 galáxias SPARC, Q = 1, alta qualidade

Com a fórmula BeeTheory modificada Kd = K0/Rd e ℓd = c – Rd, duas constantes universais se ajustam a todas as 20 galáxias simultaneamente.

A densidade de massa escura em cada galáxia é prevista somente a partir dos parâmetros do disco bariônico – sem ajuste por galáxia.

Parâmetros de melhor ajuste: K0 = 0,3759, sem dimensão, e c = 6,40, sem dimensão. Resultado: 18/20 galáxias dentro de 20% do Vf observado, χ²/dof = 0,93. Duas exceções, CamB e NGC 3741, são anãs dominadas por gás, onde a modelagem do disco estelar não funciona.

18/20

Dentro de 20% de Vf

6.8%

Erro mediano

0.93

χ²/dof

2

Constantes universais

0.91

Pearson r, Tully-Fisher

Fórmula BeeTheory modificada – ajustada para a população de galáxias \(\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. As 20 galáxias SPARC – Dados e previsões

A amostra SPARC abrange galáxias que abrangem cinco décadas de luminosidade, de irregulares anãs a espirais maciças. Para cada galáxia, os parâmetros de entrada são obtidos diretamente da Tabela 1 de Lelli et al. 2016: raio de escala do disco Rd, brilho da superfície central Σd, massa de gás HI MHI e velocidade de rotação plana Vf.

A massa estelar é calculada como M★ = Υ★ × L3,6, com Υ★ = 0,5 M⊙/L⊙. A massa de gás é calculada como Mgas = 1,33 × MHI.

Galáxia Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Erro Status
Carregando a tabela de galáxias…

Todas as velocidades em km/s. Erro = (VBT – Vf)/Vf. Parâmetros: Kd = 0,3759/Rd, ℓd = 6,40 × Rd. Previsão de BeeTheory avaliada emReval = 5Rd.

VBT vs Vf observada – 20 galáxias SPARC, BeeTheory ajustada
Dentro de 20%, 18 galáxias Excedentes: CamB, NGC 3741 Previsão perfeita, 1:1 ±20% de envelope

2. A fórmula modificada – Por que K ∝ 1/Rd

O ajuste original do BeeTheory para a Via Láctea usou uma única constante de acoplamento K = 0,02365 kpc-¹ com comprimento de coerência ℓ = 3,17Rd. Quando aplicado cegamente a 20 galáxias SPARC, esse ajuste subestimou sistematicamente a Vf em cerca de 50%.

A análise por galáxia revelou um padrão claro: a constante de acoplamento necessária varia conforme K ∝ 1/Rd.

2.1 De uma constante para uma lei de escala

O principal insight é dimensional. A densidade escura BeeTheory no raio r de um disco de escala Rd e densidade de superfície Σ0 é, no regime de rotação plana assintótica, r ≪ ℓ:

Densidade escura assintó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\)

A velocidade de rotação plana é então escalonada como:

Velocidade de rotação plana da massa escura da 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}\)

A relação Tully-Fisher bariônica observada indica Vf4 ∝ Mbar, o que significa Vf ∝ Mbar1/4. Para que isso seja reproduzido pela BeeTheory, precisamos de Vf2 ∝ M★/Rd, a densidade média da superfície do disco. Isso requer:

Escalonamento necessário para reproduzir a inclinação 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}\)

Essa escala não é um ajuste ad hoc – é o que a relação Tully-Fisher exige. Um acoplamento K ∝ 1/Rd significa que os discos mais compactos geram um campo escuro mais forte por unidade de massa.

2.2 O comprimento de coerência – Por que c = 6,40 ≠ 3,17

O ajuste da Via Láctea deu cMW = ℓd/Rd = 3,17. A amostra SPARC fornece cSPARC = 6,40, cerca de duas vezes maior. Duas explicações são possíveis:

  • Viés de seleção: as 20 galáxias do SPARC foram escolhidas por apresentarem curvas de rotação estendidas de alta qualidade, o que leva a um viés em direção a galáxias com discos HI mais estendidos.
  • Contribuição do disco de gás: em muitas galáxias SPARC, o disco de gás HI tem um raio de escalaRHI ≈ 1,7Rd. A inclusão do gás como uma fonte de disco separada aumentaria o tamanho efetivo da fonte.

Ambos os efeitos são reais e mensuráveis. O valor definitivo de c requer a modelagem dos discos de gás e estelares separadamente.

Fórmula completa da BeeTheory modificada – duas constantes universais \(\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. O cálculo – passo a passo

Para cada galáxia SPARC, a previsão do BeeTheory prossegue em cinco etapas. Nenhum parâmetro livre é ajustado por galáxia.

1
Leia as entradas bariônicas da Tabela 1 do SPARC

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

2
Calcule os parâmetros de BeeTheory a partir de Rd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd e αd = 1/ℓd. Sem ajuste.

3
Integre a densidade escura em 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’\)

Integração numérica com 60 anéis, R′ ∈ [0, 8Rd]. Em seguida, integre esfericamente para obter a massa escura fechada Mdark(<5Rd).

4
Calcule a velocidade circular bariônica
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Prever a velocidade 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}}\)

Compare com a Vf observada. Erro = (VBT – Vf)/Vf.

Relação Vdark/Vbar – Quão dominante é o componente escuro?

4. Por que o teste cego é o único teste honesto

Um modelo que reproduz os dados com os quais foi calibrado não prova nada. Todo modelo, mesmo um modelo errado, pode ser ajustado para se adequar aos dados de treinamento. O único teste cientificamente significativo é uma previsão cega: aplicar o modelo a dados que ele nunca viu, com parâmetros congelados da calibração, e relatar o resultado, seja ele qual for.

4.1 O que significa “cego” aqui

Os parâmetros do BeeTheory K0 = 0,3759 e c = 6,40 foram determinados pelo ajuste simultâneo das 20 galáxias SPARC. Agora eles estão fixados.

O teste cego seria: aplicar esses parâmetros às 155 galáxias SPARC restantes, que não foram usadas no ajuste, e relatar o resultado antes de analisar suas curvas de rotação observadas. Esse teste ainda não foi realizado – é a próxima etapa.

Os parâmetros originais da Via Láctea, Kd = 0,02365 e ℓd = 3,17Rd, foram determinados em uma única galáxia. Ao aplicá-los ao SPARC sem ajustes, 0/20 galáxias foram corrigidas – uma falha honesta e importante. Essa falha revelou o escalonamento K ∝ 1/Rd.

4.2 Significado estatístico da qualidade do ajuste

Com χ²/dof = 0,93 em 20 galáxias, o modelo se ajusta aproximadamente ao nível esperado das incertezas de velocidade de 15% assumidas.

Interpretação do qui-quadrado \(\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{(excluindo CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

Um valor de 0,93 está muito próximo do valor ideal de 1,0. O modelo é responsável pela dispersão no nível de incerteza de medição.

4.3 Os dois outliers

CamB – anã de gás puro, Vf = 2,0 km/s

CamB quase não tem massa estelar, M★ ≈ 2×10⁷ M⊙. A fórmula BeeTheory usa Σ0e-R/Rd como a fonte – mas em CamB, os bárions são quase inteiramente gás HI, não estrelas. O modelo de disco estelar é inaplicável.

NGC 3741 – superestimada em 47%

A NGC 3741 é uma pequena anã de baixo brilho superficial com um disco HI muito extenso. A fonte BeeTheory, o disco estelar, subestima a extensão bariônica real. A inclusão do disco de gás como um componente de fonte separado com raio de escala maior reduziria a massa escura prevista e corrigiria a superestimativa.

As outras 18 – previsões genuínas

Para as 18 galáxias dentro de 20%, o erro médio é de 6,8%, bem dentro das incertezas observacionais. Elas abrangem Rd de 1,3 a 5,8 kpc e Vf de 76 a 278 km/s. A BeeTheory prevê corretamente esse fator de variação de 3,7 na velocidade – a inclinação de Tully-Fisher – com duas constantes universais.

5. Significado físico – O que a escala revela

5.1 O acoplamento universal sem dimensões

Com Kd = K0/Rd e ℓd = cRd, o acoplamento BeeTheory sem dimensão é:

Acoplamento efetivo – escala com o tamanho da galáxia \(\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 cresce com Rd. Galáxias maiores geram proporcionalmente mais massa escura. Essa é a previsão da BeeTheory para explicar por que as espirais maciças são mais dominadas pela matéria escura do que as anãs.

5.2 Conexão com a relação de aceleração radial

McGaugh et al. descobriram que a aceleração centrípeta observadagobs = Vc2/R é uma função universal da contribuição bariônica gbar = GMbar/R². Na BeeTheory, essa relação surge porque:

Previsão da BeeTheory para a RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

O escalonamento gdark ∝ gbar1/2 produz a forma RAR observada. Derivar a curva RAR exata da BeeTheory é a próxima tarefa teórica imediata.

Dados: Lelli, F., McGaugh, S. S., Schombert, J. M., SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves (Modelos de massa para 175 galáxias de disco com fotometria Spitzer e curvas de rotação 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 em 2025.

Massa em relação à luz: Υ★ = 0,5 M⊙/L⊙ a 3,6 μm, McGaugh & Schombert 2014.

BeeTheory.com – Gravidade quântica baseada em ondas – © Technoplane S.A.S. 2025