跳到主要内容
Sandbox Physics
系外行星的白天与夜晚

AC21 · L4

顺着原文读化学、云与相位地图

Sandbox Physics ·

把三份固定版本原文放在旁边：Bell 与 Cowan 的氢输热模型 arXiv 1802.07725v2，Bell 等人的 WASP-43b 研究 arXiv 2401.13027v1 及其正式发表版，以及 Luger 等人的 starry 论文 arXiv 1810.06559v1。它们依次回答气体带着什么能量、哪一层的光能逃出来，以及望远镜能区分哪些空间图样。

先看氢在什么地方改变形态

阅读 Bell 与 Cowan 时，先读第 2.1 至 2.3 节，再看相位曲线。图 2 很适合停下来细读：在图示压力下，氢的解离比例会在一段超热温度范围里快速变化。白天分子氢拆成原子，储存化学能；流到较冷处重新结合时，能量可以释放回来。

每千克氢的解离能尺度约为 2.14×108 J kg−12.14\times10^8\,\mathrm{J\,kg^{-1}}。一层气体能存多少化学能，要看其质量柱和解离比例。原文相应项的结构是：

Echem=ΣqχH,dEchemdT=ΣqdχHdT.E_{\rm chem}=\Sigma q\chi_{\rm H},\qquad \frac{dE_{\rm chem}}{dT}=\Sigma q\frac{d\chi_{\rm H}}{dT}.
Σ\Sigma 采用原模型的质量柱定义，qq 是单位质量解离能，χH\chi_{\rm H} 为解离比例；要连着原文的组成假设一起使用。

解离比例对温度很敏感时，升温或降温就会同时存入或放出可观的化学能。把第 1 至第 8 式连起来读，也要保留随温度变化的显热项。这里增加的是随物质搬运的储能，能量仍来自吸收与释放的过程。

图 1 至图 4 比较的是指定环流、化学平衡等条件下的理论预测。前面课程的常热容量单层模型没有加入这套机制。选用的 WASP-43b 观测，也没有在这里被解释成探测到了超热氢解离。

对应阅读：Taylor J. Bell and Nicolas B. Cowan, Increased Heat Transport in Ultra-Hot Jupiter Atmospheres Through H2 Dissociation/Recombination · Nicolas B. Cowan and Eric Agol, A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

先读 WASP-43b 的光，再读大气解释

打开 Bell 等人的观测部分与主图 1，先找相位如何定义、食发生在哪里、实测曲线怎样变化。随文 Zenodo 产品保留时间序列及多套处理。教程从这些文件出发，说明时间数组、裁切和共用光子为什么会影响后面的解释。

接着把图 2 中的观测与大气环流预测比较。彩色模型曲线来自算出来的大气，云与环流的处理改变了它们的相位亮度。如果模型夜侧与实测不符，差异就指向大气计算还需要解释的过程。

再读图 3 的分相位光谱。云可以挡住更深、更热处出来的辐射；较弱的垂直温度梯度，也可能让光谱特征变弱。研究因此让多套大气模型同时面对曲线与光谱，逐步追问哪种结构更合适。

看光谱时还要记下波段和误差的定义。十点五微米之外的阴影区，不进入这里采用的主要白光分析。光谱来自同一条时间序列，若再把它当作独立观测乘入同一个似然，就会重复使用光子。

对应阅读：Taylor J. Bell et al., Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b · Taylor James Bell and the 17 other named Zenodo creators, Products and Models for Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

把亮度图连着波长一起读

某个红外波段推得的地图，表示哪里在该波段较亮。要转成气体温度，还需要发射源函数、不透明度、垂直结构和仪器响应。不同波长可能从不同深度出来，因此多波段地图可以提供单条宽带曲线分不开的大气信息。

相位曲线与光谱各自承担一部分工作：前者随观察经度改变，量一段波长里的总光如何起伏；后者比较这份光在波长上怎样分布。论文的亮度温度用黑体来描述光谱发射，并没有给某半球的每一压力层都指定同一个温度。

反照率还涉及另一组波长。Bond 反照率汇总所有方向、所有波长上反射的恒星功率，MIRI 则主要看行星热辐射的一部分。连接两者，需要带误差的大气模型或全波段推断，把缺失的光谱信息纳入计算。

对应阅读：Nicolas B. Cowan and Eric Agol, The Statistics of Albedo and Heat Recirculation on Hot Exoplanets · Taylor J. Bell et al., Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

看遮挡边缘多问了什么

阅读 Luger 等人的论文时，从第 2.4 节的旋转进入第 3.1 至 3.3 节的投影与遮挡。第 16 式转动亮度表示；第 17 至第 20 式把表示连到投影圆盘及未遮挡部分的积分。原文的解析工具比本课程的小规模数值模型更一般。

图 5 可以直观看出不同遮挡轨迹的作用。整盘相位观测会把许多区域加在一起，恒星边缘扫过时，被加起来的子区域却不断改变。这种额外空间权重，可以让相位曲线消去的某些模式重新出现响应。

本指南旁边的独立有限基例子，取到六阶、纬度无关的经度亮度模式，一共十三个系数：一个常数和六对正余弦。在实际保留的观测窗口中，相位矩阵即使采用约八十二度倾角，仍有四个消失方向。对这套基，倾斜本身并未消掉它们。

加入二次食采样后，十三个奇异方向都有非零响应。较细网格的最小归一化奇异值约为 0.006470.00647，而相位观测缺失的模式接近浮点零。两组分辨率帮助判断这个数值响应的稳定程度。这些量描述几何敏感性，不能当成后验探测概率。

任意二维温度图还可以改变纬度结构，亮度到温度的转换也另有未知量。所选有限基满秩回答的是一个有限维线性问题。实际噪声中，弱方向仍可能很难测到。

只考察几何的有限亮度基：一个常数与六对正余弦，均与纬度无关。结果对应十三列，不代表任意二维温度图。
算子纬度与经度格最小归一化响应
侧视相位24,  25624,\;2561.43×10−171.43\times10^{-17}
实际倾角，仅相位24,  25624,\;2561.33×10−171.33\times10^{-17}
实际倾角与食24,  25624,\;2560.006420.00642
侧视相位48,  51248,\;5121.86×10−171.86\times10^{-17}
实际倾角，仅相位48,  51248,\;5121.88×10−171.88\times10^{-17}
实际倾角与食48,  51248,\;5120.006470.00647

对应阅读：Rodrigo Luger et al., STARRY: Analytic Occultation Light Curves · Nicolas B. Cowan and Eric Agol, Inverting Phase Functions to Map Exoplanets

物理先验与看不见的模式

输运模型用微分方程、热容量、气流和加热条件限制地图；平滑先验则偏好相邻亮度变化缓慢。两种限制提供的信息不同，得到的推断也分别依赖它们。

教程里的正三阶谐波地图族给出一个具体例子：图样变了，理想侧视相位曲线仍然相同。热模型还给出另一种多解：沿稳态热容量与流速的退化线，温度场不变。前一种来自观测算子，后一种来自控制方程；仅仅把优化器做得更精确，并不会增加一次对缺失量的观测。

换一篇论文时，可以先写下发光波段、采用的辐亮度规律、观测相位与遮食路径，以及推断预设的地图或大气。再回到预测方程，找到化学能和云不透明度究竟进入哪一步。沿着这些具体选择，就能把不可见的大气同量到的曲线连接起来。

对应阅读：Nicolas B. Cowan and Eric Agol, Inverting Phase Functions to Map Exoplanets · Nicolas B. Cowan and Eric Agol, A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets · Rodrigo Luger et al., STARRY: Analytic Occultation Light Curves

练习与解释性答案

  1. 加入食采样后，选定的十三列地图算子满秩，是否任意二维温度图都能识别？

    展开解释

    这里只考察纬度无关的有限亮度基。任意纬度结构、亮度到气体温度的转换都会再增加未知量。

  2. 夜侧光谱较平，云模型和无云模型都拿来比较。还希望增加什么观测？

    展开解释

    可以扩大光谱覆盖、跨越不同吸收带，或重复不同相位，同时检验温度梯度与云的不透明度；垂直结构的多种变化都可能让光谱变平。

动手检验预测

指出原文增加的物理

  1. 比较冷得更慢的预设与普通向东预设。
  2. 找出两者分别采用的常热容量。
  3. 说明氢化学论文中的哪些储能项，需要进入原有热容量描述。

来源与计算范围

  1. Nicolas B. Cowan and Eric Agol, A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

    2010 年 11 月 1 日第一版。第 2–4 节、第 1–10 与 13 式，图 1 与图 6；第 6 节讨论近似。使用给定平流与辐射冷却，不求解风的动量。

  2. Nicolas B. Cowan and Eric Agol, Inverting Phase Functions to Map Exoplanets

    2008 年 3 月 25 日第一版。第 2 节第 1–4 式，第 3.2 节第 7 式与图 1。侧视、各向同性、无食时的经度卷积与 Fourier 零空间。

  3. Nicolas B. Cowan and Eric Agol, The Statistics of Albedo and Heat Recirculation on Hot Exoplanets

    2010 年 12 月 20 日第二版。第 1 节、第 2.1–2.2 节与第 3 式。有限波段亮度温度、全波段能量与反照率的关系。

  4. Taylor J. Bell and Nicolas B. Cowan, Increased Heat Transport in Ultra-Hot Jupiter Atmospheres Through H2 Dissociation/Recombination

    2018 年 4 月 9 日第二版。第 2.1–2.3 节第 1–8 式，图 1–4。氢解离与复合的理论输热模型，未作为常热容量演示里的化学求解。

  5. Taylor J. Bell et al., Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

    2024 年 1 月 23 日预印本第一版；另核读 2024 年 4 月 30 日正式发表版。主图 1–4；观测、裁切、阴影区与 Eureka v1 拟合方法及数据来源。公开产品固定为 Zenodo v1。

  6. Rodrigo Luger et al., STARRY: Analytic Occultation Light Curves

    2018 年 10 月 15 日第一版。第 2.4 节第 16 式，第 3.1 节第 17–20 式，第 3.2–3.3 节与第 4.2 节图 5。旋转、投影及不同遮挡轨迹的地图敏感性。

  7. Taylor James Bell and the 17 other named Zenodo creators, Products and Models for Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

    第 10525170 号记录第一版，2024 年 4 月 30 日发布。核读原 API 许可与十八位创建者、实际 HDF5 数组属性和处理配置，保留选定原文件字节及逐文件 SHA-256。

反馈这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策