Simulazione numerica del modello di massa oscura della Teoria delle Api

Un resoconto completo e riproducibile dell’integrazione numerica, della decomposizione della velocità barionica, della minimizzazione del χ² e delle scelte di implementazione alla base della simulazione della massa oscura BeeTheory.

Questa pagina tecnica spiega come riprodurre la simulazione numerica BeeTheory della massa nascosta della Via Lattea. Descrive i dati osservativi, il modello barionico, l’equazione di densità basata sulle onde, l’integrazione numerica, il metodo di adattamento, i test di convergenza e il codice di riferimento.

L’obiettivo è semplice: partire da un disco visibile della Via Lattea, applicare il modello di densità d’onda BeeTheory, calcolare la massa nascosta effettiva e confrontare la curva di velocità circolare risultante con le osservazioni dell’era Gaia.

0. Cosa stiamo elaborando e perché

La curva di rotazione della Via Lattea è la velocità circolare Vc(R) delle stelle in funzione della loro distanza R dal Centro Galattico. Oggi viene misurata in modo molto più preciso rispetto alla distribuzione della massa totale che possiamo vedere direttamente.

Il deficit tra la velocità osservata e ciò che la materia barionica visibile prevede è il problema della massa nascosta. I modelli standard invocano un alone di particelle invisibili, di solito la materia oscura fredda. La Teoria delle api propone un’interpretazione diversa: ogni elemento di massa visibile irradia un campo d’onda che decade esponenzialmente in tre dimensioni, e il campo accumulato si comporta come una massa nascosta.

La simulazione fa tre cose:

1. Calcola Vcbar(R) dal disco visibile e dal bulge utilizzando la formula analitica del disco di Freeman.
2. Integra numericamente la densità BeeTheory _ρdark(r; ℓ, λ) ad ogni raggio, poi la converte in ^VcDM(R).
3. Minimizza χ²(ℓ, λ) rispetto alla curva di rotazione della Via Lattea dell’era Gaia per trovare i parametri rappresentativi di miglior adattamento.

La simulazione è progettata per essere riproducibile. Può essere eseguita in JavaScript o Python senza una libreria astrofisica specializzata.

Notazione utilizzata in tutto il testo

• R è il raggio cilindrico galattocentrico nel piano del disco, in kpc.
• r è il raggio sferico galattocentrico.
• z è l’altezza sopra il disco.
• ℓ è la lunghezza di coerenza BeeTheory, in kpc.
• λ è l’accoppiamento onda-massa senza dimensione.
• G è utilizzato in unità di kpc km² s-² M⊙-¹.

1. Dati osservativi – Gaia 2024

Il set di dati si basa sulla curva di rotazione della Via Lattea dell’era Gaia. Utilizza raggi R, velocità circolari _Vc e incertezze σ. La versione tecnica originale utilizzava 16 punti dati che coprivano R = 4-27,3 kpc.

I fatti osservativi importanti sono:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, il punto di ancoraggio vicino all’orbita del Sole.
• _Vc è approssimativamente piatto da circa 5 a 20 kpc.
• _Vc diminuisce nel disco esterno misurato, raggiungendo circa 173 ± 17 km/s vicino a 27,3 kpc nei dati di riferimento.

Ciò richiede una distribuzione della massa nascosta effettiva che aumenta fortemente a raggi intermedi, per poi indebolirsi a raggi maggiori.

Escludiamo la Galassia più interna perché lì dominano la barra centrale e i moti non circolari. Un modello assialsimmetrico semplificato non è affidabile all’interno di questa regione.

3. L’equazione della densità oscura della Teoria delle Api

3.1 Postulato fisico

Ogni elemento di massa visibile dV in posizione r′ nel disco galattico genera un campo di onde gravitazionali che crea una densità di massa effettiva in un punto di campo r:

La densità oscura totale in qualsiasi punto (R,z) è la sovrapposizione di tutti gli elementi del disco. Poiché la sorgente è un disco, il termine di densità di volume diventa un termine di densità di superficie:

La doppia integrale completa è:

3.3 Verifica analitica del limite asintotico

Per_Rd ≪ r ≪ ℓ, l’espressione si semplifica. Il disco è molto più piccolo del raggio e la lunghezza di coerenza è ancora molto più grande del raggio.

Perché:

Questa densità asintotica si comporta come ρ ∝ r-², che dà M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}}\approx\mathrm{constant}\)

Inserendo i valori rappresentativi si ottiene una densità locale vicina al valore osservato della Via Lattea:

5. Minimizzazione del χ² e adattamento dei parametri

5.1 Funzione obiettivo

L’adattamento minimizza il chi-quadro ridotto:

con N = 16 punti dati e p = 2 parametri liberi, ℓ e λ. Questo dà 14 gradi di libertà.

5.2 Ricerca della griglia a due passaggi

Viene utilizzata una ricerca a griglia piuttosto che una discesa per gradi, perché il paesaggio presenta una lunga valle di degenerazione curva tra λ e ℓ.

• Passaggio 1: griglia grossolana su ℓ e λ.
• Passaggio 2: affinamento locale intorno al minimo grossolano.

La regione rappresentativa che si adatta meglio è:

5.3 La valle della degenerazione

Il paesaggio χ² non è una coppa simmetrica. Forma una valle allungata perché la normalizzazione della densità principale dipende fortemente dalla forza di accoppiamento e solo debolmente dalla lunghezza di coerenza all’interno del regime di rotazione piatta.

Il declino esterno della curva di rotazione rompe questa degenerazione, perché i valori più piccoli di ℓ sopprimono prima la densità effettiva.

I dati oltre i 30 kpc, compresi gli ammassi globulari, i flussi stellari, le stelle dell’alone e le galassie satellite, sono essenziali per vincolare ℓ in modo più stretto.

7. Codice di riferimento completo

La seguente implementazione di riferimento JavaScript riproduce la logica di simulazione principale. È destinata alla convalida tecnica e deve essere inserita in un ambiente di script appropriato, non direttamente all’interno di un blocco di contenuti standard di WordPress, a meno che non siano consentiti script personalizzati.

// Costanti fisiche
const G = 4.302e-3; // kpc km² s-² M☉-¹
const Sig0 = 800,0; // M☉ pc-²
const Rd = 2,6; // kpc
const Mdisk = 3,5e10; // M☉
const Mbulge = 1,2e10; // M☉
const CONV_RHO = (1,989e30 * 5,61e26) / (3,086e21)**3;

// Dati di rotazione dell'era Gaia
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];

// Segnaposto velocità barionica
funzione vBaryonic(R) {
  // Nell'implementazione completa, questa utilizza la formula del disco di Freeman
  // con le funzioni di Bessel modificate I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  restituisce Math.sqrt(Math.max(0, vb2));
}

// Densità scura di BeeTheory
funzione rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30.0 / N;
  lasciare che sum = 0;

  per (let i = 0; i < N; i++) {
    const Rp = (i + 0,5) * dRp;
    const Sig = Sig0 * Math.exp(-Rp / Rd);
    const K = (2 * Math.PI * ell / r)
      * Math.sinh(Math.min(r / ell, 30))
      * Math.exp(-Math.min((r + Rp) / ell, 30));

    somma += Sig * Rp * K * dRp;
  }

  restituisce (lam / ell) * somma;
}

// Massa oscura racchiusa
funzione enclosedDarkMass(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  lasciare M = 0;

  per (let j = 0; j < Nr; j++) {
    const rj = (j + 0,5) * dr;
    M += 4 * Math.PI * rj * rj * rhoDark(rj, ell, lam) * dr;
  }

  restituire M;
}

// Velocità oscura e totale
funzione vDM(R, ell, lam) {
  restituisce Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

funzione vTotal(R, ell, lam) {
  const vb = vBaryonic(R);
  vd = vDM(R, ell, lam);
  restituisce Math.sqrt(vb * vb + vd * vd);
}

// Chi-quadrato
funzione chiSquared(ell, lam) {
  lasciare che s = 0;

  per (let i = 0; i < OBS_R.length; i++) {
    const dv = (vTotal(OBS_R[i], ell, lam) - OBS_V[i]) / OBS_ERR[i];
    s += dv * dv;
  }

  restituisce s / (OBS_R.length - 2);
}

Una versione completa dovrebbe includere implementazioni accurate delle funzioni di Bessel modificate _I0, _I1, _K0 e _K1 per la velocità del disco di Freeman.

Riferimenti

• Abramowitz, M., Stegun, I. A. - Manuale delle funzioni matematiche, Dover, 1972.
• Bovy, J., Rix, H.-W. - Una misurazione dinamica diretta del profilo di densità superficiale del disco della Via Lattea, ApJ 779, 115, 2013.
• Freeman, K. C. - Sui dischi delle galassie a spirale e S0, ApJ 160, 811, 1970.
• McMillan, P. J. - La distribuzione di massa e il potenziale gravitazionale della Via Lattea, MNRAS 465, 76, 2017.
• 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, 2024.
• Pato, M., Iocco, F., Bertone, G. - Vincoli dinamici sulla distribuzione della materia oscura nella Via Lattea, JCAP 12, 001, 2015.
• Portail, M. et al. - Modellazione dinamica del bulge e della barra galattica, MNRAS 465, 1621, 2017.