Numerieke simulatie van het BeeTheory Donkere Massa Model

Een volledige, reproduceerbare beschrijving van de numerieke integratie, baryonische snelheidsdecompositie, χ²-minimalisatie en implementatiekeuzes achter de BeeTheory-simulatie van donkere massa.

Deze technische pagina legt uit hoe u de BeeTheory numerieke simulatie van de verborgen massa van de Melkweg kunt reproduceren. Het beschrijft de observatiegegevens, het baryonische model, de op golven gebaseerde dichtheidsvergelijking, de numerieke integratie, de aanpasmethode, de convergentietests en de referentiecode.

Het doel is eenvoudig: ga uit van een zichtbare Melkwegschijf, pas het BeeTheory-golfdichtheidsmodel toe, bereken de effectieve verborgen massa en vergelijk de resulterende cirkelsnelheidscurve met waarnemingen uit het Gaia-tijdperk.

0. Wat we aan het computeren zijn en waarom

De rotatiecurve van de Melkweg is de cirkelsnelheid Vc(R) van sterren als functie van hun afstand R van het galactisch centrum. Deze wordt tegenwoordig veel nauwkeuriger gemeten dan de totale massaverdeling die we direct kunnen zien.

Het verschil tussen de waargenomen snelheid en wat zichtbare baryonische materie voorspelt, is het verborgen-massa-probleem. Standaardmodellen beroepen zich op een onzichtbare deeltjeshalo, meestal koude donkere materie. De BeeTheory stelt een andere interpretatie voor: elk zichtbaar massa-element straalt een golfveld uit dat exponentieel vervalt in drie dimensies, en het geaccumuleerde veld gedraagt zich als verborgen massa.

De simulatie doet drie dingen:

1. Berekent Vcbar(R) van de zichtbare schijf en uitstulping met behulp van de analytische schijfformule van Freeman.
2. Integreert numeriek de BeeTheory-dichtheid _ρdark(r; ℓ, λ) bij elke straal en zet deze vervolgens om in ^VcDM(R).
3. Minimaliseert χ²(ℓ, λ) tegen de Gaia-tijdperk Melkwegrotatiecurve om representatieve best passende parameters te vinden.

De simulatie is ontworpen om reproduceerbaar te zijn. Het kan worden uitgevoerd in JavaScript of Python zonder een gespecialiseerde astrofysische bibliotheek.

Overal gebruikte notatie

• R is de cilindrische galactocentrische straal in het schijfvlak, in kpc.
• r is de sferische galactocentrische straal.
• z is de hoogte boven de schijf.
• ℓ is de BeeTheory coherentielengte, in kpc.
• λ de dimensieloze golf-massakoppeling is.
• G wordt gebruikt in eenheden van kpc km² s-² M⊙-¹.

1. Waarnemingsgegevens – Gaia 2024

De dataset is gebaseerd op de Gaia-tijdperk Melkwegrotatiecurve. Er worden stralen R, cirkelsnelheden _Vc en onzekerheden σ gebruikt. De oorspronkelijke technische versie gebruikte 16 datapunten voor R = 4-27,3 kpc.

De belangrijke waarnemingsfeiten zijn:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, het ankerpunt nabij de baan van de Zon.
• _Vc is ongeveer vlak van ongeveer 5 tot 20 kpc.
• _Vc neemt af in de buitenste gemeten schijf en bereikt ongeveer 173 ± 17 km/s bij 27,3 kpc in de gegevens waarnaar verwezen wordt.

Dit vereist een effectieve verborgen massaverdeling die sterk toeneemt bij tussenliggende stralen, en vervolgens afzwakt bij grotere stralen.

We sluiten het binnenste Melkwegstelsel uit omdat de centrale balk en niet-cirkelvormige bewegingen daar domineren. Een vereenvoudigd axisymmetrisch model is niet betrouwbaar in dat gebied.

3. De BeeTheory Donkere Dichtheidsvergelijking

3.1 Fysisch postulaat

Elk zichtbaar massa-element dV op positie r′ in de galactische schijf genereert een gravitatiegolfveld dat een effectieve massadichtheid creëert in een veldpunt r:

De totale donkere dichtheid op elk punt (R,z) is de superpositie over alle schijfelementen. Aangezien de bron een schijf is, wordt de volumedichtheidsterm een oppervlaktedichtheidsterm:

De volledige dubbele integraal is:

3.3 Analytische verificatie van de asymptotische limiet

Voor _Rd ≪ r ≪ ℓ wordt de uitdrukking eenvoudiger. De schijf is veel kleiner dan de straal, en de coherentielengte is nog steeds veel groter dan de straal.

Omdat:

Deze asymptotische dichtheid gedraagt zich als ρ ∝ r-², wat M( \(\rho(r)\propto r^{-2} \Longrightarrow M(<r)\propto r \Longrightarrow V_c=\sqrt{{frac{GM(<r)}{r}}}approx{mathrm{constant}\).

Door representatieve waarden in te voeren, krijgt u een lokale dichtheid die dicht bij de waargenomen Melkwegwaarde ligt:

5. χ²-minimalisatie en parametrering

5.1 Doelstelling

De fit minimaliseert de gereduceerde chi-kwadraat:

met N = 16 gegevenspunten en p = 2 vrije parameters, ℓ en λ. Dit geeft 14 vrijheidsgraden.

5.2 Zoeken in een raster in twee stappen

Er wordt een rasterzoekactie gebruikt in plaats van een gradiëntafdaling omdat het landschap een lange, gebogen ontaardingsvallei heeft tussen λ en ℓ.

• Stap 1: Grof rooster over ℓ en λ.
• Stap 2: lokale verfijning rond het grove minimum.

Het representatieve best passende gebied is:

5.3 De Degeneratievallei

Het χ²-landschap is geen symmetrische kom. Het vormt een langgerekte vallei omdat de leidende dichtheidsnormalisatie sterk afhangt van de koppelingssterkte en slechts zwak van de coherentielengte binnen het vlakke-rotatieregime.

De buitenste daling van de rotatiecurve doorbreekt deze ontaarding, omdat kleinere ℓ-waarden de effectieve dichtheid eerder onderdrukken.

Gegevens voorbij 30 kpc, waaronder bolvormige sterrenhopen, stellaire stromen, halosterren en satellietstelsels, zijn essentieel om ℓ nauwkeuriger te kunnen bepalen.

7. Volledige referentiecode

De volgende JavaScript referentie-implementatie reproduceert de belangrijkste simulatielogica. Het is bedoeld voor technische validatie en moet in een eigen scriptomgeving worden geplaatst, niet direct in een standaard WordPress inhoudsblok, tenzij aangepaste scripts zijn toegestaan.

// Fysische constanten
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;

// Gaia-tijdperk rotatiegegevens
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,7,8,9,10,11,13,17];

// Plaatshouder voor baryonische snelheid
functie vBaryonic(R) {
  // In de volledige implementatie gebruikt dit de schijfformule van Freeman
  // met gemodificeerde Bessel-functies I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  return Math.sqrt(Math.max(0, vb2));
}

// BeeTheory donkere dichtheid
functie rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30,0 / N;
  laat som = 0;

  for (laat 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));

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

  return (lam / ell) * som;
}

// Ingesloten donkere massa
functie ingeslotenDonkereMassa(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  laat M = 0;

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

  return M;
}

// Donkere en totale snelheid
functie vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * ingeslotenDonkereMassa(R, ell, lam) / R));
}

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

// Chi-kwadraat
functie chiSquared(ell, lam) {
  laat s = 0;

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

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

Een complete versie zou nauwkeurige implementaties van de gemodificeerde Bessel-functies _I0, _I1, _K0 en _K1 voor de Freeman-schijfsnelheid moeten bevatten.

Referenties

• Abramowitz, M., Stegun, I. A. - Handboek van wiskundige functies, Dover, 1972.
• Bovy, J., Rix, H.-W. - A Direct Dynamical Measurement of the Disk Surface Density Profile of the Milky Way, ApJ 779, 115, 2013.
• Freeman, K. C. - On the disks of spiral and S0 galaxies, ApJ 160, 811, 1970.
• McMillan, P. J. - The mass distribution and gravitational potential of the Milky Way, MNRAS 465, 76, 2017.
• Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693, 2024.
• Pato, M., Iocco, F., Bertone, G. - Dynamische beperkingen van de verdeling van donkere materie in de Melkweg, JCAP 12, 001, 2015.
• Portail, M. et al. - Dynamische modellering van de galactische uitstulping en balk, MNRAS 465, 1621, 2017.