2.4 辐射传输的数值旅程

2.4 辐射传输：光子的数值旅程

本节摘要：辐射转移方程活在位置、方向、频率的相空间里，维数本身就是墙。恒星大气用平面平行和两流近似先把墙拆窄；矩方法把方向积掉，换成闭合假设；蒙特卡罗让光子包去随机走；短特征射线在网格上沿特征积分。再电离、尘埃红外、吸积盘光深，各自挑选过不同的旅行方式。辐射从来不是“再加一个冷却曲线”那么便宜。

上手前先明确

阅读完本节，你应当能够：

1. 说明辐射转移为何比流体多出方向与频率两个负担
2. 区分两流、矩闭合、蒙特卡罗、射线追踪各自让出了什么
3. 解释光时标与流体时标的分裂，为何常把辐射做成隐式
4. 指出再电离模拟里电离锋分辨率与光子守恒为何比漂亮的气泡图更硬
5. 判断一套星系码的“辐射反馈”是在解转移方程，还是在用逃逸分数参数化

相空间这堵墙

流体场在每个空间点上是几个守恒量。辐射场在每个空间点上还要记得：光子从哪来、频率是多少。完整的转移方程因此是七维的（三维空间、二维方向、频率、时间）。解析解出现在平面平行大气、灰近似、局域热动平衡这些被削过的巷子里。恒星光谱学靠这些巷子活了一个世纪。一旦对象变成三维星风、弯曲电离锋、尘埃红外再处理，巷子不够走。

历史选择不是“上最完整的方程”，而是一次次问：我能不能少记一个维度。两流近似把方向收成向外和向内两束，恒星结构的辐射输运长期靠它。扩散近似在光深深的地方把辐射当导热。光学薄则干脆改成冷却函数，光子当场消失，不再有方向。每一种省略都在买速度，卖 nonlocal 的影子：影子里藏着电离锋的形状、喷流空腔的光照、分子云表面的加热。

完整负担                     常见省略                     卖出去的东西
  位置                       平面平行、球对称             三维阴影
  方向                       两流、矩、各向同性扩散       光束与阴影
  频率                       灰、多群、单能               线冷却与电离结构
  时间                       稳态转移                     光走时间与闪变
​

💡 关键直觉：冷却函数是辐射的骨灰盒。有用，但里面已经没有方向，也就没有影子。

四种旅行方式，以及它们占领过的领地

矩方法对方向积分，留下能量密度、通量、压力张量，再靠闭合把张量写成低阶矩的函数。Eddington 因子、M1 闭合，都是在猜光子的各向异性有多强。M1 能走光束的大致方向，却会在两束光交叉时犯傻，把交叉点合成错误的方向。再电离和点源照明因此要非常小心。矩方法便宜，能隐式耦合到流体，是许多宇宙学辐射流体的工作马。

蒙特卡罗让有限个光子包去抽样发射、散射、吸收。实现直观，阴影自然，频率线也能做。噪声按包数的根号下降，稀薄区或高光深区都可能贵得要命。方差降低技巧、离散扩散、隐式蒙特卡罗，都是在拯救根号。超新星光谱、尘埃辐射转移、某些吸积盘反射，蒙特卡罗是主力。

短特征或长特征射线在网格上沿方向角积分源函数。比完整离散纵标更省，仍能抓住阴影。恒星三维大气、星风、某些电离模拟用它。射线的方向离散本身会留下条纹伪像，加密角度和加密空间一样是分辨率声明。

射线追踪的另一种极端是把源当成点，沿光线做光学深度积分，适合数量不多的星或类星体。宇宙学再电离早期常用这种“源少、介质多”的设定。源一旦变成遍布的恒星粒子，射线数量爆炸，又得回到矩或动量耦合的近似。

与流体结婚：光时标、光子守恒、逃逸分数

辐射流体的时间分裂是方法史上反复改道的河段。光穿过一个格子可以远快于流体声速，显式把两者锁在同一 CFL 上，等于让流体为光子的时钟买单。常见做法是：流体显式，辐射隐式或稳态；或者对辐射能量做隐式扩散，动量沉积再反馈给气体。隐式会引入求解器和线性系统的新墙。稳态假设则放弃闪变和光走时间——对再电离的宇宙学时标通常可接受，对伽马射线暴火球不可接受。

光子守恒是再电离的硬账。电离一个氢原子需要一个大于阈能的光子，复合会吐回。数值上若光子在格子间泄漏、或源的注入与吸收对不上，电离分数会系统地偏快或偏慢。气泡图可以仍然好看。好的再电离码把光子计数当成和流体质量同等的诊断。读者若只看到彩色的电离区，看不到光子预算，应默认那是可视化而不是闭合的实验。

星系形成里还有一种更便宜、也更容易被说成辐射的东西：逃逸分数。它不把转移方程解完，只规定有多少电离光子离开星形成区。参数可以拟合观测的再电离历史，却把星团内部的吸收、尘埃、空洞几何全部压缩成一个数。这不是作弊，是亚网格。作弊发生在把它写成“我们包含了辐射传输”。

⚠️ 常见坑：把多群扩散或单能矩方法的结果，直接当成光谱预测。频率被压缩之后，线特征和电离结构可以整体平移。

尘埃把问题再抬一维。吸收后的再发射改频率，红外光学深度可以很大，紫外阴影和红外各向同性加热同时存在。必须允许光子改能量，蒙特卡罗和某些多群矩才够格。只在紫外上跑电离、把红外当冷却曲线，会漏掉盘和分子云里真正决定温度的那一段旅程。

概念账本：

选省略：灰 / 多群 / 线；两流 / 矩 / 射线 / 光子包
每个流体步：
    更新不透明度
    隐式或稳态求辐射场
    把加热、电离、辐射力写回气体
诊断：光子预算、电离分数、辐射能量密度
对照：光学薄极限、扩散极限、已知阴影测试
​

辐射主线至此与流体、引力并列。下一节不再发明新的物理场，而问零件如何焊接：算子分裂把模块养活，也把误差安顿在接口上。

恒星大气留下的遗产，以及三维如何把它花光

经典恒星大气把转移方程训练成一门精密手艺：线表、模型大气、光谱合成。平面平行和局域热动平衡在那里常常够用。三维对流大气把这门手艺拖进短特征射线和随时间变化的源函数。谱线不对称、粒状结构、中心到边缘的变化，成为三维辐射流体的观测接口。遗产是：频率不能随便灰掉，若你的科学问题是某条线。

再电离与星系风把遗产花在另一头：源多、介质块、需要的是电离分数和辐射压，而不是逐线轮廓。多群或单能矩方法成为工作马。两种传统互相轻视过——光谱学家觉得宇宙学辐射太粗，宇宙学家觉得逐线买不起。轻视没有帮助。正确的读法是：光谱学规定频率何时不能省，宇宙学规定方向何时可以积掉。跨过各自的红线，产品就不能再冒充对方的观测。

尘埃红外再处理是当代最容易被低估的账单。紫外光子在星团附近被吸收，红外在更大体积里以接近扩散的方式加热。只做紫外电离的模拟，分子云温度可以整体偏冷或偏热，星形成阈值跟着走。若论文的科学问题包含尘埃温度或远红外光度，冷却函数骨灰盒不够用，必须允许光子改频率。

电离锋、阴影测试与光子预算的三种死法

标准阴影测试：点源、光厚障碍、后面应有影。扩散近似会把影填亮；M1 可能在影边给出错误的光束合并；蒙特卡罗在包数不足时让影闪烁。任何声称能处理再电离几何的码，都应公开这类测试。电离锋的厚度若只有一个格子，锋的速度可以是数值的，而不是复合与光子到达率的竞速。加密格子直到锋厚收敛，才谈气泡统计。

光子预算的三种死法：源注入与模块内吸收对不上；流体平流把电离物质带走却不带光子账；频率群过粗导致阈能附近的光子被错误归类。三种死法都能留下好看的气泡。把预算表与电离分数时间史一起画，比多一张体积渲染更有说服力。逃逸分数若被用作自由参数去拟合再电离光学深度，应进入校准名单，不能再当预测。

问题：星系形成论文写了辐射反馈，是不是已经解了转移方程？

多半没有。常见的是逃逸分数、辐射压单散射、或单能矩闭合。这些可以正当，但它们省略了方向或频率，因而省略了影子或线。真正解转移的标志是：有方向离散或光子包、有光子预算表、有阴影测试。若科学问题是再电离几何或尘埃温度，冷却函数和逃逸分数不够。若科学问题只是粗略抑制冷却，参数化可以，但必须进入校准名单，不能再当预测。

要点串联

• 相空间维数是墙：方向与频率不是装饰
• 冷却函数没有影子：省略方向就是省略 nonlocal
• 四种旅行方式占领不同领地：闭合误差与根号噪声是不同的税
• 光时标通常强迫辐射隐式或稳态：否则流体为光子付时钟
• 光子守恒是再电离的硬账：气泡图不能代替预算
• 逃逸分数是亚网格：可以正当，不能冒充已解转移方程

下一节看焊接：重力、流体、冷却、辐射、磁场如何用算子分裂拼成一台能跑、也会在接口处漏油的机器。