1.3 能标、单位与天文量纲
1.3 能标、单位与天文量纲
本节摘要：粒子天体物理同时使用粒子物理单位（eV、nb）和天文单位（pc、erg），两套体系在论文里混排。本节给出最小够用的换算工具箱，并用一段可运行的代码把常见换算固化下来，读完能无障碍读谱线和流量表。
一、先看一组真实数字
- 费米实验室加速器：质子能量约 1 TeV。
- IceCube 记录的最高能中微子（2023 年公告）：约 2.6 PeV（后续研究估计可能更高，常被引作约 100 PeV 量级的 KM3NeT 事件另计）。
- Oh-My-God 宇宙线（1991，Fly's Eye）：约 3×10^20 eV = 300 EeV。
从 TeV 到 EeV 跨六个数量级，天体物理的能标谱比粒子实验宽得多，单位换算是每天的日常。
二、两套单位体系的对照
| 量 | 粒子物理习惯 | 天文习惯 | 换算 |
|---|---|---|---|
| 能量 | GeV / TeV / PeV / EeV | erg | 1 eV ≈ 1.602×10^-12 erg |
| 长度 | 飞米 fm（1 fm=10^-15 m） | 秒差距 pc | 1 pc ≈ 3.086×10^16 m |
| 截面 | cm² 或靶恩 bn | — | 1 bn = 10^-24 cm² |
| 流量 | ph cm^-2 s^-1 | erg cm^-2 s^-1 | 需乘以单光子能量 |
自然单位制里取 c = ħ = 1，于是"长度 = 时间 = 1/能量"：1 GeV^-1 ≈ 0.197 fm，1 GeV^-1 ≈ 6.58×10^-25 s。记住这两个数，推导传播方程时会快很多。
三、把换算写成代码
一个极简换算工具，展示自然单位制的实际用法：
HBARC = 0.1973269804 # GeV·fm，即 1/GeV 换算成飞米 def gev_to_fm(inv_gev): """把自然单位下的长度（GeV^-1）换成飞米""" return inv_gev * HBARC def eV_ladder(ev): """把 eV 数值折叠成带词头的可读字符串""" units = [(1e18, "EeV"), (1e15, "PeV"), (1e12, "TeV"), (1e9, "GeV"), (1e6, "MeV"), (1e3, "keV"), (1.0, "eV")] for factor, name in units: if ev >= factor: return f"{ev/factor:.3g} {name}" return f"{ev:.3g} eV" print(eV_ladder(3e20)) # 300 EeV，Oh-My-God 粒子 print(gev_to_fm(1/0.14)) # 约 1.4 fm，π介子康普顿波长量级
💡 关键直觉：能谱图横轴常用 log10(E/eV)，"膝盖"（knee）在 10^15.5 附近、"踝"（ankle）在 10^18.5 附近——记住这两个坐标，读宇宙线谱就不迷路。
本节要点回顾
- 两套单位必须双语：eV 族与 erg/pc 族在论文里共存，1 eV ≈ 1.6×10^-12 erg 是最常用桥。
- 自然单位制省墨水：c = ħ = 1 把长度、时间、能量统一，0.197 GeV·fm 是核心常数。
- 能标坐标系：knee≈10^15.5 eV、ankle≈10^18.5 eV，是宇宙线谱的两条地标。
- 代码固化换算比死记公式可靠，上面十行脚本覆盖了日常 90% 的需求。
动手估算：从观测流量反推源光度
能标与单位不是背几个词，而是能实际换算。练习：IceCube 报告一个源的中微子能流约 1.2×10^-12 TeV cm^-2 s^-1，源距离 1.5 Gpc，求中微子光度。步骤：能流乘 4πd² 再换算成 erg/s。这个量级（约 10^45 erg/s）立刻给出判断——只有喷流型活动星系核供得起，普通恒星过程直接排除。
# 能流 -> 光度 的完整换算链 F = 1.2e-12 # TeV cm^-2 s^-1 d = 1.5 # Gpc d_cm = d * 3.086e27 # Gpc -> cm tev_to_erg = 1.602e-3 # 1 TeV 的 erg 数 L = F * tev_to_erg * 4 * 3.14159 * d_cm**2 print(f"中微子光度 L ~ {L:.2e} erg/s") # 输出: 中微子光度 L ~ 1.72e+45 erg/s
常用单位速查（天文 ↔ 粒子物理）
| 量 | 天文习惯 | 粒子物理习惯 | 换算要点 |
|---|---|---|---|
| 能量 | erg | eV/GeV/TeV | 1 eV = 1.602×10^-12 erg |
| 长度 | pc、kpc、Mpc | 飞米、厘米 | 1 pc = 3.086×10^18 cm |
| 磁场 | 高斯（G） | 特斯拉 | 1 G = 10^-4 T |
| 能量密度 | GeV/cm^3 | J/m^3 | 宇宙临界密度约 5×10^-6 GeV/cm^3 |
| 流量 | Jy | ph cm^-2 s^-1 | 能段不同不可直接乘 |
# 自然单位制下"1 电子伏特有多重/多长" hbar_c = 197.327 # MeV·fm E = 1.0 # 以 1 MeV 为例 length_fm = hbar_c / E time_s = 6.582e-22 / E # hbar = 6.582e-22 MeV·s print(f"1 MeV 对应约 {length_fm:.0f} fm 的长度、{time_s:.2e} 秒的时间") # 这套换算在讨论早期宇宙粒子退耦时每天都会用到
延伸讨论：单位制背后的思维差异
天文习惯与粒子物理习惯的差异，表面是换算系数，深层是两种"提问方式"。天文的量以观测可达性为轴组织：星等、流量、视差，全都围绕"从地球看出去"的视角；粒子物理以过程截面为轴组织：能量、截面、亮度，围绕"两粒子相遇会发生什么"。粒子天体物理学家的日常就是在这两种视角间翻译——把望远镜测得的每平方厘米每秒光子数，翻译成源区每秒每立方厘米的反应数，再翻译回能被下一台望远镜检验的光子数。翻译的每一环都可能引入一个量级的系统误差，所以本领域对单位换算的较真程度在外行看来近乎苛刻：一个错误的 kpc/cm 换算足以把"超新星遗迹是宇宙线源"的结论整个翻转。养成"每个数字都重算一遍"的反射，是本章比任何换算表都重要的产出。
一个常犯的错误值得点名：把"流量"与"光度"通过距离平方直接互化时忘记红移——(1+z)因子在光度距离定义中出现两次，对z=1的源会造成4倍的系统偏差。凡涉及河外源能流换算，先写下(1+z)再动笔，是避免论文级错误的一条铁律。