Hernquist Potential

The Poisson equation ∇² Φ = 4πGρ has a finite number of solutions with potential-density pairs in which both the gravitational potential Φ and the mass density ρ can be expressed as analytic functions. These provide an important touchstone by which more general numerical codes can be checked (i.e., does your code agree with a known analytic solution?)

Use the Poisson equation to derive the density ρ(r) associated with the spherically symmetric Hernquist potential
Φ(r) = −GM/(r+a)
where M is a characteristic mass and a is a length scale.
Try this by hand with pencil & paper before resorting to codes like Mathematica, as over-reliance on such tools often leads to a failure to understand their output.

Show enough steps to follow how you got from Φ to ρ.
What are the asymptotic limits of ρ at small and large radii?
I.e., what happens as r → 0 and r → ∞?
Think about what your answers mean physically. How does the density vary with radius? Where does it increase/decrease/reach a maximum?

Milky Way Scales

• the proper motion of Sgr A* is Ω_A* = 6.379 milliarcseconds/year (Reid & Brunthaler 2004), and
• our distance from the Galactic Center is R₀ = 8.122 kpc (GRAVITY collaboration 2018),
• the rotation curve of the Galaxy can be roughly approximated as V(R) = 229−1.7(R−R₀) km/s for 5 < R < 25 kpc (Eilers et al. 2019).

What are the values of

1. Our speed relative to Sgr A* (in km/s)?
2. The solar motion? (the peculiar speed of the sun relative to the local circular speed)
3. The Oort A & B constants?
4. The orbital frequency and the corresponing period of our orbit around the Milky Way?
5. The epicyclcic frequency and corresponding period?
6. Are stellar orbits closed ellipses? Use the orbital and epicyclic periods to explain why or why not.

Appropriate Galactic units for frequencies are km/s/kpc and orbital time periods in Myr.
See the the useful numbers page for help with unit conversions.

Local Dark Matter Mass and Density

1. The baryonic mass enclosed by the solar circle (r < R₀ = 8.122 kpc) is about 6 x 10¹⁰ M_☉. If that were all of the mass, what would the circular speed at the solar circle V(R₀) be?
2. The observed circular speed is V(R₀) = 229 km/s. What mass of dark matter is required to make up the difference from part (i)?
   [Hint: masses add linearly, so velocities add in quadrature.]
3. Is the Milky Way baryon dominated or dark matter dominated within the solar circle?
4. What is the average density of dark matter (in solar mass per cubic parsec)?
5. Use your result to estimate the mass of the dark matter within the solar system (a sphere out to Neptune's orbit of 30 AU).
   Do you expect the dark matter have a noticeable impact on planetary dynamics?
   [It may be helpful to express the dark mass in terms of planetary objects.]

[See the the useful numbers page for Newton's constant in convenient units.]

NFW halos

The dark matter halos formed in numerical structure formation simulations are often approximated as the "NFW halo" with density profile

ρ(r) = ρ_sr_s³ [r (r+r_s)²]⁻¹

where ρ_s is a characteristic density and r_s is a scale radius.
In the following, you may assume spherical symmetry.

1. Obtain an expression for the mass profile M(r).
2. Find the corresponding rotation curve V(r).
   It is conventional to define a "virial" radius R₂₀₀ inside of which the average density is 200 times the critical density of the universe, ρ_crit = 3H₀²/(8πG). The mass enclosed by this radius, M₂₀₀, is often taken as a working definition of the halo mass. The circular velocity of a test particle at the virial radius is V₂₀₀.
3. Show that V₂₀₀ = R₂₀₀h when V₂₀₀ is in km/s and R₂₀₀ is in kpc and h = H₀/(100 km/s/Mpc).
   
4. The relations above mean that the mass, radius, and circular speed of an NFW halo are interchangeable quantities: if you know one, you know the others. But the over-density of 200 seems arbitrary. Why not simply obtain the mass by integrating the density profile to infinity?

Exponential Disks

Integrate Σ(r) from 0 to 2π and r = 0 to x scale lengths (r = xR_d) to obtain an expression for the luminosity enclosed by x scale lengths.

• Is the total luminosity finite as x ⇒ ∞ infinity? If not, what is it?
• How many scale lengths contain half the total light? (This is known as the half-light radius, R_e.)
• Plot the cumulative enclosed luminosity L(< x).

Disk rotation curve

• Assuming M/L is constant with radius and making the fudge of a "spherical" disk, use the expression for L(x) to compute and plot V_d(r) for an exponential disk.
• ASTR 433 only: The rotation curve of a thin exponential disk is given by Binney & Tremaine eqn. 2-169 (1st edition) or 2.165 (2nd edition) and can be written
  
  Solve this numerically for a reasonable sampling of radii and plot it on the same diagram as the spherical approximation.
  Why are the two different? Which rotates faster?
• At what radius (in scale lengths) does the rotation curve of an exponential disk peak?
• ASTR 433 only: do this for the thin disk as well as the spherical approximation.

Maximum Disk

• What is the central mass density (in solar masses per square parsec) of a maximum disk with R_d = 5 kpc and V_c = 200 km/s?
  
• ASTR 433 only: do this for the thin exponential disk as well as the spherical approximation.