[论文解读] The Atmosphere and Interior Structure of HAT-P-13b from Spitzer Secondary Eclipses

R. A. Hardy, Joseph Harrington+7 more|arXiv (Cornell University)|Jan 3, 2017

Stellar, planetary, and galactic studies参考文献 39被引用 5

一句话总结

本研究利用斯皮策空间望远镜在3.6 μm和4.5 μm波段对热木星HAT-P-13b的二次凌日观测，精确测定其轨道偏心率，并通过第二阶洛夫数$k_2$约束其内部结构。结果得到$k_2 = 0.81 \pm 0.10$，表明核心质量较低甚至可能无核心，挑战了现有模型，并对双行星系统中近日点定向一致性的假设提出疑问。

ABSTRACT

We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a two-planet system with a configuration that enables constraints on the planet's second Love number, \math{k\sb{2}}, from precise eccentricity measurements, which in turn constrains models of the planet's interior structure. We exploit the direct measurements of \math{e \cos ω} from our secondary-eclipse data and combine them with previously published radial velocity data to generate a refined model of the planet's orbit and thus an improved estimate on the possible interval for \math{k\sb{2}}. We report eclipse phases of \math{0.49154 \pm 0.00080} and \math{0.49711 \pm 0.00083} and corresponding \math{e \cos ω} estimates of \math{-0.0136 \pm 0.0013} and \math{-0.0048 \pm 0.0013}. Under the assumptions of previous work, our estimate of \math{k\sb{2}} of 0.81 {\pm} 0.10 is consistent with the lower extremes of possible core masses found by previous models, including models with no solid core. This anomalous result challenges both interior models and the dynamical assumptions that enable them, including the essential assumption of apsidal alignment. We also report eclipse depths of 0.081\% {\pm} 0.008\% in the 3.6 {\micron} channel and 0.088 \% {\pm} 0.028 \% in the 4.5 {\micron} channel. These photometric results are non-uniquely consistent with solar-abundance composition without any thermal inversion.

研究动机与目标

• 利用斯皮策望远镜的精确二次凌日时间观测，改进对热木星HAT-P-13b内部结构的约束。
• 从凌日相位测量行星的偏心率和近日点角($e\cos\omega$)，以推断其第二阶洛夫数$k_2$。
• 检验在双行星系统中估算$k_2$所依赖的近日点定向一致性假设的有效性。
• 从全波段凌日深度推断大气性质，包括成分和热反转层结构。
• 解决斯皮策望远镜3.6 μm与4.5 μm通道之间凌日时间差别的问题，该现象可能源于相关噪声或亮度分布不对称。

提出的方法

• 斯皮策空间望远镜使用子阵列模式，在IRAC通道1（3.6 μm）和通道2（4.5 μm）观测了HAT-P-13b的两次二次凌日，每组观测包含68,608帧。
• 利用POET数据处理流程提取光度，该流程建模了凌日引起的通量变化、像素内灵敏度差异以及随时间变化的灵敏度（即“阶梯”效应）。
• 凌日相位测量结果为：3.6 μm波段为0.49154 ± 0.00080，4.5 μm波段为0.49711 ± 0.00083，由此得出$e\cos\omega = -0.0136 \pm 0.0013$和$-0.0048 \pm 0.0013$。
• 将测得的$e\cos\omega$值与径向速度数据结合，基于潮汐平衡和近日点定向一致的假设，优化轨道解并估算$k_2$。
• 通过拟合大气模型到凌日深度，假设太阳丰度和高效热量再分配，推断大气热结构与成分。
• 采用贝叶斯分析比较定向一致与非一致的近日点构型，评估$k_2$估算的稳健性。

实验结果

研究问题

• RQ1 基于斯皮策二次凌日时间观测，HAT-P-13b的轨道偏心率和近日点角的精确值是什么？
• RQ2 在双行星系统中，假设近日点定向一致，能否可靠估算HAT-P-13b的第二阶洛夫数$k_2$？
• RQ3 3.6 μm与4.5 μm波段之间观测到的23分钟凌日时间差是否由相关噪声、亮度分布不对称或物理进动引起？
• RQ4 测得的凌日深度对HAT-P-13b的大气成分和热结构有何启示？
• RQ5 这些结果如何挑战或支持当前低核心质量热木星内部结构模型？

主要发现

• 二次凌日相位在3.6 μm波段为0.49154 ± 0.00080，在4.5 μm波段为0.49711 ± 0.00083，表明存在非零但极小的偏心率。
• $e\cos\omega$值在3.6 μm波段为-0.0136 ± 0.0013，在4.5 μm波段为-0.0048 ± 0.0013，与低偏心率和非零近日点角一致。
• 第二阶洛夫数$k_2$估算为0.81 ± 0.10，与低核心质量或无核心模型一致，挑战了标准内部结构模型。
• 凌日深度在3.6 μm波段为0.081% ± 0.008%，在4.5 μm波段为0.088% ± 0.028%，与太阳丰度成分一致，且无热反转层存在迹象。
• 两通道间23 ± 4分钟的凌日时间差无法由长期近日点进动解释，可能源于相关噪声或亮度分布不对称。
• 对定向一致与非一致近日点构型的贝叶斯比较结果不明确，但数据不排除显著的摆动现象，对标准假设下$k_2$测量的有效性提出质疑。