从天琴座β星的周期增长算起

先确定今天要复算哪件事

天琴座β星，也叫渐台二，是一个有复杂吸积结构的相互作用双星。周期变化提供了很好的入手点，但这个“时钟”并不是单凭一张示意光变图得到的。Harmanec与Scholz在1993年综合历史供体视向速度，在处理仪器和谱线系统差异后，给出二次历式；食时结果则作为外部比较。

本教程采用后续Mennickent与Djurašević 2013年v1文本的表1拟合质量，以及第4节采用的周期增长，复算其保守转移估计。盘的光变拟合、视向速度原始拟合和演化轨迹拟合都不在这一步里。把这些来源层次分开，才能知道结果究竟依赖哪组输入。

先下载本页的独立科学ZIP，在一个新目录解压。包中是自写科学代码、自有合成检验及模型聚合材料。原文、原生结构和上游软件文件由具名下载器留在你的电脑上；原始观测表保留在具名上游来源。

对应阅读：Harmanec & Scholz (1993), A&A279,131–147 · Mennickent & Djurašević (2013), arXiv1303.5812v1

取固定文本，而不是网页上随手抄一组质量

python verify_package.py
python fetch_sources.py --output inputs
python reproduce.py --inputs inputs --output results
python validate_science.py

先运行包内校验，再取清单中的具名上游响应。下载器检查精确大小和SHA-256，只解出两个指定的原生结构成员；若缓存损坏，就停止，不静默换成另一个版本。2013年原文使用固定v1的完整HTML，文件哈希如下，复算程序会再次检查这个值。

7c218201125132c16ea8265ac1469728561f54d40fad047d57787dd3638741bf

原文表1给出供体 $M_d=2.97M_odot$ 、受体 $M_a=13.16M_odot$ ；第4节采用 $P=12.9,mathrm d$ 和 $dot P=19,mathrm{s,yr^{-1}}$ 。程序在新取的文本中定位这些字段，而不是读取本页预先算好的速率。摘要中的圆整数字与具体表格不是同一组输入，不能混着用。

MESA固定版本的原生算例写明了混合角动量流失关系；学校剖面则用于后面的重力积分。这里分别检查源码版本和材料许可，不把“可下载”当成所有原文和图表都可再发布的许可。原生结构派生材料保留各自CC BY或CC BY-SA条款。

对应阅读：Mennickent & Djurašević (2013), arXiv1303.5812v1 · MESA 25.12.1, jdot_ml_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899 · Thoul & Chidester (2022), MESA Summer School; Zenodo 7082536 · Jared Goldberg (2025), MESA School Leuven Friday model

让周期变化变成一个速率

令 $dot x>0$ 表示供体每年送出的质量。保守、圆轨道、忽略其他轨道力矩时，总质量和轨道角动量不变。把上一章的间距响应代入开普勒周期，得到：

{\dot P\over P}=3\dot x\left({1\over M_d}-{1\over M_a}\right),\qquad\dot x={{\dot P}/P\over3(1/M_d-1/M_a)}

供体比受体轻，所以括号为正；正的周期增长对应正的供体流失率。

{\dot P\over P}={19\over12.9\times86400}\,\mathrm{yr}^{-1},\qquad\dot x\simeq2.17955\times10^{-5}M_\odot\,\mathrm{yr}^{-1}

周期先换成秒，才能与每年增加的秒数相除。结果与原文圆整到两位的估计相容；末尾数字只是这组定参的计算精度。

结果中的渐台二汇总应给出上述数值。沿输入字段与单位换算读一遍计算，再看流失通道情景。不要把两份速率接近误读成又独立测量了一次质量转移；二者共享同一组已发表的质量与周期趋势。

对应阅读：Mennickent & Djurašević (2013), arXiv1303.5812v1 · Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1

换掉保守假设，保留无正解的情况

{\dot P\over P}=\dot x\left[{3\over M_d}-{3b\over M_a}-{(1-b)(3\gamma+1)\over M_d+M_a}\right]

比例

b

由受体留下，逃逸比角动量采用

j_{ m esc}=gamma J_o/M

。受体再发射通道的gamma取当前质量比；其他给定倍数是另一组物理假设。

结果中既比较不同保留比例，也比较不同流失杠杆臂。方括号可以变小、变大，甚至变成非正数。对于已经采用的正周期趋势，非正系数不再支持正的供体流失解。程序会保留“无正解”，而不是悄悄输出一个负流失率，或剪掉不符合预期的情景。

表1给出的单项质量误差不能替代两颗质量的联合协方差，所选周期趋势也没有配套的趋势不确定性。完整的联合后验仍需要这些不确定性与相关信息。若要估计现在的转移率，还需要新的观测、系统差异处理以及其他周期变化机制的检验；这里的数值属于这组历史定参推算。

对应阅读：Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1 · Mennickent & Djurašević (2013), arXiv1303.5812v1

回到原来的视向速度历式手算一次

把具名下载的1993年扫描PDF翻到印刷页139、式5。先在原页读清二次项，再做下面的手算：OCR容易把小指数和周期平方读错。这一步明确使用你从原页抄下的系数；四条公开命令没有重新拟合1532个历史平均视向速度。

T(E)=T_0+P_0E+cE^2,\quad P_0=12.913780\,\mathrm d,\quad c=3.87196\times10^{-6}\,\mathrm d

周期编号

E

是无量纲整数；历式时标是原文的HJD。局部周期为

P(E)=P_0+2cE

，因而在参考周期处

dot P=2c/P_0

，单位为天每天。

period_days = 12.913780
quadratic_days = 3.87196e-6
seconds_per_julian_year = 2 * quadratic_days / period_days * 86400 * 365.25
print(seconds_per_julian_year)

采用一年365.25天的换算，得到约每年18.924秒，与原文的历史圆整值一致。它与2013年采用的19秒不是两份互不相关的新测量。这组手算沿用原HJD历式，结果是参考周期附近的局部趋势；更换时标或向任意未来外推都需要独立处理。

对应阅读：Harmanec & Scholz (1993), A&A279,131–147

延伸到包层：重新积分原生结构

同一复算流程还会从两个精确原生ZIP成员解出质量、半径与氢丰度列，重新选择原生核切点并积分重力成本。结果记录实际切点、相邻区和不同氢阈值的差异，再用直接Roche体积计算能量条件终态。它没有从一个预存lambda常数倒推出“新积分”，也没有运行MESA来演化新的恒星。

独立验证另查缓慢风极限、不同通道的ODE收支、体积响应和保留的拟合符号反例。自有JS合成目标由Python用独立坐标和积分重算作对照；公开Python并不重新生成受保护的浏览器Engine。两者的证据边界应该一直保留。

最后把两种事情分开写进你的记录：渐台二部分用了真实论文的历史拟合参数；包层部分用了两颗学校模型的原生结构与给定能量假设。二者没有共同测量出渐台二的公共包层效率，也没有把一颗现在看到的双星直接倒推成唯一的历史终态。

对应阅读：Thoul & Chidester (2022), MESA Summer School; Zenodo 7082536 · Jared Goldberg (2025), MESA School Leuven Friday model · MESA 25.12.1, jdot_ml_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899 · Ivanova et al. (2013), Common envelope evolution; arXiv:1209.4302v2

练习与解释性答案

把19秒每年直接除以12.9，而忘记一天有86400秒，会怎样？
展开解释
得到的分数周期变化会大了86400倍，条件转移率也会跟着大同样倍数。先让分子与周期使用相同时间单位；“每年”保留为速率的分母。
某个流失情景的方括号系数为负，而采用的周期变化为正，应该怎样报告？
展开解释
在正供体流失的定义下，这组假设没有正解。保留失败情景，说明它与采用趋势不相容。不能把代数上的负数速率改名成正转移，也不能据此宣布某个真实流失机制已经被唯一排除。

动手检验预测

把定参推算与可玩的模型放在各自的位置

在交互场景中比较相同流失量下的供体风与受体再发射，注意轨道曲线可以不同。
用独立ZIP重算历史渐台二速率和流失情景，说明为什么不能从示踪点回放速度读出这个真实单位的速率。
再比较两个原生包层结构的核切点与能量终态。哪一部分重新读取了真实来源，哪一部分仍取决于你指定的模型假设？

独立固定源复算材料

材料包运行原文参数提取、原生重力积分、条件终态与独立检验。原始响应留在读者电脑上。1993年扫描历式另作人工阅读核对；公开Python不重拟合视向速度，也不重新生成受保护Engine。

下载完整独立复算ZIP

SHA-256: dc8fd19127787fb319ac04dcf2395f72260df7a523c07e663978a35cd4c46f51

逐项文件与精确大小

LICENSE.txt · 860 字节
README.md · 3,820 字节
binary-mass-transfer-science.zip · 377,227 字节
envelope_reference.py · 4,315 字节
fetch_sources.py · 2,502 字节
independent_volume.py · 2,412 字节
native-members.json · 1,768 字节
native-production-injections.json · 2,810,072 字节
observed_reference.py · 3,737 字节
orbit_reference.py · 3,384 字节
package-manifest.json · 3,133 字节
prepare_profiles.py · 1,747 字节
produce_figures.py · 4,391 字节
reproduce.py · 7,074 字节
requirements.txt · 46 字节
roche_reference.py · 3,449 字节
source-manifest.json · 2,849 字节
stability_reference.py · 6,838 字节
validate_limits.py · 4,058 字节
validate_science.py · 507 字节
verify_native_cases.py · 6,734 字节
verify_package.py · 550 字节

来源与计算范围

Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1
第2、3节、角动量流失通道、响应指数及附录A。使用质量比前，先对照论文中供体与受体的命名。
MESA 25.12.1, jdot_ml_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899
原生inlist与run_binary_extras.f90提供混合流失的解析参照。这里用独立ODE检查数学关系，没有编译或通过MESA演化测试。原生代码为LGPL3或更新版本，本Lab未复制它。
Ivanova et al. (2013), Common envelope evolution; arXiv:1209.4302v2
能量近似、能量来源与去向、时标和核/遗骸边界讨论。标量终态不预测流动过程或抛射形态。
Harmanec & Scholz (1993), A&A279,131–147
印刷页139、式5及视向速度历式讨论。二次项系数需在扫描原页辨清；这是历史视向速度分析，不是这里新拟合的食时。
Mennickent & Djurašević (2013), arXiv1303.5812v1
表1的拟合质量，以及第4节采用12.9天周期、每年增加19秒的计算。这里不复现盘的光变拟合、演化拟合或协方差。
Thoul & Chidester (2022), MESA Summer School; Zenodo 7082536
CC BY 4.0。r22.05.1的精确成员lab3_4M_solutions/LOGS_cheb_case1/profile1.data。这里的重力/核切点汇总及能量图是具名派生结果，不是新的MESA运行或总年龄测量。
Jared Goldberg (2025), MESA School Leuven Friday model
CC BY-SA 4.0；仓库提交ade5bca00e5caae97ac0b6f880f1de7ccead1553，成员work_minimini3/ALPHA_TEST_LOGS/LOGS_mlt2mesh1time1/profile8.data。结构派生汇总保留该许可；这份15倍太阳质量模型不同于前沿论文的12倍太阳质量供体。