跳转到内容
← 返回知识库
星系天文学约 22 分钟阅读

星系光度函数

Schechter函数光度函数星系质量函数暗端斜率星系演化反馈暗物质晕

一个常见误解：星系不是"按一个平均亮度均匀分布"

人们直觉上常以为，宇宙里的星系大体差不多亮，多数挤在某个"平均值"附近。事实正相反。星系的亮度跨越极大的范围：从一个矮椭球星系（仅约 106 L⊙10^6\ L_\odot）到中央亮椭圆星系（超过 1012 L⊙10^{12}\ L_\odot），相差六个数量级。而且，暗弱星系的数目远远多于明亮星系——只是它们太暗，早期巡天容易漏掉。

星系光度函数（galaxy luminosity function，LF）就是回答这个问题的工具：单位体积内，每个光度区间里有多少个星系？它不是一条无聊的统计曲线，而是星系形成与演化理论必须复现的定量基准。任何一个宇宙学模拟，如果造不出和观测一致的光度函数，就说明它对气体冷却、恒星形成和反馈的处理有问题。

漏检不是偶然的疏忽，而是选择效应。1976 年，M. J. Disney 在 Nature 上指出：照相底片对表面亮度有硬门槛，弥散而暗的星系即使总光度不低，也会从星表里消失。后来的电荷耦合器件巡天改善了这件事，但没有取消它。任何有流量极限的样本，都必须用选择函数把「看得见的」还原成「体积里真实有的」。

因此，光度函数测量的难点从来不是画出一条曲线。难点是回答：这份星表为什么收录这个星系、为什么漏掉那个。红移巡天卖的不是数据量，是已知的完备度。没有它，暗端斜率测到的只是筛子的形状。

Schechter 函数：一个公式装下整个星系族群

1976 年，普林斯顿天文学家 Paul Schechter 提出了一个至今仍是标准的解析形式（Schechter 1976，ApJ 203，297）：

ϕ(L) dL=ϕ∗(LL∗)αexp⁡ ⁣(−LL∗)dLL∗\phi(L)\,\mathrm{d}L = \phi^{*}\left(\frac{L}{L^{*}}\right)^{\alpha} \exp\!\left(-\frac{L}{L^{*}}\right)\frac{\mathrm{d}L}{L^{*}}

它只有三个参数，却抓住了星系族群的全部主干特征：

  • ϕ∗\phi^{*}：归一化密度，决定单位体积内星系的总数标度，典型量级约 10−2 h3 Mpc−310^{-2}\ h^{3}\ \mathrm{Mpc}^{-3}。
  • L∗L^{*}：特征光度，曲线从幂律转为指数截断的"拐点"。比 L∗L^{*} 更亮的星系数目按指数急剧减少。银河系的光度大致就在 L∗L^{*} 量级附近，所以 L∗L^{*} 星系常被当作"典型大星系"的参照。
  • α\alpha：暗端斜率（faint-end slope）。在 L≪L∗L \ll L^{*} 处函数趋近幂律 ϕ∝Lα\phi \propto L^{\alpha}。观测上 α\alpha 多在 $-1.0$ 到 $-1.3$ 之间——这个负值意味着越暗的星系越多，矮星系在数目上完全主导。

换成星等形式，L∗L^{*} 对应特征绝对星等 M∗M^{*}。在 bJb_J 波段，2dF 星系红移巡天给出 M∗−5log⁡h≈−19.7M^{*} - 5\log h \approx -19.7、α≈−1.21\alpha \approx -1.21（Norberg et al. 2002，MNRAS 336，907）。

注意符号约定：Schechter 1976 原文把幂律指数写作正的 α=5/4\alpha = 5/4（即 L1/4L^{1/4} 标度的暗端），现代文献普遍改用 LαL^{\alpha} 形式，对应 α≈−1.25\alpha \approx -1.25。读不同论文时要看清楚作者用的是哪种约定。

同一套三个参数，放到不同环境里会给出不同的数。场星系生活在相对稀疏的纤维与空洞里。星系团里则多了潮汐剥离、并合，以及气体被热晕「窒息」的过程。团内光度函数往往在矮端更丰，早型星系的比例更高。因此不能把室女座星系团的计数直接当成宇宙平均。Binggeli、Sandage 与 Tammann 在 1988 年的综述里已经把场与团的差别写成光度函数研究的核心课题（Binggeli, Sandage & Tammann 1988，ARA&A 26，509）。

场样本回答的是宇宙平均密度里「随机抽到一个星系有多亮」。团样本回答的是密环境里成员如何被改造。两者的 L∗L^{*} 可以接近，暗端却常常分叉。把团内过剩的矮椭圆星系写进宇宙平均 α\alpha，会高估场中的矮星系密度。这正是为什么现代拟合要按环境、颜色或恒星质量分开，而不是把所有星系倒进一条 Schechter 曲线。

本地宇宙的定量基石来自广角光谱巡天。斯隆数字巡天自 2000 年起用新墨西哥州阿帕奇角的 2.5 米望远镜同时做成像与红移。Blanton 等人用约 14.8 万个星系，在静止系 $z=0.1$ 的 $r$ 波段给出 Schechter 拟合：ϕ∗=(1.49±0.04)×10−2 h3 Mpc−3\phi^{*}=(1.49\pm 0.04)\times 10^{-2}\,h^{3}\,\mathrm{Mpc}^{-3}，M∗−5log⁡10h=−20.44±0.01M^{*}-5\log_{10}h=-20.44\pm 0.01，α=−1.05±0.01\alpha=-1.05\pm 0.01（Blanton et al. 2003，ApJ 592，819）。这个 α\alpha 比 2dF 的 bJb_J 结果更浅。差别首先来自波段、星等系统和是否计入演化，而不是某一次测量「更正确」。

要看见更暗的星系，需要比 SDSS 主样本更深的光谱极限。星系与质量汇编巡天（GAMA）用英澳望远镜的 AAOmega 光纤光谱仪，在选定天区把极限推到 $r<19.4$ 或 $r<19.8$。Loveday 等人用第一期数据给出 ugriz 光度函数，并发现红星系在 M∗M^{*} 暗端出现凹陷后再变陡，单一 Schechter 不够用，需要双幂律形式（Loveday et al. 2012，MNRAS 420，1239）。这里刻意不把某一期的 α\alpha 写成「宇宙真值」。GAMA 的贡献是把完备度推进到 SDSS 主样本够不着的暗端，并让环境与颜色可以分开拟合。

SDSS 本身是一台专用巡天望远镜加上一套事先规定的选源规则（York et al. 2000，AJ 120，1579）。主星系样本按光纤孔径和星等极限取样，所以暗端斜率对极限附近的完备度极其敏感。GAMA 不追求全天，而追求在较小天区里把光谱做到更暗。两种策略互补：一个给体积和亮端截断，一个给暗端形状与颜色分解。把某一期拟合的 α\alpha 在不同波段之间直接相减，没有物理意义。

从光度函数到质量函数

光度反映的是"现在亮不亮"，会随恒星年龄、尘埃和近期恒星形成而变。要刻画星系本身的多少，天文学家更关心星系恒星质量函数（galaxy stellar mass function，GSMF），即单位体积内各恒星质量区间的星系数目。GSMF 同样能用 Schechter 形式拟合，特征质量 M∗≈1010.6M^{*} \approx 10^{10.6}–1011 M⊙10^{11}\ M_\odot（依红移与巡天而变）。基于 SDSS 的测量给出本地宇宙的 GSMF 在 M∗M^{*} 附近转折，低质量端斜率约 α≈−1.2\alpha \approx -1.2（Baldry et al. 2012，MNRAS 421，621）。GSMF 比 LF 更接近理论量，因为模拟直接预言的是恒星质量，而非随波段变化的瞬时光度。

同一颗星系，在紫外里可能很亮，在近红外里却只是普通成员。紫外追踪年轻大质量恒星，近红外更接近已形成的恒星质量。尘埃消光会把光学光度压低，却把能量再辐射到红外。所以把光度函数译成质量函数，必须经过恒星居族模型与质光比。Baldry 等人的 GAMA 质量函数正是在这条翻译链上工作，而不是把星等直接换算成太阳质量。

质量函数也继承了光度函数的形式问题。低质量端若被低表面亮度星系主导，漏检会让 α\alpha 变浅。高质量端若被亮团中央星系主导，并合与 AGN 反馈会决定截断有多陡。一条 Schechter 曲线同时拟合两端，常常只是一阶近似。

把光度换成质量，还要把「现在这一波段看见的光」翻译成「已经形成的恒星总量」。星族合成模型依赖初始质量函数、金属丰度与星形成史；同一套测光，换一套模型就会平移 M∗M^{*}。因此质量函数的误差预算里，常常有一半来自模型假设，而不是测光噪声。比较不同红移的 GSMF 时，若模型库不统一，演化信号会被系统误差吞掉。这不是反对用质量函数，而是提醒：它比光度函数更接近理论量，也更依赖你看不见的假设。巡天团队如今常同时公布多套质光比假设下的结果，正是为了让读者看清这条翻译链有多脆。

Baldry 等人把 GAMA 质量函数限制在 $z<0.06$，是在用较近的体积换取更暗的质量端。红移一高，同一流量极限对应的绝对星等变亮，暗矮星系立刻掉出样本。所以「本地质量函数」不是一张宇宙平均快照，而是一张被体积与完备度裁过的近邻图。比较不同论文的 M∗M^{*} 时，必须先对齐红移上限、星等系统和质光比假设。

为何是星系演化的定量基石

这里是光度函数真正有深度的地方：它的形状本身就编码了星系形成中两大反馈过程。宇宙学的出发点是暗物质晕质量函数——在冷暗物质（CDM）框架下，由 Press-Schechter 理论（Press & Schechter 1974，ApJ 187，425）及其后续模拟精确预言。晕质量函数同样是低质量端陡、高质量端指数截断。问题在于：如果星系亮度简单正比于晕质量，那么观测到的星系光度函数应当远比实际更陡、更亮端更多。两者对不上。

差异恰恰揭示了反馈：

  • 暗端没有想象中那么陡。低质量晕里，年轻恒星的紫外辐射、超新星爆发驱动的星系风把气体吹散，压制了矮星系的恒星形成效率。这就是为什么观测 α≈−1.2\alpha \approx -1.2，而非暗物质晕那样接近 $-1.9$ 的陡度。
  • 亮端呈指数截断。在大质量晕里，中央超大质量黑洞的活动星系核（AGN）反馈加热气体、阻止其冷却成恒星，给最亮星系的形成踩了刹车。

因此，"光度函数为什么是这个形状"几乎等价于问"反馈如何调节恒星形成"。半解析模型和流体力学模拟（如 EAGLE、IllustrisTNG）把复现观测 LF/GSMF 当作首要检验标准（Schaye et al. 2015，MNRAS 446，521）。

暗端的物理图像可以更具体一点。大质量恒星很快以超新星和恒星风把能量注入星际介质。在浅势阱里，气体一旦被加热就容易逃逸，下一代恒星形成被切断。观测到的矮星系于是比同质量暗物质晕「少亮很多」。这不是计数错了，而是恒星形成效率随质量急剧下降。

亮端则换了一套引擎。中央黑洞吸积时，活动星系核以辐射、喷流和星风加热晕气。气体不能冷却，就不能塌缩成恒星，L∗L^{*} 以上的指数截断由此出现。EAGLE 把复现本地 GSMF 写成项目成功标准之一。IllustrisTNG 给出的较低质量星系丰度，与近期测得的恒星质量函数相符（Pillepich et al. 2018，MNRAS 475，648）。模拟若只匹配大尺度两点相关、却造不出这条曲线，就不能算过关。

复现的意思很具体。不是「看起来有一条下降的尾巴」就算数，而是在同一套宇宙学参数下，同时对上 L∗L^{*} 附近的密度、暗端斜率，以及亮端截断。反馈开太强，亮星系太少；开太弱，矮星系过亮、过密。光度函数因此是一条沿质量切开的诊断刀，而不是拟合报告末尾的附图。大尺度结构形成可以在更大尺度上对，局部的气体物理仍可能错。两条检验缺一不可。

争议与前沿

  • 暗端斜率的真值至今仍有分歧。极暗的矮星系表面亮度低、易漏检，不同巡天给出的 α\alpha 从约 $-1.0$ 到 $-1.5$ 不等。这直接关系到"卫星星系数目偏少问题"（missing satellites）等小尺度 CDM 难题。
  • 高红移演化：JWST 测得 $z > 10$ 的紫外光度函数暗端意外地多、亮端星系也比预期亮，挑战了反馈强度与恒星形成效率的标准假设（Harikane et al. 2023，ApJS 265，5）。这是当前最活跃的争论之一。
  • 环境依赖：星系团内的光度函数和场星系不同——潮汐剥离、并合与"窒息"（气体被剥夺）会改变矮星系的丰度。单一 Schechter 形式只是一阶近似，常需双 Schechter 函数才能同时拟合矮端与巨端。

暗端争议和「失踪卫星」是同一件事的两个面孔。冷暗物质模拟预言，像银河系这样的晕应拥有远多于观测到的明亮卫星（Klypin et al. 1999，ApJ 522，82）。若那些子晕里的恒星系统极其暗弱、表面亮度极低，它们会从光度函数的暗端漏掉，问题就部分变成了完备度。若它们根本没有形成足够的恒星，缺口就指向再电离与超新星反馈。两种解释都要求把低表面亮度星系写进选择函数，而不是把 α\alpha 当成可以任意调节的拟合旋钮。

后来的广角成像确实找到了一批超暗卫星，说明「看不见」曾经是缺口的一部分。这并没有把问题结束。仍须区分：子晕是否发光，以及发光后能否进入星表。银河系结构研究把这些卫星放在我们自己的晕里计数；宇宙学光度函数则在更大体积里问同样的完备度。两条线共用一个物理：低表面亮度星系决定暗端，也决定「失踪」有多少只是没被点亮或没被看见。

把「失踪卫星」只怪到观测完备度，或只怪到反馈，都会把问题说窄。真实缺口往往是两者的卷积：子晕可以存在却几乎不发光，发光了也可以因为表面亮度太低而进不了星表。宇宙学流体模拟若只把卫星的恒星质量推到观测下限附近，却不把选择函数写进比较，仍会高估「看得见的」数目。反过来，只加深成像而不改再电离与超新星处方，也解释不了为何同质量晕的恒星含量可以差一个数量级。暗端斜率因此是一条联合约束，而不是单一旋钮。

詹姆斯韦伯空间望远镜把争论推到了宇宙黎明。Harikane 等人用早期 JWST 数据给出 z∼9z\sim 9–$16$ 的紫外光度函数，暗端比许多前 JWST 模型更丰，亮端也更亮。若这一形状在更大样本里站住，就意味着极早期的恒星形成效率被低估，或者尘埃与初始质量函数的假设需要改。它还没有推翻 Schechter 形式本身。它问的是：高红移的 L∗L^{*} 与 α\alpha 是否仍由本地那套反馈故事外推而来。

高红移紫外光度函数用的是静止系紫外，不是本地的 $r$ 或 bJb_J。直接把 $z>10$ 的 α\alpha 拿来和 Blanton 的 $z=0.1$ 相比，会把波段演化误当成反馈演化。比较必须先把滤光片和星等系统对齐，再谈物理。

环境依赖同样不是边角。把场样本与密环境样本分开后，双 Schechter 函数经常出现：一条分量描述星族序列上的巨星系，另一条托住矮端。形式不够用，本身就是物理信息。它说明「只有一种反馈、一种居族」的假设已经被数据拒绝。下一步不是再加一个自由参数，而是指出哪一种过程在哪一段光度上主导。

跨域连接

  • 暗物质晕：晕质量函数的低质量端接近 $-1.9$，星系光度函数只有约 $-1.2$；两条曲线的落差就是反馈强度沿质量的剖面。推论是暗端斜率并非拟合参数，而是矮星系反馈效率的直接读数——斜率一变，反馈模型必须跟着改；这也意味着晕质量函数不能直接拿来当星系计数用。
  • 统计学：三个参数装下六个数量级，靠的是把标度自由的幂律与标度截断的指数分开写。推论是一旦要同时拟合矮端与巨端就得换成双 Schechter 形式，这说明单一形式隐含了"只有一种反馈"的假设——形式不够用本身就是物理信息，所以函数形式的选择要当成待检验的假设，而不是拟合的技术细节。
  • 不平等经济学：斜率约 $-1.2$ 时，数目积分被暗端主导而光度积分被特征光度附近主导——"数量多"与"总量大"分属两端。推论是问"典型星系有多亮"必须先说明按数目还是按光度加权，两个答案相差几个数量级。
  • 抽样调查：任何巡天都是一次抽样，暗弱星系的入选概率天然更低；不校正入选概率，估出的分布必然偏向"亮的太多"。推论是暗端斜率在不同巡天间从约 $-1.0$ 到 $-1.5$ 的分歧，首先该怀疑选择函数而不是物理，因此这项误差应当以巡天为单位分别报告。
  • 相变与临界现象：幂律加指数截断这一形式在临界现象里也出现，截断尺度对应关联长度。这一步必须写明只是形式类比：它的用处是提示去追问"是什么物理量设定了特征光度"，而不是暗示两者机制相同。

参考文献

  1. Schechter, P. (1976). An analytic expression for the luminosity function for galaxies. Astrophysical Journal, 203, 297–306.
  2. Press, W. H. & Schechter, P. (1974). Formation of galaxies and clusters of galaxies by self-similar gravitational condensation. Astrophysical Journal, 187, 425–438.
  3. Norberg, P. et al. (2002). The 2dF Galaxy Redshift Survey: the bJ-band galaxy luminosity function and survey selection function. MNRAS, 336, 907–931.
  4. Baldry, I. K. et al. (2012). Galaxy And Mass Assembly (GAMA): the galaxy stellar mass function at z < 0.06. MNRAS, 421, 621–634.
  5. Schaye, J. et al. (2015). The EAGLE project: simulating the evolution and assembly of galaxies and their environments. MNRAS, 446, 521–554.
  6. Disney, M. J. (1976). Visibility of galaxies. Nature, 263, 573–575.
  7. Binggeli, B., Sandage, A. & Tammann, G. A. (1988). The luminosity function of galaxies. Annual Review of Astronomy and Astrophysics, 26, 509–560.
  8. Blanton, M. R. et al. (2003). The galaxy luminosity function and luminosity density at redshift z = 0.1. Astrophysical Journal, 592, 819–838.
  9. Loveday, J. et al. (2012). Galaxy and Mass Assembly (GAMA): ugriz galaxy luminosity functions. MNRAS, 420, 1239–1262.
  10. Klypin, A., Kravtsov, A. V., Valenzuela, O. & Prada, F. (1999). Where are the missing Galactic satellites? Astrophysical Journal, 522, 82–92.
  11. Pillepich, A. et al. (2018). First results from the IllustrisTNG simulations: the stellar mass content of groups and clusters of galaxies. MNRAS, 475, 648–675.
  12. York, D. G. et al. (2000). The Sloan Digital Sky Survey: technical summary. Astronomical Journal, 120, 1579–1587.

延伸阅读

  • 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.

正文中引用本条 · 2