Simulação numérica do modelo de massa escura BeeTheory

Um relato completo e reproduzível da integração numérica, decomposição da velocidade bariônica, minimização de χ² e escolhas de implementação por trás da simulação de massa escura BeeTheory.

Esta página técnica explica como reproduzir a simulação numérica BeeTheory da massa oculta da Via Láctea. Ela descreve os dados observacionais, o modelo bariônico, a equação de densidade baseada em ondas, a integração numérica, o método de ajuste, os testes de convergência e o código de referência.

O objetivo é simples: partir de um disco visível da Via Láctea, aplicar o modelo de densidade de onda BeeTheory, calcular a massa oculta efetiva e comparar a curva de velocidade circular resultante com as observações da era Gaia.

0. O que estamos computando e por quê

A curva de rotação da Via Láctea é a velocidade circular Vc(R) das estrelas em função de sua distância R do centro galáctico. Atualmente, ela é medida com muito mais precisão do que a distribuição total de massa que podemos ver diretamente.

O déficit entre a velocidade observada e o que a matéria bariônica visível prevê é o problema da massa oculta. Os modelos padrão invocam um halo de partículas invisíveis, geralmente a matéria escura fria. A BeeTheory propõe uma interpretação diferente: cada elemento de massa visível irradia um campo de onda que decai exponencialmente em três dimensões, e o campo acumulado se comporta como massa oculta.

A simulação faz três coisas:

1. Calcula ^Vcbar(R) a partir do disco visível e do bojo usando a fórmula analítica de Freeman para o disco.
2. Integra numericamente a densidade BeeTheory _ρdark(r; ℓ, λ) em cada raio e, em seguida, converte-a em ^VcDM(R).
3. Minimiza χ²(ℓ, λ) em relação à curva de rotação da Via Láctea da era Gaia para encontrar parâmetros representativos de melhor ajuste.

A simulação foi projetada para ser reproduzível. Ela pode ser executada em JavaScript ou Python sem uma biblioteca astrofísica especializada.

Notação usada em todo o texto

• R é o raio cilíndrico galactocêntrico no plano do disco, em kpc.
• r é o raio esférico galactocêntrico.
• z é a altura acima do disco.
• ℓ é o comprimento de coerência da BeeTheory, em kpc.
• λ é o acoplamento de massa de onda sem dimensão.
• G é usado em unidades de kpc km² s-² M⊙-¹.

1. Dados observacionais – Gaia 2024

O conjunto de dados é baseado na curva de rotação da Via Láctea da era Gaia. Ele usa raios R, velocidades circulares _Vc e incertezas σ. A versão técnica original usava 16 pontos de dados cobrindo R = 4-27,3 kpc.

Os fatos observacionais importantes são:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, o ponto de ancoragem próximo à órbita do Sol.
• O _Vc é aproximadamente plano de cerca de 5 a 20 kpc.
• _Vc diminui no disco externo medido, atingindo cerca de 173 ± 17 km/s perto de 27,3 kpc nos dados referenciados.

Isso requer uma distribuição efetiva de massa oculta que aumenta fortemente em raios intermediários e depois enfraquece em raios maiores.

Excluímos a galáxia mais interna porque a barra central e os movimentos não circulares predominam ali. Um modelo axissimétrico simplificado não é confiável dentro dessa região.

3. A equação de densidade escura da BeeTheory

3.1 Postulado físico

Cada elemento de massa visível dV na posição r′ no disco galáctico gera um campo de ondas gravitacionais que cria uma densidade de massa efetiva em um ponto de campo r:

A densidade escura total em qualquer ponto (R,z) é a superposição de todos os elementos do disco. Como a fonte é um disco, o termo de densidade de volume se torna um termo de densidade de superfície:

A integral dupla completa é:

3.3 Verificação analítica do limite assintótico

Para _Rd ≪ r ≪ ℓ, a expressão se simplifica. O disco é muito menor do que o raio, e o comprimento de coerência ainda é muito maior do que o raio.

Porque:

Essa densidade assintótica se comporta como ρ ∝ r-², o que dá M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}}}\approx\mathrm{constant}\)

Ao inserir valores representativos, obtém-se uma densidade local próxima ao valor observado da Via Láctea:

5. χ² Minimização e ajuste de parâmetros

5.1 Função objetiva

O ajuste minimiza o qui-quadrado reduzido:

com N = 16 pontos de dados e p = 2 parâmetros livres, ℓ e λ. Isso dá 14 graus de liberdade.

5.2 Pesquisa de grade de duas passagens

É usada uma pesquisa de grade em vez de descida de gradiente porque o cenário tem um vale de degeneração longo e curvo entre λ e ℓ.

• Passo 1: grade grossa sobre ℓ e λ.
• Passo 2: refinamento local em torno do mínimo grosseiro.

A região representativa de melhor ajuste é:

5.3 O Vale da Degeneração

A paisagem χ² não é uma tigela simétrica. Ele forma um vale alongado porque a normalização da densidade principal depende muito da força de acoplamento e apenas fracamente do comprimento de coerência dentro do regime de rotação plana.

O declínio externo da curva de rotação quebra essa degeneração porque valores menores de ℓ suprimem a densidade efetiva mais cedo.

Dados além de 30 kpc, incluindo aglomerados globulares, fluxos estelares, estrelas de halo e galáxias satélites, são essenciais para restringir ℓ com mais rigor.

7. Código de referência completo

A implementação de referência JavaScript a seguir reproduz a lógica principal da simulação. Ela se destina à validação técnica e deve ser colocada em um ambiente de script adequado, não diretamente dentro de um bloco de conteúdo padrão do WordPress, a menos que scripts personalizados sejam permitidos.

// Constantes físicas
const G = 4.302e-3; // kpc km² s-² M☉-¹
const Sig0 = 800,0; // M☉ pc-²
const Rd = 2,6; // kpc
const Mdisk = 3,5e10; // M☉
const Mbulge = 1,2e10; // M☉
const CONV_RHO = (1,989e30 * 5,61e26) / (3,086e21)**3;

// Dados de rotação da era Gaia
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];

// Espaço reservado para velocidade bariônica
function vBaryonic(R) {
  // Na implementação completa, isso usa a fórmula de disco de Freeman
  // com funções de Bessel modificadas I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  return Math.sqrt(Math.max(0, vb2));
}

// Densidade escura da BeeTheory
function rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30,0 / N;
  let sum = 0;

  for (let i = 0; i < N; i++) {
    const Rp = (i + 0,5) * dRp;
    const Sig = Sig0 * Math.exp(-Rp / Rd);
    const K = (2 * Math.PI * ell / r)
      * Math.sinh(Math.min(r / ell, 30))
      * Math.exp(-Math.min((r + Rp) / ell, 30));

    sum += Sig * Rp * K * dRp;
  }

  return (lam / ell) * sum;
}

// Massa escura fechada
function enclosedDarkMass(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  let M = 0;

  for (let j = 0; j < Nr; j++) {
    const rj = (j + 0,5) * dr;
    M += 4 * Math.PI * rj * rj * rhoDark(rj, ell, lam) * dr;
  }

  return M;
}

// Velocidade escura e total
function vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

function vTotal(R, ell, lam) {
  const vb = vBaryonic(R);
  const vd = vDM(R, ell, lam);
  return Math.sqrt(vb * vb + vd * vd);
}

// Qui-quadrado
function chiSquared(ell, lam) {
  let s = 0;

  for (let i = 0; i < OBS_R.length; i++) {
    const dv = (vTotal(OBS_R[i], ell, lam) - OBS_V[i]) / OBS_ERR[i];
    s += dv * dv;
  }

  return s / (OBS_R.length - 2);
}

Uma versão completa deve incluir implementações precisas das funções de Bessel modificadas _I0, _I1, _K0 e _K1 para a velocidade do disco de Freeman.

Referências

• Abramowitz, M., Stegun, I. A. - Handbook of Mathematical Functions, Dover, 1972.
• Bovy, J., Rix, H.-W. - A Direct Dynamical Measurement of the Milky Way's Disk Surface Density Profile (Medição dinâmica direta do perfil de densidade da superfície do disco da Via Láctea), ApJ 779, 115, 2013.
• Freeman, K. C. - On the disks of spiral and S0 galaxies (Sobre os discos de galáxias espirais e S0), ApJ 160, 811, 1970.
• 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.
• 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, 2024.
• Pato, M., Iocco, F., Bertone, G. - Dynamical constraints on the dark matter distribution in the Milky Way, JCAP 12, 001, 2015.
• Portail, M. et al. - Dynamical modelling of the Galactic bulge and bar (Modelagem dinâmica do bojo e da barra galáctica), MNRAS 465, 1621, 2017.