Data Analysis problem set
OlympiadIOAA 2017 | Phuket, Thailand | Data Analysis D2 | 75 pt

2017 D2 | Dark matter in a low surface brightness galaxy

A low surface brightness (LSB) galaxy is a diffuse galaxy whose surface brightness, seen from Earth, is at least one magnitude fainter than the night sky. Some of its matter is baryonic (neutral hydrogen gas and stars), but most of its mass is invisible dark matter.

The table gives data for the LSB galaxy UGC 4325, assumed to be edge-on. At each distance rr from the centre we have λobs\lambda_{\rm obs}, the observed wavelength of the Hα line (the Hubble expansion has already been removed); VgasV_{\rm gas}, the rotation speed the gas alone would produce (from HI surface densities); and V∗V_*, the rotation speed the stars alone would produce (from R-band photometry).

The dark matter mass inside radius rr is defined from the dark matter rotation speed VDMV_{\rm DM}:

MDM(r)=rVDM2G,Vobs2=Vgas2+V∗2+VDM2.M_{\rm DM}(r) = \frac{rV_{\rm DM}^2}{G}, \qquad V_{\rm obs}^2 = V_{\rm gas}^2 + V_*^2 + V_{\rm DM}^2.

The dark matter density is modelled by a cored profile, which gives

ρDM(r)=ρ01+(r/rC)2,MDM(r)=4πρ0rC2[r−rCarctan⁡rrC],\rho_{\rm DM}(r) = \frac{\rho_0}{1 + (r/r_C)^2}, \qquad M_{\rm DM}(r) = 4\pi\rho_0r_C^2\left[r - r_C\arctan\frac{r}{r_C}\right],

where ρ0\rho_0 is the central density and rCr_C the core radius.

(D2.1) The laboratory wavelength of Hα is λemit=656.281\lambda_{\rm emit} = 656.281 nm. Compute VobsV_{\rm obs} and VDMV_{\rm DM} (km/s) at each rr, and the dynamical mass M(r)M(r) and dark matter mass MDM(r)M_{\rm DM}(r) in solar masses. (21)

(D2.2) Plot VobsV_{\rm obs}, VDMV_{\rm DM}, VgasV_{\rm gas} and V∗V_* against rr with smooth curves. Order the contributions to the observed velocity from largest to smallest. (16)

(D2.3) Use a data point at small rr and one at large rr to estimate ρ0\rho_0 and rCr_C. For large xx, arctan⁡x≈π/2\arctan x\approx\pi/2; for small xx, arctan⁡x≈x−x3/3\arctan x\approx x - x^3/3. (7)

(D2.4) Linearise the mass equation so that ρ0\rho_0 can be found from a straight-line fit. Plot the graph and find ρ0\rho_0 in M⊙ kpc−3M_\odot\,{\rm kpc^{-3}}. If you could not find rCr_C, use rC=3.2r_C = 3.2 kpc. (19)

(D2.5) Compute ln⁡ρDM(r)\ln\rho_{\rm DM}(r) and plot it against rr. (12)

UGC 4325 rotation data
r (kpc)λ_obs (nm)V_gas (km/s)V_* (km/s)
0.70656.3712.8720.97
1.40656.4316.7532.22
2.09656.46414.1440.91
2.79656.47520.1846.75
3.49656.47824.0850.10
4.89656.48428.0847.94
6.25656.48129.2545.47
7.10656.48127.0347.78
9.03656.48225.9045.32
12.05656.48221.0342.30

Saved in this browser only, no account needed.

Reading resources

Material that covers the physics behind this problem.