Simulation numérique du modèle de masse noire de la théorie de l’abeille

Un compte-rendu complet et reproductible de l’intégration numérique, de la décomposition de la vitesse baryonique, de la minimisation du χ² et des choix d’implémentation derrière la simulation de la masse sombre de BeeTheory.

Cette page technique explique comment reproduire la simulation numérique BeeTheory de la masse cachée de la Voie Lactée. Elle décrit les données d’observation, le modèle baryonique, l’équation de densité basée sur les ondes, l’intégration numérique, la méthode d’ajustement, les tests de convergence et le code de référence.

L’objectif est simple : partir d’un disque visible de la Voie lactée, appliquer le modèle de densité d’ondes de BeeTheory, calculer la masse cachée effective et comparer la courbe de vitesse circulaire résultante avec les observations de l’ère Gaia.

0. Ce que nous calculons et pourquoi

La courbe de rotation de la Voie lactée est la vitesse circulaire Vc(R) des étoiles en fonction de leur distance R par rapport au centre galactique. Elle est mesurée aujourd’hui avec beaucoup plus de précision que la distribution de la masse totale que nous pouvons voir directement.

Le déficit entre la vitesse observée et ce que la matière baryonique visible prédit est le problème de la masse cachée. Les modèles standard invoquent un halo de particules invisibles, généralement de la matière noire froide. La théorie des abeilles propose une interprétation différente : chaque élément de masse visible émet un champ d’ondes qui se désintègre exponentiellement en trois dimensions, et le champ accumulé se comporte comme une masse cachée.

La simulation a trois fonctions :

1. Calcule ^Vcbar(R) à partir du disque visible et du bulbe en utilisant la formule analytique de Freeman.
2. Intègre numériquement la densité BeeTheory _ρdark(r ; ℓ, λ) à chaque rayon, puis la convertit en ^VcDM(R).
3. Minimise χ²(ℓ, λ) par rapport à la courbe de rotation de la Voie lactée de l’ère Gaia pour trouver les paramètres représentatifs les mieux adaptés.

La simulation est conçue pour être reproductible. Elle peut être exécutée en JavaScript ou en Python sans bibliothèque spécialisée en astrophysique.

Notation utilisée dans l’ensemble du document

• R est le rayon galactocentrique cylindrique dans le plan du disque, en kpc.
• r est le rayon galactocentrique sphérique.
• z est la hauteur au-dessus du disque.
• ℓ est la longueur de cohérence de la théorie de Bee, en kpc.
• λ est le couplage onde-masse sans dimension.
• G est utilisé en unités de kpc km² s-² M⊙-¹.

1. Données d’observation – Gaia 2024

Le jeu de données est basé sur la courbe de rotation de la Voie Lactée de l’ère Gaia. Il utilise les rayons R, les vitesses circulaires _Vc et les incertitudes σ. La version technique originale utilisait 16 points de données couvrant R = 4-27.3 kpc.

Les faits d’observation importants sont les suivants :

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, le point d’ancrage près de l’orbite du Soleil.
• _Vc est approximativement plat entre 5 et 20 kpc.
• _Vc diminue dans le disque externe mesuré, atteignant environ 173 ± 17 km/s près de 27,3 kpc dans les données référencées.

Cela nécessite une distribution effective de la masse cachée qui augmente fortement à des rayons intermédiaires, puis s’affaiblit à des rayons plus importants.

Nous excluons la galaxie la plus intérieure parce que la barre centrale et les mouvements non circulaires y dominent. Un modèle axisymétrique simplifié n’est pas fiable dans cette région.

3. L’équation de densité obscure de la théorie de Bee

3.1 Postulat physique

Chaque élément de masse visible dV à la position r′ dans le disque galactique génère un champ d’ondes gravitationnelles qui crée une densité de masse effective en un point du champ r :

La densité d’obscurité totale en tout point (R,z) est la superposition de tous les éléments du disque. La source étant un disque, le terme de densité volumique devient un terme de densité surfacique :

La double intégrale complète est :

3.3 Vérification analytique de la limite asymptotique

Pour _Rd ≪ r ≪ ℓ, l’expression se simplifie. Le disque est beaucoup plus petit que le rayon, et la longueur de cohérence est encore beaucoup plus grande que le rayon.

Parce que :

Cette densité asymptotique se comporte comme ρ ∝ r-², ce qui donne M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}\approx\mathrm{constante}\)

En introduisant des valeurs représentatives, on obtient une densité locale proche de la valeur observée dans la Voie lactée:

5. Minimisation du χ² et ajustement des paramètres

5.1 Fonction objective

L’ajustement minimise le chi-carré réduit :

avec N = 16 points de données et p = 2 paramètres libres, ℓ et λ, ce qui donne 14 degrés de liberté.

5.2 Recherche de grille à deux passages

Une recherche sur grille est utilisée plutôt qu’une descente de gradient car le paysage présente une longue vallée de dégénérescence incurvée entre λ et ℓ.

• Passe 1 : grille grossière sur ℓ et λ.
• Passe 2 : raffinement local autour du minimum grossier.

La région représentative la mieux adaptée est la suivante :

5.3 La vallée de la dégénérescence

Le paysage χ² n’est pas un bol symétrique. Il forme une vallée allongée parce que la normalisation de la densité de tête dépend fortement de la force de couplage et faiblement de la longueur de cohérence dans le régime de rotation plate.

Le déclin extérieur de la courbe de rotation rompt cette dégénérescence car les valeurs plus petites de ℓ suppriment la densité effective plus tôt.

Les données au-delà de 30 kpc, y compris les amas globulaires, les courants stellaires, les étoiles du halo et les galaxies satellites, sont essentielles pour contraindre plus étroitement ℓ.

7. Code de référence complet

La mise en œuvre de référence JavaScript suivante reproduit la logique principale de la simulation. Elle est destinée à la validation technique et doit être placée dans un environnement de script approprié, et non directement dans un bloc de contenu WordPress standard, à moins que les scripts personnalisés ne soient autorisés.

// Constantes physiques
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 ;

// Données de rotation de l'ère 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,6,6,6,6,7,7,8,9,10,11,13,17] ;

// Paramètre de la vitesse baryonique
function vBaryonic(R) {
  // Dans l'implémentation complète, cette fonction utilise la formule du disque de Freeman
  // avec les fonctions de Bessel modifiées I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2) ;
  return Math.sqrt(Math.max(0, vb2)) ;
}

// Densité d'obscurité de la théorie des abeilles
function rhoDark(r, ell, lam) {
  const N = 40 ;
  const dRp = 30.0 / N ;
  laissez 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 ;
}

// Masse sombre enfermée
function enclosedDarkMass(R, ell, lam) {
  const Nr = 30 ;
  const dr = R / Nr ;
  laissez 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 ;
  }

  retour M ;
}

// Vitesse sombre et totale
function vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R)) ;
}

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

// Khi-deux
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) ;
}

Une version complète devrait inclure des implémentations précises des fonctions de Bessel modifiées _I0, _I1, _K0 et_K1 pour la vitesse du disque de Freeman.

Références

• 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, ApJ 779, 115, 2013.
• Freeman, K. C. - On the disks of spiral and S0 galaxies, ApJ 160, 811, 1970.
• McMillan, P. J. - The mass distribution and gravitational potential of the Milky Way, 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, MNRAS 465, 1621, 2017.