BeeTheory – Galactische simulatie v2 – eerste generatie 2025 mei 17 met claude

Melkweg Verborgen Massa: Bijentheorie 3D Yukawa met fysische schijfverkleining

De gecorrigeerde simulatie: de baryonische schijfsnelheid valt Kepleriaans buiten de fysieke rand, en de 3D Yukawa kernel van BeeTheory vult de hele ruimte. Twee parameters, rotatiegegevens uit het Gaia-tijdperk en een afgekapt schijfmodel.

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

K = 0,040 kpc-¹

Golfkoppeling

α = 0,087 kpc-¹

Inverse coherentie

ℓ = 11,5 kpc

Coherentielengte

χ²/dof ≈ 0,31

Uitstekende vereenvoudigde pasvorm

0. Resultaat – vergelijkingen en parameters

Elke ringvormige ring van de galactische schijf met straal R′ genereert een 3D effectief donker massaveld door de BeeTheory Yukawa kernel. De totale donkere dichtheid bij bolstraal r is:

\(\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}\)

De kernel is afgeleid van de gecorrigeerde BeeTheory-krachtwet:

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

Het reduceert tot de Newtoniaanse omgekeerd-kwadratische vorm voor D die veel kleiner is dan de coherentielengte ℓ.

\(F(D)\propto Frac{1}{D^2}\).

De baryonische schijfsnelheid gebruikt de Freeman-formule binnen zijn fysische rand _Rtrunc ≈ _4Rd = 10,4 kpc, en gaat dan soepel over naar de Kepleriaanse val die verwacht wordt van een eindige massadistributie.

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

Overzicht pasvorm

Waarneembaar	Gaia-tijdperk waarde	Bijentheorie	Trek
_Vc(4 kpc)	220 ± 10 km/s	219,8 km/s	-0.02σ
_Vc(8 kpc)	230 ± 6 km/s	233,2 km/s	+0.53σ
_Vc(12 kpc)	226 ± 7 km/s	223,8 km/s	-0.31σ
_Vc(20 kpc)	215 ± 10 km/s	211,2 km/s	-0.38σ
_Vc(27.3 kpc)	173 ± 17 km/s	199,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⊙	sluiten
_Mtot(<200 kpc)	5-9 × 10¹¹ M⊙	3.3 × 10¹¹ M⊙	laag segment

De vereenvoudigde fit geeft χ²/dof ≈ 0,31. Het moeilijkste punt blijft de buitenste Gaia-tijdperk waarde op 27,3 kpc, waar de waargenomen afname scherper is dan dit twee-parametermodel voorspelt.

1. De schijfonderbreking – Waarom en hoe

1.1 Het probleem met een oneindige exponentiële schijf

De Freeman-schijfformule gaat uit van een exponentiële oppervlaktedichtheid die zich tot in het oneindige uitstrekt. Wiskundig bereikt deze nooit nul, maar fysisch heeft de stellaire schijf van de Melkweg een eindige omvang. Voorbij de effectieve stellaire rand is de ingesloten baryonische massa in essentie constant, en moet de snelheidsbijdrage ongeveer als een Kepleriaans puntmassa-veld vallen.

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

Voorbij de schijfrand neigt de baryonische snelheid naar:

\(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\)

Voorbeeldwaarden zijn:

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

1.2 Vloeiende Truncatieformule

De simulatie gebruikt een vloeiende overgang tussen de Freeman-schijfformule en de Kepleriaanse waarde. De overgang is gecentreerd op _Rtrunc = _4Rd = 10,4 kpc met breedte σ = 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}\)

De minimumfunctie voorkomt dat de baryonische schijf de fysieke Kepleriaanse limiet buiten de schijfrand overschrijdt.

R	_VFreeman	_VKepleriaan	Vbar_,afgekapt	Dominant regime
5 kpc	174,5 km/s	201,1 km/s	174,5 km/s	Freeman
8 kpc	161,5 km/s	159,0 km/s	161,5 km/s	Freeman ≈ Kepler
10,4 kpc	143,0 km/s	139,3 km/s	141,2 km/s	Overgang
16 kpc	112,4 km/s	112,4 km/s	112,4 km/s	Kepleriaans
25 kpc	89,9 km/s	89,9 km/s	89,9 km/s	Kepleriaans
50 kpc	63,6 km/s	63,6 km/s	63,6 km/s	Kepleriaans

2. De BeeTheory 3D Donkere Massadichtheid

2.1 Schijfringen die uitstralen in 3D

Elke ring van de galactische schijf bij straal R′ met breedte dR′ heeft massa:

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

In de BeeTheory genereert deze ring een gravitatiegolfveld dat zich in alle drie de ruimtelijke dimensies voortplant. In de monopoolbenadering is de afstand tot een 3D-veldpunt bij bolstraal r:

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

De numerieke vorm van de donkere dichtheid is:

\(\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 Ingesloten donkere massa en cirkelsnelheid

\(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 Asymptotisch gedrag

\(\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}\)

Voor α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{constant}\)

3. Simulatieresultaten – Interactieve grafieken

De simulatie hieronder behoudt het numerieke model, de schuifregelaars, de rotatiecurve, het massaprofiel, het dichtheidsprofiel en de live χ²-update. Plak deze pagina in WordPress met scriptuitvoering ingeschakeld.

Rotatiecurve – Gaia-era data vs BeeTheory met afgeknotte schijf

Alleen baryonen, afgekapt BeeTheorie totaal Donkere component Gaia-tijdperk gegevens

Parameterverkenner – pas K, α en _Rtrunc aan

K (kpc-¹) 0.040

α (kpc-¹) 0.087

_Rtrunc (kpc) 10.4

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

Massaprofiel: zichtbare schijf vs 3D donkere massa vs totaal

Zichtbare schijf + uitstulping Bijentheorie donkere massa Totale massa

r (kpc)	_Mbar (10¹⁰ M⊙)	_Mdark (10¹⁰ M⊙)	_Mtot (10¹⁰ M⊙)	DM/bar	_ρdark (GeV/cm³)
Laden…

Donker dichtheidsprofiel _ρdark(r) – log-schaal

BeeTheory Isothermische r-² referentie NFW referentie

4. Fysische interpretatie en universaliteit

4.1 coherentielengte

Binnen de coherentielengte gedraagt de Yukawa kern zich bijna als een Newtoniaanse 1/D² kernel. De donkere dichtheid volgt ongeveer r-² en de rotatiecurve is vlak. Voorbij ℓ veroorzaakt de exponentiële onderdrukking de afname die in de buitenste schijf wordt waargenomen.

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

4.2 Dimensieloze koppeling

Een dimensieloze BeeTheory-koppeling kan als volgt worden gedefinieerd:

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

Dit is in orde van grootte vergelijkbaar met de koppeling afgeleid uit de H₂-kalibratie, waarbij λ ongeveer 3-4 is. De mogelijke schaaluniversaliteit van dit getal blijft een centrale open vraag.

4.3 Vergelijking met standaardmodellen

Model	Parameters	Typische pasvorm	Schaal	Mechanisme
Isotherme halo	2	Matig	kernradius	Fenomenologische vlakke curve
NFW profiel	2	Sterk	_rs	N-lichaam simulatieprofiel
Einasto	2-3	Sterk	_r-2	Flexibel empirisch profiel
Bijentheorie 3D Yukawa	2	Veelbelovend	ℓ	Golfmassakoppeling van de schijf

Het buitenste Gaia-tijdperk punt blijft de moeilijkste beperking. Een scherpere daling kan worden verkregen met een kleinere coherentielengte, maar dat verslechtert de binnenste fit. Toekomstige gegevens van Gaia DR4, bolvormige sterrenhopen en stellaire stromen zullen belangrijke tests zijn.

Referenties

Ou, X. et al. - The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 2024.
Dutertre, X. - Bee Theory™: Op golven gebaseerde modellering van zwaartekracht, 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 - Op golven gebaseerde kwantumzwaartekracht