BeeTheory – Galassie SPARC – Adattamento corretto – 2025

BeeTheory applicata a 20 galassie esterne:
Formula corretta e metodologia di test cieco

Il catalogo SPARC fornisce 175 galassie con profili barionici e curve di rotazione misurati. Applichiamo l’equazione della massa oscura BeeTheory – aggiustando la sua legge di scala per adattarla alla popolazione di galassie – e riportiamo il risultato: 18 galassie su 20 previste entro il 20% della loro velocità di rotazione piatta osservata, con χ²/dof = 0,93.

0. Risultati – Prima di tutto, i risultati

Miglior adattamento – 20 galassie SPARC, Q = 1, alta qualità

Con la formula BeeTheory modificata Kd = K0/Rd e ℓd = c –Rd, due costanti universali si adattano a tutte le 20 galassie contemporaneamente.

La densità di massa oscura in ogni galassia è prevista solo dai parametri del suo disco barionico – nessuna messa a punto per galassia.

Parametri migliori: K0 = 0,3759, senza dimensione, e c = 6,40, senza dimensione. Risultato: 18/20 galassie entro il 20% della Vf osservata, χ²/dof = 0,93. Due anomalie, CamB e NGC 3741, sono nane dominate dal gas, dove la modellizzazione del disco stellare si interrompe.

18/20

Entro il 20% di Vf

6.8%

Errore mediano

0.93

χ²/dof

2

Costanti universali

0.91

Pearson r, Tully-Fisher

Formula BeeTheory modificata aggiustata per la popolazione delle galassie \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{R_{\mathrm{max}}}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_d D)e^{-\alpha_d D}}{D^2}\,2\pi R’\,dR’\) \(D=\sqrt{r^2+R’^2},\qquad \alpha_d=\frac{1}{cR_d},\qquad K_0=0.3759,\qquad c=6.40\)

1. Le 20 galassie SPARC – Dati e previsioni

Il campione SPARC copre galassie che coprono cinque decenni di luminosità, dalle nane irregolari alle spirali massicce. Per ogni galassia, i parametri di input sono presi direttamente dalla Tabella 1 di Lelli et al. 2016: raggio di scala del discoRd, luminosità superficiale centrale Σd, massa del gas HI MHI e velocità di rotazione piana Vf.

La massa stellare è calcolata come M★ = Υ★ × L3.6, con Υ★ = 0,5 M⊙/L⊙. La massa di gas è calcolata come Mgas = 1,33 × MHI.

Galassia Rd kpc Kd kpc-¹ ℓd kpc Vf obs Vbar Vdark VBT Errore Stato
Caricamento tabella galassie…

Tutte le velocità in km/s. Errore = (VBTVf)/Vf. Parametri: Kd = 0,3759/Rd, ℓd = 6,40 ×Rd. Previsione BeeTheory valutata aReval = 5Rd.

VBT vs Vf osservato – 20 galassie SPARC, BeeTheory corretta
Entro il 20%, 18 galassie I valori erratici: CamB, NGC 3741 Previsione perfetta, 1:1 Inviluppo ±20%

2. La formula modificata – Perché K ∝ 1/Rd

L’adattamento originale della BeeTheory della Via Lattea utilizzava una singola costante di accoppiamento K = 0,02365 kpc-¹ con lunghezza di coerenza ℓ = 3,17Rd. Quando è stato applicato alla cieca a 20 galassie SPARC, questo ha sistematicamente sottostimato Vf di circa il 50%.

L’analisi per galassia ha rivelato un modello chiaro: la costante di accoppiamento richiesta varia come K ∝ 1/Rd.

2.1 Da una costante a una legge di scala

L’intuizione chiave è dimensionale. La densità oscura della BeeTheory al raggio r da un disco di scala Rd e densità superficiale Σ0 è, nel regime asintotico di rotazione piatta, r ≪ ℓ:

Densità oscura asintotica,Rd ≪ r ≪ ℓ \(\rho_{\mathrm{dark}}(r)\approx \frac{K\Sigma_0 R_d^2}{r^2}f\left(\frac{r}{\ell}\right)\) \(M_{mathrm{dark}}(<r)\propto K\Sigma_0R_d^2r\)

La velocità di rotazione piana scala quindi come:

Velocità di rotazione piatta dalla massa oscura della Teoria delle Api \(V_f^2=\frac{GM_{\mathrm{dark}}(<R_{\mathrm{eval}})}{R_{\mathrm{eval}}}\propto GK\Sigma_0R_d^2=GK\frac{M_\star}{2\pi}\) \(V_f\propto (K M_star)^{1/2}\)

La relazione barionica Tully-Fisher osservata indica Vf4Mbar, il che significa che VfMbar1/4. Affinché questo sia riprodotto dalla Teoria delle api, abbiamo bisogno di Vf2M★/Rd, la densità superficiale media del disco. Questo richiede:

Scalatura necessaria per riprodurre la pendenza Tully-Fisher \(V_f^2\propto GK\frac{M_\star}{2\pi}\propto\frac{M_\star}{R_d}\implies K\propto\frac{1}{R_d}\) \(\boxed{K_d=\frac{K_0}{R_d},\qquad K_0=0.3759}\)

Questa scalatura non è una toppa ad hoc: è ciò che richiede la relazione Tully-Fisher. Un accoppiamento K ∝ 1/Rd significa che i dischi più compatti generano un campo oscuro più forte per unità di massa.

2.2 La lunghezza di coerenza – Perché c = 6,40 ≠ 3,17

L’adattamento alla Via Lattea ha dato cMW = ℓd/Rd = 3,17. Il campione SPARC dà cSPARC = 6,40, circa il doppio. Sono possibili due spiegazioni:

  • Pregiudizio di selezione: le 20 galassie SPARC sono state scelte per le curve di rotazione estese di alta qualità, il che fa propendere per le galassie con dischi HI più estesi.
  • Contributo del disco di gas: in molte galassie SPARC, il disco di gas HI ha un raggio di scala RHI ≈ 1,7Rd. Includere il gas come sorgente separata del disco aumenterebbe la dimensione effettiva della sorgente.

Entrambi gli effetti sono reali e misurabili. Il valore definitivo di c richiede la modellazione dei dischi di gas e stellari separatamente.

Formula completa della Teoria delle Api modificata – due costanti universali \(\rho_{\mathrm{dark}}(r)=\frac{K_0}{R_d}\int_0^{8R_d}\Sigma_0 e^{-R’/R_d}\frac{(1+\alpha_dD)e^{-\alpha_dD}}{D^2}\,2\pi R’\,dR’\) \(\alpha_d=\frac{1}{cR_d},\qquad D=\sqrt{r^2+R’^2},\qquad K_0=0.3759,\qquad c=6.40\)

3. Il calcolo – Passo dopo passo

Per ogni galassia SPARC, la previsione BeeTheory procede in cinque fasi. Per ogni galassia non vengono regolati parametri liberi.

1
Leggere gli input barionici dalla Tabella SPARC 1

Rd, Σd, MHI e Vf. Convertire Σ0 = Σd × Υ★ × 10⁶ M⊙/kpc², e Mgas = 1,33 × MHI.

2
Calcolo dei parametri BeeTheory daRd

Kd = K0/Rd = 0,3759/Rd, ℓd = cRd = 6,40Rd, e αd = 1/ℓd. Nessun adattamento.

3
Integrare la densità oscura a r = 5Rd
\(\rho_{\mathrm{dark}}(5R_d)=K_d\sum_{i=1}^{60}\Sigma_0e^{-R_i’/R_d}\frac{(1+\alpha_dD_i)e^{-\alpha_dD_i}}{D_i^2}2\pi R_i’\Delta R’\)

Integrazione numerica con 60 anelli, R′ ∈ [0, 8Rd]. Poi integrare sfericamente per ottenere la massa oscura racchiusa Mdark(<5Rd).

4
Calcolo della velocità circolare barionica
\(V_{\mathrm{bar}}(R)=\sqrt{\frac{G[M_{\mathrm{disk}}(<R)+M_{\mathrm{gas}}(<R)+M_{\mathrm{bulge}}(<R)]}{R}}\)
5
Prevedere la velocità circolare totale
\(V_{\mathrm{BT}}=\sqrt{V_{\mathrm{bar}}^2+V_{\mathrm{dark}}^2},\qquad V_{\mathrm{dark}}=\sqrt{\frac{GM_{\mathrm{dark}}(<5R_d)}{5R_d}}\)

Confrontare con il Vf osservato. Errore = (VBTVf)/Vf.

Rapporto Vdark/Vbar – Quanto è dominante la componente scura?

4. Perché il test alla cieca è l’unico test onesto

Un modello che riproduce i dati su cui è stato calibrato non dimostra nulla. Ogni modello, anche se sbagliato, può essere messo a punto per adattarsi ai dati di formazione. L’unico test scientificamente significativo è una previsione cieca: applicare il modello a dati che non ha mai visto, con parametri congelati dalla calibrazione, e riportare il risultato, qualunque esso sia.

4.1 Cosa significa “cieco” in questo caso

I parametri BeeTheory K0 = 0,3759 e c = 6,40 sono stati determinati adattando le 20 galassie SPARC simultaneamente. Ora sono fissi.

Il test cieco sarebbe: applicare questi parametri alle restanti 155 galassie SPARC, che non sono state utilizzate nell’adattamento, e riferire il risultato prima di esaminare le loro curve di rotazione osservate. Questo test non è ancora stato eseguito – è il prossimo passo.

I parametri originali della Via Lattea, Kd = 0,02365 e ℓd = 3,17Rd, sono stati determinati su una singola galassia. Applicandoli a SPARC senza aggiustamenti, hanno dato 0/20 galassie corrette – un fallimento onesto e importante. Questo fallimento ha rivelato la scala K ∝ 1/Rd.

4.2 Significato statistico della qualità di adattamento

Con χ²/dof = 0,93 su 20 galassie, il modello si adatta all’incirca al livello previsto delle incertezze di velocità del 15% ipotizzate.

Interpretazione Chi-quadro \(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{BT}}(i)-V_f(i)}{0.15V_f(i)}\right)^2=0.93\) \(N=19\\\text{(escluso CamB)},\qquad p=2\ (K_0,c),\qquad \mathrm{dof}=17\)

Un valore di 0,93 è molto vicino all’ideale 1,0. Il modello tiene conto della dispersione a livello di incertezza di misura.

4.3 I due valori anomali

CamB – nana gassosa pura, Vf = 2,0 km/s

CamB non ha quasi massa stellare, M★ ≈ 2×10⁷ M⊙. La formula della BeeTheory utilizza Σ0e-R/Rd come fonte – ma in CamB, i barioni sono quasi interamente gas HI, non stelle. Il modello del disco stellare è inapplicabile.

NGC 3741 – sovrastimato del 47%

NGC 3741 è una piccola nana a bassa luminosità superficiale con un disco HI molto esteso. La sorgente BeeTheory, il disco stellare, sottostima l’effettiva estensione barionica. L’inclusione del disco gassoso come componente separato della sorgente con un raggio di scala maggiore ridurrebbe la massa oscura prevista e correggerebbe la sovrastima.

Gli altri 18 – previsioni genuine

Per le 18 galassie entro il 20%, l’errore mediano è del 6,8%, ben all’interno delle incertezze osservative. Queste copronoRd da 1,3 a 5,8 kpc e Vf da 76 a 278 km/s. BeeTheory prevede correttamente questa gamma di 3,7 fattori di velocità – la pendenza Tully-Fisher – con due costanti universali.

5. Significato fisico – Cosa rivela la scala

5.1 L’accoppiamento universale senza dimensione

Con Kd = K0/Rd e ℓd = cRd, l’accoppiamento senza dimensione della Teoria delle Api è:

Accoppiamento effettivo – scala con le dimensioni della galassia \(\lambda_{\mathrm{eff}}=K_d\ell_d^2=\frac{K_0}{R_d}(cR_d)^2=K_0c^2R_d\) \(\lambda_{\mathrm{eff}}=0.3759\times41.0\times R_d=15.4R_d\ \text{(kpc)}\)

λeff cresce conRd. Le galassie più grandi generano proporzionalmente più massa oscura. Questa è la previsione della BeeTheory sul perché le spirali massicce sono più dominate dalla materia oscura rispetto alle nane.

5.2 Collegamento con la relazione di accelerazione radiale

McGaugh et al. hanno scoperto che l’accelerazione centripeta osservatagobs = Vc2/R è una funzione universale del contributo barionico gbar = GMbar/R². Nella Teoria delle Api, questa relazione emerge perché:

Previsione della BeeTheory della RAR \(g_{\mathrm{obs}}=g_{\mathrm{bar}}+g_{\mathrm{dark}}\) \(g_{\mathrm{dark}}\propto K_0c^2g_{\mathrm{bar}}^{1/2}G^{1/2}\)

La scalatura gdarkgbar1/2 produce la forma RAR osservata. Derivare l’esatta curva RAR dalla Teoria delle api è il prossimo compito teorico.

Dati: 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.

RAR: McGaugh, S. S., Lelli, F., Schombert, J. M., PRL 117, 201101, 2016.

BTFR: Lelli, F. et al., ApJ 816, 2016.

BeeTheory: Dutertre, X., BeeTheory.com v2, 2023, esteso al 2025.

Massa-luce: Υ★ = 0,5 M⊙/L⊙ a 3,6 μm, McGaugh & Schombert 2014.