--- title: "顺着原文读化学、云与相位地图" description: "逐式逐图阅读固定版本原文，比较氢化学输热、WASP-43b 的云与遮食地图。" canonical_url: "https://sandboxphysics.com/zh/stories/planet-phase-curves/clouds-chemistry-and-maps/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/planet-phase-curves/clouds-chemistry-and-maps/" last_modified: "2026-10-28" --- [系外行星的白天与夜晚](https://sandboxphysics.com/zh/stories/planet-phase-curves/) AC21 · L4 # 顺着原文读化学、云与相位地图 Sandbox Physics · 2026-10-28 把三份固定版本原文放在旁边：Bell 与 Cowan 的氢输热模型 arXiv 1802.07725v2，Bell 等人的 WASP-43b 研究 arXiv 2401.13027v1 及其正式发表版，以及 Luger 等人的 starry 论文 arXiv 1810.06559v1。它们依次回答气体带着什么能量、哪一层的光能逃出来，以及望远镜能区分哪些空间图样。 ## 开始前需要什么 先读气团与投影两章，把链接的原文 PDF 放在旁边。 ## 先看氢在什么地方改变形态 阅读 Bell 与 Cowan 时，先读第 2.1 至 2.3 节，再看相位曲线。图 2 很适合停下来细读：在图示压力下，氢的解离比例会在一段超热温度范围里快速变化。白天分子氢拆成原子，储存化学能；流到较冷处重新结合时，能量可以释放回来。 每千克氢的解离能尺度约为 $2.14\times10^8\,\mathrm{J\,kg^{-1}}$。一层气体能存多少化学能，要看其质量柱和解离比例。原文相应项的结构是： $$ E_{\rm chem}=\Sigma q\chi_{\rm H},\qquad \frac{dE_{\rm chem}}{dT}=\Sigma q\frac{d\chi_{\rm H}}{dT}. $$ $\Sigma$ 采用原模型的质量柱定义，$q$ 是单位质量解离能，$\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.00647$，而相位观测缺失的模式接近浮点零。两组分辨率帮助判断这个数值响应的稳定程度。这些量描述几何敏感性，不能当成后验探测概率。 任意二维温度图还可以改变纬度结构，亮度到温度的转换也另有未知量。所选有限基满秩回答的是一个有限维线性问题。实际噪声中，弱方向仍可能很难测到。 只考察几何的有限亮度基：一个常数与六对正余弦，均与纬度无关。结果对应十三列，不代表任意二维温度图。 | 算子 | 纬度与经度格 | 最小归一化响应 | | --- | --- | --- | | 侧视相位 | $24,\;256$ | $1.43\times10^{-17}$ | | 实际倾角，仅相位 | $24,\;256$ | $1.33\times10^{-17}$ | | 实际倾角与食 | $24,\;256$ | $0.00642$ | | 侧视相位 | $48,\;512$ | $1.86\times10^{-17}$ | | 实际倾角，仅相位 | $48,\;512$ | $1.88\times10^{-17}$ | | 实际倾角与食 | $48,\;512$ | $0.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. 夜侧光谱较平，云模型和无云模型都拿来比较。还希望增加什么观测？ 展开解释 可以扩大光谱覆盖、跨越不同吸收带，或重复不同相位，同时检验温度梯度与云的不透明度；垂直结构的多种变化都可能让光谱变平。 ## 动手检验预测 ### [指出原文增加的物理](https://sandboxphysics.com/zh/labs/exoplanet-thermal-phase-curve/) 1. 比较冷得更慢的预设与普通向东预设。 2. 找出两者分别采用的常热容量。 3. 说明氢化学论文中的哪些储能项，需要进入原有热容量描述。 [上一篇：重跑 WASP-43b 的 MIRI 相位曲线](https://sandboxphysics.com/zh/stories/planet-phase-curves/replay-wasp43-miri/) ## 来源与计算范围 1. [Nicolas B. Cowan and Eric Agol, A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets](https://arxiv.org/pdf/1011.0428v1) 2010 年 11 月 1 日第一版。第 2–4 节、第 1–10 与 13 式，图 1 与图 6；第 6 节讨论近似。使用给定平流与辐射冷却，不求解风的动量。 2. [Nicolas B. Cowan and Eric Agol, Inverting Phase Functions to Map Exoplanets](https://arxiv.org/pdf/0803.3622v1) 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](https://arxiv.org/pdf/1001.0012v2) 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](https://arxiv.org/pdf/1802.07725v2) 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](https://arxiv.org/pdf/2401.13027v1) 2024 年 1 月 23 日预印本第一版；另核读 2024 年 4 月 30 日正式发表版。主图 1–4；观测、裁切、阴影区与 Eureka v1 拟合方法及数据来源。公开产品固定为 Zenodo v1。 6. [Rodrigo Luger et al., STARRY: Analytic Occultation Light Curves](https://arxiv.org/pdf/1810.06559v1) 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](https://zenodo.org/records/10525170) 第 10525170 号记录第一版，2024 年 4 月 30 日发布。核读原 API 许可与十八位创建者、实际 HDF5 数组属性和处理配置，保留选定原文件字节及逐文件 SHA-256。