BeeTheory – Fondements – Note technique XII
Formalisation :
Le calcul de la théorie des abeilles à l’échelle galactique
Cette note formalise le cadre de la Théorie de l’abeille tel qu’il est appliqué à une galaxie à disque. Elle précise les données d’observation, la décomposition géométrique de la distribution baryonique, les équations intégrales définissant le champ d’ondes pour chaque composante, et la chaîne d’opérations qui permet d’obtenir la courbe de rotation prédite. La procédure est strictement unidirectionnelle : la structure baryonique observée détermine le champ d’ondes, qui détermine la courbe de rotation – jamais l’inverse.
1. Le calcul dans un diagramme
Une chaîne unidirectionnelle
Photométrie observée $\;\longrightarrow\;$ Décomposition baryonique $(\rho_\text{bar})$
Flèche descendante $\big\downarrow$
Convolution du champ d’onde $\;\longrightarrow;$ Densité d’onde $\rho_\text{wave})$
$\big\downarrow$ (flèche descendante)
Intégration de la masse $\;\longrightarrow;$ Masse d’onde enfermée $(M_\text{wave})$
$\big\downarrow$ (flèche de bas de page)
Relation newtonienne $\;\longrightarrow;$ Courbe de rotation prédite $(V_c)$
Aucune étape n’est inversée. La courbe de rotation $V_c(R)$ n’est jamais utilisée comme entrée.
2. Données d’observation
Pour chaque galaxie, le calcul nécessite cinq observables publiés. Ce sont les seules quantités spécifiques aux galaxies ; tout le reste est calculé à partir d’elles. Aucun ajustement par rapport à la courbe de rotation n’est effectué à ce stade.
| Symbole | Quantité | Source |
|---|---|---|
| $T$ | Type morphologique de Hubble | Catalogue (de Vaucouleurs et al. 1991, SPARC) |
| R_d$ | Longueur d’échelle du disque stellaire (kpc) | Photométrie Spitzer 3.6 µm (SPARC) |
| $\Sigma_d$ | Luminosité de surface du disque central ($L_\odot/\text{pc}^2$) | Photométrie Spitzer 3.6 µm (SPARC) |
| $M_\text{HI}$ | Masse totale d’hydrogène atomique ($M_odot$) | Observations radio à 21 cm (SPARC) |
| $\NUpsilon_\Nstar$ (étoile) | Rapport masse/lumière stellaire à 3,6 µm | Universel fixe : 0,5\N$, M_\odot/L_\odot$ (McGaugh 2014) |
Pour la Voie Lactée, $R_d$, $Sigma_d$, et $M_text{HI}$ sont remplacés par les valeurs analogues déterminées à partir des relevés stellaires internes (Bovy & Rix 2013) et des cartes à 21 cm. Le même vecteur d’entrée à cinq quantités est utilisé.
3. Décomposition baryonique – cinq composantes géométriques
À partir des cinq données d’observation, la masse baryonique est divisée en cinq composantes géométriques distinctes. Chaque composante possède son propre profil de densité et sa propre échelle caractéristique.
3.1 Masse totale des étoiles et du gaz
$$M_\Nstar \N;=\N ; 2\Npi\N;R_d^2\N;\NSigma_d\N;\NUpsilon_\Nstar$$$$
$$M_\text{gas} \N;=\N;1.33\N;M_\text{HI} \qquad \text{(He correction ; Arnett 1996)}$$$
3.2 Masses et échelles des composants
| Composant | Masse | Échelle | Activation |
|---|---|---|---|
| Le bourrelet | $M_b = 0,20\N- M_\Nétoile$ | $r_b = \max(0.5\N,R_d,\N,0.3\Ntext{ kpc})$ | Si $T \leq 4$ |
| Disque mince | M_\text{thin} = 0,75\,(M_\star – M_b)$ | R_d$ | Toujours |
| Disque épais | M_\text{thick} = 0,25\,(M_\star – M_b)$ | $1.5\,R_d$ | Toujours |
| Anneau de gaz | $M_\text{gas} = 1.33\\NM_\text{HI}$ | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997) | Toujours |
| Bras en spirale | $M_\text{arm} = 0,10\\NM_\text{thin}$ (effectif) | $R_d$ (suit le disque mince) | Toujours |
3.3 Profils de densité
Renflement (Hernquist 3D)
$$\rho_b(r) \;=\ ; \frac{M_b\,r_b}{2\pi\r,r\,(r + r_b)^3}$$$.
Disques stellaires minces et épais (exponentielle 2D)
$$\Sigma_\text{thin}(R) \;=\; \frac{M_\text{thin}}{2\pi\,R_d^2}\,e^{-R/R_d}$$
$$\Sigma_\text{thick}(R) \;=\; \frac{M_\text{thick}}{2\pi\,(1.5\,R_d)^2}\,e^{-R/(1.5R_d)}$$
Anneau de gaz (exponentielle 2D avec trou central)
$$\Sigma_\text{gas}(R) \;=\ ; \frac{M_\text{gas}}{2\pi\ R_g^2}\,\exp!\à gauche(-\frac{R_\text{gaz}{R} – \frac{R}{R_g}\à droite), \quad R_\text{gaz} = 0.5\ R_g$$$.
Excès du bras spiralé (2D, suit un disque mince)
$$\Sigma_\text{bras}(R) \;=\ ; 0.10\\Sigma_\text{mince}(R)$$
4. Le noyau d’onde
Chaque élément de masse baryonique génère un champ d’ondes BeeTheory. Le champ en un point $vec{r}$ produit par un élément source à $vec{r},’$ séparé par $D = |vec{r} – vec{r},’|$ est régi par le noyau de type Yukawa dérivé de la fonction d’onde régularisée de la Note I :
Noyau d’onde de la théorie des abeilles
$$mathcal{K}_i(D) \;=\ ; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\ ; \frac{1}{\ell_i}$$$$
Ici, $K_0$ est l’amplitude universelle de la masse d’onde (un seul nombre sans dimension) et $\ell_i$ est la longueur de cohérence de la composante $i$. Le noyau encode un comportement quasi-newtonien $1/D^2$ à de courtes séparations, modulé par une coupure exponentielle à des échelles au-delà de $\ell_i$. La forme $(1 + \alpha D)\,e^{-\alpha D}$ assure la continuité et une masse totale enfermée finie à l’infini.
4.1 Longueurs de cohérence des composantes
La longueur de cohérence de chaque composant est déterminée par son échelle géométrique naturelle, multipliée par une constante sans dimension propre à sa dimensionnalité :
| Composant | Longueur de cohérence | Constante géométrique |
|---|---|---|
| Bulbe (sphère 3D) | $\ell_b = c_\text{sph}\,r_b$ | $c_\text{sph}$ |
| Disque mince (2D) | $\ell_\text{thin} = c_\text{disk}\,R_d$ | $c_\texte{disque}$ |
| Disque épais (2D) | $\N-\N-\N-\N-{épaisseur} = c_\N-\N-{disque}\N-\N(1.5\N,R_d)$ | $c_\texte{disque}$ |
| Anneau de gaz (2D) | $\ell_\text{gaz} = c_\text{disk}\,R_g$ | $c_\texte{disque}$ |
| Bras en spirale (2D, concentration azimutale) | $\ell_\text{arm} = c_\text{arm}\,R_d$ | $c_\text{arm}$ |
Les trois constantes géométriques $(c_\text{sph},\\text{disk},\text{arm})$ sont universelles – elles ne varient pas d’une galaxie à l’autre. Avec l’amplitude globale de la masse d’onde $K_0$ et le couplage champ d’onde $\lambda$, ils constituent l’ensemble des paramètres du niveau théorique.
5. Convolution du champ d’ondes – équations intégrales par composante
La densité du champ d’ondes à une position $\vec{r}$ est la convolution de la distribution des sources baryoniques avec le noyau d’ondes. Pour un système à symétrie galactique (symétrie axiale, approximation monopolaire), chaque composante baryonique contribue de manière additive :
Densité totale du champ d’ondes à un rayon de $r$.
$$\rho_\text{wave}(r) \;=\ ; \lambda \;\sum_{i \in \{\text{thin, thick, gas, arm, bulge}\}} \rho_\text{wave}^{(i)}(r)$$
Les cinq intégrales sont écrites ci-dessous, une par composante. Chaque intégrale convertit une distribution de masse baryonique en une distribution de masse du champ d’ondes au même point spatial.
5.1 Bulge – Intégration de la coque en 3D
$$\rho_\text{wave}^{(b)}(r) \;=\ ; \int_0^{r_\text{max}} \Nrho_b(r’)\N;\Nmathcal{K}_b\N!\Nà gauche(\Nsqrt{r^2 + r’^2}\Nà droite)\N;4\Npi r’^2\N,dr’$$
L’intégration se fait sur des coques sphériques concentriques de rayon $r’$. Le point de champ au rayon $r$ du centre voit chaque coquille à une séparation effective $D = \sqrt{r^2 + r’^2}$ dans l’approximation monopolaire. L’intégration s’étend jusqu’à $r_\text{max} = 6\,r_b$, au-delà de laquelle la densité du renflement est numériquement négligeable.
5.2 Intégration disque mince – anneau 2D
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{thin}(R’)\N;\Mathcal{K}_\text{thin}\!\Nà gauche(\sqrt{r^2 + R’^2}\Nà droite)\N;2\pi R’\N,dR’$$
Le disque est décomposé en anneaux concentriques de rayon $R’$ et de largeur infinitésimale $dR’$, chacun portant une masse surfacique $\Sigma_\text{thin}(R’)\,2\pi R’\,dR’$. La même approximation monopolaire s’applique : le champ d’ondes au rayon $r$ du centre reçoit des contributions de chaque anneau à la séparation effective $D = \sqrt{r^2 + R’^2}$. Le domaine d’intégration est $R_\text{max} = 8\,R_d$.
5.3 Intégration disque épais – anneau 2D
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{épaisseur}(R’)\N;\Mathcal{K}_\text{épaisseur}\!\Nà gauche(\sqrt{r^2 + R’^2}\Nà droite)\N;2\pi R’\N,dR’$$
Identique à l’intégration du disque mince, avec $\Sigma_\text{thick}(R’)$ comme densité de source et un paramètre de noyau $\alpha_\text{thick} = 1/(c_\text{disk}\\N,\cdot 1.5\N,R_d)$. L’étendue radiale plus large du disque épais se traduit par une gamme de cohérence d’onde légèrement plus large.
5.4 Anneau de gaz – intégration d’un anneau en 2D avec appauvrissement central
$$\rho_\text{wave}^{(\text{gas})}(r) \;=\; \int_{R_\text{hole}}^{R_\text{max}} \Sigma_\text{gas}(R’)\;\mathcal{K}_\text{gas}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
La distribution du gaz a un trou central, capturé par la coupure radiale à $R_\text{hole} = 0.5\,R_g$ dans la limite inférieure de l’intégration. En dehors de cette coupure, le gaz s’étend plus loin que le disque stellaire, ce qui se traduit par une échelle caractéristique plus grande $R_g = 1.7\,R_d$, qui alimente la longueur de cohérence $\ell_\text{gas} = c_\text{disk}\,R_g$.
5.5 Excès de bras en spirale – intégration d’un anneau en 2D avec une amplitude réduite
$$\rho_\text{wave}^{(\text{arm})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_\text{arm}(R’)\;\mathcal{K}_\text{arm}\!\left(\sqrt{r^2 + R’^2}\right)\;2\pi R’\,dR’$$
Les bras spiraux sont traités comme un renforcement moyen axial de la densité de surface du disque mince au niveau de $10\%$, avec leur propre longueur de cohérence $\ell_\text{arm} = c_\text{arm}\,R_d$. Le noyau est donc plus étroit que celui du disque mince, ce qui reflète la concentration azimutale de la structure en spirale.
6. Masse d’onde enfermée et courbe de rotation prévue
Une fois connue la densité totale du champ d’ondes $\rho_\text{wave}(r)$, la masse du champ d’ondes enfermée dans une sphère de rayon $R$ est obtenue par intégration radiale :
Masse du champ d’ondes fermé
$$M_\text{wave}(R) \;=\ ; \int_0^{R} 4\pi\,r^2\,\rho_\text{wave}(r)\nbsp;dr$$$.
La vitesse circulaire prédite au rayon $R$ découle alors de la relation newtonienne, combinant les contributions baryoniques et du champ d’ondes en quadrature :
Vitesse circulaire prévue
$$V_c^2(R) \;=\ ; V_\text{bar}^2(R) \;+\ ; \frac{G\\N,M_\text{wave}(R)}{R}$$$
La vitesse baryonique $V_\text{bar}(R)$ est elle-même la somme quadratique des contributions des quatre composantes du disque (formule de Freeman 1970 pour chaque profil exponentiel) et du bulbe (formule de masse fermée de Hernquist) :
$$V_\text{bar}^2(R) \;=\ ; V_\text{bulge}^2 + V_\text{thin}^2 + V_\text{thick}^2 + V_\text{gas}^2$$.
où chaque $V_i(R)$ est la vitesse circulaire newtonienne standard de la distribution de masse correspondante.
7. Paramètres au niveau de la théorie
Le cadre complet de la Théorie de l’abeille, tel qu’il est appliqué aux galaxies, contient cinq paramètres au niveau de la théorie. Ces paramètres sont universels : ils ne varient pas d’une galaxie à l’autre.
| Symbole | Signification | Rôle |
|---|---|---|
| $K_0$ | Amplitude de la masse d’onde | Définit l’échelle sans dimension du noyau d’onde |
| $c_\text{sph}$ | Constante géométrique 3D | Rapport $\ell/r_\text{scale}$ pour les sources sphériques (bulbe) |
| $c_\texte{disque}$ | Constante géométrique 2D | Rapport $\ell/R_\text{scale}$ pour les sources en disque et en anneau |
| $c_\text{arm}$ | Constante géométrique de la spirale | Rapport $\ell/R_d$ pour l’excès de bras concentré azimutalement |
| $\lambda$ | Couplage global champ d’ondes | Échelle de la densité totale du champ d’ondes |
Universalité des paramètres
Les cinq paramètres sont globaux. Les mêmes valeurs numériques s’appliquent à la Voie lactée, aux irrégulières naines et aux spirales massives. L’information spécifique à la galaxie n’entre que par les cinq entrées observationnelles $(T,\N- R_d,\N- Sigma_d,\N- M_text{HI},\N- Upsilon_\Nstar)$. Le modèle ne contient aucun paramètre accordable par galaxie.
8. La nature unidirectionnelle du calcul
Une chaîne ouverte – pas de rétroaction
L’ensemble du calcul va des entrées aux sorties, dans une seule direction. Les observations photométriques et à 21 cm déterminent la décomposition baryonique. La décomposition baryonique détermine la densité du champ d’ondes. La densité du champ d’ondes détermine la masse d’ondes incluse. La masse d’onde enfermée détermine la courbe de rotation prédite. La courbe de rotation n’influence à aucun moment une étape antérieure du calcul.
Cette unidirectionnalité a trois conséquences importantes.
(a) Une fois les cinq paramètres théoriques fixés, la courbe de rotation est une prédiction stricte et non un ajustement. La comparaison avec la courbe de rotation observée est un test et non un étalonnage.
(b) Le modèle ne comporte aucun mécanisme d’ajustement galaxie par galaxie. Toute modification de la prédiction de la courbe de rotation doit provenir d’une modification du vecteur d’entrée $(T,\N-R_d,\N-Sigma_d,\N-M_text{HI},\N-Upsilon_\Nstar)$ ou d’un changement des paramètres théoriques universels $(K_0,\N-c_text{sph},\N-c_text{disk},\N-c_text{arm},\N-lambda)$.
(c) Calibrer $lambda$ sur une galaxie de référence n’est pas la même chose que l’ajuster à la courbe de rotation de cette galaxie. L’étalonnage détermine un nombre global unique ; la courbe de rotation à tous les autres rayons de la galaxie de référence, et les courbes de rotation de toutes les autres galaxies, sont alors des prédictions strictes du cadre étalonné.
9. Le rôle de la densité de la surface centrale (Note XI révision)
Le diagnostic de la Note XI a identifié que l’erreur de prédiction résiduelle est fortement corrélée avec la densité baryonique centrale de surface $\Sigma_d$, indépendamment de la longueur d’échelle du disque $R_d$. La formalisation présentée ci-dessus est la version du modèle avant que cette découverte ne soit incorporée – elle n’utilise que $R_d$ dans les expressions de la longueur de cohérence $\ell_i = c_i\,R_d$.
Où le raffinement entrera en jeu
Dans le modèle affiné, les longueurs de cohérence $\ell_i$ dépendront à la fois de $R_d$ et de $\Sigma_d$, remplaçant la relation linéaire stricte $\ell_i = c_i\,R_d$ par une fonction $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_\text{ref})$ qui absorbe le résidu identifié dans la note XI. La forme fonctionnelle de $\phi$ et ses paramètres seront déterminés dans les notes suivantes, d’abord sur l’ensemble d’étalonnage de 22 galaxies, puis validés par prédiction aveugle sur le reste de l’échantillon SPARC.
La structure unidirectionnelle du calcul est préservée par ce raffinement : $\Sigma_d$ est une entrée observationnelle, les longueurs de cohérence modifiées alimentent les mêmes intégrales de convolution, et la courbe de rotation émerge comme auparavant. Seul un lien opérationnel est ajouté – la dépendance de $\ell_i$ à l’égard d’une seconde observable.
10. Résumé de la méthodologie
1. Données d’entrée. Cinq observables par galaxie : le type de Hubble $T$, l’échelle du disque $R_d$, la luminosité de surface $\Sigma_d$, la masse HI $M_\text{HI}$, et le rapport universel masse stellaire/lumière $\Upsilon_\star$.
2. Décomposition baryonique. Cinq composantes : bulbe (si $T \leq 4$), disque mince, disque épais, anneau de gaz, excès du bras spiral. Chacune porte un profil analytique de densité.
3. Noyau d’onde. Forme universelle de type Yukawa $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ avec une longueur de cohérence $\ell_i = c_i\,R_\text{scale}$ déterminée par l’étendue géométrique de chaque composant.
4. Convolution. Chaque composante génère une densité de champ d’ondes par le biais d’une intégrale unidimensionnelle sur les anneaux (composantes 2D) ou les coquilles (renflement 3D). La densité totale du champ d’ondes est la somme des cinq composantes, mise à l’échelle par le couplage global $\lambda$.
5. Résultat. La masse d’onde incluse $M_\text{wave}(R)$ est intégrée et combinée avec la vitesse baryonique $V_\text{bar}(R)$ pour obtenir la courbe de rotation prédite $V_c(R)$.
6. Paramètres au niveau de la théorie. (K_0,\Nc_texte{sph},\Nc_texte{disque},\Nc_texte{bras},\Nc_texte{bras},\Nc_lambda)$ – universel, pas de réglage par galaxie. Un raffinement en cours d’étude ajoutera une dépendance à $\Sigma_d$.
7. Direction. Entrées → baryons → champ d’ondes → courbe de rotation. Pas de rétroaction. La courbe de rotation est une prédiction, pas un ajustement.
Références. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC : Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves, AJ 152, 157 (2016). – Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970). – Hernquist, L. – Un modèle analytique pour les galaxies sphériques et les bulbes, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – Short 21-cm WSRT observations of spiral and irregular galaxies, A&A 324, 877 (1997). – McGaugh, S. S. – The third law of galactic rotation, Galaxies 2, 601 (2014). – Bovy, J., Rix, H.-W. – A direct dynamical measurement of the Milky Way’s disk surface density profile, ApJ 779, 115 (2013). – Arnett, D. – Supernovae and Nucleosynthesis, Princeton (1996). – Dutertre, X. – Bee Theory™ : Wave-Based Modeling of Gravity, v2, BeeTheory.com (2023).
BeeTheory.com – Gravité quantique basée sur les ondes – Méthodologie galactique – © Technoplane S.A.S. 2026