Nota técnica XXXVIII

Primeiro teste da BeeTheory em relação à Relação de Aceleração Radial

BeeTheory.com – Xavier Dutertre – Technoplane S.A.S. – 20 de maio de 2026

Resultado. Confrontamos o modelo de 3 parâmetros da BeeTheory (\(\lambda = 12,696\), \(c = 0,163\), \(\ell_{\rm floor} = 3,00\) kpc) com a Radial Acceleration Relation (RAR) de McGaugh, Lelli & Schombert (2016). Avaliando a previsão de \(g_{\rm obs}(g_{\rm bar})\) em 909 pontos amostrados em 101 galáxias SPARC sem bulbo, a nuvem BeeTheory tem a forma geral correta – segue a curva de McGaugh em \(g_{\rm bar}\) baixo – mas é tendenciosamente alta em +0,25 dex em média, com uma dispersão de 0,39 dex (vs. 0,13 dex empiricamente). O viés não é uniforme: em baixa aceleração (\(g_{\rm bar} < 10^{-12}\) m/s²), o BeeTheory concorda com o RAR, enquanto em alta aceleração (em torno e acima de \(g_\dagger\)) ele supera a previsão em 0,4 a 0,5 dex. O modelo captura o comportamento assintótico semelhante ao MOND, mas não consegue recuperar o limite de Newton \(g_{\rm obs} \to g_{\rm bar}\) em um raio pequeno. Esse é um sucesso parcial e um indicador claro para o próximo refinamento.

1. O teste

A Radial Acceleration Relation (Relação de Aceleração Radial) expressa a correlação empírica extremamente estreita entre duas quantidades medidas em curvas de rotação: \(g_{\rm bar}(R)\), a aceleração newtoniana que seria produzida apenas pela matéria visível, e \(g_{\rm obs}(R)\), a aceleração centrípeta total inferida a partir da velocidade de rotação. McGaugh, Lelli e Schombert (2016) mostraram que ~150 galáxias SPARC caem em uma única curva bem ajustada por:

\[ g_{\rm obs}(g_{\rm bar}) \;=\; \frac{g_{\rm bar}}{1 – \exp\!\left(-\sqrt{g_{\rm bar}/g_\dagger}\right)},\qquad g_\dagger = (1,20 \pm 0,02)\times 10^{-10}\;{\rm m/s}^2 \]

A dispersão em torno dessa curva é de apenas ~0,13 dex – pouco maior do que as incertezas observacionais. Qualquer teoria de gravidade modificada que pretenda substituir a matéria escura deve reproduzir essa relação ponto a ponto. O RAR é o filtro empírico mais rigoroso do mercado.

A previsão da BeeTheory é direta. Em cada raio \(R\) dentro de uma galáxia, a matéria visível gera o campo de ondas cuja massa contida se soma à massa bariônica contida com acoplamento \(\lambda\):

\[ \begin{aligned} g_{\rm bar}(R) &= G\,M_{\rm bar}(<\!R)\,/\,R^2 \\ g_{\rm obs}^{\rm BT}(R) &= g_{\rm bar}(R) + G\,\lambda\,M_{\rm wave}(<\!R)\,/\,R^2 \\ M_{\rm wave}(<\!R) &= M\,\bigl[\,1 - (1 + x + x^2/2)\,e^{-x}\bigr],\quad x = R/\ell_{\rm wave} \\ \ell_{\rm wave} &= c\,R_d + \ell_{\rm floor} \end{aligned} \]

Usamos a mesma decomposição bariônica que a calibração: \(M_{\rm disk} = \Upsilon\cdot 2\pi\,\Sigma_d\,R_d^2\) com \(\Upsilon = 0,5\), \(M_{\rm gas} = 1,33\,M_{\rm HI}\), \(R_{d,{\rm gas}} = 2,5\,R_{d,\star}\). Todos os três parâmetros universais são mantidos em seus valores estáveis. A nuvem é gerada pela amostragem de cada galáxia em nove raios de \(0,5\,R_d\) a \(10\,R_d\), produzindo 909 pontos nos 101 discos.

2. O diagrama

Nuvem BeeTheory e curva RAR empírica McGaugh no plano g_bar - g_obs
Figura 1 – Nuvem do modelo BeeTheory (pontos verdes) e RAR empírico de McGaugh (curva vermelha) no plano \(g_{\rm bar}\)-\(g_{\rm obs}\). Os pontos observados em \(R = 5\,R_d\) calculados a partir do \(V_f\) medido são sobrepostos para verificação de sanidade: calibração (círculos dourados, 20 galáxias) e cego (triângulos azuis, 81 galáxias).

Três recursos se destacam.

(a) Os pontos observados ficam próximos à curva de McGaugh. A dispersão dos valores SPARC \(V_f\) em torno do RAR empírico é de 0,17 dex com um pequeno viés de +0,10 dex – comparável aos 0,13 dex relatados por McGaugh+2016. Nossa amostra de 101 galáxias é saudável e consistente com a RAR publicada. Isso não era garantido: isso confirma que os pares \((g_{\rm bar}, g_{\rm obs})\) derivados do SPARC que usamos aqui traçam a mesma relação empírica.

(b) A nuvem BeeTheory acompanha globalmente a curva McGaugh. Ela não é plana (não entra em colapso em Newton) e se curva na direção certa com \(g_{\rm bar}\) baixo. A forma funcional é qualitativamente correta.

(c) A nuvem está sistematicamente acima da curva, com deslocamento mediano de +0,25 dex e dispersão total de 0,39 dex – três vezes a dispersão empírica. A BeeTheory nessa forma prevê excessivamente \(g_{\rm obs}\).

3. Onde a BeeTheory se distancia da RAR

BeeTheory residuals vs g_bar mostrando um aumento sistemático em g_bar alto
Figura 2 – Resíduos \(\log_{10}(g_{\rm obs}^{\rm BT} / g_{\rm obs}^{\rm McGaugh})\) em função de \(g_{\rm bar}\). A banda de dispersão de McGaugh (\(\pm 0,13\) dex) está sombreada em vermelho. Pontos azuis: mediana em compartimentos dos resíduos de BeeTheory; faixa azul: percentil 25-75.

A estrutura não é aleatória. No menor \(g_{\rm bar}\) (em torno de \(3\times 10^{-13}\) m/s², ou seja, as regiões mais externas dos discos LSB), os resíduos do BeeTheory se agrupam dentro de \(\pm 0,15\) dex de zero – o modelo cai na curva McGaugh dentro da dispersão empírica. Esse é o regime em que a calibração foi realizada (\(V_f\) é atingido em \(R \approx 5\,R_d\) onde \(g_{\rm bar}\) é baixo), portanto, a consistência aqui era esperada.

A tendência, então, aumenta constantemente com \(g_{\rm bar}\). Em torno de \(g_\dagger = 1,2\times 10^{-10}\) m/s² (a aceleração de transição no RAR empírico), a mediana binned fica em +0,45 a +0,50 dex. Em unidades lineares, \(g_{\rm obs}^{\rm BT}\) é aproximadamente 3 vezes o valor de McGaugh no \(g_{\rm bar}\) correspondente. Essas são as regiões internas das galáxias, onde \(R\) é pequeno em comparação com \(R_d\).

Por que ele falha internamente. O RAR empírico exige que em alto \(g_{\rm bar}\) (profundamente dentro do disco), \(g_{\rm obs} \to g_{\rm bar}\) – a matéria visível domina, e a contribuição do campo de ondas deve se tornar subdominante. Na parametrização atual da BeeTheory, \(\ell_{\rm wave} = c\,R_d + \ell_{\rm floor}\) é de aproximadamente 3 kpc para quase todas as galáxias (porque \(c = 0,16\) é pequeno). Em \(R \ll \ell_{\rm wave}\), \(x = R/\ell_{\rm wave}\) é pequeno e \(M_{\rm wave}(<\!R) \approx M\,x^3/6\) - a massa da onda cresce com \(R\), mas cresce a partir de um acoplamento \(\lambda = 12,7\). O produto \(\lambda \cdot M_{\rm wave}(<\!R)\) não desaparece com rapidez suficiente para recuperar Newton em \(R\) pequeno. O modelo "vaza" a contribuição da onda para dentro.

4. Síntese e próximas etapas

MétricaNuvem BeeTheoryObservado (\(V_f\) pontos)McGaugh+2016
Deslocamento mediano (dex)+0.25+0.100,00 (definição)
Dispersão \(\sigma\) (dex)0.390.17~0.13
Número de pontos909 (101 gal.)101~2700 (~150 gal.)
Forma qualitativamente correta?SimSim
Limite de Newton em alto \(g_{\rm bar}\)?Não (+0,5 dex)Sim (integrado)

Leitura. A BeeTheory passa o RAR assintoticamente em baixa aceleração – o regime que gera curvas de rotação planas, onde a calibração foi ancorada. Ele falha no regime de alta aceleração, em que o RAR converge para Newton e o BeeTheory não. Esse é um defeito limpo e localizado, não uma incompatibilidade total. A forma de 3 parâmetros nunca foi necessária para reproduzir o comportamento do disco interno porque \(V_f\) está definido em \(R \sim 5\,R_d\). O encontro com o RAR força o modelo a estar correto em todos os raios, não apenas em \(V_f\).

Caminhos para explorar

(i) Saturando \(\ell_{\rm wave}\) em \(R\) pequeno. Um refinamento natural é fazer com que a extensão efetiva do campo de ondas seja uma função das condições locais, por exemplo, \(\ell_{\rm wave}(R)\) crescendo com \(R\) em vez de ser fixada em \(c\,R_d + \ell_{\rm floor}\). Isso suprimiria o acoplamento de onda em \(R\) pequeno e deixaria Newton dominar onde \(g_{\rm bar}\) é alto.

(ii) Acoplamento dependente de \(g\)\. Se o próprio \(\lambda\) depender de \(g_{\rm bar}\) – por exemplo, \(\lambda \to \lambda \cdot f(g_{\rm bar}/g_\dagger)\) com \(f \to 0\) em \(g_{\rm bar}\) alto e \(f \to 1\) em baixo – o modelo poderia reproduzir exatamente a curva de McGaugh. O desafio é motivar essa dependência a partir da microfísica da função de onda, e não impô-la fenomenologicamente.

(iii) Limite em \(\ell_{\rm floor}\) para galáxias maciças. O mesmo defeito que produz o viés high-\(g_{\rm bar}\) aqui pode estar relacionado à previsão excessiva de +43% na NGC3198 na calibração e às previsões excessivas de +25-45% em várias galáxias cegas Sc/Sbc. Um \(\ell_{\rm floor}\) saturado (ou um que dependa de \(M_{\rm visible}\)) poderia corrigir ambos simultaneamente.

Nota sobre a metodologia. Todos os cálculos: 101 galáxias SPARC bulgeless (\(T \geq 4\)), 909 \((R, g_{\rm bar}, g_{\rm obs})\) tripletos em \(R/R_d \in \{0.5, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0, 7.0, 10.0\}\). Parâmetros estáveis \((\lambda, c, \ell_{\rm floor}) = (12,696,\, 0,163,\, 3,00\ {\rm kpc})\). RAR empírico de McGaugh com \(g_\dagger = 1,20\times 10^{-10}\) m/s². A comparação com as medições SPARC individuais publicadas de McGaugh+2016 (não reproduzidas aqui a partir das amostras originais de curva de rotação) é uma próxima etapa natural.


BeeTheory.com – Primeiro teste contra a Relação de Aceleração Radial – Geração inicial: 2026-05-20 com Claude.ai – © Technoplane S.A.S. 2026