Extragalactic Astronomy: Galaxies, AGN & Clusters
This unit broadens our cosmological perspective to extragalactic scales, investigating the morphology, physics, and evolution of external galaxies, active galactic nuclei (AGN), and rich galaxy clusters. We examine galaxy classification via the Hubble tuning fork, surface brightness laws (de Vaucouleurs and exponential disks), and empirical scaling relations (Tully-Fisher and Faber-Jackson). We construct the unified model of AGN, calculating the Eddington luminosity limit, accretion efficiency of supermassive black holes, and the kinematic illusion of superluminal jet motion. Finally, we explore the physics of galaxy clusters, the hot X-ray emitting intra-cluster medium (ICM), and apply the Virial Theorem to reproduce Fritz Zwicky's 1933 discovery of dark matter.
§6.1 Galaxy Classification & Empirical Scaling Relations
Edwin Hubble (1926, 1936) classified galaxies based on their optical morphologies, constructing the famous Hubble Tuning Fork Diagram.
1. Elliptical Galaxies (E0 - E7)
Elliptical galaxies appear as smooth, featureless spheroids with little cool interstellar gas and negligible ongoing star formation, populated primarily by old Population II stars. The ellipticity class $n$ is defined by apparent major axis $a$ and minor axis $b$:
ranging from E0 (circular, $b/a = 1$) to E7 (highly flattened, $b/a = 0.3$). Ellipticals are supported dynamically not by bulk rotation, but by random anisotropic velocity dispersion $\sigma$. Their radial surface brightness profile obeys the de Vaucouleurs $R^{1/4}$ Law:
where $R_e$ is the effective half-light radius enclosing $50\%$ of total galaxy luminosity.
2. Spiral Galaxies (S and SB)
Spirals feature a central spheroidal bulge and a thin rotating disk with spiral arms. They divide into normal spirals (S) and barred spirals (SB), sub-classified as a, b, c:
- Sa / SBa: Large dominant central bulge, tightly wound smooth spiral arms, low gas fraction.
- Sb / SBb: Intermediate bulge and arm openness (e.g., Milky Way is SBbc, Andromeda is SA(s)b).
- Sc / SBc: Small central bulge, loosely wound knotty arms rich in luminous H II star-forming regions.
The disk surface brightness profile obeys an exponential law:
Empirical Scaling Relations
- The Tully-Fisher Relation (Spirals, 1977): Connects total stellar luminosity $L$ to maximum rotation velocity $V_{\text{max}}$:
$$L \propto V_{\text{max}}^4 \implies M_B = -10 \log_{10}(V_{\text{max}}) + \text{const}$$Derived physically from $M \propto V^2 R/G$ and constant mean surface brightness $I_0 \propto L/R^2 \implies R \propto \sqrt{L}$, giving $M \propto V^2 \sqrt{L} \propto L \implies L \propto V^4$. This relation serves as a potent extragalactic distance indicator out to $\sim 100\text{ Mpc}$.
- The Faber-Jackson Relation (Ellipticals, 1976): Relates elliptical galaxy luminosity to central 1D velocity dispersion $\sigma$:
$$L \propto \sigma^4$$Generalized in 3D parameter space as the Fundamental Plane: $\log R_e = \alpha \log \sigma + \beta \langle \mu_e \rangle + \gamma$.
§6.2 Galaxy Formation, Hierarchical Mergers & Feedback Physics
Modern extragalactic astrophysics understands galaxy assembly within the framework of $\Lambda\text{CDM}$ hierarchical bottom-up structure formation.
Hierarchical Merging vs Monolithic Collapse
Early models (Eggen, Lynden-Bell, & Sandage 1962) proposed that galaxies formed via rapid monolithic gravitational collapse of a single giant protogalactic gas cloud. Today, cosmological simulations (e.g., Illustris, EAGLE) establish that structure forms hierarchically:
- Cold dark matter clumps collapse first on sub-galactic scales ($M \sim 10^6 - 10^8 M_\odot$).
- Baryonic gas cools radiatively via atomic hydrogen line transitions ($T > 10^4\text{ K}$) and sinks to the centers of dark matter potential wells, spinning up to form rotationally supported gas disks.
- Repeated minor mergers build stellar halos and thick disks, while major mergers (mass ratio $> 1:4$) violently disrupt disks, scrambling stellar orbits into pressure-supported elliptical galaxies (Toomre merger hypothesis).
Stellar & AGN Feedback Mechanisms
Without energetic feedback, numerical simulations predict that all gas would rapidly cool and convert into stars, producing an overabundance of hyper-luminous galaxies, contradicting the observed Schechter luminosity function. Two feedback mechanisms regulate galaxy growth:
- Supernova & Stellar Feedback: In low-mass dwarf galaxies ($M < 10^{10} M_\odot$), collective supernova explosions and stellar winds drive supersonic galactic superwinds ($v \sim 500\text{ km/s}$), blowing gas completely out of shallow potential wells and suppressing dwarf galaxy formation.
- Active Galactic Nucleus (AGN) Feedback: In massive galaxies ($M > 10^{11} M_\odot$), the central supermassive black hole releases colossal energy via relativistic radio jets ('maintenance/radio mode') and radiation pressure ('quasar mode'), heating intra-cluster gas and preventing gas from cooling and collapsing, thereby 'quenching' star formation and capping maximum galaxy masses.
§6.3 Active Galactic Nuclei (AGN), Quasars & The Eddington Limit
Active Galactic Nuclei (AGN) are the most luminous steady sources in the Universe, emitting up to $10^{41}\text{ Watts}$ ($10^{14} L_\odot$) from a compact region no larger than the Solar System ($< 10^{-4}\text{ pc}$).
The Eddington Luminosity Limit
Consider an ionized gas consisting of free electrons and protons surrounding an accreting black hole. Outward radiation pressure acts primarily on electrons via Thomson scattering (cross section $\sigma_T = 6.652 \times 10^{-29}\text{ m}^2$), while inward gravitational attraction acts primarily on protons ($m_p \gg m_e$). Electrostatic coupling prevents charge separation. The outward radiative radiation force on an electron-proton pair at radius $r$ is:
The inward gravitational force is:
Setting $F_{\text{rad}} = F_{\text{grav}}$ defines the maximum steady luminosity an object can radiate without blowing away its accreting material—the Eddington Limit ($L_{\text{Edd}}$):
Eddington Accretion Rate & Salpeter Timescale
The radiant luminosity is produced by gravitational accretion at mass rate $\dot{M}$ with radiative efficiency $\eta$: $L = \eta \dot{M} c^2$ (where $\eta \approx 0.10$ for standard thin accretion disks). The Eddington accretion rate is:
Because the growth rate $\dot{M} \propto M$ is exponential, a seed black hole growing at the Eddington limit increases its mass as $M(t) = M_0 e^{t / \tau_S}$, where the Salpeter growth timescale is:
To grow a billion solar mass ($10^9 M_\odot$) quasar black hole at $z \sim 7$ (only $800\text{ Myr}$ after the Big Bang) requires continuous, near-uninterrupted Eddington accretion from early stellar seed remnants.
§6.4 The Unified Model of AGN & Relativistic Superluminal Jets
The observational diversity of AGN—Seyfert 1 and 2 galaxies, radio galaxies (FR I and FR II), quasars, and blazars—is explained by the Unified Model of AGN (Antonucci 1993, Urry & Padovani 1995).
Anatomy of the Unified AGN Engine
- Central Supermassive Black Hole: Mass $M_{\text{BH}} \sim 10^6 - 10^{10} M_\odot$.
- Shakura-Sunyaev Accretion Disk: Geometrically thin, optically thick plasma disk ($r \sim 10^{-4} - 10^{-2}\text{ pc}$) releasing intense thermal UV/optical radiation ('Big Blue Bump').
- Broad Line Region (BLR): High-density gas clouds ($n_e > 10^9\text{ cm}^{-3}$) orbiting close to the black hole ($r \sim 0.01 - 0.1\text{ pc}$) with high Keplerian velocities ($v \sim 1{,}000 - 10{,}000\text{ km/s}$), producing Doppler-broadened permitted emission lines (H$\alpha$, H$\beta$, C IV).
- Dusty Molecular Torus: Thick obscuring donut of gas and dust ($r \sim 1 - 10\text{ pc}$) aligned with the accretion disk plane.
- Narrow Line Region (NLR): Low-density gas clouds ($n_e \sim 10^3 - 10^6\text{ cm}^{-3}$) at large distances ($r \sim 100 - 1000\text{ pc}$) producing narrow forbidden lines ([O III], [N II]) with velocities $v \sim 300-500\text{ km/s}$.
- Relativistic Bipolar Jets: Magnetically collimated plasma beams launched perpendicular to the disk along the black hole spin axis.
Orientation-Based Unification
- Face-On / Unobscured View ($\theta < \theta_{\text{torus}}$): The observer looks directly into the inner core, viewing both the accretion disk and BLR $\implies$ Type 1 AGN (Seyfert 1, Quasars) showing broad + narrow lines.
- Edge-On / Obscured View ($\theta > \theta_{\text{torus}}$): The dusty torus completely blocks the line of sight to the central engine and BLR. Only the extended NLR is visible $\implies$ Type 2 AGN (Seyfert 2, Narrow-line radio galaxies) showing only narrow lines. Spectropolarimetry (Antonucci & Miller 1985) revealed hidden broad lines in polarized scattered light of NGC 1068, proving the model!
- Looking Directly Down the Jet ($\theta \approx 0^\circ$): The emission is dominated by violently variable, Doppler-boosted synchrotron radiation $\implies$ Blazar / BL Lac object.
The Kinematics of Superluminal Motion
Very Long Baseline Interferometry (VLBI) radio observations of quasar jets (e.g., 3C 273, M87) track plasma blobs apparently moving across the sky at velocities $v_{\text{app}} = 5c - 10c$, seemingly violating special relativity. This is a purely geometric time-compression effect.
Consider a blob launched from the nucleus at $t=0$ moving at true relativistic speed $v = \beta c$ at an angle $\theta$ relative to the observer's line of sight. At time $t$, the blob has traveled distance $v t$.
The physical transverse displacement across the sky plane is $\Delta x = v t \sin\theta$. Meanwhile, the blob has moved closer to the observer by $\Delta z = v t \cos\theta$. Photons emitted at time $t$ have a shorter distance to travel, arriving at the observer at time:
The apparent transverse speed measured by the observer is:
Differentiating with respect to $\theta$ to find the maximum apparent speed: $\frac{d}{d\theta}\left(\frac{\sin\theta}{1 - \beta\cos\theta}\right) = 0 \implies \cos\theta_{\text{max}} = \beta$. Substituting into the equation gives:
For an ultra-relativistic jet with $\beta = 0.995$ ($\gamma \approx 10$), $v_{\text{app},\text{max}} \approx 10 c$, creating the illusion of superluminal expansion without violating relativity.
§6.5 Galaxy Clusters, Intra-Cluster Medium & Large-Scale Cosmic Web
Galaxy clusters are the largest gravitationally virialized structures in the universe, containing hundreds to thousands of galaxies bound within a common dark matter potential well spanning $R \sim 1-3\text{ Mpc}$ with total masses $M \sim 10^{14} - 10^{15} M_\odot$.
The Intra-Cluster Medium (ICM)
Contrary to optical appearances, galaxies contain only $\sim 2-5\%$ of the total baryonic mass of a cluster. The vast majority of cluster baryons ($\sim 12-15\%$) reside in the Intra-Cluster Medium (ICM)—a tenuous ($n_e \sim 10^{-4} - 10^{-2}\text{ cm}^{-3}$), shock-heated plasma at temperatures $T \sim 10^7 - 10^8\text{ K}$ ($k T \sim 1 - 10\text{ keV}$).
At these extreme temperatures, the ICM is completely ionized and radiates intensely in X-rays via thermal bremsstrahlung (free-free emission) with emissivity:
Observatories such as Chandra and XMM-Newton map cluster X-ray emission to determine plasma temperature $T(r)$ and density $n_e(r)$, enabling hydrostatic mass estimation.
The Sunyaev-Zel'dovich (SZ) Effect
When low-energy photons of the Cosmic Microwave Background (CMB) pass through the hot ICM of a galaxy cluster, they undergo inverse Compton scattering off relativistic thermal electrons, gaining a small energy boost ($h\nu' > h\nu$). This distorts the Planck CMB blackbody spectrum, producing a distinctive temperature decrement at $\nu < 217\text{ GHz}$ and an increment at $\nu > 217\text{ GHz}$. The thermal SZ temperature shift is:
Because the SZ effect is a spectral distortion independent of distance (redshift), instruments like the Planck satellite and South Pole Telescope detect distant galaxy clusters out to $z > 1.5$.
§6.6 Virial Mass Estimation & Dark Matter Proof in Galaxy Clusters
The existence of dark matter was first discovered not in individual galaxies, but in galaxy clusters by the Swiss astrophysicist Fritz Zwicky in 1933.
Zwicky's 1933 Discovery of Dark Matter
Zwicky measured the radial velocities of galaxies in the Coma Cluster using the Mount Wilson 100-inch telescope. He observed an unexpectedly large line-of-sight velocity dispersion: $\sigma_r \approx 1000\text{ km s}^{-1}$.
Applying the Virial Theorem ($2K + \Omega = 0$) to a cluster of $N$ galaxies with total mass $M$ and virial radius $R_{\text{vir}}$:
Zwicky computed the total dynamical mass $M_{\text{vir}}$ required to gravitationally bind the cluster against its galaxy velocities. Comparing this to the luminous mass inferred from galaxy counts and stellar mass-to-light ratios ($M_{\text{lum}} \approx N_{\text{gal}} \times 10^{11} M_\odot$), Zwicky discovered:
Zwicky concluded that the Coma Cluster must be dominated by invisible mass, which he termed dunkle Materie (dark matter).
Gravitational Lensing & The Bullet Cluster (1E 0657-56)
Einstein's General Relativity predicts that mass bends light rays by deflection angle $\hat{\alpha} = \frac{4GM}{c^2 b}$ ($b$ is impact parameter). Rich galaxy clusters act as colossal cosmic gravitational lenses, warping background galaxies into magnified arcs and multiple images.
The Bullet Cluster (Clowe et al. 2006) provides the definitive proof of particle dark matter over modified Newtonian dynamics (MOND):
- Two galaxy clusters recently collided at high speed ($v \sim 4500\text{ km s}^{-1}$).
- The collisional intra-cluster gas (observed in X-rays by Chandra, representing most of the baryonic mass) experienced ram-pressure hydrodynamic drag and stalled in the center.
- The galaxies (collisionless) passed through unhindered.
- Weak gravitational lensing maps (measuring total gravitational potential) revealed that the dominant gravitational mass peaks coincide directly with the collisionless galaxies, distinctly separated from the baryonic gas!
This physical spatial separation between the baryonic mass and the gravitational potential proves that the majority of matter in the Universe is non-baryonic and collisionless.
Honors Examination Worked Problems & Solutions
Rigorous step-by-step mathematical proofs and solutions to university degree examination problems.
The prominent quasar 3C 273 has a measured bolometric luminosity $L_{\text{bol}} = 4.0 \times 10^{39}\text{ Watts}$ ($4.0 \times 10^{46}\text{ erg s}^{-1}$).\n\n(a) Calculate the minimum mass of the central supermassive black hole $M_{\text{BH}}$ assuming the quasar radiates at or below the Eddington limit.\n(b) Assuming standard accretion disk radiative efficiency $\eta = 0.10$, compute the mass accretion rate $\dot{M}$ in $\text{kg s}^{-1}$ and in solar masses per year ($M_\odot\text{ yr}^{-1}$).\n(c) Calculate the Schwarzschild radius $R_s$ in AU and in light-hours.\n(d) Calculate the minimum Salpeter e-folding growth timescale $\tau_S$.
(a) Minimum Black Hole Mass from Eddington Limit
The Eddington limit is $L_{\text{Edd}} = \frac{4\pi G M_{\text{BH}} m_p c}{\sigma_T}$. Setting $L_{\text{bol}} \le L_{\text{Edd}}$:
In solar masses ($M_\odot = 1.989 \times 10^{30}\text{ kg}$):
The central engine contains a black hole of at least $318\text{ million solar masses}$.
(b) Mass Accretion Rate
From $L = \eta \dot{M} c^2$ with $\eta = 0.10$:
Converting to solar masses per year ($1\text{ yr} = 3.1558 \times 10^7\text{ s}$):
The black hole devours roughly 7 full stars worth of mass every year.
(c) Schwarzschild Radius
The event horizon radius is:
In Astronomical Units ($1\text{ AU} = 1.496 \times 10^{11}\text{ m}$):
In light-hours ($c \times 3600\text{ s} = 1.079 \times 10^{12}\text{ m}$):
(d) Salpeter Growth Timescale
The Salpeter timescale is:
Complete rigorous derivation and proof detailed above.
Complete rigorous derivation and proof detailed above.
VLBI radio observations of a knot in the relativistic jet of quasar 3C 279 show an apparent proper motion $\mu = 0.52\text{ mas yr}^{-1}$. The quasar is at redshift $z = 0.536$, corresponding to an angular diameter distance $d_A = 1320\text{ Mpc}$.\n\n(a) Compute the apparent transverse velocity $v_{\text{app}}$ in units of $c$ (superluminal parameter $\beta_{\text{app}} = v_{\text{app}}/c$).\n(b) Using the superluminal equation $\beta_{\text{app}} = \frac{\beta \sin\theta}{1 - \beta \cos\theta}$, derive the minimum possible true physical velocity $\beta_{\text{min}} = v_{\text{min}}/c$ and minimum Lorentz factor $\gamma_{\text{min}}$.\n(c) For $\beta = 0.995$ ($\gamma = 10.0$), find the maximum viewing angle $\theta_{\text{max}}$ that still produces the observed superluminal speed.
(a) Apparent Transverse Velocity Calculation
The angular displacement rate is $\mu = 0.52\text{ mas yr}^{-1} = 0.52 \times 10^{-3} \times \left(\frac{\pi}{180 \times 3600}\right) \approx 2.521 \times 10^{-15}\text{ rad yr}^{-1}$.
The linear transverse apparent speed is $v_{\text{app}} = d_A \mu$:
Dividing by $c = 2.998 \times 10^8\text{ m/s}$:
The radio knot appears to move across the sky at nearly $11$ times the speed of light!
(b) Minimum True Jet Velocity & Lorentz Factor
The maximum apparent speed achievable for a given true velocity $\beta$ occurs at $\cos\theta = \beta$, giving $\beta_{\text{app},\text{max}} = \beta \gamma = \frac{\beta}{\sqrt{1 - \beta^2}}$. Setting this equal to the observed $\beta_{\text{app}} = 10.85$:
The minimum Lorentz factor is:
The jet must be moving at least $99.58\%$ the speed of light ($\gamma \ge 10.9$).
(c) Permissible Viewing Angle Range
If $\beta = 0.995$, the equation $\beta_{\text{app}} = \frac{\beta \sin\theta}{1 - \beta \cos\theta} = 10.85$ rearranges to:
Squaring both sides using $\sin^2\theta = 1 - \cos^2\theta$ and solving the quadratic equation yields the permissible range:
The jet must be pointed within $\approx 7^\circ$ of our direct line of sight.
Complete rigorous derivation and proof detailed above.
Complete rigorous derivation and proof detailed above.
Spectroscopic surveys of the Coma Galaxy Cluster (Abell 1656) measure a line-of-sight velocity dispersion $\sigma_r = 1008\text{ km s}^{-1}$ for its member galaxies. The virial radius enclosing the bound cluster is $R_{\text{vir}} = 2.0\text{ Mpc}$.\n\n(a) Use the Virial Theorem $M_{\text{vir}} = \frac{5 \sigma_r^2 R_{\text{vir}}}{G}$ to calculate the total cluster dynamical mass $M_{\text{vir}}$ in kg and in solar masses ($M_\odot$).\n(b) Optical counts detect $\approx 1000$ bright galaxies with mean stellar mass $M_ \approx 3.0 \times 10^{10} M_\odot$. Calculate the stellar mass $M_{\text{stars}}$ and the stellar mass fraction $f_ = M_{\text{stars}}/M_{\text{vir}}$.\n(c) X-ray observations by the Chandra satellite reveal an intra-cluster gas mass $M_{\text{gas}} = 1.40 \times 10^{14} M_\odot$. Calculate the total baryonic mass $M_{\text{baryon}} = M_{\text{stars}} + M_{\text{gas}}$ and the total dark matter mass $M_{\text{DM}}$.\n(d) Calculate the dark matter fraction $f_{\text{DM}} = M_{\text{DM}}/M_{\text{vir}}$ and compare with cosmic $\Omega_{\text{DM}}/\Omega_m$.
(a) Virial Dynamical Mass
Given $\sigma_r = 1008\text{ km s}^{-1} = 1.008 \times 10^6\text{ m s}^{-1}$ and $R_{\text{vir}} = 2.0\text{ Mpc} = 2.0 \times 3.0857 \times 10^{22}\text{ m} = 6.1714 \times 10^{22}\text{ m}$:
In solar masses ($M_\odot = 1.989 \times 10^{30}\text{ kg}$):
The Coma cluster contains over $2.3\text{ quadrillion solar masses}$.
(b) Stellar Mass & Stellar Fraction
The total stellar mass in member galaxies is:
Visible stars account for barely $1.3\%$ of the total mass of the cluster!
(c) Baryon Mass vs Dark Matter Mass
Total baryonic mass:
Notice that the hot X-ray gas contains $\frac{1.40 \times 10^{14}}{3.0 \times 10^{13}} \approx 4.7$ times more mass than all the stars combined!
The dark matter mass is:
(d) Dark Matter Fraction & Cosmic Concordance
The dark matter fraction in the cluster is:
Over $92\%$ of the cluster is composed of invisible dark matter. The ratio of baryonic to total matter in the cluster ($f_b \approx 7.2-15\%$) closely mirrors the universal cosmological baryon fraction $\Omega_b / \Omega_m = \frac{0.049}{0.315} \approx 15.6\%$, demonstrating that galaxy clusters represent fair, representative cosmological samples of the Universe.
Complete rigorous derivation and proof detailed above.
Complete rigorous derivation and proof detailed above.