BeeTheory – Foundations – Note technique XIV
Étape 1 – Voie lactée :
Application du noyau de Yukawa de BeeTheory
La méthodologie de la Note XII est appliquée à la Voie Lactée en utilisant le noyau explicite de la forme d’onde de Yukawa $\mathcal{K}(D) = K_0\,(1+\alpha D)\,e^{-\alpha D}/D^2$, l’intégration étant effectuée séparément pour chaque composant baryonique en fonction de sa géométrie. Le résultat est comparé point par point à la courbe de rotation de Gaia 2024 et à la « masse manquante » du modèle standard. Cette note établit la base du cadre tel qu’il est appliqué avec la nouvelle formulation entièrement géométrique.
1. Le résultat d’abord
Résultat de base – Voie lactée avec noyau Yukawa explicite
Courbe de rotation. Le modèle reproduit la vitesse de Gaia 2024 à 2 km/s près à $R = 4$ kpc, mais sur-prédit de $+33$ km/s au rayon solaire ($R = 8$ kpc) et de $+64$ km/s à $R = 27.3$ kpc. $\chi^2/\text{dof} = 1.27$.
Masse manquante. La masse du champ d’onde deBeeTheory correspond à la « masse manquante » du modèle standard à 5% près à $R = 4$ kpc, mais la dépasse d’un facteur 2.2 à $R = 27.3$ kpc. Le modèle produit trop de masse de champ sombre aux grands rayons.
Densité locale. $\rho_\text{wave}(R_\odot) = 0.72$ GeV/cm³, comparé à l’intervalle observé $0.39$-$0.45$ GeV/cm³. Sur-prédiction d’un facteur d’environ 1,7.
Ce que cela signifie
La formulation explicite de Yukawa produit une courbe de rotation trop plate aux grands rayons. La longueur de décroissance du champ d’ondes $\ell$ est trop longue, ce qui permet au champ d’ondes de continuer à contribuer à la masse au-delà du disque visible. Il s’agit de la ligne de base structurelle avant l’incorporation du raffinement de la densité de surface identifié dans la Note XI.
2. Ce que nous voulons calculer
La Voie Lactée est le cas de test naturel car c’est la galaxie sur laquelle le couplage global $lambda$ a été calibré à l’origine, et parce qu’il existe deux observations indépendantes auxquelles se comparer :
(a) La courbe de rotation $V_c(R)$ de Gaia 2024 (Ou et al., MNRAS 528), qui mesure la vitesse circulaire à dix rayons de $R = 2$ kpc à $R = 27.3$ kpc avec des incertitudes statistiques de $7$-$17$ km/s. C’est cette vitesse que BeeTheory doit reproduire, en combinant la contribution baryonique $V_\text{bar}(R)$ avec la contribution du champ d’ondes $V_\text{wave}(R)$.
(b) La « masse manquante » $M_text{missing}(. L’interprétation standard invoque la matière noire particulaire pour fournir cette masse. La théorie des abeilles prédit au contraire que cette masse est le champ d’ondes $M_\text{wave}(
(c) La densité locale de matière noire au Soleil, mesurée à $\rho \approx 0.39$-$0.45$ GeV/cm³ à partir d’études cinématiques du voisinage solaire. BeeTheory prédit une valeur pour $\rho_\text{wave}(R_\odot)$ à partir du même calcul du champ d’onde.
L’accord (ou le désaccord) sur ces trois observables permet de tester trois aspects différents du modèle : la forme de la courbe de rotation, le profil de la masse enfermée et la normalisation de la densité locale.
3. Le noyau d’onde – forme explicite
Chaque élément de masse baryonique génère un champ d’ondes BeeTheory, dont l’intensité en un point du champ séparé par la distance $D$ est donnée par :
Noyau d’onde de la forme de Yukawa
$$\mathcal{K}(D) \;=\ ; K_0 \cdot \frac{(1 + \alpha D)\,e^{-\alpha D}{D^2}, \qquad \alpha = \frac{1}{\ell}$$$.
L’amortissement exponentiel $e^{-\alpha D}$ garantit que le champ d’ondes a une masse totale finie – sans lui, l’intégrale du champ d’ondes divergerait à l’infini. Le préfacteur $(1 + alpha D)$ provient de la fonction d’onde régularisée de la Note I ; avec l’exponentielle, il rend la structure spatiale du noyau quasi-newtonienne à $D ll ell$ et exponentiellement supprimée à $D gg ell$.
La longueur caractéristique $\ell$ dépend de la composante générant le champ :
| Composant | Longueur de cohérence $\ell$ (kpc) | Échelle géométrique |
|---|---|---|
| Bulge (3D Hernquist) | $\_ell_b = c_\text{sph}\,r_b = 0.41 \Nfois 0.61$ | $0.25$ |
| Disque mince | $\ell_\text{thin} = c_\text{disk}\,R_d = 3.17 \times 2.6$ | $8.24$ |
| Disque épais | $\N-\N-\N-\N-{épaisseur} = c_\N-\N-{disque}\N,(1.5\N,R_d)$ | $12.36$ |
| Anneau de gaz | $\N-\N-\N{gaz} = c_\N-\N{disque}\N,(1.7\N,R_d)$\N-\N-\N-\N-\N-\N-\N-\N | $14.01$ |
| Bras en spirale | $\ell_\text{arm} = c_\text{arm}\,R_d = 2.0 \\Nfois 2.6$ | $5.20$ |
4. Géométrie d’intégration, composant par composant
Pour chaque composante, nous intégrons la densité de la source convoluée avec le noyau, en utilisant l’élément de volume approprié pour la géométrie. Le résultat est la densité du champ d’ondes $\rho_\text{wave}^{(i)}(r)$ au point de champ.
4.1 Intégration bulbe – coque sphérique
Le bulbe est une distribution sphérique tridimensionnelle. Chaque coquille fine de rayon $r’$ contient la masse $\rho_b(r’)\,4\pi r’^2\,dr’$ et contribue au champ à une distance radiale $r$ du centre. Dans l’approximation monopolaire, la séparation effective entre le point du champ et un point générique de la coquille est $D = \sqrt{r^2 + r’^2}$ :
$$\rho_\text{wave}^{(b)}(r) \;=\ ; \int_0^{6r_b} \rho_b(r’) \cdot K_0\,\frac{(1 + \alpha_b D)\,e^{-\alpha_b D}}{D^2} \cdot 4\pi r’^2 \, dr’, \quad D = \sqrt{r^2 + r’^2}$$.
Avec $\rho_b(r’) = M_b r_b / [2\pi r'(r’+r_b)^3]$ (Hernquist) et $\alpha_b = 1/0,25 = 4,0$ kpc$^{-1}$. L’intégration est interrompue à $6,r_b$, au-delà de laquelle la densité est numériquement négligeable.
4.2 Disques minces et épais – intégration des anneaux concentriques
Chaque disque est une fine distribution axisymétrique. Le disque est décomposé en anneaux concentriques de rayon galactocentrique $R’$ et de largeur $dR’$, chacun portant la masse $\Sigma(R’)\,2\pi R’\,dR’$. La contribution au point de champ au rayon $r$ (dans le plan du disque) nécessite la séparation effective $D = \sqrt{r^2 + R’^2}$ dans la même approximation monopolaire :
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\ ; \int_0^{8R_d} \Sigma_\text{thin}(R’) \cdot K_0\,\frac{(1 + \alpha_\text{thin} D)\,e^{-\alpha_\text{thin} D}}{D^2} \cdot 2\pi R’ \c, dR’, \cquad D = \sqrt{r^2 + R’^2}$$.
avec $\Sigma_\text{thin}(R’) = (M_\text{thin}/2\pi R_d^2)\,e^{-R’/R_d}$ et $\alpha_\text{thin} = 1/8.24 = 0.121$ kpc$^{-1}$. Le disque épais est identique avec sa propre échelle : $\Sigma_\text{thick}$, $R_\text{thick} = 3.9$ kpc, $\alpha_\text{thick} = 0.081$ kpc$^{-1}$.
4.3 Anneau de gaz – intégration de l’anneau avec appauvrissement central
La distribution de gaz présente un trou central (HI négligeable à $R \lesssim 2$ kpc) et s’étend plus loin que le disque stellaire. Le profil $\Sigma_\text{gas}(R’) = \Sigma_0\,\exp(-R_\text{hole}/R’ – R’/R_g)$ rend compte de ces deux caractéristiques :
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{8R_g} \Sigma_\text{gas}(R’) \cdot K_0\,\frac{(1 + \alpha_\text{gas} D)\,e^{-\alpha_\text{gas} D}}{D^2} \cdot 2\pi R’ \, dR’$$
avec $R_g = 4.42$ kpc, $R_\text{hole} = 2.21$ kpc, $\alpha_\text{gas} = 0.071$ kpc$^{-1}$. La longueur de cohérence plus importante reflète la distribution plus étendue du gaz.
4.4 Bras en spirale – intégration de l’anneau avec une amplitude réduite et un noyau plus étroit
Les bras spiraux portent $10\%$ de la densité de surface du disque fin et ont leur propre longueur de cohérence $\ell_\text{arm} = 5.2$ kpc, plus étroite que celle du disque pour refléter la concentration azimutale des bras :
$$\rho_\text{wave}^{(\text{arm})}(r) \;=\ ; \int_0^{8R_d} 0.10\\N-\text{thin}(R’) \cdot K_0\N-\cfrac{(1 + \alpha_\text{arm} D)\N,e^{-\alpha_\text{arm} D}}{D^2} \cdot 2\pi R’ \, dR’$$
4.5 Densité totale du champ d’ondes
$$\rho_\text{wave}(r) \;=\ ; \lambda \,\sum_{i} \rho_\text{wave}^{(i)}(r), \quad \lambda = 0.189$$
5. De la densité d’ondes à la courbe de rotation
Une fois que $\rho_\text{wave}(r)$ est connu à chaque rayon, la masse totale du champ d’ondes enfermé est obtenue par intégration radiale :
$$M_\text{wave}(R) \;=\ ; \int_0^{R} 4\pi r^2 \rho_\text{wave}(r) \, dr$$$
La vitesse circulaire prédite combine alors les contributions baryoniques et du champ d’ondes en quadrature, selon la relation newtonienne :
$$V_c^2(R) \;=\ ; V_\text{bar}^2(R) \;+\ ; \frac{G\\NM_text{wave}(R)}{R}$$$.
6. Courbe de rotation – résultats point par point
| $R$ (kpc) | $V_\text{bar}$ (km/s) | $M_\text{wave}/10^{10}$ | $V_\text{wave}$ (km/s) | V_\text{tot}$ (km/s) | $V_\text{obs}$ (Gaia 2024) | $\Delta$ |
|---|---|---|---|---|---|---|
| 2.0 | 157.8 | 0.67 | 120.0 | 198.3 | 250 ± 12 | -51.7 |
| 4.0 | 164.1 | 2.54 | 165.4 | 233.0 | 235 ± 10 | -2.0 |
| 6.0 | 163.5 | 5.17 | 192.5 | 252.6 | 230 ± 8 | +22.6 |
| 8.0 | 157.3 | 8.13 | 209.1 | 261.7 | 229 ± 7 | +32.7 |
| 10.0 | 148.7 | 11.18 | 219.3 | 265.0 | 224 ± 8 | +41.0 |
| 12.0 | 139.6 | 14.15 | 225.2 | 265.0 | 217 ± 9 | +48.0 |
| 15.0 | 126.7 | 18.29 | 229.0 | 261.8 | 208 ± 10 | +53.8 |
| 20.0 | 109.3 | 24.10 | 227.7 | 252.6 | 195 ± 12 | +57.6 |
| 25.0 | 96.7 | 28.54 | 221.6 | 241.8 | 180 ± 15 | +61.8 |
| 27.3 | 92.0 | 30.18 | 218.1 | 236.7 | 173 ± 17 | +63.7 |
Le modèle correspond parfaitement à l’observation à $R = 4$ kpc (Δ = -2 km/s) mais la sur-prédiction augmente avec le rayon. Au rayon solaire, la sur-prédiction est de +33 km/s (4,7σ au-dessus de l’incertitude Gaia). À la limite extérieure $R = 27,3$ kpc, la sur-prédiction atteint +64 km/s (3,8σ). La courbe prédite est trop plate – le champ d’ondes continue à contribuer à la masse au-delà du disque visible, parce que la coupure exponentielle à $D \sim \ell$ le permet.
7. La masse d’onde par rapport à la « masse manquante » du modèle standard
Pour chaque rayon, nous comparons trois quantités : la masse baryonique incluse (matière visible uniquement), la masse dynamique requise par la vitesse observée (loi de Newton appliquée à $V_\text{obs}$), et la masse du champ d’ondes de la théorie de Bee. La différence entre la seconde et la première est ce que le modèle standard appelle la « masse manquante »:
$$M_\text{missing}(<R) \;=\ ; \frac{R\,V_\text{obs}^2(R)}{G} \;-\ ; M_\text{bar}(<R)$$
Le rapport $M_text{wave}/M_text{missing}$ nous indique dans quelle mesure le champ ondulatoire de la théorie de l’abeille remplace la matière noire particulaire, rayon par rayon :
| $R$ (kpc) | $M_\text{bar}(| $M_\text{dyn}( | $M_\text{missing}$ | $M_\text{wave}$ (BT) | Rapport $M_\text{wave}/M_\text{miss}$ (BT) | |
|---|---|---|---|---|---|
| 2.0 | 1.16e+10 | 2.91e+10 | 1.75e+10 | 6.69e+09 | 0.38 |
| 4.0 | 2.51e+10 | 5.13e+10 | 2.63e+10 | 2.54e+10 | 0.97 |
| 6.0 | 3.73e+10 | 7.38e+10 | 3.65e+10 | 5.17e+10 | 1.42 |
| 8.0 | 4.60e+10 | 9.75e+10 | 5.15e+10 | 8.13e+10 | 1.58 |
| 10.0 | 5.14e+10 | 1.17e+11 | 6.52e+10 | 1.12e+11 | 1.71 |
| 12.0 | 5.44e+10 | 1.31e+11 | 7.70e+10 | 1.41e+11 | 1.84 |
| 15.0 | 5.60e+10 | 1.51e+11 | 9.49e+10 | 1.83e+11 | 1.93 |
| 20.0 | 5.56e+10 | 1.77e+11 | 1.21e+11 | 2.41e+11 | 1.99 |
| 25.0 | 5.43e+10 | 1.88e+11 | 1.34e+11 | 2.85e+11 | 2.13 |
| 27.3 | 5.37e+10 | 1.90e+11 | 1.36e+11 | 3.02e+11 | 2.22 |
Lecture quantitative
A $R = 4$ kpc, le champ d’ondes correspond essentiellement à la masse manquante (ratio $0.97$). Entre $R = 6$ et $R = 8$ kpc, le modèle dépasse déjà la masse manquante de 40-60%. Au-delà de $R = 15$ kpc, la masse du champ d’ondes est à peu près le double de ce que le modèle standard invoque comme matière noire. Le modèle produit une masse supplémentaire aux grands rayons – exactement le symptôme d’une longueur de cohérence $\ell$ qui est trop longue pour le disque stellaire visible.
8. Contributions des composantes au champ d’ondes au rayon solaire
L’évaluation de la contribution de chaque composant à $\rho_\text{wave}(R_\odot = 8\,\text{kpc})$ montre quelle source baryonique domine le champ d’ondes à cet endroit :
| Composant | $\rho_\text{wave}(R_\odot)$ ($M_\odot$/kpc³) | Fraction du total |
|---|---|---|
| Disque stellaire mince | 6,05 $ \Npar 10^7$ | 60.6% |
| Disque stellaire épais | 1,91 $ x 10^7$ | 19.1% |
| Anneau de gaz | 1,62 $ x 10^7$ | 16.2% |
| Bras en spirale | 4,15 \Npar 10^6$ | 4.1% |
| Le bourrelet | 1,55 $ \N- fois 10^{-5}$ | $\sim$0% |
Le disque stellaire mince domine le champ d’ondes à la position du Soleil (60%), le disque épais et l’anneau de gaz contribuant à peu près également (16-19%). La contribution du bulbe est négligeable car $\ell_b = 0.25$ kpc est beaucoup plus petit que $R_\odot = 8$ kpc – la suppression exponentielle dans le noyau tue la contribution du bulbe à cette distance.
La conversion de la densité totale en unités de physique des particules donne $\rho_\text{wave}(R_\odot) = 0,717$ GeV/cm³, à comparer avec la mesure cinématique de $0,39$-$0,45$ GeV/cm³ de Read 2014 et des analyses ultérieures. La prédiction dépasse la densité locale observée d’un facteur de $1.6$-$1.8$ – ce qui est cohérent avec la sur-prédiction de la courbe de rotation au même rayon.
9. Ce que cette ligne de base établit
Le mécanisme fonctionne en principe
Autour de $R = 4$ kpc – le corps central du disque – le champ d’ondes intégré est égal à la masse manquante du modèle standard à 5% près, et la courbe de rotation est reproduite à 2 km/s près. Le noyau ondulatoire, appliqué aux baryons visibles, produit une masse gravitationnelle quantitativement comparable à la matière noire particulaire à ce rayon. Aucune nouvelle particule n’est nécessaire ; le champ d’ondes de la matière visible rend compte de la gravité manquante.
Mais la courbe est trop plate aux grands rayons
Au-delà du disque central, le modèle surestime la vitesse de rotation d’une quantité qui croît de façon monotone avec le rayon. Le champ d’ondes continue d’accumuler de la masse au-delà du disque visible, car la longueur de cohérence du noyau $\N_\text{thin} = 8,24$ kpc est comparable à la taille du disque lui-même, ce qui permet des contributions significatives à $D = 15$-$25$ kpc. La courbe de rotation de Gaia, en revanche, décline légèrement au-delà de $R \sim 10$ kpc – une caractéristique que la formulation actuelle ne reproduit pas.
Une base, pas une réponse définitive
Ce calcul établit la base du modèle avec $\ell_i$ dépendant linéairement de $R_d$ seul. Le diagnostic de la Note XI a identifié que $\Sigma_d$ – la densité de la surface centrale – doit entrer dans la détermination de $\ell_i$ pour corriger la courbe aux grands rayons. Plus le disque est dense, plus la réponse des ondes doit être localisée. L’intégration de ce raffinement fait l’objet de notes ultérieures. La ligne de base de la Voie Lactée rapportée ici est ce que ces notes doivent améliorer.
10. Résumé
1. La courbe de rotation de la Voie Lactée est calculée en intégrant chaque composante baryonique par rapport au noyau d’onde de Yukawa $mathcal{K}(D) = K_0(1 + alpha D),e^{-alpha D}/D^2$, avec une géométrie appropriée : coquilles sphériques pour le bulbe, anneaux concentriques pour les disques, le gaz et les bras spiraux.
2. A $R = 4$ kpc, la masse du champ d’ondes de BeeTheory est en accord avec la « masse manquante » du modèle standard à 5% près (ratio 0.97), et la vitesse prédite correspond à Gaia 2024 à 2 km/s près.
3. Au rayon solaire ($R = 8$ kpc), le modèle surestime la vitesse de rotation de $+33$ km/s et la densité locale de matière noire d’un facteur de 1,6 – les deux étant cohérents entre eux.
4. Au-delà de $R = 15$ kpc, la masse prédite du champ d’ondes dépasse la masse manquante du modèle standard d’un facteur 2 ou plus. La courbe de rotation prédite ne diminue pas comme l’ exigent les données Gaia.
5. Le disque stellaire mince domine le champ d’ondes à la position du Soleil (60% de $\rho_\text{wave}$). La contribution du bulbe est négligeable. La décomposition est cohérente avec la géométrie d’intégration décrite.
6. La sur-prédiction à grand $R$ est la signature structurelle de $\ell_i$ qui est trop long. La note XI identifie que $\Sigma_d$ doit entrer dans la formule de la longueur de cohérence. Le raffinement de $\ell_i$ via $\Sigma_d$ est la prochaine étape.
Références. Ou, X. et al – The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693 (2024). Courbe de rotation Gaia 2024. – Bland-Hawthorn, J., Gerhard, O. – The Galaxy in Context, ARA&A 54, 529 (2016). Décomposition structurelle de MW. – Read, J. I. – The Local Dark Matter Density, J. Phys. G 41, 063101 (2014). Local DM density measurements (Mesures de la densité locale de la matière noire). – Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970). – Hernquist, L. – An analytical model for spherical galaxies and bulges, ApJ 356, 359 (1990). – Dutertre, X. – Bee Theory™ : Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Gravité quantique basée sur les ondes – Application de l’étape 1 – © Technoplane S.A.S. 2026