Part III: Galaxies & Cosmology | Chapter 3

Dark Matter Halos

The invisible scaffolding of the Universe: the NFW density profile, galaxy rotation curves, virial mass, and the concentration-mass relation

Overview

Dark matter constitutes approximately 85% of the matter in the Universe ($\Omega_{\text{DM}} \approx 0.25$), yet interacts with ordinary matter only through gravity. N-body simulations of structure formation in the CDM paradigm reveal that dark matter organizes into gravitationally bound halos with a universal density profile. These halos serve as the gravitational hosts for galaxies and galaxy clusters, determining their dynamics, lensing properties, and assembly history.

In this chapter we derive the NFW density profile and its properties, analyze galaxy rotation curves as evidence for dark matter, define the virial mass and radius, derive the concentration-mass relation, and discuss the small-scale challenges to CDM.

1. The Navarro-Frenk-White Density Profile

Navarro, Frenk, and White (1996, 1997) discovered from N-body simulations that CDM halos of all masses follow a universal density profile characterized by two parameters.

1.1 The NFW Profile

The NFW density profile is:

$$\boxed{\rho(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2}}$$

where $r_s$ is the scale radius and $\rho_s$ is the characteristic density. The profile has two distinct regimes:

Inner region ($r \ll r_s$):$\rho \propto r^{-1}$, a cuspy inner slope.

Outer region ($r \gg r_s$):$\rho \propto r^{-3}$, a steep outer falloff.

The logarithmic slope varies continuously: $d\ln\rho/d\ln r = -(1 + 2x)/(1+x)$where $x = r/r_s$, transitioning from $-1$ at the center to $-3$ at large radii. This is steeper than an isothermal sphere ($\rho \propto r^{-2}$) at large radii and shallower at the center.

1.2 Enclosed Mass

The mass enclosed within radius $r$ is obtained by integrating the density profile:

$$M(r) = 4\pi\rho_s r_s^3\left[\ln\left(\frac{r_s + r}{r_s}\right) - \frac{r}{r_s + r}\right] = 4\pi\rho_s r_s^3\left[\ln(1+x) - \frac{x}{1+x}\right]$$

Note that the total mass diverges logarithmically — the profile must be truncated at some outer boundary, conventionally taken to be the virial radius.

2. Virial Mass and the Concentration Parameter

The NFW profile is fully specified by two parameters, conveniently chosen as the virial mass and concentration.

2.1 Definition of Virial Quantities

The virial radius $R_{\text{vir}}$ is defined as the radius within which the mean density equals $\Delta_{\text{vir}}$ times the critical density of the Universe:

$$\boxed{M_{\text{vir}} = \frac{4}{3}\pi\Delta_{\text{vir}}\rho_{\text{crit}}R_{\text{vir}}^3}$$

The overdensity $\Delta_{\text{vir}}$ is given by the spherical collapse model; for the standard cosmology at $z = 0$,$\Delta_{\text{vir}} \approx 200$ (often written $M_{200}$and $R_{200}$). For the Milky Way,$M_{200} \approx 1.0 \times 10^{12}\,M_\odot$and $R_{200} \approx 200$ kpc.

2.2 The Concentration Parameter

The concentration is defined as the ratio of virial radius to scale radius:

$$\boxed{c = \frac{R_{\text{vir}}}{r_s}}$$

Given $c$ and $M_{\text{vir}}$, all other parameters are determined. The characteristic density is:

$$\rho_s = \frac{\Delta_{\text{vir}}}{3}\frac{c^3}{\ln(1+c) - c/(1+c)}\,\rho_{\text{crit}}$$

2.3 The Concentration-Mass Relation

N-body simulations reveal a systematic decrease in concentration with increasing halo mass:

$$\boxed{c(M,z) \approx 10\left(\frac{M}{10^{12}\,M_\odot}\right)^{-0.1}(1+z)^{-0.5}}$$

This reflects the assembly history: lower-mass halos collapse earlier when the Universe was denser, producing more concentrated profiles. Galaxy-mass halos have $c \approx 10\text{--}15$, while cluster-mass halos have $c \approx 4\text{--}7$.

3. Galaxy Rotation Curves

The most direct evidence for dark matter in galaxies comes from the flat rotation curves of spiral galaxies, first systematically measured by Vera Rubin, Kent Ford, and collaborators in the 1970s and 1980s.

3.1 Circular Velocity from Mass Distribution

For a spherically symmetric mass distribution, the circular velocity at radius $r$ is:

$$v_c(r) = \sqrt{\frac{GM(r)}{r}}$$

For the NFW profile:

$$v_c^2(r) = \frac{4\pi G\rho_s r_s^3}{r}\left[\ln(1 + r/r_s) - \frac{r/r_s}{1 + r/r_s}\right]$$

The maximum circular velocity occurs at $r \approx 2.16\,r_s$ and is related to the virial velocity by:

$$\frac{V_{\text{max}}}{V_{\text{vir}}} = \left(\frac{0.216\,c}{\ln(1+c) - c/(1+c)}\right)^{1/2}$$

3.2 The Rotation Curve Decomposition

The observed rotation curve of a galaxy is the quadrature sum of contributions from multiple mass components:

$$v_c^2(r) = v_{\text{bulge}}^2(r) + v_{\text{disk}}^2(r) + v_{\text{gas}}^2(r) + v_{\text{halo}}^2(r)$$

The key observation is that beyond the optical disk ($r > 3\text{--}5\,R_d$), where the baryonic contributions decline, the rotation curve remains flat or even rises slightly. This requires a dark matter halo with $M(r) \propto r$(i.e., $\rho \propto r^{-2}$) over the observed range.

4. Gravitational Lensing by Dark Matter Halos

Gravitational lensing provides a direct probe of the total mass distribution (luminous plus dark) of halos, independent of the dynamical state of the system.

4.1 Weak Lensing Shear

A mass distribution produces a convergence $\kappa$ (magnification) and shear $\gamma$ (distortion) on background galaxy images. For an NFW halo, the projected surface mass density (convergence) is:

$$\kappa(R) = \frac{\Sigma(R)}{\Sigma_{\text{crit}}}, \qquad \Sigma_{\text{crit}} = \frac{c^2}{4\pi G}\frac{D_s}{D_l D_{ls}}$$

where $D_s$, $D_l$, $D_{ls}$ are the angular diameter distances to the source, lens, and between lens and source. The tangential shear profile around an NFW halo declines as roughly $\gamma_t \propto R^{-1}$at large radii, enabling mass measurement from the shapes of background galaxies.

4.2 Galaxy-Galaxy Lensing

By stacking the shear signal around many foreground galaxies, galaxy-galaxy lensing measures the average halo mass profile. This technique has established that the total mass-to-light ratio increases with radius, confirming that galaxies are embedded in extended dark matter halos with $M/L \sim 100\text{--}300$ at the virial radius.

5. Small-Scale Challenges to CDM

While CDM is remarkably successful on large scales, several tensions exist at the scale of individual galaxies.

5.1 The Cusp-Core Problem

The NFW profile predicts a central cusp ($\rho \propto r^{-1}$), but observations of dwarf and low-surface-brightness galaxies often show constant-density cores ($\rho \propto r^0$). The Burkert profile provides a better fit to many observed rotation curves:

$$\rho_{\text{Burkert}}(r) = \frac{\rho_0}{(1 + r/r_0)(1 + (r/r_0)^2)}$$

Baryonic feedback (supernova-driven outflows that rapidly fluctuate the gravitational potential) can transform cusps into cores in dwarf galaxies, though whether this mechanism works for all systems remains debated.

5.2 Missing Satellites Problem

CDM simulations predict thousands of subhalos around a Milky Way-mass galaxy, but only$\sim 60$ satellite galaxies are observed. This discrepancy is largely resolved by understanding that most low-mass halos are too small to form stars (below the atomic hydrogen cooling limit of $\sim 10^8\,M_\odot$, and further suppressed by reionization). The discovery of ultra-faint dwarf galaxies by SDSS and DES has closed much of the gap.

5.3 The Diversity Problem

Observed rotation curves of galaxies with the same maximum circular velocity show surprisingly diverse inner profiles. Some rise steeply (cuspy), while others rise slowly (cored). CDM with baryonic feedback can produce some diversity, but whether it fully accounts for the observed range remains an active research question.

Applications

Dark Matter Direct Detection

The local dark matter density at the Sun's position is $\rho_\odot \approx 0.3\text{--}0.4$ GeV cm$^{-3}$($\approx 0.01\,M_\odot$ pc$^{-3}$). Experiments such as LZ, XENONnT, and PandaX-4T attempt to detect the nuclear recoils produced by dark matter particles (assumed to be WIMPs) passing through underground detectors. The local velocity distribution, modeled by the NFW halo profile with corrections for the disk and substructure, directly determines the expected recoil spectrum.

Cluster Mass Calibration

Galaxy cluster masses, essential for cosmological constraints, are calibrated using NFW profile fits to weak lensing data, X-ray observations of the hot intracluster medium, and the Sunyaev-Zel'dovich effect. The mass-observable scaling relations are a primary systematic uncertainty in cluster cosmology.

Substructure and Tidal Streams

CDM halos contain a rich population of subhalos — the remnants of smaller halos that were accreted during the hierarchical assembly process. The subhalo mass function follows $dN/dM_{\text{sub}} \propto M_{\text{sub}}^{-1.9}$, meaning most of the substructure is in small subhalos. Tidal interactions with the host halo strip matter from subhalos, creating long, thin stellar streams that trace the gravitational potential. The Gaia satellite has revealed dozens of stellar streams in the Milky Way halo. Gaps in these streams, caused by encounters with dark subhalos, provide a method to detect completely dark substructure with no baryonic content — a direct test of CDM predictions on the smallest scales.

Alternative Dark Matter Models

While CDM remains the standard paradigm, several alternative models address the small-scale challenges. Warm dark matter (WDM) with particle mass $m_{\text{WDM}} \sim 1\text{--}10$ keV suppresses structure below a free-streaming scale, reducing the number of small subhalos. Self-interacting dark matter (SIDM) with cross-section $\sigma/m \sim 1$ cm$^2$/g produces constant-density cores through thermalization in halo centers. Fuzzy dark matter (ultra-light axions with $m \sim 10^{-22}$ eV) produces wave-like behavior on galactic scales. Each model makes distinct predictions for halo structure that can be tested observationally.

Strong Gravitational Lensing

Galaxy clusters with masses $\sim 10^{15}\,M_\odot$ produce dramatic strong gravitational lensing, distorting background galaxies into giant arcs and multiple images. The Einstein radius for a cluster at typical redshift is $\theta_E \sim 20''\text{--}40''$. Modeling the positions and shapes of multiple images constrains the cluster mass distribution with exquisite precision, revealing the NFW profile and substructure. The Hubble Frontier Fields program used massive clusters as cosmic telescopes, magnifying galaxies at $z > 6$ by factors of 10–100, enabling the study of the earliest galaxies.

Historical Notes

The dark matter problem originated with Fritz Zwicky's 1933 observation that galaxies in the Coma Cluster move too fast to be gravitationally bound by visible matter alone. He coined the term "dunkle Materie" (dark matter). The rotation curve evidence was established by Vera Rubin and Kent Ford (1970) for M31 and extended to many galaxies in the following decade. Jeremiah Ostriker, James Peebles, and Amos Yahil (1974) provided theoretical arguments for massive dark halos. The NFW profile was published in 1996 and 1997, based on cosmological N-body simulations. The CDM paradigm, with dark matter composed of cold, collisionless particles, remains the standard model despite decades of experimental searches that have not yet identified the dark matter particle.

The Bullet Cluster (1E 0657-558), observed in 2006, provided some of the most compelling evidence for dark matter as a particle rather than a modification of gravity. In this system of two colliding galaxy clusters, the gravitational mass (mapped by weak lensing) is offset from the visible baryonic mass (hot gas mapped by X-ray emission). The dark matter component passed through the collision largely unimpeded, while the gas was slowed by ram pressure — precisely the behavior expected for weakly interacting particles and inconsistent with modified gravity theories like MOND. This observation, combined with the CMB power spectrum, galaxy cluster abundances, and BBN constraints, firmly establishes that dark matter constitutes $\sim 85\%$ of the matter content of the Universe.

The question of what dark matter is made of remains one of the deepest open problems in physics. The leading particle candidate for decades has been the Weakly Interacting Massive Particle (WIMP), which naturally produces the observed dark matter abundance through thermal freeze-out in the early Universe (the "WIMP miracle"). However, decades of increasingly sensitive direct detection experiments (LUX, XENON, PandaX, and now LZ and XENONnT) have not found WIMPs, pushing the allowed parameter space to ever-smaller cross-sections. Alternative candidates include axions (searched for by ADMX, ABRACADABRA, and CASPEr), sterile neutrinos (constrained by X-ray observations), and primordial black holes (constrained by microlensing and gravitational wave observations). The identification of the dark matter particle would represent one of the greatest discoveries in the history of science, simultaneously solving problems in cosmology, particle physics, and astrophysics.

The Milky Way's dark matter halo has been extensively studied through stellar kinematics, satellite galaxy orbits, and stellar stream modeling. Current estimates give a total halo mass of $M_{200} \approx 1.0 \times 10^{12}\,M_\odot$ with a concentration $c \approx 10$. The local dark matter density at the Sun's position ($R = 8$ kpc) is measured to be $\rho_\odot \approx 0.3\text{--}0.4$ GeV cm$^{-3}$, a value critical for interpreting direct detection experiments. The shape of the inner halo is constrained by the precession of stellar streams and the orbits of the Sagittarius dwarf galaxy. The Gaia satellite has transformed this field by providing 6D phase-space information for over a billion stars, enabling dynamical mass modeling of unprecedented precision and revealing the rich substructure of the stellar halo as a fossil record of the Milky Way's hierarchical assembly.

Halo Shapes and Dynamical Friction

CDM halos are not spherical. N-body simulations show they are triaxial ellipsoids with typical axis ratios $c/a \approx 0.5\text{--}0.7$ (minor-to-major axis). More massive halos tend to be more prolate (elongated), while the inner regions are rounder due to the dissipative infall of baryons. The halo shape affects gravitational lensing predictions and the orbits of satellite galaxies.

Satellite galaxies and subhalos experience dynamical friction as they move through the dark matter background. The Chandrasekhar dynamical friction formula gives the deceleration:

$$\frac{d\mathbf{v}}{dt} = -\frac{4\pi G^2 M_{\text{sat}} \rho \ln\Lambda}{v^3}\left[\text{erf}(X) - \frac{2X}{\sqrt{\pi}}e^{-X^2}\right]\mathbf{v}$$

where $X = v/(\sqrt{2}\sigma)$, $\ln\Lambda$ is the Coulomb logarithm, and $\rho$ is the local dark matter density. Dynamical friction causes massive satellites to spiral inward and merge with the central galaxy on a timescale $t_{\text{df}} \propto M_{\text{host}}/M_{\text{sat}}$. For the Magellanic Clouds ($M_{\text{LMC}} \sim 10^{11}\,M_\odot$), the dynamical friction timescale is comparable to the Hubble time, consistent with their being on their first infall into the Milky Way.

The Splashback Radius and Halo Boundaries

The splashback radius $R_{\text{sp}}$ marks the physical boundary of a dark matter halo, defined as the radius where recently accreted particles reach the apocentre of their first orbit after infall. At $R_{\text{sp}}$, the density profile exhibits a sharp steepening (a caustic), producing a characteristic feature in the logarithmic slope:

$$\frac{d\ln\rho}{d\ln r}\bigg|_{R_{\text{sp}}} \lesssim -3$$

The splashback radius is typically $R_{\text{sp}} \approx 1.5\text{--}2.5\,R_{200}$and depends on the mass accretion rate: rapidly accreting halos have smaller$R_{\text{sp}}/R_{200}$ because infalling material has less time to orbit outward. This feature has been detected observationally through galaxy number density profiles around clusters (using SDSS and DES data) and through weak lensing, providing a physically motivated definition of halo boundaries that is independent of the arbitrary overdensity threshold used in $R_{200}$ or $R_{500}$.

Computational Exploration

The following simulation computes NFW density profiles for different halo masses, generates rotation curve decompositions, illustrates the concentration-mass relation, and compares cuspy (NFW) vs cored (Burkert) profiles.

NFW Profiles, Rotation Curves, and the Concentration-Mass Relation

Python

script.py230 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a0a')
for ax in axes.flat:
    ax.set_facecolor('#0f0f0f')
    ax.tick_params(colors='white')
    ax.xaxis.label.set_color('white')
    ax.yaxis.label.set_color('white')
    ax.title.set_color('white')
    for spine in ax.spines.values():
        spine.set_color('#333')

# Constants
G = 4.302e-3  # pc (km/s)^2 / M_sun
G_kpc = G * 1e3  # kpc (km/s)^2 / M_sun
H_0 = 70.0  # km/s/Mpc
rho_crit = 3 * H_0**2 / (8 * np.pi * G * 1e6)  # M_sun / kpc^3
# rho_crit ~ 1.36e-7 M_sun / pc^3 = 1.36e2 M_sun / kpc^3
rho_crit = 277.5 * (H_0/100)**2  # M_sun / kpc^3 / h^2 ... use simplified
rho_crit = 1.4e-7 * 1e9  # M_sun / kpc^3 ~ 140 M_sun/kpc^3

# --- Plot 1: NFW density profiles for different halo masses ---
ax1 = axes[0, 0]

Delta_vir = 200

def nfw_params(M_vir, c):
    """Compute NFW parameters from virial mass and concentration"""
    R_vir = (3 * M_vir / (4 * np.pi * Delta_vir * rho_crit))**(1.0/3.0)  # kpc
    r_s = R_vir / c
    rho_s = M_vir / (4 * np.pi * r_s**3 * (np.log(1+c) - c/(1+c)))
    return R_vir, r_s, rho_s

halo_configs = [
    (1e10, 15, 'Dwarf galaxy', '#44bbff'),
    (1e12, 10, 'Milky Way', '#ffa500'),
    (1e14, 5, 'Galaxy cluster', '#ff4444'),
]

for M_vir, c, label, col in halo_configs:
    R_vir, r_s, rho_s = nfw_params(M_vir, c)
    r_range = np.logspace(np.log10(r_s * 0.01), np.log10(R_vir * 2), 200)
    x = r_range / r_s
    rho = rho_s / (x * (1 + x)**2)

ax1.loglog(r_range, rho, color=col, linewidth=2.5, label=f'{label} (M={M_vir:.0e})')
    ax1.axvline(x=R_vir, color=col, linestyle=':', linewidth=1, alpha=0.4)
    ax1.axvline(x=r_s, color=col, linestyle='--', linewidth=1, alpha=0.3)

# Reference slopes
r_ref = np.logspace(-1, 1, 50)
ax1.loglog(r_ref, 1e7 * r_ref**(-1), 'w--', linewidth=1, alpha=0.3, label=r'$\rho propto r^{-1}$')
r_ref2 = np.logspace(2, 3, 50)
ax1.loglog(r_ref2, 5e13 * r_ref2**(-3), 'w:', linewidth=1, alpha=0.3, label=r'$\rho propto r^{-3}$')

ax1.set_xlabel('Radius (kpc)', fontsize=12)
ax1.set_ylabel('Density (M_sun / kpc^3)', fontsize=12)
ax1.set_title('NFW Density Profiles', fontsize=14, fontweight='bold')
ax1.legend(fontsize=8, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax1.set_xlim(0.1, 3000)
ax1.set_ylim(1e-1, 1e10)
ax1.grid(True, alpha=0.15)

# --- Plot 2: Rotation curve decomposition (MW-like) ---
ax2 = axes[0, 1]

r_kpc = np.linspace(0.5, 50, 200)

# Bulge: Hernquist profile
M_bulge = 1e10  # M_sun
a_bulge = 0.5   # kpc
v_bulge = np.sqrt(G_kpc * M_bulge * r_kpc / (r_kpc + a_bulge)**2)

# Disk: exponential disk (Freeman approximation)
M_disk = 5e10  # M_sun
R_d = 3.0      # kpc
y = r_kpc / (2 * R_d)
# Approximate: v_disk^2 ~ GM_disk/(2R_d) * y^2 * (I0*K0 - I1*K1)
# Simplified: use empirical fit
v_disk = np.sqrt(0.5 * G_kpc * M_disk / R_d * y**2 * 3.36 * np.exp(-y) * (1 + 0.5*y))
v_disk = np.minimum(v_disk, 180)  # cap

# Gas disk
M_gas = 8e9
v_gas = np.sqrt(G_kpc * M_gas * r_kpc / (r_kpc + 5.0)**2) * 0.5

# Dark matter halo (NFW)
M_halo = 1e12
c_halo = 10
R_vir_h, r_s_h, rho_s_h = nfw_params(M_halo, c_halo)
x_h = r_kpc / r_s_h
M_enc_halo = 4 * np.pi * rho_s_h * r_s_h**3 * (np.log(1 + x_h) - x_h/(1 + x_h))
v_halo = np.sqrt(G_kpc * M_enc_halo / r_kpc)

# Total
v_total = np.sqrt(v_bulge**2 + v_disk**2 + v_gas**2 + v_halo**2)

ax2.plot(r_kpc, v_bulge, '--', color='#ff4444', linewidth=1.5, alpha=0.7, label='Bulge')
ax2.plot(r_kpc, v_disk, '--', color='#ffa500', linewidth=1.5, alpha=0.7, label='Disk')
ax2.plot(r_kpc, v_gas, '--', color='#44ff44', linewidth=1.5, alpha=0.7, label='Gas')
ax2.plot(r_kpc, v_halo, '--', color='#44bbff', linewidth=1.5, alpha=0.7, label='DM Halo')
ax2.plot(r_kpc, v_total, '-', color='#ffa500', linewidth=3, label='Total')

# Simulated observations
rng = np.random.default_rng(42)
r_obs = np.linspace(2, 45, 25)
v_obs_interp = np.interp(r_obs, r_kpc, v_total)
v_obs = v_obs_interp + rng.normal(0, 5, len(r_obs))
ax2.errorbar(r_obs, v_obs, yerr=8, fmt='o', color='#ffdd00', markersize=4, capsize=2,
             ecolor='#ffdd00', alpha=0.7, label='Observed')

ax2.set_xlabel('Radius (kpc)', fontsize=12)
ax2.set_ylabel('Circular Velocity (km/s)', fontsize=12)
ax2.set_title('Rotation Curve Decomposition (MW-like)', fontsize=14, fontweight='bold')
ax2.legend(fontsize=8, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white', ncol=2)
ax2.set_xlim(0, 50)
ax2.set_ylim(0, 300)
ax2.grid(True, alpha=0.15)

# Keplerian reference
r_kep = np.linspace(10, 50, 50)
v_kep = v_total[100] * np.sqrt(10 / r_kep)
ax2.plot(r_kep, v_kep, ':', color='#888', linewidth=1.5, alpha=0.5)
ax2.text(40, 100, 'Keplerian', color='#888', fontsize=9, rotation=-15)

# --- Plot 3: Concentration-mass relation ---
ax3 = axes[1, 0]

M_range = np.logspace(8, 16, 200)

redshifts_cm = [0, 0.5, 1, 2]
colors_cm = ['#ffa500', '#ffdd00', '#44bbff', '#ff4444']

for z_val, col in zip(redshifts_cm, colors_cm):
    c_val = 10 * (M_range / 1e12)**(-0.1) * (1 + z_val)**(-0.5)
    ax3.semilogx(M_range, c_val, color=col, linewidth=2.5, label=f'z = {z_val}')

# Scatter
M_scatter = np.logspace(9, 15, 50)
c_scatter = 10 * (M_scatter / 1e12)**(-0.1)
c_noise = c_scatter * np.exp(rng.normal(0, 0.15, len(M_scatter)))
ax3.scatter(M_scatter, c_noise, s=15, color='#ffa500', alpha=0.4, zorder=0)

# Mark MW
ax3.plot(1e12, 10, '*', color='#ffdd00', markersize=15, zorder=10)
ax3.text(1.5e12, 11, 'Milky Way', color='#ffdd00', fontsize=10)

ax3.set_xlabel('Halo Mass (M_sun)', fontsize=12)
ax3.set_ylabel('Concentration c = R_vir / r_s', fontsize=12)
ax3.set_title('Concentration-Mass Relation', fontsize=14, fontweight='bold')
ax3.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax3.set_xlim(1e8, 1e16)
ax3.set_ylim(0, 30)
ax3.grid(True, alpha=0.15)

# --- Plot 4: NFW vs Burkert (cusp vs core) ---
ax4 = axes[1, 1]

r_inner = np.logspace(-2, 1.5, 200)  # r/r_s

# NFW
rho_NFW = 1 / (r_inner * (1 + r_inner)**2)

# Burkert (with same rho_0 and r_0 = r_s)
rho_Burkert = 1 / ((1 + r_inner) * (1 + r_inner**2))

# Isothermal sphere
rho_ISO = 1 / r_inner**2

# Einasto (with alpha = 0.17)
alpha_ein = 0.17
rho_Einasto = np.exp(-2/alpha_ein * (r_inner**alpha_ein - 1))

ax4.loglog(r_inner, rho_NFW, color='#ffa500', linewidth=3, label='NFW (cusp)')
ax4.loglog(r_inner, rho_Burkert, color='#44bbff', linewidth=3, label='Burkert (core)')
ax4.loglog(r_inner, rho_Einasto, '--', color='#ff4444', linewidth=2, label='Einasto')
ax4.loglog(r_inner, rho_ISO, ':', color='#888', linewidth=1.5, alpha=0.5, label='Isothermal')

# Highlight core region
ax4.axvspan(0.01, 0.3, alpha=0.08, color='#ffa500')
ax4.text(0.05, 2e-3, 'Core/cusp\nregion', color='#aaa', fontsize=10)

ax4.set_xlabel('r / r_s', fontsize=12)
ax4.set_ylabel('Density / rho_s', fontsize=12)
ax4.set_title('Dark Matter Density Profiles Compared', fontsize=14, fontweight='bold')
ax4.legend(fontsize=10, facecolor='#1a1a1a', edgecolor='#444', labelcolor='white')
ax4.set_xlim(0.01, 30)
ax4.set_ylim(1e-5, 1e3)
ax4.grid(True, alpha=0.15)

# Print summary
print("=" * 65)
print("DARK MATTER HALOS: KEY RESULTS")
print("=" * 65)

print("\nNFW halo parameters:")
for M_vir, c, label, _ in halo_configs:
    R_vir, r_s, rho_s = nfw_params(M_vir, c)
    V_vir = np.sqrt(G_kpc * M_vir / R_vir)
    print(f"  {label:20s}: M_vir = {M_vir:.0e} M_sun, R_vir = {R_vir:.0f} kpc, "
          f"r_s = {r_s:.1f} kpc, V_vir = {V_vir:.0f} km/s")

print(f"\nCritical density: rho_crit = {rho_crit:.0f} M_sun/kpc^3")

print("\nConcentration-mass relation (z=0):")
for M in [1e10, 1e12, 1e14, 1e15]:
    c = 10 * (M / 1e12)**(-0.1)
    print(f"  M = {M:.0e} M_sun: c = {c:.1f}")

print("\nRotation curve at solar position (R=8 kpc):")
idx_sun = np.argmin(np.abs(r_kpc - 8))
print(f"  v_bulge = {v_bulge[idx_sun]:.0f} km/s")
print(f"  v_disk  = {v_disk[idx_sun]:.0f} km/s")
print(f"  v_halo  = {v_halo[idx_sun]:.0f} km/s")
print(f"  v_total = {v_total[idx_sun]:.0f} km/s")

print("\nLocal DM density (r=8 kpc, MW halo):")
x_sun = 8.0 / r_s_h
rho_local = rho_s_h / (x_sun * (1 + x_sun)**2)
print(f"  rho_DM = {rho_local:.2e} M_sun/kpc^3")
print(f"  rho_DM ~ {rho_local * 1e-9 * 1.989e30 / (3.086e16)**3 / 1.783e-27:.2f} GeV/cm^3")

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
print("\nPlot saved to output.png")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Chapter Summary

Dark matter halos follow the universal NFW density profile $\rho \propto r^{-1}(1+r/r_s)^{-2}$, transitioning from a $r^{-1}$ cusp to a $r^{-3}$ outer falloff. The profile is specified by two parameters: the virial mass $M_{\text{vir}}$ and concentration $c = R_{\text{vir}}/r_s$.

Galaxy rotation curves provide the most direct evidence for dark matter, remaining flat far beyond the optical disk where baryonic matter is negligible. The rotation curve decomposition reveals that dark matter dominates the mass budget at large radii.

The concentration-mass relation $c \propto M^{-0.1}$ reflects hierarchical assembly. Small-scale challenges (cusp-core, missing satellites, diversity) can be partially addressed by baryonic feedback effects but remain active areas of research.

← Active Galactic Nuclei Large-Scale Structure →