BeeTheory – Fondazioni – Nota tecnica XII
Formalizzazione:
Il calcolo della teoria delle api su scala galattica
Questa nota formalizza il quadro BeeTheory applicato a una galassia a disco. Specifica gli input osservativi, la decomposizione geometrica della distribuzione barionica, le equazioni integrali che definiscono il campo d’onda per ogni componente e la catena di operazioni che produce la curva di rotazione prevista. La procedura è strettamente unidirezionale: la struttura barionica osservata determina il campo d’onda, che determina la curva di rotazione – mai il contrario.
1. Il calcolo in un diagramma
Una catena unidirezionale
Fotometria osservata $\;\longrightarrow;$ Decomposizione barionica $(\rho_testo{bar})$
$molto grande$
Convoluzione del campo d’onda $\;\longrightarrow};$ Densità d’onda $(\rho_testo{wave})$
$molto grande
Integrazione di massa $\;\longrightarrow;$ Massa d’onda racchiusa $(M_testo{onda})$
$molto in basso
Relazione newtoniana $\;\longrightarrow\;$ Curva di rotazione prevista $(V_c)$
Nessun passo viene invertito. La curva di rotazione $V_c(R)$ non viene mai utilizzata come input.
2. Input osservativi
Per ogni galassia, il calcolo richiede cinque osservabili pubblicati. Queste sono le uniche quantità specifiche della galassia; tutto il resto viene calcolato a partire da esse. In questa fase non viene eseguito alcun adattamento alla curva di rotazione.
| Simbolo | Quantità | Fonte |
|---|---|---|
| $T$ | Tipo morfologico di Hubble | Catalogo (de Vaucouleurs et al. 1991, SPARC) |
| $R_d$ | Lunghezza di scala del disco stellare (kpc) | Fotometria Spitzer 3,6 µm (SPARC) |
| $Sigma_d$ | Luminosità della superficie del disco centrale ($L_\odot/\text{pc}^2$) | Fotometria Spitzer 3,6 µm (SPARC) |
| $M_testo{HI}$ | Massa atomica totale dell ‘idrogeno ($M_odot$) | Osservazioni radio a 21 cm (SPARC) |
| $Upsilon_Stella$ | Rapporto massa-luce stellare a 3,6 µm | Universale fisso: $0,5\,M_\odot/L_\odot$ (McGaugh 2014) |
Per la Via Lattea, $R_d$, $Sigma_d$ e $M_text{HI}$ sono sostituiti dai valori analoghi determinati dalle indagini stellari interne (Bovy & Rix 2013) e dalle mappe a 21 cm. Viene utilizzato lo stesso vettore di input a cinque quantità.
3. Decomposizione barionica – cinque componenti geometriche
Dai cinque input osservativi, la massa barionica viene suddivisa in cinque componenti geometriche distinte. Ogni componente ha un proprio profilo di densità e una scala caratteristica.
3.1 Masse stellari e di gas totali
$$M_star \;=\; 2\pi\,R_d^2\,\Sigma_d\,\Upsilon_star$$.
$$M_testo{gas} \;=\; 1,33\,M_testo{HI} \qquad \text{(correzione di He; Arnett 1996)}$$
3.2 Masse e scale dei componenti
| Componente | Massa | Scala | Attivazione |
|---|---|---|---|
| Sporgenza | $M_b = 0,20\, M_star$ | $r_b = \max(0,5\,R_d,\,0,3\text{ kpc})$ | Se $T \leq 4$ |
| Disco sottile | $M_testo{sottile} = 0,75\,(M_stella – M_b)$ | $R_d$ | Sempre |
| Disco spesso | $M_testo{spessore} = 0,25\,(M_stella – M_b)$ | $1.5\,R_d$ | Sempre |
| Anello di gas | $M_testo{gas} = 1,33\, M_testo{HI}$. | $R_g = 1,7\,R_d$ (Broeils & Rhee 1997). | Sempre |
| Bracci a spirale | $M_testo{braccio} = 0,10\, M_testo{sottile}$ (effettivo) | $R_d$ (segue disco sottile) | Sempre |
3.3 Profili di densità
Rigonfiamento (3D Hernquist)
$$\rho_b(r) \;=\; \frac{M_b\,r_b}{2\pi\,r\, (r + r_b)^3}$$
Dischi stellari sottili e spessi (esponenziale 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)}$$
Anello di gas (esponenziale 2D con foro centrale)
$$Sigma_testo{gas}(R) \;=\; \frac{M_testo{gas}}{2\pi\,R_g^2}\,\exp\!\left(-\frac{R_testo{foro}}{R} – \frac{R}{R_g}\right), \quad R_testo{foro} = 0,5\,R_g$$
Eccesso di braccio a spirale (2D, segue disco sottile)
$$\Sigma_testo_{braccio}(R) \;=\; 0,10\;\Sigma_testo_{sottile}(R)$$
4. Il kernel d’onda
Ogni elemento di massa barionica genera un campo d’onda in BeeTheory. Il campo in un punto $vec{r}$ prodotto da un elemento sorgente a $vec{r},’$ separato da $D = |vec{r} – vec{r},’|$ è governato dal kernel di tipo Yukawa derivato dalla funzione d’onda regolarizzata della Nota I:
Kernel d’onda della Teoria delle Api
$$\mathcal{K}_i(D) \;=\; K_0\,\frac{(1 + \alpha_i\,D)\,e^{-\alpha_i\,D}}{D^2}, \qquad \alpha_i \;=\; \frac{1}{ell_i}$$.
Qui $K_0$ è l’ampiezza universale della massa d’onda (un singolo numero senza dimensione) e $\ell_i$ è la lunghezza di coerenza del componente $i$. Il kernel codifica un comportamento quasi newtoniano $1/D^2$ a brevi separazioni, modulato da un cutoff esponenziale a scale superiori a $\ell_i$. La forma $(1 + \alpha D)\,e^{-\alpha D}$ assicura la continuità e la massa totale racchiusa finita all’infinito.
4.1 Lunghezze di coerenza dei componenti
La lunghezza di coerenza di ogni componente è stabilita dalla sua scala geometrica naturale, moltiplicata per una costante adimensionale specifica della sua dimensionalità:
| Componente | Lunghezza di coerenza | Costante geometrica |
|---|---|---|
| Bulge (sfera 3D) | $\ell_b = c_\text{sph}\,r_b$ | $c_testo{sph}$ |
| Disco sottile (2D) | $$$ell_testo{sottile} = c_testo{disco}},R_d$ | $c_testo{disco}$ |
| Disco spesso (2D) | $$$ell_testo{spessore} = c_testo{disco}\,(1.5\,R_d)$ | $c_testo{disco}$ |
| Anello di gas (2D) | $\ell_\text{gas} = c_\text{disk}\,R_g$ | $c_testo{disco}$ |
| Bracci a spirale (2D, concentrati azimutalmente) | $$$ell_testo{arm} = c_testo{arm}\,R_d$ | $c_testo{arm}$ |
Le tre costanti geometriche $(c_testo{sph},\,c_testo{disco},\,c_testo{braccio})$ sono universali – non variano da galassia a galassia. Insieme all’ampiezza globale della massa d’onda $K_0$ e all’accoppiamento del campo d’onda $\lambda$, costituiscono l’insieme completo dei parametri a livello teorico.
5. Convoluzione del campo d’onda – equazioni integrali per componente
La densità del campo d’onda in una posizione ${vec{r}$ è la convoluzione della distribuzione della sorgente barionica con il kernel d’onda. Per un sistema galatticamente simmetrico (simmetria assiale, approssimazione monopolare), ogni componente barionica contribuisce in modo additivo:
Densità totale del campo d’onde al raggio $r
$$\rho_testo{onda}(r) \;=\; \lambda \;\somma_{i \in \{{sottile, denso, gas, braccio, rigonfiamento}}} \rho_testo{onda}^{(i)}(r)$$
I cinque integrali sono scritti di seguito, uno per componente. Ogni integrale converte una distribuzione di massa barionica in una distribuzione di massa del campo d’onda nello stesso punto spaziale.
5.1 Bulge – Integrazione del guscio 3D
$$\rho_testo{onda}^{(b)}(r) \;=\; \int_0^{r_testo{max}} \rho_b(r’)\;\mathcal{K}_b\!\sinistra(\sqrt{r^2 + r’^2}\destra)\;4\pi r’^2\,dr’$$
L’integrazione avviene su gusci sferici concentrici di raggio $r’$. Il punto di campo a raggio $r$ dal centro vede ogni guscio a una separazione effettiva $D = \sqrt{r^2 + r’^2}$ nell’approssimazione monopolare. L’integrazione si estende fino a $r_testo{max} = 6\,r_b$, oltre il quale la densità del rigonfiamento è numericamente trascurabile.
5.2 Disco sottile – integrazione ad anello 2D
$$\rho_\text{wave}^{(\text{thin})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_testo{sottile}(R’)\;\mathcal{K}_testo{sottile}\!\sinistra(\sqrt{r^2 + R’^2}\destra)\;2\pi R’\,dR’$$
Il disco viene scomposto in anelli concentrici di raggio $R’$ e larghezza infinitesimale $dR’$, ciascuno dei quali trasporta la massa superficiale $\Sigma_testo{thin}(R’)\,2\pi R’\,dR’$. Si applica la stessa approssimazione monopolare: il campo d’onda al raggio $r$ dal centro riceve contributi da ciascun anello alla separazione effettiva $D = \sqrt{r^2 + R’^2}$. L’intervallo di integrazione è $R_testo{max} = 8\,R_d$.
5.3 Disco spesso – integrazione ad anello 2D
$$\rho_\text{wave}^{(\text{thick})}(r) \;=\; \int_0^{R_\text{max}} \Sigma_testo{spessore}(R’)\;\mathcal{K}_testo{spessore}\!\sinistra(\sqrt{r^2 + R’^2}\destra)\;2\pi R’\,dR’$$
Identico all’integrazione del disco sottile, con $\Sigma_{text{thick}(R’)$ come densità della sorgente e un parametro kernel $\alpha_{text{thick} = 1/(c_{text{disk}\,\cdot 1,5\,R_d)$. L’estensione radiale più ampia del disco spesso comporta un intervallo di coerenza d’onda leggermente più ampio.
5.4 Anello di gas – Integrazione dell’anello 2D con esaurimento centrale
$$\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 distribuzione del gas ha un buco centrale, catturato dal cutoff radiale a $R_{text{hole} = 0,5\,R_g$ nel limite inferiore dell’integrazione. Al di fuori di questo cutoff, il gas si estende più del disco stellare; ciò si riflette nella scala caratteristica più grande $R_g = 1,7\,R_d$, che alimenta la lunghezza di coerenza $\ell_{text{gas} = c_{text{disk}\,R_g$.
5.5 Eccesso del braccio a spirale – Integrazione dell’anello 2D con ampiezza ridotta
$$\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’$$
I bracci a spirale sono trattati come un potenziamento mediato assialmente della densità di superficie del disco sottile al livello di $10\\\i}, con la propria lunghezza di coerenza $\ell_testo{arm} = c_testo{arm}\,R_d$. Il kernel è quindi più stretto del kernel del disco sottile, riflettendo la concentrazione azimutale della struttura a spirale.
6. Massa d’onda chiusa e curva di rotazione prevista
Una volta conosciuta la densità totale del campo d’onda $\rho_testo{onda}(r)$, la massa del campo d’onda racchiusa all’interno di una sfera di raggio $R$ si ottiene mediante integrazione radiale:
Massa del campo d’onda chiuso
$$M_testo{onda}(R) \;=\; \int_0^{R} 4\pi\,r^2\,\rho_testo{onda}(r)\,dr$$
La velocità circolare prevista al raggio $R$ deriva quindi dalla relazione newtoniana, combinando i contributi barionici e del campo d’onda in quadratura:
Velocità circolare prevista
$$V_c^2(R) \;=\; V_testo{bar}^2(R) \;+\; \frac{G\,M_testo{onda}(R)}{R}$$.
La velocità barionica $V_testo{bar}(R)$ è di per sé la somma quadratica dei contributi delle quattro componenti discoidali (formula di Freeman 1970 per ogni profilo esponenziale) e del bulge (formula della massa racchiusa di Hernquist):
$$V_testo{bar}^2(R) \;=\; V_testo{bulge}^2 + V_testo{thin}^2 + V_testo{thick}^2 + V_testo{gas}^2$.
dove ogni $V_i(R)$ è la velocità circolare newtoniana standard della distribuzione di massa corrispondente.
7. Parametri a livello teorico
Il quadro completo della BeeTheory, applicato alle galassie, contiene cinque parametri a livello di teoria. Questi sono universali: non variano da galassia a galassia.
| Simbolo | Significato | Ruolo |
|---|---|---|
| $K_0$ | Ampiezza della massa d’onda | Imposta la scala adimensionale del kernel d’onda |
| $c_testo{sph}$ | Costante geometrica 3D | Rapporto $\ell/r_testo{scala}$ per sorgenti sferiche (rigonfiamento) |
| $c_testo{disco}$ | Costante geometrica 2D | Rapporto $\ell/R_testo{scala}$ per le sorgenti a disco e ad anello |
| $c_testo{arm}$ | Costante geometrica della spirale | Rapporto $\ell/R_d$ per l’eccesso di braccio concentrato azimutalmente |
| $\lambda$ | Accoppiamento globale del campo d’onda | Scala la densità totale del campo d’onda |
Universalità dei parametri
Tutti e cinque i parametri sono globali. Gli stessi valori numerici si applicano alla Via Lattea, alle nane irregolari, alle spirali massicce. Le informazioni specifiche della galassia entrano solo attraverso i cinque input osservativi $(T,\,R_d,\,\Sigma_d,\,M_text{HI},\,\Upsilon_\star)$. Il modello non contiene alcun parametro sintonizzabile per galassia.
8. La natura unidirezionale del calcolo
Una catena aperta – nessun feedback
L’intero calcolo scorre dagli ingressi alle uscite, in una sola direzione. Le osservazioni fotometriche e a 21 cm determinano la decomposizione barionica. La decomposizione barionica determina la densità del campo d’onda. La densità del campo d’onda determina la massa d’onda racchiusa. La massa d’onda racchiusa determina la curva di rotazione prevista. In nessun momento la curva di rotazione influenza una fase precedente del calcolo.
Questa unidirezionalità ha tre conseguenze importanti.
(a) Una volta fissati i cinque parametri a livello di teoria, la curva di rotazione è una previsione rigorosa, non un adattamento. Il confronto con la curva di rotazione osservata è un test, non una calibrazione.
(b) Il modello non ha un meccanismo di regolazione galassia per galassia. Ogni modifica della previsione della curva di rotazione deve derivare da una modifica del vettore di ingresso $(T,\,R_d,\,\Sigma_d,\,M_testo{HI},\,\Upsilon_Stella)$ o da una modifica dei parametri universali a livello di teoria $(K_0,\,c_testo{sph},\,c_testo{disco},\,c_testo{braccio},\,\lambda)$.
(c) Calibrare $lambda$ su una galassia di riferimento non significa adattarla alla curva di rotazione di quella galassia. La calibrazione determina un singolo numero globale; la curva di rotazione a tutti gli altri raggi della galassia di riferimento e le curve di rotazione di tutte le altre galassie sono poi previsioni rigorose del quadro calibrato.
9. Il ruolo della densità della superficie centrale (Nota XI revisione)
La diagnostica della Nota XI ha identificato che l’errore di previsione residuo è fortemente correlato con la densità di superficie barionica centrale $\Sigma_d$, indipendentemente dalla lunghezza di scala del disco $R_d$. La formalizzazione presentata qui sopra è la versione del modello prima dell’ incorporazione di questa scoperta – utilizza solo $R_d$ nelle espressioni della lunghezza di coerenza $\ell_i = c_i\,R_d$.
Dove entrerà la raffinatezza
Nel modello perfezionato, le lunghezze di coerenza $\ell_i$ dipenderanno sia da $R_d$ che da $\Sigma_d$, sostituendo la stretta relazione lineare $\ell_i = c_i\,R_d$ con una funzione $\ell_i = c_i\,R_d\,\phi(\Sigma_d/\Sigma_testo{ref})$ che assorbe il residuo identificato nella Nota XI. La forma funzionale di $phi$ e i suoi parametri saranno determinati nelle note successive, prima sul set di calibrazione di 22 galassie, poi convalidati dalla previsione cieca sul campione SPARC rimanente.
La struttura unidirezionale del calcolo viene preservata da questo perfezionamento: $\Sigma_d$ è un input osservativo, le lunghezze di coerenza modificate alimentano gli stessi integrali di convoluzione e la curva di rotazione emerge come prima. Viene aggiunto solo un collegamento operativo: la dipendenza di $\ell_i$ da un secondo osservabile.
10. Riassunto della metodologia
1. Ingressi. Cinque osservabili per ogni galassia: tipo di Hubble $T$, scala del disco $R_d$, luminosità della superficie $Sigma_d$, massa HI $M_testo{HI}$ e il rapporto universale massa stellare/luce $Upsilon_\star$.
2. Decomposizione barionica. Cinque componenti: bulge (se $T \leq 4$), disco sottile, disco spesso, anello di gas, eccesso del braccio di spirale. Ciascuna porta con sé un profilo analitico di densità.
3. Kernel d’onda. Forma universale di tipo Yukawa $\mathcal{K}_i(D) = K_0\,(1 + \alpha_i D)\,e^{-\alpha_i D}/D^2$ con lunghezza di coerenza $\ell_i = c_i\,R_testo{scala}$ determinata dall’estensione geometrica di ogni componente.
4. Convoluzione. Ogni componente genera una densità di campo d’onda tramite un integrale monodimensionale su anelli (componenti 2D) o gusci (rigonfiamento 3D). La densità totale del campo d’onda è la somma dei cinque componenti, scalata dall’accoppiamento globale $\lambda$.
5. Risultato. La massa d’onda racchiusa $M_testo{onda}(R)$ viene integrata e combinata con la velocità barionica $V_testo{bar}(R)$ per ottenere la curva di rotazione prevista $V_c(R)$.
6. Parametri a livello di teoria. $(K_0,\,c_testo{sph},\,c_testo{disco},\,c_testo{braccio},\,\lambda)$ – universale, nessuna regolazione per galassia. Una raffinatezza in fase di studio aggiungerà una dipendenza da $Sigma_d$.
7. Direzione. Ingressi → barioni → campo d’onda → curva di rotazione. Nessun feedback. La curva di rotazione è una previsione, non un adattamento.
Riferimenti. Lelli, F., McGaugh, S. S., Schombert, J. M. – SPARC: Modelli di massa per 175 galassie a disco con fotometria Spitzer e curve di rotazione accurate, AJ 152, 157 (2016). – Freeman, K. C. – Sui dischi delle galassie a spirale e S0, ApJ 160, 811 (1970). – Hernquist, L. – Un modello analitico per galassie sferiche e bulge, ApJ 356, 359 (1990). – Broeils, A. H., Rhee, M.-H. – Osservazioni WSRT brevi a 21 cm di galassie a spirale e irregolari, A&A 324, 877 (1997). – McGaugh, S. S. – La terza legge della rotazione galattica, Galassie 2, 601 (2014). – Bovy, J., Rix, H.-W. – Una misura dinamica diretta del profilo di densità superficiale del disco della Via Lattea, ApJ 779, 115 (2013). – Arnett, D. – Supernovae e Nucleosintesi, Princeton (1996). – Dutertre, X. – Teoria delle api™: Modellazione della gravità basata sulle onde, v2, BeeTheory.com (2023).
BeeTheory.com – Gravità quantistica basata sulle onde – Metodologia galattica – © Technoplane S.A.S. 2026