Numerisk simulering af BeeTheory-modellen for mørk masse

En komplet, reproducerbar redegørelse for den numeriske integration, nedbrydning af baryonisk hastighed, χ²-minimering og implementeringsvalg bag BeeTheory-simuleringen af mørk masse.

Denne tekniske side forklarer, hvordan man reproducerer BeeTheorys numeriske simulering af Mælkevejens skjulte masse. Den beskriver observationsdataene, den baryoniske model, den bølgebaserede tæthedsligning, den numeriske integration, tilpasningsmetoden, konvergenstests og referencekoden.

Målet er enkelt: Tag udgangspunkt i en synlig Mælkevejsskive, anvend BeeTheorys bølgetæthedsmodel, beregn den effektive skjulte masse, og sammenlign den resulterende cirkulære hastighedskurve med observationer fra Gaia-æraen.

0. Hvad vi beregner og hvorfor

Mælkevejens rotationskurve er stjernernes cirkulære hastighed Vc(R) som en funktion af deres afstand R fra det galaktiske centrum. Den måles i dag langt mere præcist end den samlede massefordeling, vi direkte kan se.

Underskuddet mellem den observerede hastighed og det, som synligt baryonisk stof forudsiger, er problemet med den skjulte masse. Standardmodeller påberåber sig en usynlig partikelhalo, som regel koldt mørkt stof. BeeTheory foreslår en anden fortolkning: Hvert synligt masseelement udstråler et bølgefelt, der henfalder eksponentielt i tre dimensioner, og det akkumulerede felt opfører sig som skjult masse.

Simuleringen gør tre ting:

1. Beregner Vcbar(R) fra den synlige skive og bulge ved hjælp af Freemans analytiske skiveformel.
2. Integrerer numerisk BeeTheory-tætheden _ρdark(r; ℓ, λ) ved hver radius og konverterer den derefter til ^VcDM(R).
3. Minimerer χ²(ℓ, λ) mod Mælkevejens rotationskurve fra Gaia-æraen for at finde repræsentative best-fit-parametre.

Simuleringen er designet til at være reproducerbar. Den kan køre i JavaScript eller Python uden et specialiseret astrofysisk bibliotek.

Notation brugt hele vejen igennem

• R er den cylindriske galaktocentriske radius i skivens plan, i kpc.
• r er den sfæriske galaktocentriske radius.
• z er højden over disken.
• ℓ er BeeTheory-kohærenslængden i kpc.
• λ er den dimensionsløse bølge-masse-kobling.
• G bruges i enheder af kpc km² s-² M⊙-¹.

1. Observationsdata – Gaia 2024

Datasættet er baseret på Mælkevejens rotationskurve fra Gaia-æraen. Det bruger radier R, cirkulære hastigheder _Vc og usikkerheder σ. Den oprindelige tekniske version brugte 16 datapunkter, der dækkede R = 4-27,3 kpc.

De vigtige observationelle fakta er:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, ankerpunktet nær solens bane.
• _Vc er omtrent flad fra ca. 5 til 20 kpc.
• _Vc falder i den ydre målte skive og når ca. 173 ± 17 km/s nær 27,3 kpc i de refererede data.

Det kræver en effektiv skjult massefordeling, som stiger kraftigt ved mellemliggende radier og derefter aftager ved større radier.

Vi udelukker den inderste galakse, fordi den centrale bjælke og de ikke-cirkulære bevægelser dominerer der. En forenklet aksesymmetrisk model er ikke pålidelig i det område.

3. Bi-teoriens ligning for mørk tæthed

3.1 Fysisk postulat

Hvert synligt masseelement dV på position r′ i den galaktiske skive genererer et tyngdebølgefelt, der skaber en effektiv massetæthed i et feltpunkt r:

Den samlede mørke tæthed i ethvert punkt (R,z) er superpositionen over alle diskens elementer. Da kilden er en skive, bliver udtrykket for volumendensitet til et udtryk for overfladetæthed:

Det fulde dobbelte integral er:

3.3 Analytisk verificering af den asymptotiske grænse

For_Rd ≪ r ≪ ℓ forenkles udtrykket. Disken er meget mindre end radius, og kohærenslængden er stadig meget større end radius.

Fordi:

Denne asymptotiske tæthed opfører sig som ρ ∝ r-², hvilket giver M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}}\approx\mathrm{konstant}\).

Ved at indsætte repræsentative værdier får man en lokal tæthed tæt på den observerede værdi for Mælkevejen:

5. χ²-minimering og parametertilpasning

5.1 Målfunktion

Tilpasningen minimerer den reducerede chi-kvadrat:

med N = 16 datapunkter og p = 2 frie parametre, ℓ og λ. Det giver 14 frihedsgrader.

5.2 Gittersøgning med to gennemløb

Der bruges en gittersøgning i stedet for gradientnedstigning, fordi landskabet har en lang, buet degenerationsdal mellem λ og ℓ.

• Pass 1: groft gitter over ℓ og λ.
• Pass 2: lokal forfining omkring det grove minimum.

Den repræsentative best-fit region er:

5.3 Degenerationsdalen

χ²-landskabet er ikke en symmetrisk skål. Det danner en langstrakt dal, fordi den førende tæthedsnormalisering afhænger stærkt af koblingsstyrken og kun svagt af kohærenslængden inden for det flade rotationsregime.

Det ydre fald i rotationskurven bryder denne degeneration, fordi mindre ℓ-værdier undertrykker den effektive tæthed tidligere.

Data ud over 30 kpc, herunder kuglehobe, stjernestrømme, halostjerner og satellitgalakser, er afgørende for at begrænse ℓ mere præcist.

7. Fuld referencekode

Følgende JavaScript-referenceimplementering gengiver den vigtigste simuleringslogik. Den er beregnet til teknisk validering og bør placeres i et ordentligt scriptmiljø, ikke direkte i en standard WordPress-indholdsblok, medmindre brugerdefinerede scripts er tilladt.

// Fysiske konstanter
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;

// Rotationsdata fra Gaia-æraen
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];

// Baryonisk hastigheds pladsholder
funktion vBaryonic(R) {
  // I den fulde implementering bruger dette Freemans diskformel
  // med modificerede Bessel-funktioner I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  return Math.sqrt(Math.max(0, vb2));
}

// BeeTheory mørk tæthed
funktion rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30,0 / N;
  lad sum = 0;

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

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

  return (lam / ell) * sum;
}

// Indesluttet mørk masse
funktion lukketMørkMasse(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  lad M = 0;

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

  return M;
}

// Mørk og samlet hastighed
funktion vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

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

// Chi-kvadrat
function chiSquared(ell, lam) {
  lad s = 0;

  for (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;
  }

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

En komplet version bør indeholde nøjagtige implementeringer af de modificerede Bessel-funktioner _I0, _I1, _K0 og_K1 for Freeman-skivens hastighed.

Referencer

• Abramowitz, M., Stegun, I. A. - Handbook of Mathematical Functions, Dover, 1972.
• Bovy, J., Rix, H.-W. - A Direct Dynamical Measurement of the Milky Way's Disk Surface Density Profile, 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. - Dynamiske begrænsninger på fordelingen af mørkt stof i Mælkevejen, JCAP 12, 001, 2015.
• Portail, M. et al - Dynamical modelling of the Galactic bulge and bar, MNRAS 465, 1621, 2017.