The Missing Mass of the Milky Way: Visible Matter, Rotation Curves, and Dark Matter
TL;DR: The visible matter in the Milky Way—stars, gas, and dust—does not provide enough gravity to explain the observed orbital speeds of stars and gas. From the rotation curve, astronomers infer a larger dynamical mass. The difference between this dynamical mass and the visible mass is called the missing mass, usually modeled as a dark matter halo.
1. The basic problem
In a galaxy, the circular orbital speed v(r) at distance r from the Galactic center depends on the mass enclosed inside that radius. If gravity were produced only by the visible disk, the rotation speed should decline at large radius. Instead, the Milky Way’s rotation curve remains broadly flat over a large radial range, which implies more mass than we directly observe. Rotation-curve studies commonly use the relation between circular speed and enclosed mass to reconstruct the Galaxy’s mass distribution. :contentReference[oaicite:0]{index=0}
2. Dynamical mass from the rotation curve
For an approximately circular orbit, Newtonian dynamics gives:
\[ M_{\rm dyn}(r)=\frac{v(r)^2\,r}{G} \]
where:
- Mdyn(r) is the dynamical mass enclosed within radius r,
- v(r) is the observed circular velocity,
- G is Newton’s gravitational constant.
If the rotation curve is approximately flat, so that:
\[ v(r)\approx v_0 \]
then:
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}\,r \]
This is the central mathematical reason why a flat rotation curve implies a mass that continues to grow roughly linearly with radius.
3. Visible mass of the Milky Way disk
The visible disk of the Milky Way is often approximated by an exponential surface-density profile:
\[ \Sigma_{\rm vis}(r)=\Sigma_0 e^{-r/R_d} \]
where:
- Σ0 is the central surface density,
- Rd is the disk scale length,
- r is the distance from the Galactic center.
The visible mass inside radius r is obtained by adding circular annuli of the disk:
\[ dM_{\rm vis}=2\pi r\,\Sigma_{\rm vis}(r)\,dr \]
which gives:
\[ M_{\rm vis}(r)=2\pi \Sigma_0 R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Equivalently, defining the total disk mass as:
\[ M_d=2\pi\Sigma_0R_d^2 \]
we can write:
\[ M_{\rm vis}(r)=M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
This visible mass saturates at large radius:
\[ M_{\rm vis}(r)\rightarrow M_d \quad \text{for} \quad r\gg R_d \]
This saturation is crucial: the visible disk does not keep adding enough mass to explain the nearly flat rotation curve at large radius.
4. Definition of the missing mass
The missing mass is defined as the difference between the dynamical mass required by the rotation curve and the visible mass actually observed:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Using the equations above:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi \Sigma_0 R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
For a flat rotation curve, v(r) ≈ v0:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
At large radius, because the visible disk mass approaches a constant:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
And asymptotically:
\[ M_{\rm miss}(r)\propto r \]
5. Density of the missing mass
If the missing mass is modeled as a roughly spherical halo, then the corresponding volume density is obtained from:
\[ \rho_{\rm miss}(r)= \frac{1}{4\pi r^2} \frac{dM_{\rm miss}}{dr} \]
In the outer region, where the rotation curve is approximately flat and the visible mass changes slowly:
\[ \frac{dM_{\rm miss}}{dr}\approx \frac{v_0^2}{G} \]
therefore:
\[ \rho_{\rm miss}(r)\approx \frac{v_0^2}{4\pi G r^2} \]
This means that the missing mass density decreases as 1/r², while the enclosed missing mass grows approximately as r.
6. Standard dark matter halo interpretation
In the standard cosmological interpretation, the missing mass is modeled as a dark matter halo surrounding the visible galaxy. A commonly used halo profile is the Navarro–Frenk–White profile, or NFW profile. Milky Way mass models often combine baryonic components—bulge, stellar disk, gas disk—with a dark halo component to fit the observed rotation curve and other dynamical constraints. :contentReference[oaicite:1]{index=1}
The NFW density profile is:
\[ \rho_{\rm NFW}(r)= \frac{\rho_s}{(r/r_s)(1+r/r_s)^2} \]
where:
- ρs is a characteristic density,
- rs is a scale radius.
The enclosed NFW mass is:
\[ M_{\rm NFW}(r)= 4\pi\rho_s r_s^3 \left[ \ln\left(1+\frac{r}{r_s}\right) – \frac{r/r_s}{1+r/r_s} \right] \]
This profile does not produce a perfectly linear mass growth at all radii, but it can reproduce approximately flat rotation curves over the radial range where galaxies are observed.
7. Simplified Milky Way picture
The Milky Way can therefore be summarized with three mass functions:
\[ M_{\rm vis}(r)= M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
The core issue is that:
\[ M_{\rm vis}(r)\rightarrow \text{constant} \]
\[ M_{\rm dyn}(r)\propto r \]
\[ M_{\rm miss}(r)\propto r \]
That mismatch is the mathematical signature of the missing mass problem.
8. Physical interpretation
The visible disk is concentrated: most of its mass lies within a few scale lengths. But the gravitational field inferred from orbital speeds behaves as if additional mass continues to exist far beyond the bright disk. This is why the Milky Way is modeled as a visible baryonic disk embedded in a much larger dark matter halo. Sofue’s Milky Way rotation-curve work, for example, fits bulge, disk, and dark halo components and reports halo parameters using an NFW-type profile. :contentReference[oaicite:2]{index=2}
9. Key equations summary
Visible surface density:
\[ \Sigma_{\rm vis}(r)=\Sigma_0e^{-r/R_d} \]
Visible disk mass:
\[ M_{\rm vis}(r)= 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Dynamical mass:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Missing mass:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Outer-halo approximation:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r \]
Missing-mass density:
\[ \rho_{\rm miss}(r)\approx \frac{v_0^2}{4\pi G r^2} \]
10. Limitations of this simplified model
- The Milky Way is not a perfect exponential disk; it also contains a bulge, bar, gas layers, spiral structure, and stellar halo.
- The relation \(M(r)=v(r)^2r/G\) is exact only for ideal spherical mass distributions; for a thin disk, the gravitational field is more geometrically complex.
- The rotation curve is not perfectly flat at all radii.
- The NFW profile is a model for a dark matter halo, not a direct observation of invisible matter.
- Mass estimates depend on tracer populations, distance measurements, solar position, and assumptions about equilibrium.
Conclusion
The missing mass problem of the Milky Way can be stated mathematically: the observed rotation curve implies a dynamical mass that keeps growing with radius, while the visible disk mass approaches a finite value. This leads to the standard inference of an extended dark matter halo. The essential equations are the visible exponential disk mass, the dynamical mass inferred from rotation, and the missing mass defined as their difference.
Further reading
- McMillan, P. J. “The mass distribution and gravitational potential of the Milky Way.” :contentReference[oaicite:3]{index=3}
- Sofue, Y. “A Grand Rotation Curve and Dark Matter Halo in the Milky Way Galaxy.” :contentReference[oaicite:4]{index=4}
- McMillan, P. J. “Mass models of the Milky Way.” :contentReference[oaicite:5]{index=5}
The visible matter in the Milky Way—stars, gas, and dust—does not provide enough gravity to explain the observed orbital speeds of stars and gas. From the rotation curve, astronomers infer a larger dynamical mass. The difference between this dynamical mass and the visible mass is called the missing mass, usually modeled as a dark matter halo.
1. The basic problem
In a galaxy, the circular orbital speed v(r) at distance r from the Galactic center depends on the mass enclosed inside that radius. If gravity were produced only by the visible disk, the rotation speed should decline at large radius. Instead, the Milky Way’s rotation curve remains broadly flat over a large radial range, which implies more mass than we directly observe. Rotation-curve studies commonly use the relation between circular speed and enclosed mass to reconstruct the Galaxy’s mass distribution. :contentReference[oaicite:0]{index=0}
2. Dynamical mass from the rotation curve
For an approximately circular orbit, Newtonian dynamics gives:
\[ M_{\rm dyn}(r)=\frac{v(r)^2\,r}{G} \]
where:
- Mdyn(r) is the dynamical mass enclosed within radius r,
- v(r) is the observed circular velocity,
- G is Newton’s gravitational constant.
If the rotation curve is approximately flat, so that:
\[ v(r)\approx v_0 \]
then:
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}\,r \]
This is the central mathematical reason why a flat rotation curve implies a mass that continues to grow roughly linearly with radius.
3. Visible mass of the Milky Way disk
The visible disk of the Milky Way is often approximated by an exponential surface-density profile:
\[ \Sigma_{\rm vis}(r)=\Sigma_0 e^{-r/R_d} \]
where:
- Σ0 is the central surface density,
- Rd is the disk scale length,
- r is the distance from the Galactic center.
The visible mass inside radius r is obtained by adding circular annuli of the disk:
\[ dM_{\rm vis}=2\pi r\,\Sigma_{\rm vis}(r)\,dr \]
which gives:
\[ M_{\rm vis}(r)=2\pi \Sigma_0 R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Equivalently, defining the total disk mass as:
\[ M_d=2\pi\Sigma_0R_d^2 \]
we can write:
\[ M_{\rm vis}(r)=M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
This visible mass saturates at large radius:
\[ M_{\rm vis}(r)\rightarrow M_d \quad \text{for} \quad r\gg R_d \]
This saturation is crucial: the visible disk does not keep adding enough mass to explain the nearly flat rotation curve at large radius.
4. Definition of the missing mass
The missing mass is defined as the difference between the dynamical mass required by the rotation curve and the visible mass actually observed:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Using the equations above:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi \Sigma_0 R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
For a flat rotation curve, v(r) ≈ v0:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
At large radius, because the visible disk mass approaches a constant:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
And asymptotically:
\[ M_{\rm miss}(r)\propto r \]
5. Density of the missing mass
If the missing mass is modeled as a roughly spherical halo, then the corresponding volume density is obtained from:
\[ \rho_{\rm miss}(r)= \frac{1}{4\pi r^2} \frac{dM_{\rm miss}}{dr} \]
In the outer region, where the rotation curve is approximately flat and the visible mass changes slowly:
\[ \frac{dM_{\rm miss}}{dr}\approx \frac{v_0^2}{G} \]
therefore:
\[ \rho_{\rm miss}(r)\approx \frac{v_0^2}{4\pi G r^2} \]
This means that the missing mass density decreases as 1/r², while the enclosed missing mass grows approximately as r.
6. Standard dark matter halo interpretation
In the standard cosmological interpretation, the missing mass is modeled as a dark matter halo surrounding the visible galaxy. A commonly used halo profile is the Navarro–Frenk–White profile, or NFW profile. Milky Way mass models often combine baryonic components—bulge, stellar disk, gas disk—with a dark halo component to fit the observed rotation curve and other dynamical constraints. :contentReference[oaicite:1]{index=1}
The NFW density profile is:
\[ \rho_{\rm NFW}(r)= \frac{\rho_s}{(r/r_s)(1+r/r_s)^2} \]
where:
- ρs is a characteristic density,
- rs is a scale radius.
The enclosed NFW mass is:
\[ M_{\rm NFW}(r)= 4\pi\rho_s r_s^3 \left[ \ln\left(1+\frac{r}{r_s}\right) – \frac{r/r_s}{1+r/r_s} \right] \]
This profile does not produce a perfectly linear mass growth at all radii, but it can reproduce approximately flat rotation curves over the radial range where galaxies are observed.
7. Simplified Milky Way picture
The Milky Way can therefore be summarized with three mass functions:
\[ M_{\rm vis}(r)= M_d \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
The core issue is that:
\[ M_{\rm vis}(r)\rightarrow \text{constant} \]
\[ M_{\rm dyn}(r)\propto r \]
\[ M_{\rm miss}(r)\propto r \]
That mismatch is the mathematical signature of the missing mass problem.
8. Physical interpretation
The visible disk is concentrated: most of its mass lies within a few scale lengths. But the gravitational field inferred from orbital speeds behaves as if additional mass continues to exist far beyond the bright disk. This is why the Milky Way is modeled as a visible baryonic disk embedded in a much larger dark matter halo. Sofue’s Milky Way rotation-curve work, for example, fits bulge, disk, and dark halo components and reports halo parameters using an NFW-type profile. :contentReference[oaicite:2]{index=2}
9. Key equations summary
Visible surface density:
\[ \Sigma_{\rm vis}(r)=\Sigma_0e^{-r/R_d} \]
Visible disk mass:
\[ M_{\rm vis}(r)= 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Dynamical mass:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Missing mass:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Outer-halo approximation:
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r \]
Missing-mass density:
\[ \rho_{\rm miss}(r)\approx \frac{v_0^2}{4\pi G r^2} \]
10. Limitations of this simplified model
- The Milky Way is not a perfect exponential disk; it also contains a bulge, bar, gas layers, spiral structure, and stellar halo.
- The relation \(M(r)=v(r)^2r/G\) is exact only for ideal spherical mass distributions; for a thin disk, the gravitational field is more geometrically complex.
- The rotation curve is not perfectly flat at all radii.
- The NFW profile is a model for a dark matter halo, not a direct observation of invisible matter.
- Mass estimates depend on tracer populations, distance measurements, solar position, and assumptions about equilibrium.
Conclusion
The missing mass problem of the Milky Way can be stated mathematically: the observed rotation curve implies a dynamical mass that keeps growing with radius, while the visible disk mass approaches a finite value. This leads to the standard inference of an extended dark matter halo. The essential equations are the visible exponential disk mass, the dynamical mass inferred from rotation, and the missing mass defined as their difference.
Further reading
- McMillan, P. J. “The mass distribution and gravitational potential of the Milky Way.” :contentReference[oaicite:3]{index=3}
- Sofue, Y. “A Grand Rotation Curve and Dark Matter Halo in the Milky Way Galaxy.” :contentReference[oaicite:4]{index=4}
- McMillan, P. J. “Mass models of the Milky Way.” :contentReference[oaicite:5]{index=5}