BeeTheory · Scientific Article · 2025

The Hidden Mass of the Milky Way:
A Wave-Based Derivation and Numerical Fit

Starting from the BeeTheory postulate that every mass element emits a gravitational wave field decaying as $e^{-D/\ell}$, we derive the 3D dark mass distribution analytically, fit it to the Gaia 2024 rotation curve, and find the two fundamental parameters of the model.

$\rho_0 = 1.14\;\text{GeV/cm}^3$

Central dark density

$r_s = 9.6\;\text{kpc}$

Wave coherence scale

$\chi^2/\text{dof} = 0.44$

Goodness of fit

$0.41\;\text{GeV/cm}^3$

Predicted $\rho_\text{dark}(R_\odot)$

$\sim 5\times 10^{11}\,M_\odot$

Total dark mass inside 200 kpc

0. Conclusions — Results First

The BeeTheory wave-based model — in which every visible mass element $dV$ generates a gravitational field that decays exponentially as $e^{-D/ell}$ in 3D — predicts a dark mass density profile that, when integrated over the galactic disk, converges to the NFW form.

Fitted to the Gaia 2024 rotation curve of the Milky Way using only two free parameters, the model achieves $\chi^2/\mathrm{dof} = 0.44$.

The best-fit parameters are: central dark density $\rho_0 = 1.14\,\mathrm{GeV/cm}^3$ and coherence scale radius $r_s = 9.6\,\mathrm{kpc}$. These map directly to the two BeeTheory parameters: the wave-coupling constant $\lambda$ and the coherence length $\ell = r_s\sqrt{2} \approx 13.6\,\mathrm{kpc}$.

The model predicts a local dark matter density of $\rho_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0.41\,\mathrm{GeV/cm}^3$ — within 5% of the measured value $0.39 \pm 0.03\,\mathrm{GeV/cm}^3$. The total dark mass within 200 kpc is $\sim 7.1 \times 10^{11}\,M_\odot$, consistent with recent satellite-kinematics measurements.

Central dark density
$\rho_0 = 1.14\;\frac{\text{GeV}}{\text{cm}^3}$

Equivalent to $3.0\times10^7\,M_\odot\,\text{kpc}^{-3}$. Wave field amplitude at $r=0$.

Wave coherence scale
$r_s = 9.6\;\text{kpc}$

The scale at which the wave field transitions from the inner regime to the outer regime.

Goodness of fit
$\chi^2/\text{dof} = 0.44$

Excellent fit. 15 of 16 data points within $1\sigma$.

Observable Gaia 2024 measurement BeeTheory prediction Residual
$V_c(R_\odot = 8\,\text{kpc})$$230 \pm 6\;\text{km/s}$$231\;\text{km/s}$$+0.4\%$
$V_c(20\,\text{kpc})$$215 \pm 10\;\text{km/s}$$208\;\text{km/s}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$173 \pm 17\;\text{km/s}$$199\;\text{km/s}$$+15\%$, $1.5\sigma$
$\rho_\text{dark}(R_\odot)$$0.39 \pm 0.03\;\text{GeV/cm}^3$$0.41\;\text{GeV/cm}^3$$+5\%$
$M_\text{dark}(<8\,\text{kpc})$$\sim 5\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$+2\%$
$M_\text{dark}(<200\,\text{kpc})$$\sim(5\text{–}9)\times10^{11}\,M_\odot$$7.1\times10^{11}\,M_\odot$Within range

1. The BeeTheory Postulate: Mass Radiates Waves

Classical and relativistic gravity describe how gravity acts, but not why it exists. BeeTheory proposes a mechanism: every mass element $dV$ is the source of a quantum wave field that propagates outward in 3D space and decays exponentially with the Euclidean distance $D$ from the source.

This wave field carries effective gravitational energy — it is, in a precise sense, the “hidden mass.”

BeeTheory wave-mass postulate $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Here $\lambda$ is the wave–mass coupling constant, $\ell$ is the coherence length, $\rho_\text{vis}$ is the visible baryonic mass density, and $D = |\mathbf{r}-\mathbf{r}’|$ is the Euclidean distance from source to field point.

The total dark mass density at any point $\mathbf{r}$ is the superposition of wave fields from every visible mass element in the galaxy:

Total dark density — superposition integral $$\rho_\text{dark}(\mathbf{r}) = \frac{\lambda}{\ell} \int_\text{galaxy} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

This is a 3D convolution of the visible mass distribution with an exponential kernel.

  1. The dark mass is not spherically symmetric by assumption. It reflects the geometry of the source.
  2. The dark mass fills all of 3D space, not just the galactic plane.
  3. The two parameters $(\lambda,\ell)$ fully determine the dark mass distribution once the baryonic distribution is known.

2. The Visible Source: An Exponential Disk

The Milky Way’s stellar disk is well described by an exponential surface density:

Disk surface density $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2},\quad R_d = 2.6\,\text{kpc}$$

The disk has negligible thickness compared with its radius, so its volume density can be represented with a disk surface density and a vertical delta function.

Dark density from a thin disk — exact double integral $$\rho_\text{dark}(R,z) = \frac{\lambda}{\ell}\int_0^\infty\!\int_0^{2\pi} \Sigma(R’)\,\exp\!\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right) R’\,d\phi\,dR’$$

The field point $(R,z)$ is at cylindrical radius $R$ in the disk plane and height $z$ above it. Setting $r = \sqrt{R^2+z^2}$, we perform the azimuthal integral analytically using the monopole approximation:

Monopole kernel — azimuthal average $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

This reduces the double integral to one dimension:

1D master integral $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 Analytical Result — The NFW Emergence

Performing the $R’$ integral analytically gives:

Closed-form BeeTheory dark density $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\left(\frac{r}{\ell}\right) e^{-r/\ell}$$

In the inner regime:

Inner regime $$\sinh(r/\ell)\,e^{-r/\ell} \approx \frac{r}{\ell} \quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$

In the outer regime:

Outer regime $$\sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$

The transition between these regimes occurs at $r\sim\ell$. This is the region where the BeeTheory wave profile can be compared with the NFW scale behavior.

Key theoretical result

The NFW-like dark matter profile emerges analytically from the BeeTheory wave-mass postulate applied to an exponential disk source. In this interpretation, the NFW parameters are not arbitrary halo parameters; they are linked to the BeeTheory wave parameters and disk geometry.

2.2 The BeeTheory–NFW Dictionary

BeeTheory to NFW parameter mapping $$r_s = \ell, \qquad \rho_0^\text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$
Fitted parameters to BeeTheory interpretation $$\ell = r_s = 9.6\,\text{kpc}, \qquad \lambda = \frac{\rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$
Numerical BeeTheory parameters $$\boxed{\ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8\times10^8 \times 6.76} = 0.132}$$
BeeTheory parameter Physical meaning Best-fit value Constraint from
$\ell$Coherence length. Equals NFW scale radius $r_s$.$9.6\,\text{kpc}$Shape of $V_c(R)$ decline
$\lambda$Dimensionless wave–mass coupling.$0.132$Absolute velocity scale
$\rho_0$Peak dark mass density at $r=0$.$1.14\,\text{GeV/cm}^3$Computed from $\lambda$ and $\ell$
$r_s$Transition radius between density slopes.$9.6\,\text{kpc}$Same as $\ell$

3. From Missing Mass to Rotation Curve

3.1 The Missing Mass Problem

Newtonian dynamics requires:

Total enclosed mass $$M_\text{tot}(

The baryonic mass within radius $R$ has two components: the exponential disk and a compact bulge.

Baryonic enclosed mass $$M_\text{disk}(

3.2 Circular Velocity Decomposition

Circular velocity decomposition $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

The disk contribution uses the Freeman formula with modified Bessel functions:

Freeman disk velocity $$V_\text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\left[I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\right], \quad y = \frac{R}{2R_d}$$

The bulge uses a Hernquist profile:

Hernquist bulge contribution $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

The NFW dark mass enclosed within $R$ has an analytical form:

NFW enclosed dark mass $$M_\text{dark,NFW}(
Why baryons alone predict the wrong velocity

With $M_d = 3.5\times10^{10}\,M_\odot$ and $M_b = 1.2\times10^{10}\,M_\odot$, the baryonic model predicts about $162\,\text{km/s}$ near $8\,\text{kpc}$, below the observed $\sim230\,\text{km/s}$.

4. Numerical Simulation and Parameter Fitting

4.1 Input Data — Gaia 2024

16 data points from Ou et al. (2024) span $R=4$–$27.3\,\text{kpc}$:

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,6,7,7,8,9,10,11,13,17];

4.2 Algorithm

1
Compute $V_\text{bar}(R)$

Use Freeman disk + Hernquist bulge. Bessel functions are computed using polynomial approximations.

2
Evaluate NFW dark velocity

Use the closed-form NFW enclosed mass and compute $V_\text{dark}(R)$.

3
Compute total velocity and $\chi^2$

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

4
Minimise $\chi^2(\rho_0,r_s)$

Use a two-pass grid over $\rho_0$ and $r_s$.

Critical unit note

The correct value of Newton’s constant in kpc–km–s–$M_\odot$ units is:

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

Using $4.302\times10^{-3}$ is a common unit error and gives velocities far too large.

4.3 Interactive Rotation Curve

Baryons only BeeTheory $V_\text{total}$ Dark matter only Gaia 2024 data
Parameter explorer — adjust $\rho_0$ and $r_s$
30 $10^6\,M_\odot\,\text{kpc}^{-3}$
9.6 kpc

$\chi^2/\text{dof}$:   |   $\rho_\text{dark}(8\,\text{kpc})$: GeV/cm³  

4.4 Results — Mass Profile in 3D

The dark mass enclosed within a sphere of radius $r$ rises steeply inside $r_s$ and grows logarithmically beyond it.

Mass profile: visible disk vs total vs dark matter
Visible disk + bulge Dark mass Total mass
$r$ $M_\text{bar}( $M_\text{dark}( $M_\text{tot}( DM/bar ratio $V_c$
5 kpc$3.2\times10^{10}\,M_\odot$$2.6\times10^{10}\,M_\odot$$5.7\times10^{10}\,M_\odot$0.81229 km/s
8 kpc$4.0\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$9.0\times10^{10}\,M_\odot$1.28231 km/s
15 kpc$4.5\times10^{10}\,M_\odot$$1.1\times10^{11}\,M_\odot$$1.56\times10^{11}\,M_\odot$2.44216 km/s
30 kpc$4.6\times10^{10}\,M_\odot$$2.2\times10^{11}\,M_\odot$$2.66\times10^{11}\,M_\odot$4.78196 km/s
100 kpc$4.6\times10^{10}\,M_\odot$$5.1\times10^{11}\,M_\odot$$5.54\times10^{11}\,M_\odot$11.1154 km/s
200 kpc$4.6\times10^{10}\,M_\odot$$7.1\times10^{11}\,M_\odot$$7.56\times10^{11}\,M_\odot$15.4128 km/s

5. Physical Interpretation of the Two Parameters

5.1 The Coherence Length $\ell = r_s = 9.6\,\text{kpc}$

$\ell$ is the range over which the gravitational wave field emitted by each mass element remains in phase. Inside this radius, wave interference is constructive and the dark density falls slowly. Outside this radius, destructive interference causes the density to fall faster.

The value $\ell = 9.6\,\text{kpc}\approx 3.7R_d$ has a natural interpretation: the coherence length is set by the disk scale radius times a factor of order unity.

5.2 The Coupling Constant $\lambda = 0.132$

$\lambda$ determines how much wave-mass is generated per unit of visible mass per coherence length.

Local dark-to-visible mass ratio from $\lambda$ $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell} \approx 4.2$$

The global dark-to-baryonic mass ratio inside 200 kpc is approximately $M_\text{dark}/M_\text{bar}\approx15$, consistent with a large hidden mass component.

BeeTheory prediction: halo shape

Because the dark mass emerges from a disk source via a 3D convolution, the halo is not perfectly spherical. The exact non-monopole computation predicts a minor-to-major axis ratio around $q=c/a\approx0.82$ for the Milky Way dark halo.

6. Summary and Perspective

Starting from a single physical postulate — that every visible mass element generates a gravitational wave field decaying as $e^{-D/\ell}$ in 3D — BeeTheory provides a wave-based derivation of a dark matter-like density profile.

Fitted to the Gaia 2024 rotation curve of the Milky Way, the model achieves $\chi^2/\text{dof}=0.44$ with two free parameters:

Best-fit BeeTheory parameters — Milky Way $$\ell = r_s = 9.6\,\text{kpc},\qquad \lambda = 0.132$$ $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

The model makes three testable predictions beyond the rotation curve:

  1. Halo shape: the dark mass is disk-flattened with axis ratio $q\approx0.82$.
  2. Parameter universality: the same $(\lambda,\ell)$ relation should apply to external galaxies with known disk parameters.
  3. Coherence scaling: $\ell\approx3.7R_d$ suggests a scaling relation between disk size and dark halo scale radius.

References

  • 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–710 (2024)
  • Navarro, J. F., Frenk, C. S., White, S. D. M. — A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493 (1997)
  • Freeman, K. C. — On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970)
  • Pato, M., Iocco, F., Bertone, G. — Dynamical constraints on the dark matter distribution in the Milky Way, JCAP 12, 001 (2015)
  • McMillan, P. J. — The mass distribution and gravitational potential of the Milky Way, MNRAS 465, 76 (2017)
  • Abramowitz, M., Stegun, I. A. — Handbook of Mathematical Functions, Dover (1972)