--- title: "从天琴座β星的周期增长算起" description: "下载并核对固定一手来源，从表1拟合质量和历史周期增长重算条件转移率，再检查流失假设与原生包层能量。" canonical_url: "https://sandboxphysics.com/zh/stories/binary-mass-transfer/beta-lyrae-period-growth/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/binary-mass-transfer/beta-lyrae-period-growth/" last_modified: "2026-11-03" --- [双星物质转移与公共包层](https://sandboxphysics.com/zh/stories/binary-mass-transfer/) AC29 · L3 # 从天琴座β星的周期增长算起 Sandbox Physics · 2026-11-03 历史定参下，天琴座β星每年约19秒的周期增长，对应保守圆轨道模型中约每年十万分之二太阳质量的供体流失。本教程从固定原站文本重新提取所需参数，算清这个有条件的数值；目标是复核已发表的推算。 ## 开始前需要什么 读过质量与角动量一章。独立复算需要Python、NumPy和SciPy；无需安装MESA，也不需要运行流体模拟。 ## 先确定今天要复算哪件事 天琴座β星，也叫渐台二，是一个有复杂吸积结构的相互作用双星。周期变化提供了很好的入手点，但这个“时钟”并不是单凭一张示意光变图得到的。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 ## 练习与解释性答案 1. 把19秒每年直接除以12.9，而忘记一天有86400秒，会怎样？ 展开解释 得到的分数周期变化会大了86400倍，条件转移率也会跟着大同样倍数。先让分子与周期使用相同时间单位；“每年”保留为速率的分母。 2. 某个流失情景的方括号系数为负，而采用的周期变化为正，应该怎样报告？ 展开解释 在正供体流失的定义下，这组假设没有正解。保留失败情景，说明它与采用趋势不相容。不能把代数上的负数速率改名成正转移，也不能据此宣布某个真实流失机制已经被唯一排除。 ## 动手检验预测 ### [把定参推算与可玩的模型放在各自的位置](https://sandboxphysics.com/zh/labs/lagrange-channels/?analysis=mass-transfer) 1. 在交互场景中比较相同流失量下的供体风与受体再发射，注意轨道曲线可以不同。 2. 用独立ZIP重算历史渐台二速率和流失情景，说明为什么不能从示踪点回放速度读出这个真实单位的速率。 3. 再比较两个原生包层结构的核切点与能量终态。哪一部分重新读取了真实来源，哪一部分仍取决于你指定的模型假设？ [上一篇：物质转移什么时候会越流越满](https://sandboxphysics.com/zh/stories/binary-mass-transfer/when-transfer-feeds-itself/) [下一篇：公共包层释放的能量去了哪里](https://sandboxphysics.com/zh/stories/binary-mass-transfer/where-common-envelope-energy-goes/) ## 独立固定源复算材料 材料包运行原文参数提取、原生重力积分、条件终态与独立检验。原始响应留在读者电脑上。1993年扫描历式另作人工阅读核对；公开Python不重拟合视向速度，也不重新生成受保护Engine。 [下载完整独立复算ZIP](https://sandboxphysics.com/stories/binary-mass-transfer/materials/binary-mass-transfer-science.zip) SHA-256: `dc8fd19127787fb319ac04dcf2395f72260df7a523c07e663978a35cd4c46f51` 逐项文件与精确大小 - [LICENSE.txt](https://sandboxphysics.com/stories/binary-mass-transfer/materials/LICENSE.txt) · 860 字节 - [README.md](https://sandboxphysics.com/stories/binary-mass-transfer/materials/README.md) · 3,820 字节 - [binary-mass-transfer-science.zip](https://sandboxphysics.com/stories/binary-mass-transfer/materials/binary-mass-transfer-science.zip) · 377,227 字节 - [envelope\_reference.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/envelope_reference.py) · 4,315 字节 - [fetch\_sources.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/fetch_sources.py) · 2,502 字节 - [independent\_volume.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/independent_volume.py) · 2,412 字节 - [native-members.json](https://sandboxphysics.com/stories/binary-mass-transfer/materials/native-members.json) · 1,768 字节 - [native-production-injections.json](https://sandboxphysics.com/stories/binary-mass-transfer/materials/native-production-injections.json) · 2,810,072 字节 - [observed\_reference.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/observed_reference.py) · 3,737 字节 - [orbit\_reference.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/orbit_reference.py) · 3,384 字节 - [package-manifest.json](https://sandboxphysics.com/stories/binary-mass-transfer/materials/package-manifest.json) · 3,133 字节 - [prepare\_profiles.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/prepare_profiles.py) · 1,747 字节 - [produce\_figures.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/produce_figures.py) · 4,391 字节 - [reproduce.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/reproduce.py) · 7,074 字节 - [requirements.txt](https://sandboxphysics.com/stories/binary-mass-transfer/materials/requirements.txt) · 46 字节 - [roche\_reference.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/roche_reference.py) · 3,449 字节 - [source-manifest.json](https://sandboxphysics.com/stories/binary-mass-transfer/materials/source-manifest.json) · 2,849 字节 - [stability\_reference.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/stability_reference.py) · 6,838 字节 - [validate\_limits.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/validate_limits.py) · 4,058 字节 - [validate\_science.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/validate_science.py) · 507 字节 - [verify\_native\_cases.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/verify_native_cases.py) · 6,734 字节 - [verify\_package.py](https://sandboxphysics.com/stories/binary-mass-transfer/materials/verify_package.py) · 550 字节 ## 来源与计算范围 1. [Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1](https://arxiv.org/abs/astro-ph/9703016v1) 第2、3节、角动量流失通道、响应指数及附录A。使用质量比前，先对照论文中供体与受体的命名。 2. [MESA 25.12.1, jdot\_ml\_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899](https://docs.mesastar.org/en/25.12.1/test_suite/jdot_ml_check.html) 原生inlist与run\_binary\_extras.f90提供混合流失的解析参照。这里用独立ODE检查数学关系，没有编译或通过MESA演化测试。原生代码为LGPL3或更新版本，本Lab未复制它。 3. [Ivanova et al. (2013), Common envelope evolution; arXiv:1209.4302v2](https://arxiv.org/abs/1209.4302v2) 能量近似、能量来源与去向、时标和核/遗骸边界讨论。标量终态不预测流动过程或抛射形态。 4. [Harmanec & Scholz (1993), A&A279,131–147](https://articles.adsabs.harvard.edu/pdf/1993A%26A...279..131H) 印刷页139、式5及视向速度历式讨论。二次项系数需在扫描原页辨清；这是历史视向速度分析，不是这里新拟合的食时。 5. [Mennickent & Djurašević (2013), arXiv1303.5812v1](https://arxiv.org/abs/1303.5812v1) 表1的拟合质量，以及第4节采用12.9天周期、每年增加19秒的计算。这里不复现盘的光变拟合、演化拟合或协方差。 6. [Thoul & Chidester (2022), MESA Summer School; Zenodo 7082536](https://zenodo.org/records/7082536) CC BY 4.0。r22.05.1的精确成员lab3\_4M\_solutions/LOGS\_cheb\_case1/profile1.data。这里的重力/核切点汇总及能量图是具名派生结果，不是新的MESA运行或总年龄测量。 7. [Jared Goldberg (2025), MESA School Leuven Friday model](https://vincentvanlaer.github.io/mesa-school-labs/friday/morning-session/) CC BY-SA 4.0；仓库提交ade5bca00e5caae97ac0b6f880f1de7ccead1553，成员work\_minimini3/ALPHA\_TEST\_LOGS/LOGS\_mlt2mesh1time1/profile8.data。结构派生汇总保留该许可；这份15倍太阳质量模型不同于前沿论文的12倍太阳质量供体。