BeeTheory · Gravity & Wave Physics
The Radial Equations of the Milky Way’s Hidden Mass
From density profiles to ring integrals and rotation curves — a mathematical treatment of hidden mass as a function of galactic radius R.
This page introduces the radial equations used to describe the Milky Way’s hidden mass. It compares classical dark matter density profiles, ring and shell integrals, enclosed mass equations, rotation curves, and the BeeTheory interpretation of missing mass as a possible wave-interference effect.
~10¹² M⊙
Classical estimate for the total Milky Way halo mass.
0.39 GeV/cm³
Typical local dark matter density near the Sun.
R⊙ ≈ 8 kpc
Approximate distance of the Sun from the Galactic Center.
~200 kpc
Approximate outer scale used for Milky Way halo estimates.
Contents
- Why is mass missing?
- Density profiles ρ(R)
- Ring and annulus mass dM/dR
- Enclosed dark matter mass M(<R)
- The rotation curve V(R)
- Current observational estimates
- Competing hypotheses
- The BeeTheory perspective
1. Why Is Mass Missing?
In 1933, Fritz Zwicky noticed that galaxies in the Coma Cluster moved far too fast to be held together by their visible mass alone. In the 1970s, Vera Rubin and Kent Ford measured the rotation curves of spiral galaxies and found something equally striking: stars at large radii orbit almost as fast as those near the center, while Newtonian gravity from visible matter predicts that they should slow down.
For a simple Keplerian orbit around a central mass, we expect:
\(V(R)\propto \frac{1}{\sqrt{R}}\)What is observed instead is an approximately flat, or only slowly declining, rotation curve:
\(V(R)\approx V_{\infty}=\mathrm{const}\qquad \mathrm{for}\ R\gtrsim 5\,\mathrm{kpc}\)Reconciling these facts with Newtonian gravity requires an additional invisible mass component whose density falls off approximately as:
\(\rho(r)\propto r^{-2}\)This produces a total enclosed mass proportional to radius:
\(M(Key quantitative puzzle
The luminous baryonic mass of the Milky Way is approximately 5 × 10¹⁰ M⊙. The total dynamical mass inferred from kinematics out to roughly 200 kpc is around 10¹² M⊙. This implies a dark-to-luminous mass ratio of roughly 10 to 1.
2. Density Profiles ρ(R)
A density profile is a mathematical function that describes how the dark matter density ρ varies with galactocentric radius r, or with cylindrical radius R in the galactic plane.
2.1 NFW Profile
The NFW profile, introduced by Navarro, Frenk, and White, is derived from N-body cosmological simulations. It produces a characteristic double power law with a central cusp.
\(\rho_{\mathrm{NFW}}(r)=\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}\)| Parameter | Symbol | Milky Way estimate | Role |
|---|---|---|---|
| Scale radius | rs | 15–25 kpc | Transition between inner and outer slopes |
| Characteristic density | ρ0 | Calibrated to local dark matter density | Overall normalization |
| Inner slope | γ | −1 | Cuspy behavior |
| Outer slope | — | −3 | Rapid decline at large radius |
2.2 Einasto Profile
The Einasto profile avoids a strict central divergence and uses a shape parameter α that allows the density slope to change smoothly with radius.
\(\rho_{\mathrm{Ein}}(r)=\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)| Parameter | Symbol | Milky Way estimate | Role |
|---|---|---|---|
| Shape index | α | Model-dependent | Controls how fast the slope changes |
| Scale radius | r−2 | ~18–22 kpc | Radius where logarithmic slope equals −2 |
| Density at r−2 | ρ−2 | Calibrated to local density | Normalization |
Recent observational tension
Recent Gaia-based studies suggest that the Milky Way rotation curve may decline more rapidly beyond the solar radius than a standard NFW halo would predict. This makes cored or smoothly varying profiles such as Einasto particularly important in current discussions.
2.3 Pseudo-Isothermal Profile
The pseudo-isothermal profile is often used as a simple analytic approximation for a cored halo.
\(\rho_{\mathrm{iso}}(r)=\frac{\rho_0}{1+\left(\frac{r}{r_s}\right)^2}\)At small radius, the density approaches a constant value. At large radius, it falls as r⁻² and produces a flat rotation curve.
\(V_{\infty}=\sqrt{4\pi G\rho_0 r_s^2}\)Cusp versus core problem
N-body simulations often predict cuspy NFW profiles, while many observed dwarf galaxies appear to prefer cored density profiles. This cusp-core problem remains one of the main unresolved issues in dark matter physics.
3. Ring and Annulus Mass — dM/dR
To compute how much dark matter sits in each radial slice of the Galaxy, we integrate the density over a thin shell or annulus. The geometry depends on whether the halo is treated as spherical or flattened.
3.1 Spherical Thin Shell
For a spherically symmetric halo, the mass in a shell of thickness dr at radius r is:
\(\frac{dM_{\mathrm{shell}}}{dr}=4\pi r^2\rho(r)\)3.2 Disk-Plane Annular Ring
For a ring lying in the Galactic plane, with cylindrical radius R and effective half-thickness H(R), the annulus mass is:
\(dM_{\mathrm{ann}}=2\pi R\,\rho(R,0)\,2H(R)\,dR\)For a spherical halo, this can be written as an integral over height z:
\(\frac{dM}{dR}=2\pi R\int_{-\infty}^{+\infty}\rho\left(\sqrt{R^2+z^2}\right)dz\)In the spherical approximation, this connects back to:
\(\frac{dM}{dR}\approx4\pi R^2\rho(R)\)3.3 NFW Mass Per Shell
\(\frac{dM_{\mathrm{NFW}}}{dr}=4\pi r^2\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}=\frac{4\pi\rho_0 r_s r}{\left(1+\frac{r}{r_s}\right)^2}\)This function peaks around the scale radius rs, meaning that much of the dark matter mass per shell is deposited in the intermediate halo rather than only at the center or at the outskirts.
3.4 Einasto Mass Per Shell
\(\frac{dM_{\mathrm{Ein}}}{dr}=4\pi r^2\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)The Einasto enclosed mass is usually evaluated numerically.
Physical meaning
The function dM/dr tells us which Galactic radius contributes most to the hidden mass budget. A steeper outer profile reduces the inferred total halo mass, while a shallower profile increases it.
4. Enclosed Dark Matter Mass M(<R)
Integrating the shell element from 0 to R gives the total dark matter mass enclosed within radius R:
\(M_{\mathrm{DM}}(The total enclosed dynamical mass can be decomposed into visible and hidden components:
[latex]M_{\mathrm{tot}}(~3 × 10¹⁰ M⊙
Approximate hidden mass inside the solar circle.
~1–2 × 10¹¹ M⊙
Approximate hidden mass inside 20 kpc.
5–12 × 10¹¹ M⊙
Approximate hidden mass inside the virial region, strongly model-dependent.
The mass profile remains model-dependent.
Estimates of the Milky Way halo mass depend strongly on how the outer halo is extrapolated beyond the region with strong observational constraints.
5. The Rotation Curve V(R)
The circular velocity at radius R is set by the total enclosed mass through the balance of gravitational pull and centripetal acceleration:
[latex]V_c(R)=\sqrt{\frac{G\,M_{\mathrm{tot}}(5.1 Baryonic Disk Contribution
The stellar thin disk follows an exponential surface density profile:
\(\Sigma(R)=\Sigma_0\exp\left(-\frac{R}{R_d}\right)\)The corresponding circular velocity for an exponential disk is:
\(V_{\mathrm{disk}}^2(R)=\frac{2GM_d}{R_d}y^2\left[I_0(y)K_0(y)-I_1(y)K_1(y)\right],\qquad y=\frac{R}{2R_d}\)Here In and Kn are modified Bessel functions. Typical Milky Way thin-disk parameters are Rd ≈ 2.6 kpc and Md ≈ 3.5 × 10¹⁰ M⊙.
5.2 Dark Matter Contribution
\(V_{\mathrm{DM,NFW}}(R)=\sqrt{\frac{4\pi G\rho_0r_s^3}{R}\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]}\)5.3 Baryonic Tully-Fisher Relation
The baryonic Tully-Fisher relation connects the flat rotation velocity of a galaxy to its total baryonic mass:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0,\qquad a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)~230 km/s
Circular velocity near the solar radius.
~170–180 km/s
Possible declining value in the outer disk, depending on tracer data.
~150 km/s
Approximate outer-halo velocity scale from halo tracers.
6. Current Observational Estimates
The table below summarizes representative values for dark matter density and mass at key Galactic radii. Exact values vary by dataset, tracer population, and halo model.
| Radius R | Dark matter density | Enclosed dark mass | Method |
|---|---|---|---|
| Center | Divergent in NFW, finite in cored models | Model-dependent | N-body simulations and inner-Galaxy modeling |
| R⊙ ≈ 8 kpc | ~0.39 GeV/cm³ | ~3 × 10¹⁰ M⊙ | Rotation curve and vertical kinematics |
| 20 kpc | ~0.05 GeV/cm³ | ~1–2 × 10¹¹ M⊙ | Gaia and spectroscopic tracers |
| 50 kpc | ~5 × 10⁻³ GeV/cm³ | ~3–5 × 10¹¹ M⊙ | Globular clusters and halo stars |
| 100–200 kpc | ≤10⁻³ GeV/cm³ | ~5–12 × 10¹¹ M⊙ | Satellite galaxies and escape-velocity methods |
Combining globular cluster kinematics, halo stars, satellite galaxies, and Gaia-era astrometry suggests that the Milky Way’s outer halo profile remains uncertain. This uncertainty is central to the hidden-mass problem.
7. Competing Hypotheses for the Missing Mass
Several major families of explanation remain active. None has been definitively confirmed or ruled out across all observational scales.
7.1 Cold Dark Matter Particles
Cold dark matter remains the leading paradigm. Candidate particles include WIMPs, sterile neutrinos, and other beyond-Standard-Model possibilities. These candidates form extended halos often modeled with NFW or Einasto profiles.
\(m_{\chi}\sim10\text{–}1000\,\mathrm{GeV}\)The main tension is experimental: direct detection has not yet found a confirmed dark matter particle.
7.2 Ultralight or Fuzzy Dark Matter
Fuzzy dark matter uses ultralight particles whose de Broglie wavelength can become astrophysically large, suppressing small-scale structure.
\(m_a\sim10^{-22}\,\mathrm{eV}\) \(\lambda_{\mathrm{dB}}\sim\mathrm{kpc}\)This framework naturally produces smoother inner density cores, but it is constrained by Lyman-alpha forest data and dwarf-galaxy structure.
7.3 Modified Newtonian Dynamics
MOND modifies the effective gravitational acceleration below a characteristic scale:
\(a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)In the deep-MOND regime, the effective acceleration becomes:
\(g_{\mathrm{eff}}=\sqrt{g_Na_0}\)MOND predicts the baryonic Tully-Fisher relation:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0\)It works well for many galaxy rotation curves, but galaxy clusters and cosmology remain difficult.
7.4 Self-Interacting Dark Matter
Self-interacting dark matter proposes that dark matter particles interact with each other strongly enough to reshape inner halo density profiles.
\(\frac{\sigma}{m}\sim1\text{–}100\,\mathrm{cm^2/g}\)This may help explain the diversity of halo cores, but no specific particle candidate has yet been confirmed.
7.5 Primordial Black Holes
Primordial black holes formed in the early universe could make up part of the hidden mass. Many mass windows are strongly constrained by microlensing, cosmic microwave background, and gravitational-wave observations.
\(10^{-16}\text{–}10^{-11}\,M_\odot\)They remain speculative as a full explanation for the Milky Way’s hidden mass.
8. The BeeTheory Perspective
BeeTheory proposes that gravity may be understood as an emergent effect arising from wave behavior rather than as a fundamental force carried only by a particle or produced only by spacetime curvature.
In this framework, every massive system is associated with a wave function ψ(r,t). A basic quantum starting point is the three-dimensional Schrödinger equation:
\(i\hbar\frac{\partial\psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\psi+V(\mathbf{r})\psi\)When two mass distributions approach each other, their wave functions overlap. The convolution of these wave functions can be written as:
\(\Psi_{12}(\mathbf{r})=(\psi_1*\psi_2)(\mathbf{r})=\int\psi_1(\mathbf{r}’)\psi_2(\mathbf{r}-\mathbf{r}’)\,d^3r’\)BeeTheory interprets gravitational attraction as a large-scale manifestation of structured wave overlap, resonance, and field coherence.
8.1 BeeTheory Reinterpretation of Hidden Mass
In BeeTheory, what is usually called dark matter may be interpreted as the cumulative gravitational effect of wave interference from many oscillatory systems distributed throughout the galactic halo.
\(\rho_{\mathrm{eff}}(R)=\rho_{\mathrm{bar}}(R)+\Delta\rho_{\mathrm{wave}}(R)\)Here Δρwave(R) represents an additional effective gravitational density arising from coherent wave-field structure rather than from directly visible baryonic matter.
This term would need to reproduce the radial behavior normally attributed to dark matter. In particular, it would need to generate approximately flat rotation curves over the relevant Galactic range.
\(\rho_{\mathrm{wave}}(R)\propto R^{-2}\)Open quantitative challenge
BeeTheory must show whether a wave-based interference model can reproduce the precise radial density profile required by observed rotation curves. It must also explain why the effective hidden mass is often much larger than the visible baryonic mass.
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.
- Einasto, J. — On the construction of a composite model for the Galaxy, Trudy 5, 87, 1965.
- Watkins, L. L., van der Marel, R. P. et al. — Evidence for an Anticorrelation between the Masses of the Milky Way and Andromeda Galaxies, ApJ 873, 111, 2019.
- Milgrom, M. — A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis, ApJ 270, 365, 1983.
- McGaugh, S. S. et al. — Radial Acceleration Relation in Rotationally Supported Galaxies, PRL 117, 201101, 2016.
Note: recent or future-dated references should be verified before publication if the page is used as a scientific citation source.
BeeTheory Perspective
The hidden mass problem is not only a question of how much matter is missing. It is a question of what kind of physical structure produces gravity at galactic scale.
Classical dark matter models interpret the missing mass as invisible matter. BeeTheory explores a complementary possibility: part of the hidden gravitational effect may arise from structured wave coherence.
The next step is mathematical: define the radial wave-density term, derive its rotation curve, and compare it directly with Gaia-era Milky Way data.