Foundations of Observational & Theoretical Astronomy / 观测与理论天文学基础
第 1 章 观测与理论天文学基础
Chapter 1 — Foundations of Observational & Theoretical Astronomy
> 核心问题 / Core Question — 我们如何把遥远天体从像素点转化为可测量的物理量，又如何把跨越 13 个数量级的尺度压缩到同一套数学语言之中？这一章将回答这一问题：从球面坐标系的定义出发，过渡到测光定标与星际消光模型，搭建宇宙距离阶梯的误差预算，再用牛顿两体问题与水星近日点进动打开通往广义相对论的大门，最后以 Gaia EDR3 为案例，给出一份可直接在 Python 中复现的实践流程。
How do astronomers transform raw photon counts into physically meaningful quantities, and how do we set up the mathematical frameworks that allow extrapolation across thirteen orders of magnitude in scale? This chapter answers that question through coordinate systems, photometry, the distance ladder, Newtonian two-body dynamics, and a hands-on Gaia EDR3 analysis.
---
1.1 球面天文学与坐标系统 / Spherical Astronomy and Coordinate Systems
球面天文学的首要任务是把二维探测器平面上的像元坐标映射到三维天体方向上，并规定唯一的、可重复的数值表达。所有球面坐标都由三个要素构成：原点（极轴方向）、参考平面（赤道面或类似平面）与参数化方法（角度或复合角）。天文学中最常用的是赤道坐标系、银道坐标系与超银道坐标系，三者之间的转换矩阵可以在任何需要精度的台址（VLT、Keck、ALMA、FAST）上实时计算，且全球测地学界已基本统一使用 IAU 2006/2006A 规范。
The first task of spherical astronomy is to map detector pixel coordinates to physical directions on the celestial sphere, parameterized by a reference frame, a fundamental plane, and an angular convention.
1.1.1 赤道坐标系 / The Equatorial Frame
赤道坐标系以地球自转轴方向为极轴，以地球赤道面（投影到天球上）为参考平面。位置用赤经 $\alpha$（由春分点沿黄赤交线向东量度，$0\!\le\!\alpha\!<\!2\pi$）和赤纬 $\delta$（由赤道面向北量度，$-\pi/2\!\le\!\delta\!\le\!\pi/2$）表示。在 *ICRS*（International Celestial Reference System）下，赤道坐标系被推广到以太阳系质心为原点的版本，并通过一组河外射电源（QSO、AGN）的 VLBI 位置来固定坐标轴方向。*ICRS* 坐标与历元 J2000.0 的平赤道平春分点（mean equator and equinox of J2000.0）重合到约 0.02 arcsec 的精度。
从仪器坐标系到 *ICRS* 的标准链条为：
```text detector pixel → tangent plane (ξ, η) → apparent (α, δ, J{epoch}) → proper-motion propagated → mean (α, δ) → ICRS ```
实现层面，IAU 2006/2006A 推荐使用高度近似 $d\epsilon_{\text{prec}}$ 项的 *P03* 模型，每百年误差小于 0.2 mas；若以更高精度进入 ms 域，则需在岁差之外加入章动矩阵 $N(\psi,\epsilon,\Delta\psi,\Delta\epsilon)$ 与帧偏差（frame bias）旋转。
1.1.2 银道与超银道 / Galactic & Supergalactic Frames
对恒星物理与星系天文学而言，赤道坐标不方便讨论银盘结构。因而 IAU 在 1958 年定义银道坐标系，以银心方向 $(\alpha_{G},\delta_{G}) = (266.40^\circ, -29.01^\circ)$ 作为升银点，银北极方向 $(\alpha_{P},\delta_{P}) = (192.86^\circ, 27.13^\circ)$。其与 *ICRS* 的固定角动量为：
$$ \begin{aligned} \alpha_G &= 266.40510^\circ - 0.01390^\circ\, t, \\ \delta_G &= -28.93617^\circ - 0.00089^\circ\, t, \\ \alpha_P &= 192.85950^\circ - 0.00096^\circ\, t, \\ \delta_P &= \phantom{-}27.12825^\circ + 0.00076^\circ\, t, \end{aligned} $$
其中 $t$ 为自 J2000.0 起算的儒略世纪数。超银道坐标系（de Vaucouleurs 1975）则把基平面转为包含本地星系（Virgo Cluster 等）的"超银道面"，其在宇宙学中用于描述本地宇宙的大尺度结构。
The galactic frame is defined by the position of the Galactic Center, while the supergalactic frame aligns with the Local Supercluster plane.
1.1.3 岁差、章动与帧偏差 / Precession, Nutation, Frame Bias
地球自转轴在空间中运动可分解为三部分：
1. 岁差 (luni-solar precession)：周期约 25,772 年的圆锥运动，主要由太阳与月球施加在地球赤道隆起上的非球形引力造成，量级约 50.3 arcsec/yr。 2. 章动 (nutation)：叠加在岁差之上的周期性振荡，主项 18.6 年月球章动 $\Delta\psi \approx 17.2''\sin\Omega$。 3. 帧偏差：在 CIP（Celestial Intermediate Pole）与 IAU 2000A 框架下的固定偏移 $\eta_0, \xi_0, d\alpha_0$，量级约 17 mas，需通过 IAU 2006 起算模型修正。
精化到毫角秒级别，观测归算从旧的等赤角等赤纬到天球中介原点 (*CIO-based*) 与赤道中间坐标 (*CIRS*) 的演化，已经成为 VLBI 与 Gaia 时序归算的默认规范。
1.1.4 儒略日及其变体 / Julian Date and Its Variants
时间系统对球面天文学至关重要，因为赤道与岁差均依赖历元。儒略日（JD）定义为从公元前 4713 年 1 月 1 日 12 TT 起算的连续天数：
$$ \text{JD} = 2451545.0 + \Delta T / 86400\,\text{s}. $$
当 $\Delta T$ 极小时，可由以下公式计算 JD：
$$ \text{JD} = 367\,Y - \text{int}\!\left(\tfrac{7(Y+\text{int}((M+9)/12))}{4}\right) + \text{int}\!\left(\tfrac{275\,M}{9}\right) + D + 1721013.5 + \tfrac{h + m/60 + s/3600}{24}. $$
修正变体包括 MJD $=\text{JD} - 2400000.5$、BJD（Barycentric Julian Date，太阳系质心时）与 TDB，后者与 TT 之差通常在 1.7 ms 量级。BJD 对多普勒与历元归算尤其关键：地面 1 m 望远镜相邻两端的轨道光程差在 100 s 内可达 10 μs，等价于 ~25 m 的光程差异。
The choice of time system — JD, MJD, BJD, TDB — encodes both the calendar convention and the inertial reference adopted for that measurement.
---
1.2 测光系统与星际消光 / Photometric Systems and Extinction Corrections
测光是把光子流转换成数字星等（$m$）的关键一环。任何有效的测光系统都必须先规定：光谱响应 $S(\lambda)$、标准源星表与零点。UBV、SDSS、Gaia $G$、2MASS、HSC、LSST 等多套系统并存，使得"同一源在不同过滤片下的星等"成为巡天交叉定标的基本数据。
1.2.1 视星等、绝对星等与距离模数 / Apparent, Absolute Magnitudes, and Distance Modulus
视星等 $m$ 与光子流密度 $F_\nu$ 由 Pogson 对数尺度定义：
$$ m_1 - m_2 = -2.5 \log_{10}\!\left(\frac{F_{1,\nu}}{F_{2,\nu}}\right). $$
将 $m$ 推回到 *绝对星等* $M$（光源位于 $d_0=10\,$pc 时的视星等）需要距离与消光的联合修正：
$$ \boxed{\;m - M = 5\log_{10}\!\left(\frac{d}{10\,\text{pc}}\right) + A, $$
其中 $A$ 是给定波段下的总消光，$m - M$ 称为 *距离模数*。这是度量宇宙的关键公式之一，但因为 $A$ 通常沿视线方向随路径长度而累积，所以 *单一线方向的距离模数本身并不直接给出 $d$*，必须先做消光剥离。
The distance modulus equation is the cornerstone of extragalactic distance measurement, but its dependence on extinction means de-reddening must precede distance estimation.
1.2.2 星际消光的物理机制 / Physical Origin of Interstellar Extinction
星际消光主要由 $\lesssim 0.1\,\mu\text{m}$ 尺寸的 *碳质* 与 *硅酸盐* 颗粒组成（Mathis, Rumpl & Nordsieck 1977, MRN）。在波长 $\lambda \gtrsim 0.3\,\mu\text{m}$ 的 *紫外到近红外* 区，消光曲线可参数化为幂律
$$ A(\lambda)/A(V) = a(x) + b(x)/R_V, $$
其中 $x\equiv 1/\lambda$（$\lambda$ 单位 $\mu\text{m}^{-1}$），且
$$ R_V \equiv A(V)/E(B-V) $$
反映总消光与选择消光之比，典型扩散云 $R_V\approx 3.1$，致密云团 $R_V\approx 5.5$。
1.2.3 Cardelli 律与 O'Donnell 修订 / Cardelli's Law
Cardelli, Clayton & Mathis (1989, *ApJ*, 345, 245) 给出单参数家族：
$$ \begin{aligned} a(x) &= 0.574\,x^{1.61}, \\ b(x) &= -0.527\,x^{1.61}, \end{aligned} \quad 0.3\le x\le 1.1. $$
此区域覆盖经典 UBVRIJK 主要波段，对 $R_V=3.1$ 平均误差 $\sigma(A)/A(V)\approx 0.05$。O'Donnell (1994) 修正了 $3.3<x<8.0$（即 $1150\,$Å $\,\lesssim\lambda\lesssim 3030\,$Å）的 $a(x),b(x)$ 系数，引入近紫外跳变特性：
$$ a(x)= 1.752-0.316\,x-0.104/(x^2-4.067^2) \quad (3.3\le x\le 8.0). $$
在 $\lambda>8\,\mu\text{m}^{-1}$ 区，则需红外模型（Indebetouw et al. 2005 *Spitzer* / GLIMPSE；Hensley & Draine 2023 偏振率自洽模型）。当消光过强（如 $A_V>5$）或视线穿越弥漫星际带（Diffuse Interstellar Bands, DIBs）时，建议使用 *光谱能量分布拟合*（SED-fitting, e.g., Le Phare、FAST、CIGALE）代替单色定标。
1.2.4 热改正 / Bolometric Correction
视星等通常只在某窄波段测量，而 *bolometric 流量* 涵盖全部波长：
$$ M_{\text{bol}} = m_\lambda + BC_\lambda, $$
其中 $BC_\lambda$ 是 *bolometric correction*。对主序恒星，Bessell, Castelli & Plez (1998) 给出 $-0.07 \lesssim BC_V \lesssim -5.6$，误差与 $T_{\text{eff}}$、$\log g$ 与金属丰度 $[\text{Fe/H}]$ 三维相关。对星系与活动星系核，则需以 *SED 模型积分* 计算或借助 *K-改正*（一阶 Taylor 展开 $K(z) \sim 2.5(1+z)\!\int\!d\log\nu\,(\partial L_\nu/\partial z)$）。
---
1.3 宇宙距离阶梯与系统误差预算 / The Cosmic Distance Ladder and Systematic Error Budgets
The cosmic distance ladder is the cascade of rungs that bring us from a few hundred parsecs to the Hubble flow. 每"一级"都以更远处的天体为标定物，把前一级校准 $d$ 的方法延伸到更远的尺度。
1.3.1 第一级：三角视差 / Rung 1 — Trigonometric Parallax
对于一台视差基线 $B=1\,$AU 的望远镜，周年视差 $\varpi$（单位 arcsec）通过 $d=\frac{1}{\varpi}\,$pc 给出（参见图 1.1）。
```mermaid flowchart TD A[Earth at t1 - 6 months] --> B[Sin axis toward Star] C[Earth at t0] --> B D[Earth at t2 + 6 months] --> B B --> E[Observed shift ellipse] ```
Figure 1.1 — *Trigonometric parallax geometry.* The star at distance $d$ traces a small ellipse on the sky of *angular* semimajor axis $\varpi$ during Earth's orbit. The *Gaia* spacecraft provides the same baseline from L2 with comparable amplitude but far higher angular precision.
Hipparcos (1989–1993) 给出了 117,955 颗恒星 $\sigma_{\varpi}/\varpi \le 0.2$ 的视差结果，量级约 1 mas。Gaia DR3 (2022) 将精度推到 0.02 mas 量级（$G<15$），对最近的暗弱源甚至达到 $\sim 20\,\mu$as。但即便如此，仍需注意 *Lutz–Kelker* 偏差：横截面上视差分布的 $\text{noise}$ 会把均值抬高，需做贝叶斯反向重构（如 Bailer-Jones 等 2021 的 `Apsis` 距离推断）。
1.3.2 第二级：造父变星周光关系 / Rung 2 — Cepheid Period-Luminosity
$\delta$ Cephei 型造父变星的 *周光关系*（P-L）由 Leavitt (1908) 发现，21 世纪被 HST 在河外造父、SPOLLWITR (Riess et al. 2009；2019) 通过 *HST WFC3* 的红外观测大幅收紧：
$$ M_W = \mu_0 - 2.5\log_{10}PL + \beta_W, $$
其中 $\mu_0$ 是大小麦哲伦云或 NGC 4258 的距离模数锚定，$\beta_W$ 是金属丰度修正（$-0.2\lesssim\beta_W\lesssim -0.3\,\text{dex}\,\mu\text{mag}$）。Riess et al. (2022, *ApJ* 938, L26) 给出 $H_0 = 73.04 \pm 1.04\,\text{km}\,\text{s}^{-1}\text{Mpc}^{-1}$，但与 Planck 协作的 CMB 推断 $H_0 = 67.4 \pm 0.5$ 存在 $4.99\sigma$ "Hubble tension"。