The Milky Way, Galactic Dynamics & Spiral Structure

The Milky Way is a barred spiral galaxy classified as SBbc in the Hubble-de Vaucouleurs sequence. It contains approximately $(1-4) \times 10^{11}$ stars and spans an optical diameter of roughly $30\text{ kpc}$.

Structural Components of the Galaxy

Sagittarius A*: The Central Supermassive Black Hole

At the dynamic center of the Milky Way lies Sagittarius A* (Sgr A*). Infrared speckle and adaptive optics observations (Ghez, Genzel, Nobel Prize 2020) tracked the Keplerian orbital trajectories of individual stars ('S-stars') over decades. Star S2 completes an elliptical orbit ($e = 0.884$) with period $P = 16.05\text{ years}$, approaching within $r_{\text{peri}} = 120\text{ AU} \approx 1400 R_s$ of Sgr A* at speeds exceeding $7700\text{ km s}^{-1}$ ($2.6\% c$). Applying Kepler's Third Law proves that a mass of:

is packed inside a volume smaller than the orbit of Mercury, definitively proving the existence of a central supermassive black hole, as confirmed directly by the Event Horizon Telescope (EHT) shadow image in 2022.

To analyze structures within our own galaxy from our internal vantage point, the IAU established the Galactic Coordinate System.

Galactic Coordinates: Longitude ($l$) and Latitude ($b$)

• Galactic Center: Defined as $(l = 0^\circ, b = 0^\circ)$, located in the constellation Sagittarius at Right Ascension $\alpha = 17^\text{h} 45^\text{m} 40^\text{s}$, Declination $\delta = -29^\circ 00' 28''$ (J2000).
• Galactic Longitude ($l$): Measured eastward along the Galactic plane from the Galactic Center ($0^\circ \le l < 360^\circ$). $l = 90^\circ$ points in the direction of Galactic rotation (Cygnus), $l = 180^\circ$ is the Galactic Anticenter (Auriga), and $l = 270^\circ$ is Vela.
• Galactic Latitude ($b$): Measured perpendicular to the Galactic plane ($-90^\circ \le b \le +90^\circ$), with the North Galactic Pole (NGP) at $(b = +90^\circ)$ in Coma Berenices.

Interstellar Dust Extinction & Reddening

Microscopic interstellar dust grains ($a \sim 0.01 - 0.2\ \mu\text{m}$, composed of silicates, graphite, and polycyclic aromatic hydrocarbons) absorb and scatter optical light. The scattering cross-section obeys a power law $\sigma_{\text{scat}} \propto \lambda^{-1}$ (Mie scattering regime). Consequently, blue light is scattered far more than red light, causing interstellar reddening.

Extinction along the line of sight reduces apparent stellar flux according to Beer's law: $F = F_0 e^{-\tau_\lambda} = F_0 10^{-0.4 A_\lambda}$. In the Galactic disk plane ($b \approx 0^\circ$), average visual extinction is:

Toward the Galactic Center ($d \approx 8.2\text{ kpc}$), total visual extinction is $A_V \approx 30\text{ magnitudes}$! This attenuates optical light by a factor of $10^{-0.4 \times 30} = 10^{-12}$ (one-trillionth), completely blinding optical telescopes to the inner galaxy and creating the 'Zone of Avoidance'.

The 21-cm Hyperfine Neutral Hydrogen (H I) Line

In 1944, Hendrik van de Hulst realized that neutral atomic hydrogen (H I) can be observed at radio wavelengths. Neutral hydrogen consists of one proton and one electron. In the ground state ($1s$), the nuclear magnetic dipole moment of the proton and the spin magnetic dipole moment of the electron can be aligned (parallel, higher energy) or anti-aligned (antiparallel, lower energy).

The quantum mechanical spin-flip transition releases an energy of $\Delta E = 5.874 \times 10^{-6}\text{ eV}$, corresponding to a photon wavelength and frequency of:

The transition is forbidden by electric dipole selection rules, proceeding via magnetic dipole transition with an Einstein $A$ coefficient of $A_{10} = 2.87 \times 10^{-15}\text{ s}^{-1}$ (spontaneous radiative lifetime $\tau = 1/A_{10} \approx 11\text{ million years}$).

However, because interstellar hydrogen is vastly abundant, collisions easily maintain the population. Most importantly, at $\lambda = 21\text{ cm}$, interstellar dust extinction is utterly negligible ($A_{21} \approx 10^{-4} A_V$). 21-cm radio telescopes peer directly through the entire Galactic disk, measuring Doppler velocity profiles that reveal the spiral arms of our galaxy.

Stars in the disk of the Milky Way do not rotate as a rigid body. Instead, they exhibit differential galactic rotation, orbiting the center with an angular velocity $\Omega(R)$ that depends on galactocentric radius $R$.

Derivation of Oort's Formulae

Consider the Sun at galactocentric radius $R_0 \approx 8.2\text{ kpc}$ moving in a circular orbit with speed $V_0 = \Omega_0 R_0 \approx 220\text{ km s}^{-1}$. A star located at distance $d \ll R_0$ and Galactic longitude $l$ moves at radius $R$ with circular speed $V = \Omega R$.

The observed radial velocity $v_r$ and transverse velocity $v_t$ of the star relative to the Local Standard of Rest (LSR) are:

Using the sine and cosine laws for the triangle formed by the Galactic Center, Sun, and Star: $R \sin\alpha = R_0 \cos l - d$ and $R \cos\alpha = R_0 \sin l$:

Expanding $\Omega(R)$ in a Taylor series around $R = R_0$ for small distances ($d \ll R_0$), where $R - R_0 \approx -d \cos l$:

Substituting into the velocity equations yields Oort's Kinematic Equations:

where the Oort Constants $A$ and $B$ are defined as:

Direct subtraction and addition yield fundamental kinematic quantities:

Modern astrometric determinations (Gaia DR3) yield:

Epicyclic Approximation of Stellar Orbits

Real stellar orbits in the disk are not perfectly circular. In Lindblad's epicyclic approximation, a star executes a small elliptical retrograde oscillation (an epicycle) centered on a guiding center that orbits circularly at radius $R_0$. The radial oscillation frequency is the epicyclic frequency ($\kappa$):

In the solar neighborhood, $\kappa_0 = \sqrt{-4 (-11.9)(27.2)} = \sqrt{1294.7} \approx 36.0\text{ km s}^{-1}\text{kpc}^{-1}$. The radial oscillation period is $P_\kappa = 2\pi / \kappa_0 \approx 171\text{ Myr}$. Because $\kappa_0 \ne \Omega_0$, stellar orbits do not close, tracing rosettes in the Galactic plane.

Newtonian mechanics predicts how orbital speeds must behave in any self-gravitating disk system. Measuring the actual circular velocity profile $V(R)$ of the Milky Way and external spiral galaxies provides the most decisive empirical evidence for the existence of non-baryonic dark matter.

The Keplerian Prediction vs The Observation

Consider a galaxy whose mass is entirely luminous (baryonic stars and gas concentrated in a central bulge and exponential disk with scale length $R_d \sim 3\text{ kpc}$). At radii well beyond the visible disk ($R \gg R_d$), the enclosed mass is constant: $M(R) \to M_{\text{lum}} = \text{const}$. Balancing gravitational attraction and centripetal acceleration:

Just as the orbital speeds of planets in the Solar System fall off as $V \propto r^{-1/2}$ according to Kepler's Third Law, astronomers expected galactic rotation speeds to decline sharply beyond the optical edge ($R \sim 15\text{ kpc}$).

In the 1970s, Vera Rubin and Kent Ford (optical spectroscopy of H II regions) alongside Morton Roberts and Albert Bosma (21-cm radio observations of H I gas) measured rotation curves out to several times the visible optical radius. Shockingly, they discovered that:

The rotation curve does not decline!

Inferring the Dark Matter Halo Density Profile

If $V(R) = V_0 = \text{constant}$, the enclosed mass must grow linearly with radius:

Since the enclosed mass is $M(R) = \int_0^R 4\pi r^2 \rho(r) dr$, differentiating with respect to $R$ yields:

This $R^{-2}$ density profile corresponds to an isothermal sphere. Because the visible stars and gas drop off exponentially ($\rho_{\text{lum}} \propto e^{-R/R_d}$), the mass must be dominated by a vast, non-luminous, roughly spherical dark matter halo.

Cosmological Halo Profiles: The NFW Model

High-resolution cosmological $N$-body simulations of collisionless Cold Dark Matter (Navarro, Frenk, & White 1996) reveal a universal halo density profile:

where $r_s$ is a characteristic scale radius. The profile exhibits a central cusp ($\rho \propto r^{-1}$ for $r \ll r_s$), transitions to isothermal behavior ($\rho \propto r^{-2}$ near $r \sim r_s$), and steepens to $\rho \propto r^{-3}$ at large distances ($r \gg r_s$). The dark matter halo of the Milky Way extends out to $R_{\text{vir}} \approx 200\text{ kpc}$, comprising $(1.0-1.5) \times 10^{12} M_\odot$—over $85-90\%$ of the total mass of our galaxy!

Spiral arms are the most visually stunning features of disk galaxies. However, explaining their persistent survival posed a fundamental theoretical crisis known as the winding dilemma.

The Winding Dilemma

If spiral arms were rigid, physical material structures composed of the same stars and gas over time, differential galactic rotation would destroy them. Because stars at smaller radii orbit faster than stars at larger radii ($d\Omega/dR < 0$), an initially radial material arm would be wound into a tight spiral. The number of turns $n$ wound after time $t$ is:

Over the $10\text{ Gyr}$ lifetime of the Milky Way, differential rotation would wind the arms into more than 50 tightly wound turns, obliterating the open 2- and 4-arm patterns universally observed. Spiral arms cannot be material entities!

The Lin-Shu Density Wave Theory

In 1964, C.C. Lin and Frank Shu resolved the dilemma by proposing that spiral arms are quasi-stationary density waves—gravitational perturbation waves that propagate through the stellar and gaseous disk. The spiral pattern rotates as a rigid shape with a constant angular pattern speed $\Omega_p$, while individual stars and gas clouds orbit at their local speeds $\Omega(R)$, continually entering, passing through, and exiting the wave.

Resonances in the Galactic Disk

A star oscillating with epicyclic frequency $\kappa(R)$ encounters a spiral pattern with $m$ arms at an apparent frequency $m(\Omega - \Omega_p)$. Resonances occur when this driving frequency matches the natural epicyclic frequency:

• Corotation Resonance (CR): Where the stars rotate at exactly the same speed as the spiral pattern: $\Omega(R_{\text{CR}}) = \Omega_p$. Inside corotation ($R < R_{\text{CR}}$), stars overtake the spiral wave. Outside corotation ($R > R_{\text{CR}}$), the pattern sweeps past slower-moving stars.
• Inner Lindblad Resonance (ILR): $\Omega_p = \Omega(R) - \frac{\kappa(R)}{m}$. Stars complete two epicycles per encounter with a 2-arm ($m=2$) pattern.
• Outer Lindblad Resonance (OLR): $\Omega_p = \Omega(R) + \frac{\kappa(R)}{m}$.

Self-consistent spiral patterns can only survive between the ILR and OLR.

Star Formation in Spiral Arms

As cold interstellar gas clouds orbit supersonic relative to the pattern ($v_{\perp} = [\Omega(R) - \Omega_p] R \sin i > c_s$), they slam into the gravitational potential minimum of the density wave. The sudden deceleration produces a sharp hydrodynamic galactic shock wave. The gas is compressed by factors of $5-10$, driving the local density above the Jeans threshold ($M > M_J$) and triggering the birth of massive, short-lived O and B stars. Because OB stars live only $\sim 10-30\text{ Myr}$, they die near their birth sites, lighting up the trailing edges of the density wave with glowing blue star clusters and pink H II emission nebulae.

The interstellar medium is permeated by relativistic charged particles known as cosmic rays and a pervasive, coherent galactic magnetic field.

Composition & Energy Spectrum of Cosmic Rays

Cosmic rays consist of relativistic atomic nuclei ($89\%$ protons, $10\%$ alpha particles, $1\%$ heavier elements up to uranium) and relativistic electrons ($1\%$). Their energy spectrum spans an immense range from $10^9\text{ eV}$ ($1\text{ GeV}$) to $> 10^{20}\text{ eV}$ ($100\text{ EeV}$), governed by a broken power law:

• For $10^9\text{ eV} < E < 3 \times 10^{15}\text{ eV}$, $s \approx 2.7$.
• At the 'Knee' ($E \approx 3 \times 10^{15}\text{ eV} = 3\text{ PeV}$), the spectrum steepens to $s \approx 3.1$, marking the maximum confinement and acceleration energy of Galactic supernova remnants.
• At the 'Ankle' ($E \approx 3 \times 10^{18}\text{ eV} = 3\text{ EeV}$), the spectrum flattens back to $s \approx 2.6$, signaling a transition to Ultra-High Energy Cosmic Rays (UHECRs) of extragalactic origin.

Diffusive Shock Acceleration: The First-Order Fermi Mechanism

Enrico Fermi (1949, 1954) proposed that cosmic rays gain energy by scattering off magnetized plasma clouds. In a supernova blast wave expanding into the ISM at shock velocity $v_s$, particles scatter elastically off magnetic turbulence on both sides of the shock front. Because the downstream gas moves toward the upstream gas at relative speed $u = \frac{3}{4} v_s$, every round-trip crossing across the shock is a head-on collision. The fractional energy gain per cycle is:

Because the energy gain is linear in $v_s/c$ (first-order) and particles have a finite escape probability $P_{\text{esc}}$ per crossing, repeated shock cycling naturally produces an exact power-law distribution $dN/dE \propto E^{-2}$, which steepens to $E^{-2.7}$ as cosmic rays diffuse out of the Galactic magnetic trap.

Galactic Magnetic Fields & Synchrotron Radiation

The Milky Way hosts a large-scale magnetic field of strength $B \approx 3-6\ \mu\text{G}$ ($0.3-0.6\text{ nT}$), ordered along the spiral arms. Relativistic electrons gyrating around magnetic field lines with Lorentz factor $\gamma = E/(m_e c^2)$ emit beamed synchrotron radiation at critical frequency:

A power-law electron distribution $N(E) \propto E^{-p}$ produces a synchrotron radio spectrum with flux $S_\nu \propto \nu^{-\alpha}$, where the spectral index is $\alpha = (p - 1)/2 \approx 0.7-0.8$, illuminating the Milky Way at radio frequencies.