跳到主要内容
Sandbox Physics
从行星轨道读到凌星计时

L3 · 凌星计时变化

为什么测时波很清楚，质量仍不确定

Sandbox Physics ·

一个测时信号可以很好测，却不好解释。慢振幅同时依赖质量和偏心率组合；较短的会合分量，有时能补上缺少的独立信息。

先把角位移换成时间

参考遮星附近，一点角度偏差会把事件沿时钟移动。把角位置对时间作一阶展开，未受扰角速度就乘在时刻偏差上。角度若领先，遮星时间便提前，对应负的时间延迟。

摄动不仅改变平均经度，也改变相对平均经度的偏心位移。后一项在小偏心率时也有贡献。若只把平均经度变化换成时间，就会漏掉部分质量信息。

0≃n δt+δθ,δt=−δθn0\simeq n\,\delta t+\delta\theta,\qquad \delta t=-\frac{\delta\theta}{n}
角度相对于参考遮星量起，并把观测方向定为经度零点；一阶展开要求位移较小。

对应阅读：Lithwick · Deck & Agol

分清被驱动的偏心率与自由部分

偏心率既有大小，也有指向近心点的方向。把两个分量放进一个复数，就能同时记录它们；乘一个复指数相当于旋转向量，省去反复书写两套三角函数。

接近一阶共振时，偏心率向量包含共振驱动的部分，以及由初始轨道选定的自由部分。驱动项随缓慢的会合经度转动；在这个近似里，自由部分要到更长的世俗时间尺度上才明显改变。

领先阶的慢测时只看得到两颗行星自由向量的一个加权和，无法分别测出四个实分量。即使完全没有测量噪声，这条公式也有这样的信息限制。

z=e eiϖ=zfree+zforced,Zfree=fz1,free+gz2,freez=e\,e^{i\varpi}=z_{\rm free}+z_{\rm forced},\qquad Z_{\rm free}=f z_{1,\rm free}+g z_{2,\rm free}
虚数单位记录第二个分量，近心点经度决定方向；系数来自引力摄动函数，二比一附近的领先值约为负一点一九零与零点四二八四。

对应阅读：Lithwick · Deck & Agol

从振幅里读出退化

用两个平均周期定义偏离二比一的相对量。复振幅同时储存慢测时的大小和相位；质量因子乘着一个含自由偏心率的括号。改变质量，就可能用偏心率变化补偿。

同时测两颗行星当然有帮助，但它们仍共享同一个偏心率组合，允许的质量区域可能是一条弯曲长脊。把自由偏心率固定为零，只是从这条脊上截取一片，并未证明数据测到了零偏心率。

Δ=P22P1−1,ϕ=2λ2−λ1,δti=Im⁡(Vieiϕ)\Delta=\frac{P_2}{2P_1}-1,\quad \phi=2\lambda_2-\lambda_1,\quad \delta t_i=\operatorname{Im}(V_i e^{i\phi})
周期与凌星起点共同规定相位；教程让这些历表参数联合变化，而不把它们当成无误差输入。
V1≃P1μ2π 22/3Δ(−f−3Zfree∗2Δ),V2≃P2μ12πΔ(−g+3Zfree∗2Δ)V_1\simeq\frac{P_1\mu_2}{\pi\,2^{2/3}\Delta}\left(-f-\frac{3Z_{\rm free}^{*}}{2\Delta}\right),\qquad V_2\simeq\frac{P_2\mu_1}{2\pi\Delta}\left(-g+\frac{3Z_{\rm free}^{*}}{2\Delta}\right)
这是接近共振的领先阶表达式；质量比以恒星质量归一化，上标星号表示复共轭。接近精确共振时，这个展开失效，须改用能解析共振的动力学计算。

对应阅读：Lithwick · Deck & Agol

找一个以不同方式依赖质量的分量

小偏心率下，会合引起会合角及其整数倍的谐波，振幅由伴星质量和轨道间距函数决定。若分辨出一个没有混叠的分量，就能给质量增加一条不带相同自由偏心率因子的约束。

凌星观测像频闪照片，每绕一圈才取一个点。接近二比一时，内行星的第二个会合谐波会落到采样后的慢频率上，外行星的第一个谐波也会如此。它们已经包含零自由偏心率的共振贡献；把完整共振公式再加上去，就会把这部分算两遍。

教程采用 Deck 与 Agol 的拆法：先算零自由偏心率的会合项，再只加慢频率处的自由偏心率修正。这里最干净的额外质量信息来自内行星第一会合谐波。这套近似仍须与实际积分产生的凌星时刻比较。

δt1,chop=P12πμ2∑j≥1f1(j)(α)sin⁡[j(λ1−λ2)],α=a1/a2\delta t_{1,\rm chop}=\frac{P_1}{2\pi}\mu_2\sum_{j\ge1}f_1^{(j)}(\alpha)\sin[j(\lambda_1-\lambda_2)],\qquad \alpha=a_1/a_2
无量纲系数由收敛的拉普拉斯系数积分得到；轨道过近、偏心率过大或真正进入共振，都可能让谐波解释失效。

对应阅读：Deck & Agol · Lithwick

在一条完整轨道里量出小分量

采用保存的默认模型：伴星二十地球质量、初始周期比二点零四、内行星为圆、伴星初始经度一弧度，积分初始内行星二百圈。对一百九十九个内侧事件，在拟合钟表的同时，按伴星平均频率拟合多个谐波。第一谐波的半振幅约零点四四三分钟，远小于完整残差约二十九点五分钟的极差。

独立计算的领先阶会合系数预言约零点四五六分钟，差异为百分之二点六七，包含有限质量、有限时窗和更高阶偏心效应。采用四、八、十二个谐波时，第一分量保持稳定。这是对实际合成牛顿事件的定量近似检查，不是一项新的真实行星质量测量。

保存的合成量结果
初始二百圈内的内侧事件数一百九十九
完整测时极差二十九点五分钟
拟合第一谐波半振幅零点四四三分钟
领先阶摄动预言零点四五六分钟

对应阅读：Deck & Agol · Deck et al.

完整做一个小反问题

取一个算术模型：以地球为单位的质量，乘二分之一加十倍的无量纲偏心组合，得到分钟振幅。十分钟的慢振幅可以来自质量十、偏心组合零点零五；也可以来自质量二十、偏心组合零；或质量四十、偏心组合负零点零二五。这只是连续退化上的三个点。

现在再测一个独立分量，系数为每地球质量十分之一分钟，结果为两分钟，高斯标准误为五分之一分钟。条件于这个精确系数，质量就是二十、标准误两个地球质量。代回慢振幅，才开始约束偏心组合。

这道练习说明信息方向如何相交，采用的是人为设定的算术系数，不是 Kepler-289 标定。物理拟合中，周期、恒星质量、谐波系数以及模型本身的可靠程度，也会影响答案。

A=m(0.5+10z)  min,B=0.1m  minA=m(0.5+10z)\;\mathrm{min},\qquad B=0.1m\;\mathrm{min}
质量变量取以地球为单位的数值，另一变量是实数教学参数。两次观测约束质量与偏心平面上的不同方向。

对应阅读：Lithwick · Deck & Agol

把不确定性放进似然

似然描述：给定一个模型，这组测得时刻具有多大的概率密度。若误差相互独立且服从高斯分布，每个残差平方就除以相应方差。再拟合额外白噪声散布时，它的方差与原误差方差相加；总方差的对数项也必须保留，否则把散布无限加大就能白白改善拟合。

后验把似然与参数的正规化先验合起来。先验说明，使用这组时刻之前允许哪些质量和偏心率。把不关心的参数积分掉叫边缘化；把它们换成各自最佳值，得到的则是轮廓似然，两者不同。

真实数据计算联合采样两个质量、自由向量的两个分量、四个历表参数和一个正的额外散布。教程保存先验、独立采样群诊断，以及加宽偏心率先验后的重算。很大的额外散布，或者新出现而未充分探索的模式，都要求重新审视简化模型和质量区间的稳健程度。

ln⁡L=−12∑n[(tn−tn,model)2σn2+s2+ln⁡ ⁣(2π(σn2+s2))],p(ϑ∣t)∝L p(ϑ)\ln L=-\frac12\sum_n\left[\frac{(t_n-t_{n,\rm model})^2}{\sigma_n^2+s^2}+\ln\!\bigl(2\pi(\sigma_n^2+s^2)\bigr)\right],\qquad p(\vartheta\mid t)\propto L\,p(\vartheta)
额外散布与原测时误差用相同时间单位。这个对角似然省略了未公布的时间协方差，结果明确条件于近似两行星模型。

对应阅读：Greklek-McKeon et al. · Deck & Agol

亲手观察

区分轨道敏感性与质量反演

  1. 先比较偏心轨道与圆轨道默认值，再改变质量。
  2. 改变间距，观察完整测时图案的哪些部分移动，别把大残差直接当成唯一质量。
  3. 继续读实测时刻教程，把条件后验与残差诊断一起比较。

练习与答案

拟合把自由偏心率全部固定为零，给出很窄的质量区间。先问什么？

先问数据是否独立约束了偏心率。固定它是在假设中移除退化；解释窄区间之前，应比较自由向量联合变化的结果，并检验先验敏感性。

为什么不能把完整二比一共振式加在全部会合谐波上？

凌星采样把内行星第二、外行星第一会合谐波混叠到慢频率；会合和式中已经包含它们的零自由偏心率贡献。只能再加入自由偏心率修正，并统一经度零点。

原测时误差很小，模型却需要十五分钟额外散布。窄条件质量区间能解决真实质量问题吗？

它描述的是带额外散布项的近似模型。显著差异要求更完整的动力学或光度模型，或解释散布来源；相应的物理准确度需要能够解释观测残差的模型。

来源与计算范围

  1. Lithwick, Xie & Wu, Extracting Planet Mass and Eccentricity From TTV Data (2012)

    固定 v2，第二节、公式一至十及附录 B：自由偏心率向量、近共振近似的适用条件。

  2. Deck & Agol, Measurement of Planet Masses with Synodic Chopping (2015)

    固定 v1，第二至第五点一节、公式十至十五及附录 A：凌星采样如何混叠分量，PH3 案例如何比较加入会合踢动前后的质量约束。

  3. Deck et al., TTVFast (2014)

    固定 v1，第二点一点三、第二点二至三和第三节：几何凌星中心取前景天体的星心投影距离极小值。教学积分复用本站 RK4 核并独立比较 IAS15，不移植 GPL TTVFast 代码。

  4. Greklek-McKeon et al., Kepler-289 Densities and Transit Timing (2023)

    固定 v1，第二至三节、表一及图二至四、九至十，许可为 CC BY 4.0。作者先联合拟合遮星形状，再用三行星动力学拟合测得的时刻。

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

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