Dark Matter Density Profiles
We start with the (collisionless) Boltzmann equation:
which says that the change in the distribution of velocities in a volume which is comoving with the fluid (), is related to the change in distribution of a non-comoving element (), plus a factor which represents a flux over the boundary of the comoving volume (both in position and velocity space). This equation is not equal to 0 when collisions become important.
We now begin taking moments of this equation.
where the last term of the first equation was eliminated using the divergence theorem: . Now is the number density of stars, and , so we have our first Jean’s equation:
This equation is the star analog for the continuity equation in fluids:
For stars, instead of conserving mass, we conserve number. Next, we use:
Physically, the last term is , and represents the correlation between velocities in two different directions. Instead of going through the derivation, we will jump straight to the second Jean’s equation:
where . This equation is the star-equivalent of the fluid momentum equation:
In fluids, we have the an energy conservation equation. In the stellar dynamics case, this equation is not useful because the only way to change the energy of a galaxy is to hit it with another galaxy. The one exception to this is in globular clusters, where stars can get evaporated out of the group.
Now, unfortunately, we need to talk about these equations in cylindrical coordinates. The Jeans conservation equation becomes:
The second equation becomes:
For a plane-parallel sheet in steady state, we know that , , and , so in our second equation, that leaves us with:
Which is an equation which agrees nicely with observations of disk galaxies. Combining this with Poisson’s equation:
which is an acceleration in the z direction. Thus:
Thus, we can measure the velocity and number density in the z direction and we can get an estimate of the total mass density in the midplane of a galaxy. Typically, this is done using K giants because they are old (and therefore, more representative of the total population). We make the assumption here that all stars have the same dynamical temperature (). It turns out that this assumption breaks down, as it was discovered that there are “old thin disk” and “thick disk” stellar populations.
If we integrate the above equation, we can get an equation for the total density along a column . This determination was used to show there was no excess dark matter associated with the disk compared to the halo.
Graphing , we find it grows linearly with up to 300 pc, with . Because grows linearly, we can infer that is relatively constant near the midplane. We measure and