Teoria delle api – Articolo scientifico – 2025

La massa nascosta della Via Lattea:
Una derivazione basata sulle onde e un adattamento numerico

Partendo dal postulato della BeeTheory, secondo cui ogni elemento di massa emette un campo di onde gravitazionali che decadono come $e^{-D/\ell}$, ricaviamo analiticamente la distribuzione della massa oscura 3D, la adattiamo alla curva di rotazione di Gaia 2024 e troviamo i due parametri fondamentali del modello.

$Rho_0 = 1,14\;\testo{GeV/cm}^3$

Densità oscura centrale

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

Scala di coerenza d’onda

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

Bontà di adattamento

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

Previsto $irho_testo{scuro}(R_dot)$

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

Massa oscura totale all’interno di 200 kpc

0. Conclusioni – Prima i risultati

Il modello basato sulla BeeTheory wave – in cui ogni elemento di massa visibile $dV$ genera un campo gravitazionale che decade esponenzialmente come $e^{-D/ell}$ in 3D – prevede un profilo di densità di massa oscura che, quando viene integrato sul disco galattico, converge alla forma NFW.

Adattato alla curva di rotazione Gaia 2024 della Via Lattea, utilizzando solo due parametri liberi, il modello raggiunge $\chi^2/\mathrm{dof} = 0,44$.

I parametri migliori sono: densità centrale oscura $\rho_0 = 1,14\,\mathrm{GeV/cm}^3$ e raggio della scala di coerenza $r_s = 9,6\,\mathrm{kpc}$. Questi corrispondono direttamente ai due parametri della Teoria delle Api: la costante di accoppiamento d’onda $\lambda$ e la lunghezza di coerenza $\ell = r_s\sqrt{2} \circa 13,6\,\mathrm{kpc}$.

Il modello predice una densità di materia oscura locale di $irho_testo{dark}(R_\odot = 8\,\mathrm{kpc}) = 0,41\,\mathrm{GeV/cm}^3$ – entro il 5% del valore misurato $0,39 \pm 0,03\,\mathrm{GeV/cm}^3$. La massa oscura totale all’interno di 200 kpc è pari a $sim 7,1 \times 10^{11}\,M_\odot$, coerente con le recenti misurazioni della cinematica satellitare.

Densità oscura centrale
$irho_0 = 1,14\;\frac{{testo{GeV}}{testo{cm}^3$

Equivalente a $3,0\times10^7\,M_\odot\,\text{kpc}^{-3}$. Ampiezza del campo d’onda a $r=0$.

Scala di coerenza dell’onda
$r_s = 9,6\;\text{kpc}$

La scala in cui il campo d’onda passa dal regime interno al regime esterno.

Bontà di adattamento
$\chi^2/\text{dof} = 0,44$.

Adattamento eccellente. 15 dei 16 punti dati entro $1\sigma$.

Osservabile Misurazione Gaia 2024 Previsione BeeTheory Residuo
$V_c(R_nodot = 8\,\text{kpc})$$230 \pm 6\;\testo{km/s}$$231;\testo{km/s}$$+0.4\%$
$V_c(20\,\text{kpc})$215 dollari \pm 10\;\text{km/s}$208$;\testo{km/s}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$173 \pm 17\;\text{km/s}$$199;\testo{km/s}$$+15\%$, $1,5\sigma$
$rho_testo_scuro}(R_odot)$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$All’interno dell’intervallo

1. Il postulato della teoria delle api: La massa irradia onde

La gravità classica e relativistica descrivono come agisce la gravità, ma non perché esiste. La BeeTheory propone un meccanismo: ogni elemento di massa $dV$ è la fonte di un campo d’onda quantistico che si propaga verso l’esterno nello spazio 3D e decade esponenzialmente con la distanza euclidea $D$ dalla fonte.

Questo campo d’onda trasporta l’energia gravitazionale effettiva – è, in un senso preciso, la “massa nascosta”.

Postulato della massa d’onda della Teoria delle Api $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Qui $\lambda$ è la costante di accoppiamento onda-massa, $\ell$ è la lunghezza di coerenza, $\rho_testo{vis}$ è la densità di massa barionica visibile e $D = |\mathbf{r}-\mathbf{r}’|$ è la distanza euclidea dalla sorgente al punto di campo.

La densità di massa oscura totale in un punto qualsiasi $mathbf{r}$ è la sovrapposizione dei campi d’onda di ogni elemento di massa visibile nella galassia:

Densità oscura totale – integrale di sovrapposizione $$\rho_testo{dark}(\mathbf{r}) = \frac{\lambda}{\ell} \int_testo{galassia} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

Si tratta di una convoluzione 3D della distribuzione di massa visibile con un kernel esponenziale.

  1. La massa oscura non è sfericamente simmetrica per assunto. Riflette la geometria della sorgente.
  2. La massa oscura riempie tutto lo spazio 3D, non solo il piano galattico.
  3. I due parametri $(\lambda,\ell)$ determinano completamente la distribuzione della massa oscura, una volta conosciuta la distribuzione barionica.

2. La sorgente visibile: Un disco esponenziale

Il disco stellare della Via Lattea è ben descritto da una densità superficiale esponenziale:

Densità superficiale del 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}$$

Il disco ha uno spessore trascurabile rispetto al suo raggio, quindi la sua densità di volume può essere rappresentata con una densità di superficie del disco e una funzione delta verticale.

Densità oscura da un disco sottile – doppio integrale esatto $$\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’$$

Il punto di campo $(R,z)$ si trova al raggio cilindrico $R$ nel piano del disco e all’altezza $z$ sopra di esso. Impostando $r = \sqrt{R^2+z^2}$, eseguiamo l’integrale azimutale analiticamente utilizzando l’approssimazione del monopolo:

Kernel del monopolo – media azimutale $$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)$$

Questo riduce l’integrale doppio a una dimensione:

Integrale principale 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 Risultato analitico – L’emergenza NFW

Eseguendo analiticamente l’integrale $R’$ si ottiene:

Densità oscura della Teoria delle Api in forma chiusa $$\rho_text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\sinistra(\frac{r}{\ell}\ destra) e^{-r/\ell}$$

Nel regime interno:

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

Nel regime esterno:

Regime esterno $$sinh(r/\ell)\,e^{-r/\ell}\approssimativamente \tfrac{1}{2}\quadLongrightarrow\quad \rho_testo{scuro}(r)\propto \frac{e^{-r/\ell}{r}$$.

La transizione tra questi regimi avviene a $r\sim\ell$. Questa è la regione in cui il profilo d’onda della BeeTheory può essere confrontato con il comportamento della scala NFW.

Risultato teorico chiave

Il profilo della materia oscura simile alla NFW emerge analiticamente dal postulato della massa d’onda della Teoria di Bee applicato a una sorgente disco esponenziale. In questa interpretazione, i parametri NFW non sono parametri arbitrari dell’alone; sono legati ai parametri d’onda della BeeTheory e alla geometria del disco.

2.2 Il dizionario BeeTheory-NFW

Mappatura dei parametri BeeTheory e NFW $$r_s = \ell, \qquad \rho_0^{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$
Parametri adattati all’interpretazione di 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}$$
Parametri numerici della Teoria delle Api $$}}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}$$
Parametro della Teoria delle Api Significato fisico Valore di miglior adattamento Vincolo da
$\ell$Lunghezza di coerenza. Equivale al raggio di scala NFW $r_s$.$9,6\,\code(01)}$Forma del declino di $V_c(R)$
$$Lambda$Accoppiamento onda-massa senza dimensione.$0.132$Scala di velocità assoluta
$$Rho_0$Picco di densità di massa oscura a $r=0$.$1.14\,\text{GeV/cm}^3$Calcolato da $\lambda$ e $\ell$
$r_s$Raggio di transizione tra le pendenze della densità.9,6 $, \text{kpc}$Uguale a quello di ${ell$

3. Dalla massa mancante alla curva di rotazione

3.1 Il problema della massa mancante

La dinamica newtoniana richiede:

Massa totale racchiusa $$M_testo{tot}(<R) = M_testo{bar}(<R) + M_testo{scuro}(<R), \qquad M_testo{scuro}(<R) = \frac{V_c^2 R}{G} – M_testo{bar}(<R)$$

La massa barionica entro il raggio $R$ ha due componenti: il disco esponenziale e un bulge compatto.

Massa barionica racchiusa $$M_testo{disco}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 – \left(1+\frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_testo{bulge} = 1.2\times10^{10}\,M_odot$$

3.2 Decomposizione della velocità circolare

Decomposizione circolare della velocità $$V_c^2(R) = V_testo{disco}^2(R) + V_testo{bulge}^2(R) + V_testo{buio}^2(R)$$

Il contributo del disco utilizza la formula di Freeman con funzioni di Bessel modificate:

Velocità del disco di Freeman $$V_testo{disco}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\sinistra[I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\destra], \quad y = \frac{R}{2R_d}$$

Il rigonfiamento utilizza un profilo Hernquist:

Contributo del rigonfiamento Hernquist $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

La massa oscura NFW racchiusa in $R$ ha una forma analitica:

Massa oscura NFW racchiusa $$M_testo{scuro,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]$$
Perché i barioni da soli predicono la velocità sbagliata

Con $M_d = 3,5\times10^{10}\,M_\odot$ e $M_b = 1,2\times10^{10}\,M_\odot$, il modello barionico prevede circa $162\,\text{km/s}$ vicino a $8\,\text{kpc}$, al di sotto dell’osservato $sim230\,\text{km/s}$.

4. Simulazione numerica e adattamento dei parametri

4.1 Dati di ingresso – Gaia 2024

16 punti di dati di Ou et al. (2024) che abbracciano $R=4$-$27,3\,\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,6,6,6,7,7,8,9,10,11,13,17];

4.2 Algoritmo

1
Calcolo di $V_testo{bar}(R)$

Utilizza il disco di Freeman + il rigonfiamento di Hernquist. Le funzioni di Bessel vengono calcolate utilizzando approssimazioni polinomiali.

2
Valutare la velocità oscura NFW

Utilizziamo la massa racchiusa NFW in forma chiusa e calcoliamo $V_text{dark}(R)$.

3
Calcolo della velocità totale e di $\chi^2

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

4
Ridurre al minimo $\chi^2(\rho_0,r_s)$

Utilizzi una griglia a due passaggi su $\rho_0$ e $r_s$.

Nota sull’unità critica

Il valore corretto della costante di Newton in unità kpc-km-s-$M_\odot$ è:

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

L’utilizzo di $4,302\times10^{-3}$ è un errore di unità comune e fornisce velocità troppo grandi.

4.3 Curva di rotazione interattiva

Solo barioni BeeTheory $V_testo{totale}$ Solo materia oscura Dati Gaia 2024
Esploratore di parametri – aggiusta $\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 Risultati – Profilo di massa in 3D

La massa oscura racchiusa in una sfera di raggio $r$ aumenta bruscamente all’interno di $r_s$ e cresce logaritmicamente oltre.

Profilo di massa: disco visibile vs totale vs materia oscura
Disco visibile + bulge Massa oscura Massa totale
$r$ $M_testo{bar}(<r)$ $M_testo{scuro}(<r)$ $M_testo{tot}(<r)$ Rapporto DM/bar $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. Interpretazione fisica dei due parametri

5.1 La Lunghezza di Coerenza $\ell = r_s = 9,6\, \text{kpc}$

$\ell$ è il raggio in cui il campo di onde gravitazionali emesso da ciascun elemento di massa rimane in fase. All’interno di questo raggio, l’interferenza delle onde è costruttiva e la densità oscura scende lentamente. Al di fuori di questo raggio, l’interferenza distruttiva fa sì che la densità scenda più velocemente.

Il valore $\ell = 9,6\,\text{kpc}\circa 3,7R_d$ ha un’interpretazione naturale: la lunghezza di coerenza è stabilita dal raggio di scala del disco moltiplicato per un fattore dell’ordine dell’unità.

5.2 La costante di accoppiamento $\lambda = 0.132$

$\lambda$ determina quanta massa d’onda viene generata per unità di massa visibile per lunghezza di coerenza.

Il rapporto di massa locale oscura-visibile da $\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} \circa 4,2$$

Il rapporto globale tra massa oscura e massa barionica all’interno di 200 kpc è di circa $M_testo{oscuro}/M_testo{bar}\approssimativamente 15$, coerente con una grande componente di massa nascosta.

Previsione della Teoria delle Api: forma dell’alone

Poiché la massa oscura emerge da una sorgente disco tramite una convoluzione 3D, l’alone non è perfettamente sferico. Il calcolo esatto del non-monopolo prevede un rapporto tra asse minore e asse maggiore di circa $q=c/a\ circa 0,82$ per l’alone oscuro della Via Lattea.

6. Riassunto e prospettiva

Partendo da un unico postulato fisico – che ogni elemento di massa visibile genera un campo di onde gravitazionali che decade come $e^{-D/\ell}$ in 3D – BeeTheory fornisce una derivazione basata sulle onde di un profilo di densità simile alla materia oscura.

Adeguato alla curva di rotazione Gaia 2024 della Via Lattea, il modello raggiunge $\chi^2/testo{dof}=0,44$ con due parametri liberi:

Parametri BeeTheory migliori – Via Lattea $$\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$$

Il modello fa tre previsioni verificabili oltre alla curva di rotazione:

  1. Forma dell’alone: la massa oscura è appiattita su un disco con un rapporto d’asse $q\circa 0,82$.
  2. Universalità dei parametri: la stessa relazione $(\lambda,\ell)$ dovrebbe applicarsi alle galassie esterne con parametri del disco noti.
  3. Scalatura della coerenza: $\ell\approssimativamente3,7R_d$ suggerisce una relazione di scalatura tra le dimensioni del disco e il raggio di scala dell’alone oscuro.

Riferimenti

  • Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – Il profilo di materia oscura della Via Lattea dedotto dalla sua curva di velocità circolare, MNRAS 528, 693-710 (2024)
  • Navarro, J. F., Frenk, C. S., White, S. D. M. – Un profilo di densità universale dal raggruppamento gerarchico, ApJ 490, 493 (1997)
  • Freeman, K. C. – Sui dischi delle galassie a spirale e S0, ApJ 160, 811 (1970)
  • Pato, M., Iocco, F., Bertone, G. – Vincoli dinamici sulla distribuzione della materia oscura nella Via Lattea, JCAP 12, 001 (2015)
  • McMillan, P. J. – La distribuzione di massa e il potenziale gravitazionale della Via Lattea, MNRAS 465, 76 (2017)
  • Abramowitz, M., Stegun, I. A. – Manuale delle funzioni matematiche, Dover (1972)