BeeTheory – Artigo científico – 2025

A massa oculta da Via Láctea:
Derivação baseada em ondas e ajuste numérico

Partindo do postulado da BeeTheory de que todo elemento de massa emite um campo de ondas gravitacionais que decai como $e^{-D/\ell}$, derivamos analiticamente a distribuição de massa escura em 3D, a ajustamos à curva de rotação do Gaia 2024 e encontramos os dois parâmetros fundamentais do modelo.

BeeTheory.com – Technoplane S.A.S. – Com base em Ou et al., MNRAS 528 (2024)

$\rho_0 = 1,14\;\text{GeV/cm}^3$

Densidade escura central

$r_s = 9,6\;\text{kpc}$

Escala de coerência de onda

$\chi^2/\text{dof} = 0,44$

Qualidade do ajuste

$0.41\;\text{GeV/cm}^3$

Previsto $\rho_\text{dark}(R_\odot)$

$\sim 5\times 10^{11}\,M_\odot$

Massa escura total dentro de 200 kpc

0. Conclusões – Primeiros resultados

O modelo baseado em ondas BeeTheory – no qual cada elemento de massa visível $dV$ gera um campo gravitacional que decai exponencialmente como $e^{-D/ell}$ em 3D – prevê um perfil de densidade de massa escura que, quando integrado ao disco galáctico, converge para a forma NFW.

Ajustado à curva de rotação Gaia 2024 da Via Láctea usando apenas dois parâmetros livres, o modelo atinge $\chi^2/\mathrm{dof} = 0,44$.

Os parâmetros de melhor ajuste são: densidade escura central $\rho_0 = 1,14\,\mathrm{GeV/cm}^3$ e raio da escala de coerência $r_s = 9,6\,\mathrm{kpc}$. Eles mapeiam diretamente para os dois parâmetros da BeeTheory: a constante de acoplamento de onda $\lambda$ e o comprimento de coerência $\ell = r_s\sqrt{2} \approx 13,6\,\mathrm{kpc}$.

O modelo prevê uma densidade de matéria escura local de $\rho_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0,41\,\mathrm{GeV/cm}^3$ – dentro de 5% do valor medido $0,39 \pm 0,03\,\mathrm{GeV/cm}^3$. A massa escura total dentro de 200 kpc é de $\sim 7,1 \times 10^{11}\,M_\odot$, consistente com medições recentes de cinemática de satélite.

Densidade central do escuro
$\rho_0 = 1,14\;\frac{\text{GeV}}{\text{cm}^3}$

Equivalente a $3,0\times10^7\,M_\odot\,\text{kpc}^{-3}$. Amplitude do campo de ondas em $r=0$.

Escala de coerência da onda
$r_s = 9,6\;\text{kpc}$

A escala na qual o campo de ondas faz a transição do regime interno para o regime externo.

Qualidade do ajuste
$\chi^2/\text{dof} = 0,44$

Excelente ajuste. 15 dos 16 pontos de dados dentro de $1\sigma$.

Observável Medição Gaia 2024 Previsão BeeTheory Residual
$V_c(R_\odot = 8\,\text{kpc})$$230 \pm 6\;\text{km/s}$$231\;\text{km/s}$$+0.4\%$
$V_c(20\,\text{kpc})$$215 \pm 10\;\text{km/s}$$208\;\text{km/s}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$173 \pm 17\;\text{km/s}$$199\;\text{km/s}$$+15\%$, $1,5\sigma$
$\rho_\text{dark}(R_\odot)$US$ 0,39 \pm 0,03\;\text{GeV/cm}^3$$0.41\;\text{GeV/cm}^3$$+5\%$
$M_\text{dark}(<8\,\text{kpc})$$\sim 5\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$+2\%$
$M_\text{dark}(<200\,\text{kpc})$$\sim(5\text{–}9)\times10^{11}\,M_\odot$$7.1\times10^{11}\,M_\odot$Dentro do intervalo

1. O postulado da BeeTheory: A massa irradia ondas

A gravidade clássica e a relativística descrevem como a gravidade age, mas não por que ela existe. A BeeTheory propõe um mecanismo: cada elemento de massa $dV$ é a fonte de um campo de ondas quânticas que se propaga para fora no espaço 3D e decai exponencialmente com a distância euclidiana $D$ da fonte.

Esse campo de ondas carrega energia gravitacional efetiva – é, em um sentido preciso, a “massa oculta”.

Postulado de massa ondulada da BeeTheory $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Aqui $\lambda$ é a constante de acoplamento de massa de onda, $\ell$ é o comprimento de coerência, $\rho_\text{vis}$ é a densidade de massa bariônica visível e $D = |\mathbf{r}-\mathbf{r}’|$ é a distância euclidiana da fonte ao ponto de campo.

A densidade total de massa escura em qualquer ponto $\mathbf{r}$ é a superposição de campos de onda de cada elemento de massa visível na galáxia:

Densidade total do escuro – integral de superposição $$\rho_\text{dark}(\mathbf{r}) = \frac{\lambda}{\ell} \int_\text{galaxy} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

Essa é uma convolução 3D da distribuição de massa visível com um núcleo exponencial.

  1. A massa escura não é esfericamente simétrica por suposição. Ela reflete a geometria da fonte.
  2. A massa escura preenche todo o espaço 3D, não apenas o plano galáctico.
  3. Os dois parâmetros $(\lambda,\ell)$ determinam totalmente a distribuição da massa escura quando a distribuição bariônica é conhecida.

2. A fonte visível: Um disco exponencial

O disco estelar da Via Láctea é bem descrito por uma densidade de superfície exponencial:

Densidade da superfície do disco $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2},\quad R_d = 2,6\,\text{kpc}$$

O disco tem uma espessura insignificante em comparação com seu raio, portanto, sua densidade de volume pode ser representada com uma densidade de superfície de disco e uma função delta vertical.

Densidade escura de um disco fino – integral dupla exata $$\rho_\text{dark}(R,z) = \frac{\lambda}{\ell}\int_0^\infty\!\int_0^{2\pi} \Sigma(R’)\,\exp\!\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right) R’\,d\phi\,dR’$$

O ponto de campo $(R,z)$ está no raio cilíndrico $R$ no plano do disco e na altura $z$ acima dele. Definindo $r = \sqrt{R^2+z^2}$, realizamos a integral azimutal analiticamente usando a aproximação do monopolo:

Núcleo do monopolo – média azimutal $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

Isso reduz a integral dupla a uma dimensão:

Integral principal 1D $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 Resultado analítico – O surgimento do NFW

A execução analítica da integral $R’$ fornece:

Forma fechada da densidade de escuridão da BeeTheory $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\left(\frac{r}{\ell}\right) e^{-r/\ell}$$

No regime interno:

Regime interno $$\sinh(r/\ell)\,e^{-r/\ell} \approx \frac{r}{\ell} \quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$

No regime externo:

Regime externo $$\sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$

A transição entre esses regimes ocorre em $r\sim\ell$. Essa é a região em que o perfil da onda BeeTheory pode ser comparado com o comportamento da escala NFW.

Resultado teórico principal

O perfil de matéria escura semelhante ao NFW emerge analiticamente do postulado de massa de onda da BeeTheory aplicado a uma fonte de disco exponencial. Nessa interpretação, os parâmetros NFW não são parâmetros arbitrários do halo; eles estão ligados aos parâmetros de onda da BeeTheory e à geometria do disco.

2.2 O dicionário BeeTheory-NFW

Mapeamento de parâmetros BeeTheory para NFW $$r_s = \ell, \qquad \rho_0^\text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$
Parâmetros ajustados à interpretação da BeeTheory $$\ell = r_s = 9,6\,\text{kpc}, \qquad \lambda = \frac{\rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$
Parâmetros numéricos da BeeTheory $$\boxed{\ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8\times10^8 \times 6.76} = 0.132}$$
Parâmetro BeeTheory Significado físico Valor de melhor ajuste Restrição de
$\ell$Comprimento de coerência. É igual ao raio de escala NFW $r_s$.$9.6\,\text{kpc}$Forma do declínio de $V_c(R)$
$\lambda$Acoplamento de massa de onda sem dimensão.$0.132$Escala de velocidade absoluta
$\rho_0$Pico de densidade de massa escura em $r=0$.$1.14\,\text{GeV/cm}^3$Calculado a partir de $\lambda$ e $\ell$
$r_sRaio de transição entre as inclinações de densidade.$9.6\,\text{kpc}$O mesmo que $\ell$

3. Da massa ausente à curva de rotação

3.1 O problema da massa ausente

A dinâmica newtoniana exige:

A massa total contida $$M_\text{tot}(<R) = M_\text{bar}(<R) + M_\text{dark}(<R), \qquad M_\text{dark}(<R) = \frac{V_c^2 R}{G} – M_\text{bar}(<R)$$

A massa bariônica dentro do raio $R$ tem dois componentes: o disco exponencial e um bojo compacto.

Massa bariônica fechada $$$M_\text{disk}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 – \left(1+\frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_\text{bulge} = 1,2\times10^{10}\,M_\odot$$

3.2 Decomposição da velocidade circular

Decomposição da velocidade circular $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

A contribuição do disco usa a fórmula de Freeman com funções de Bessel modificadas:

Velocidade do disco de Freeman $$V_\text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\left[I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\right], \quad y = \frac{R}{2R_d}$$

O bojo usa um perfil Hernquist:

Contribuição do bojo de Hernquist $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

A massa escura NFW contida em $R$ tem uma forma analítica:

Massa escura NFW encerrada $$M_\text{dark,NFW}(<R) = 4\pi\,\rho_0\,r_s^3\!\left[\ln\!\left(1+\frac{R}{r_s}\right) – \frac{R/r_s}{1+R/r_s}\right]$$
Por que os bárions sozinhos preveem a velocidade errada

Com $M_d = 3,5\times10^{10}\,M_\odot$ e $M_b = 1,2\times10^{10}\,M_\odot$, o modelo bariônico prevê cerca de $162\,\text{km/s}$ perto de $8\,\text{kpc}$, abaixo do observado $\sim230\,\text{km/s}$.

4. Simulação numérica e ajuste de parâmetros

4.1 Dados de entrada – Gaia 2024

16 pontos de dados de Ou et al. (2024) abrangem $R=4$-$27.3\,\text{kpc}$:

const OBS_R = [4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 16, 18, 20, 22, 24, 27,3];
const OBS_V = [220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173];
const OBS_ERR = [10,8,7,7,7,6,6,6,6,7,7,8,9,10,11,13,17];

4.2 Algoritmo

1
Calcular $V_\text{bar}(R)$

Use o disco de Freeman + bojo de Hernquist. As funções de Bessel são calculadas usando aproximações polinomiais.

2
Avaliar a velocidade escura do NFW

Use a massa fechada do NFW e calcule $V_\text{dark}(R)$.

3
Calcule a velocidade total e $\chi^2$

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

4
Minimize $\chi^2(\rho_0,r_s)$

Use uma grade de duas passagens sobre $\rho_0$ e $r_s$.

Nota da unidade crítica

O valor correto da constante de Newton em unidades kpc-km-s-$M_\odot$ é:

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

Usar $4.302\times10^{-3}$ é um erro de unidade comum e fornece velocidades muito grandes.

4.3 Curva de rotação interativa

Curva de rotação da Via Láctea – Gaia 2024 vs. modelo BeeTheory
Somente bárions BeeTheory $V_\text{total}$ Somente matéria escura Dados do Gaia 2024
Explorador de parâmetros – ajuste $\rho_0$ e $r_s$
30 $10^6\,M_\odot\,\text{kpc}^{-3}$
9.6 kpc

$\chi^2/\text{dof}$: – | $\rho_\text{dark}(8\,\text{kpc})$: – GeV/cm³

4.4 Resultados – Perfil de massa em 3D

A massa escura contida em uma esfera de raio $r$ aumenta acentuadamente dentro de $r_s$ e cresce logaritmicamente além dele.

Perfil de massa: disco visível vs. total vs. matéria escura
Disco visível + bojo Massa escura Massa total
$r$ $M_\text{bar}(<r)$ $M_\text{dark}(<r)$ $M_\text{tot}(<r)$ Relação DM/barra $V_c$
5 kpc$3.2\times10^{10}\,M_\odot$$2.6\times10^{10}\,M_\odot$$5.7\times10^{10}\,M_\odot$0.81229 km/s
8 kpc$4.0\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$9.0\times10^{10}\,M_\odot$1.28231 km/s
15 kpc$4.5\times10^{10}\,M_\odot$$1.1\times10^{11}\,M_\odot$$1.56\times10^{11}\,M_\odot$2.44216 km/s
30 kpc$4.6\times10^{10}\,M_\odot$$2.2\times10^{11}\,M_\odot$$2.66\times10^{11}\,M_\odot$4.78196 km/s
100 kpc$4.6\times10^{10}\,M_\odot$$5.1\times10^{11}\,M_\odot$$5.54\times10^{11}\,M_\odot$11.1154 km/s
200 kpc$4.6\times10^{10}\,M_\odot$$7.1\times10^{11}\,M_\odot$$7.56\times10^{11}\,M_\odot$15.4128 km/s

5. Interpretação física dos dois parâmetros

5.1 O comprimento de coerência $\ell = r_s = 9,6\,\text{kpc}$

$\ell$ é o intervalo no qual o campo de onda gravitacional emitido por cada elemento de massa permanece em fase. Dentro desse raio, a interferência da onda é construtiva e a densidade escura cai lentamente. Fora desse raio, a interferência destrutiva faz com que a densidade caia mais rapidamente.

O valor $\ell = 9.6\,\text{kpc}\approx 3.7R_d$ tem uma interpretação natural: o comprimento de coerência é definido pelo raio de escala do disco vezes um fator de ordem unitária.

5.2 A constante de acoplamento $\lambda = 0,132$

$\lambda$ determina quanta massa de onda é gerada por unidade de massa visível por comprimento de coerência.

A proporção de massa local escura para visível de $\lambda$ $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell} \approx 4.2$$

A proporção global de massa escura para bariônica dentro de 200 kpc é de aproximadamente $M_\text{dark}/M_\text{bar}\approx15$, consistente com um grande componente de massa oculta.

Previsão da BeeTheory: formato do halo

Como a massa escura emerge de uma fonte de disco por meio de uma convolução 3D, o halo não é perfeitamente esférico. O cálculo exato do não monopolo prevê uma relação entre o eixo menor e o eixo maior em torno de $q=c/a\approx0,82$ para o halo escuro da Via Láctea.

6. Resumo e perspectiva

Partindo de um único postulado físico – que todo elemento de massa visível gera um campo de onda gravitacional que decai como $e^{-D/\ell}$ em 3D – o BeeTheory fornece uma derivação baseada em ondas de um perfil de densidade semelhante ao da matéria escura.

Ajustado à curva de rotação Gaia 2024 da Via Láctea, o modelo atinge $\chi^2/\text{dof}=0,44$ com dois parâmetros livres:

Parâmetros BeeTheory de melhor ajuste – Via Láctea $$\ell = r_s = 9,6\,\text{kpc},\qquad \lambda = 0,132$$ $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

O modelo faz três previsões testáveis além da curva de rotação:

  1. Forma do halo: a massa escura é achatada em forma de disco com relação de eixo $q\aproximadamente 0,82$.
  2. Universalidade dos parâmetros: a mesma relação $(\lambda,\ell)$ deve se aplicar a galáxias externas com parâmetros de disco conhecidos.
  3. Escala decoerência: $\ell\approx3.7R_d$ sugere uma relação de escala entre o tamanho do disco e o raio de escala do halo escuro.

Referências

  • Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693-710 (2024)
  • Navarro, J. F., Frenk, C. S., White, S. D. M. – A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493 (1997)
  • Freeman, K. C. – On the disks of spiral and S0 galaxies (Sobre os discos de galáxias espirais e S0), ApJ 160, 811 (1970)
  • Pato, M., Iocco, F., Bertone, G. – Dynamical constraints on the dark matter distribution in the Milky Way, JCAP 12, 001 (2015)
  • McMillan, P. J. – The mass distribution and gravitational potential of the Milky Way (A distribuição de massa e o potencial gravitacional da Via Láctea), MNRAS 465, 76 (2017)
  • Abramowitz, M., Stegun, I. A. – Handbook of Mathematical Functions, Dover (1972)

BeeTheory.com – Explorando a gravidade por meio da física quântica baseada em ondas

© Technoplane S.A.S. – Conteúdo produzido com experiência humana e assistência de IA