BeeTheory – Galactic Simulation v2 – første generation 2025 maj 17 med claude

Mælkevejens skjulte masse: BeeTheory 3D Yukawa med fysisk diskafkortning

Den korrigerede simulering: Den baryoniske skives hastighed falder keplersk ud over dens fysiske kant, og BeeTheory 3D Yukawa-kernen fylder hele rummet. To parametre, Gaia-æraens rotationsdata og en afkortet diskmodel.

BeeTheory.com – Ou et al., MNRAS 528, 2024 – Korrigeret BeeTheory v2

0. Resultat – ligninger og parametre

Hver ring i den galaktiske skive med radius R′ genererer et 3D effektivt mørkt massefelt gennem BeeTheory Yukawa-kernen. Den samlede mørke tæthed ved den sfæriske radius r er:

\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\) \(D(r,R’)=\sqrt{r^2+R’^2}\)

Kernen er afledt af den korrigerede BeeTheory-kraftlov:

\(F(D)\propto\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)

Den reduceres til den newtonske invers-kvadratform for D, der er meget mindre end kohærenslængden ℓ.

\(D\ll\ell=\frac{1}{\alpha}\quad\Longrightarrow\quad F(D)\propto\frac{1}{D^2}\)

Den baryoniske skives hastighed bruger Freeman-formlen inden for dens fysiske kant Rtrunc4Rd = 10,4 kpc og overgår derefter jævnt til det Keplerske fald, der forventes fra en endelig massefordeling.

\(K=0.0397\,\mathrm{kpc}^{-1},\qquad \alpha=0.0868\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.5\,\mathrm{kpc}\)

Oversigt over fit

ObserverbarGaia-æraens værdiBeeTheoryTræk
Vc(4 kpc)220 ± 10 km/s219,8 km/s-0.02σ
Vc(8 kpc)230 ± 6 km/s233,2 km/s+0.53σ
Vc(12 kpc)226 ± 7 km/s223,8 km/s-0.31σ
Vc(20 kpc)215 ± 10 km/s211,2 km/s-0.38σ
Vc(27,3 kpc)173 ± 17 km/s199,0 km/s+1.53σ
ρdark(R⊙ = 8 kpc)0,39 ± 0,03 GeV/cm³0,47 GeV/cm³+2.3σ
Mdark(<8 kpc)~5 × 10¹⁰ M⊙5.3 × 10¹⁰ M⊙tæt på
Mtot(<200 kpc)5-9 × 10¹¹ M⊙3.3 × 10¹¹ M⊙lav ende

Den forenklede tilpasning giver χ²/dof ≈ 0,31. Det sværeste punkt er stadig den yderste Gaia-værdi ved 27,3 kpc, hvor det observerede fald er skarpere, end denne model med to parametre forudsiger.

1. Diskafkortning – hvorfor og hvordan

1.1 Problemet med en uendelig eksponentiel disk

Freemans skiveformel forudsætter en eksponentiel overfladetæthed, der strækker sig ud i det uendelige. Matematisk set når den aldrig nul, men fysisk set har Mælkevejens stjerneskive en begrænset udstrækning. Ud over den effektive stjernekant er den indesluttede baryoniske masse stort set konstant, og hastighedsbidraget må falde omtrent som et Keplersk punktmassefelt.

\(\Sigma(R)=\Sigma_0e^{-R/R_d}\)

Ud over diskens kant tenderer den baryoniske hastighed mod:

\(V_{\mathrm{bar}}(R)\xrightarrow{R\gg R_d}\sqrt{\frac{GM_{\mathrm{bar,tot}}}{R}}\) \(M_{\mathrm{bar,tot}}=M_{\mathrm{disk}}+M_{\mathrm{bulge}}\approx4.7\times10^{10}M_\odot\)

Eksempler på værdier er:

\(V_{\mathrm{bar}}(30\,\mathrm{kpc})\approx82\,\mathrm{km/s},\qquad V_{\mathrm{bar}}(50\,\mathrm{kpc})\approx63\,\mathrm{km/s}\)

1.2 Formel for glat trunkering

Simuleringen bruger en jævn overgang mellem Freemans diskformel og den keplerske værdi. Overgangen er centreret ved Rtrunc = 4Rd = 10,4 kpc med bredden σ = 1,5 kpc.

\(V_{\mathrm{bar}}(R)=\sqrt{(1-w)V_{\mathrm{Freeman}}^2(R)+w\,\min(V_{\mathrm{Freeman}},V_{\mathrm{Kepler}})^2}\) \(w(R)=\frac{1}{2}\left[1+\tanh\left(\frac{R-R_{\mathrm{trunc}}}{\sigma}\right)\right]\) \(R_{\mathrm{trunc}}=4R_d=10.4\,\mathrm{kpc},\qquad \sigma=1.5\,\mathrm{kpc}\)

Minimumsfunktionen forhindrer den baryoniske skive i at overskride den fysiske keplerske grænse uden for skivens kant.

RVFreemanVKeplerianVbar,afkortetDominerende regime
5 kpc174,5 km/s201,1 km/s174,5 km/sFreeman
8 kpc161,5 km/s159,0 km/s161,5 km/sFreeman ≈ Kepler
10,4 kpc143,0 km/s139,3 km/s141,2 km/sOvergang
16 kpc112,4 km/s112,4 km/s112,4 km/sKeplerian
25 kpc89,9 km/s89,9 km/s89,9 km/sKeplerian
50 kpc63,6 km/s63,6 km/s63,6 km/sKeplerian

2. BeeTheory 3D-tæthed af mørk masse

2.1 Diskringe, der udstråler i 3D

Hver ring i den galaktiske skive med radius R′ og bredden dR′ har en masse:

\(dM=\Sigma(R’)\,2\pi R’\,dR’\)

I BeeTheory genererer denne ring et gravitationsbølgefelt, der udbreder sig i alle tre rumlige dimensioner. I monopoltilnærmelsen er afstanden til et 3D-feltpunkt ved den sfæriske radius r:

\(D(r,R’)=\sqrt{r^2+R’^2}\)

Den numeriske form af den mørke tæthed er:

\(\rho_{\mathrm{dark}}(r)=K\sum_{i=1}^{N}\Sigma_0e^{-R’_i/R_d}\frac{(1+\alpha D_i)e^{-\alpha D_i}}{D_i^2}\,2\pi R’_i\Delta R’\) \(D_i=\sqrt{r^2+R_i’^2},\qquad R’_i=\left(i-\frac{1}{2}\right)\frac{R_{\mathrm{max}}}{N}\) \(N=60,\qquad R_{\mathrm{max}}=25\,\mathrm{kpc}\)

2.2 Indesluttet mørk masse og cirkulær hastighed

\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\) \(M_{\mathrm{dark}}(<r)\approx\sum_{j=1}^{30}4\pi r_j^2\rho_{\mathrm{dark}}(r_j)\Delta r\) \(V_{\mathrm{dark}}(R)=\sqrt{\frac{GM_{\mathrm{dark}}(<R)}{R}}\) \(V_c(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{dark}}^2(R)}\)

2.3 Asymptotisk opførsel

\(\rho_{\mathrm{dark}}(r)\approx\frac{2\pi K\Sigma_0R_d^2}{r^2}\left(1+\alpha r+\frac{\alpha^2r^2}{2}\right)e^{-\alpha r}\)

For αr ≪ 1:

\(\rho_{\mathrm{dark}}(r)\xrightarrow{\alpha r\ll1}\frac{2\pi K\Sigma_0R_d^2}{r^2}\) \(M_{\mathrm{dark}}(<r)\propto r\qquad\Longrightarrow\qquad V_{\mathrm{dark}}\approx\mathrm{konstant}\)

3. Simuleringsresultater – interaktive diagrammer

Simuleringen nedenfor beholder den numeriske model, skyderne, rotationskurven, masseprofilen, tæthedsprofilen og live χ²-opdateringen. Indsæt denne side i WordPress med scriptudførelse aktiveret.

Kun baryoner, afkortet BeeTheory i alt Mørk komponent Data fra Gaia-æraen
Parameterudforsker – juster K, α og Rtrunc
0.040
0.087
10.4

χ²/dof: | ℓ: kpc | ρ(R⊙): – GeV/cm³

Masseprofil: synlig skive vs 3D mørk masse vs total
Synlig skive + udbuling BeeTheory mørk masse Samlet masse
r (kpc) Mbar (10¹⁰ M⊙) Mdark (10¹⁰ M⊙) Mtot (10¹⁰ M⊙) DM/bar ρdark (GeV/cm³)
Indlæsning…
Mørk tæthedsprofil ρdark(r) – logaritmisk skala
Bi-teori Isotermisk r-²-reference NFW-reference

4. Fysisk fortolkning og universalitet

4.1 Kohærenslængde

Inden for kohærenslængden opfører Yukawa-kernen sig næsten som en newtonsk 1/D²-kerne. Den mørke tæthed følger omtrent r-², og rotationskurven er flad. Ud over ℓ giver den eksponentielle undertrykkelse det fald, der er observeret i den ydre skive.

\(\ell=\frac{1}{\alpha}\approx11.5\,\mathrm{kpc}\) \(\frac{\ell}{R_d}=\frac{11.5}{2.6}\approx4.4\)

4.2 Dimensionsløs kobling

En dimensionsløs BeeTheory-kobling kan defineres som:

\(\lambda_{\mathrm{galaxy}}=K\ell^2\) \(\lambda_{\mathrm{galaxy}}=0.040\times(11.5)^2\approx5.3\)

Dette er sammenligneligt i størrelsesorden med den kobling, der udledes af H₂-kalibreringen, hvor λ er omkring 3-4. Den mulige skalauniversalitet af dette tal er stadig et centralt åbent spørgsmål.

4.3 Sammenligning med standardmodeller

ModelParametreTypisk pasformSkalaMekanisme
Isotermisk halo2ModeratKerneradiusFænomenologisk flad kurve
NFW-profil2StærkrsProfil for N-kropssimulering
Einasto2-3Stærkr-2Fleksibel empirisk profil
BeeTheory 3D Yukawa2LovendeBølge-masse-kobling fra skiven

Det yderste punkt i Gaia-æraen er stadig den sværeste begrænsning. Et skarpere fald kan frembringes med en mindre kohærenslængde, men det forværrer den indre tilpasning. Fremtidige data fra Gaia DR4, kuglehobe og stjernestrømme vil være vigtige tests.

Referencer

  • Ou, X. et al. - Mælkevejens mørke stofprofil udledt af dens cirkulære hastighedskurve, MNRAS 528, 2024.
  • Dutertre, X. - Bee Theory™: Bølgebaseret modellering af tyngdekraften, BeeTheory.com v2, 2023.
  • 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.
  • Navarro, J. F., Frenk, C. S., White, S. D. M. - A Universal Density Profile from Hierarchical Clustering, ApJ 490, 1997.

BeeTheory.com - Bølgebaseret kvantegravitation

© Technoplane S.A.S. - 2025