星系光度函数
一个常见误解:星系不是"按一个平均亮度均匀分布"
人们直觉上常以为,宇宙里的星系大体差不多亮,多数挤在某个"平均值"附近。事实正相反。星系的亮度跨越极大的范围:从一个矮椭球星系(仅约 )到中央亮椭圆星系(超过 ),相差六个数量级。而且,暗弱星系的数目远远多于明亮星系——只是它们太暗,早期巡天容易漏掉。
星系光度函数(galaxy luminosity function,LF)就是回答这个问题的工具:单位体积内,每个光度区间里有多少个星系?它不是一条无聊的统计曲线,而是星系形成与演化理论必须复现的定量基准。任何一个宇宙学模拟,如果造不出和观测一致的光度函数,就说明它对气体冷却、恒星形成和反馈的处理有问题。
Schechter 函数:一个公式装下整个星系族群
1976 年,普林斯顿天文学家 Paul Schechter 提出了一个至今仍是标准的解析形式(Schechter 1976,ApJ 203,297):
它只有三个参数,却抓住了星系族群的全部主干特征:
- :归一化密度,决定单位体积内星系的总数标度,典型量级约 。
- :特征光度,曲线从幂律转为指数截断的"拐点"。比 更亮的星系数目按指数急剧减少。银河系的光度大致就在 量级附近,所以 星系常被当作"典型大星系"的参照。
- :暗端斜率(faint-end slope)。在 处函数趋近幂律 。观测上 多在 $-1.0$ 到 $-1.3$ 之间——这个负值意味着越暗的星系越多,矮星系在数目上完全主导。
换成星等形式, 对应特征绝对星等 。在 波段,2dF 星系红移巡天给出 、(Norberg et al. 2002,MNRAS 336,907)。
注意符号约定:Schechter 1976 原文把幂律指数写作正的 (即 标度的暗端),现代文献普遍改用 形式,对应 。读不同论文时要看清楚作者用的是哪种约定。
从光度函数到质量函数
光度反映的是"现在亮不亮",会随恒星年龄、尘埃和近期恒星形成而变。要刻画星系本身的多少,天文学家更关心星系恒星质量函数(galaxy stellar mass function,GSMF),即单位体积内各恒星质量区间的星系数目。GSMF 同样能用 Schechter 形式拟合,特征质量 –(依红移与巡天而变)。基于 SDSS 的测量给出本地宇宙的 GSMF 在 附近转折,低质量端斜率约 (Baldry et al. 2012,MNRAS 421,621)。GSMF 比 LF 更接近理论量,因为模拟直接预言的是恒星质量,而非随波段变化的瞬时光度。
为何是星系演化的定量基石
这里是光度函数真正有深度的地方:它的形状本身就编码了星系形成中两大反馈过程。宇宙学的出发点是暗物质晕质量函数——在冷暗物质(CDM)框架下,由 Press-Schechter 理论(Press & Schechter 1974,ApJ 187,425)及其后续模拟精确预言。晕质量函数同样是低质量端陡、高质量端指数截断。问题在于:如果星系亮度简单正比于晕质量,那么观测到的星系光度函数应当远比实际更陡、更亮端更多。两者对不上。
差异恰恰揭示了反馈:
- 暗端没有想象中那么陡。低质量晕里,年轻恒星的紫外辐射、超新星爆发驱动的星系风把气体吹散,压制了矮星系的恒星形成效率。这就是为什么观测 ,而非暗物质晕那样接近 $-1.9$ 的陡度。
- 亮端呈指数截断。在大质量晕里,中央超大质量黑洞的活动星系核(AGN)反馈加热气体、阻止其冷却成恒星,给最亮星系的形成踩了刹车。
因此,"光度函数为什么是这个形状"几乎等价于问"反馈如何调节恒星形成"。半解析模型和流体力学模拟(如 EAGLE、IllustrisTNG)把复现观测 LF/GSMF 当作首要检验标准(Schaye et al. 2015,MNRAS 446,521)。
争议与前沿
- 暗端斜率的真值至今仍有分歧。极暗的矮星系表面亮度低、易漏检,不同巡天给出的 从约 $-1.0$ 到 $-1.5$ 不等。这直接关系到"卫星星系数目偏少问题"(missing satellites)等小尺度 CDM 难题。
- 高红移演化:JWST 测得 $z > 10$ 的紫外光度函数暗端意外地多、亮端星系也比预期亮,挑战了反馈强度与恒星形成效率的标准假设(Harikane et al. 2023,ApJS 265,5)。这是当前最活跃的争论之一。
- 环境依赖:星系团内的光度函数和场星系不同——潮汐剥离、并合与"窒息"(气体被剥夺)会改变矮星系的丰度。单一 Schechter 形式只是一阶近似,常需双 Schechter 函数才能同时拟合矮端与巨端。
跨域连接
- 暗物质晕:晕质量函数的低质量端接近 $-1.9$,星系光度函数只有约 $-1.2$;两条曲线的落差就是反馈强度沿质量的剖面。推论是暗端斜率并非拟合参数,而是矮星系反馈效率的直接读数——斜率一变,反馈模型必须跟着改;这也意味着晕质量函数不能直接拿来当星系计数用。
- 统计学:三个参数装下六个数量级,靠的是把标度自由的幂律与标度截断的指数分开写。推论是一旦要同时拟合矮端与巨端就得换成双 Schechter 形式,这说明单一形式隐含了"只有一种反馈"的假设——形式不够用本身就是物理信息,所以函数形式的选择要当成待检验的假设,而不是拟合的技术细节。
- 不平等经济学:斜率约 $-1.2$ 时,数目积分被暗端主导而光度积分被特征光度附近主导——"数量多"与"总量大"分属两端。推论是问"典型星系有多亮"必须先说明按数目还是按光度加权,两个答案相差几个数量级。
- 抽样调查:任何巡天都是一次抽样,暗弱星系的入选概率天然更低;不校正入选概率,估出的分布必然偏向"亮的太多"。推论是暗端斜率在不同巡天间从约 $-1.0$ 到 $-1.5$ 的分歧,首先该怀疑选择函数而不是物理,因此这项误差应当以巡天为单位分别报告。
- 相变与临界现象:幂律加指数截断这一形式在临界现象里也出现,截断尺度对应关联长度。这一步必须写明只是形式类比:它的用处是提示去追问"是什么物理量设定了特征光度",而不是暗示两者机制相同。
参考文献
- Schechter, P. (1976). An analytic expression for the luminosity function for galaxies. Astrophysical Journal, 203, 297–306.
- Press, W. H. & Schechter, P. (1974). Formation of galaxies and clusters of galaxies by self-similar gravitational condensation. Astrophysical Journal, 187, 425–438.
- Norberg, P. et al. (2002). The 2dF Galaxy Redshift Survey: the bJ-band galaxy luminosity function and survey selection function. MNRAS, 336, 907–931.
- Baldry, I. K. et al. (2012). Galaxy And Mass Assembly (GAMA): the galaxy stellar mass function at z < 0.06. MNRAS, 421, 621–634.
- Schaye, J. et al. (2015). The EAGLE project: simulating the evolution and assembly of galaxies and their environments. MNRAS, 446, 521–554.
延伸阅读
- Harikane, Y. et al. (2023). A comprehensive study of galaxies at z ~ 9–16 found in early JWST data. Astrophysical Journal Supplement Series, 265, 5.
- Mo, H., van den Bosch, F. & White, S. (2010). Galaxy Formation and Evolution. Cambridge University Press.