跳到主要内容
Sandbox Physics
双星物质转移与公共包层

AC29 · L2

轻星送出物质,轨道为什么反而张开

轻的供体把物质送给重的伴星,轨道可以张开;较重的供体开始送物质时,轨道却可以收紧。关键不在于“总质量少了多少”这一句话,而在于质量怎样重新分配,以及什么物质带走了多少角动量。

把两颗星携带的角动量加起来

先做最干净的圆轨道模型:两颗星用点质量代表,暂不计算自转、盘中储存的角动量和其他外力矩。每颗星到共同质心的距离不同,轻星离得更远。把两颗星的轨道角动量相加,得到一个同时含质量分配与间距的量。

M=Md+Ma,Jo=MdMaGaMM=M_d+M_a,\qquad J_o=M_dM_a\sqrt{Ga\over M}
总质量记为MM,间距为aa。供体与受体的下标沿用上一章。
dln⁡a=2 dln⁡Jo−2 dln⁡Md−2 dln⁡Ma+dln⁡Md\ln a=2\,d\ln J_o-2\,d\ln M_d-2\,d\ln M_a+d\ln M
这是同一个圆轨道关系的对数微分。质量和角动量的变化给定后,间距响应就随之确定;不用再凭“引力更强”来猜。

这里已经包含重要的物理选择。如果盘或自转暂时存下了显著角动量,轨道角动量就不能擅自视为不变。相反,若声明这些储存与交换都可忽略,下面的保守计算就有清楚的用途。

对应阅读:Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1

保守转移的最紧轨道在等质量处

令x>0x>0表示从供体移走、又全部交给受体的质量。于是Md=Md0−xM_d=M_{d0}-x、Ma=Ma0+xM_a=M_{a0}+x,总质量与轨道角动量都保持不变。积分上式,可以直接比较前后间距。

aa0=(Md0Ma0MdMa)2,dln⁡adx=2Md−2Ma{a\over a_0}=\left({M_{d0}M_{a0}\over M_dM_a}\right)^2,\qquad {d\ln a\over dx}={2\over M_d}-{2\over M_a}
供体较重时导数为负;供体较轻时导数为正。质量相等的地方是转折,不是动画中任意选定的一帧。

例如最初两颗星分别为3Modot3M_odot和2Modot2M_odot。移走0.5Modot0.5M_odot以后,两颗星都为2.5Modot2.5M_odot,质量乘积最大,因此轨道最紧。

amin⁡a0=(3×22.5×2.5)2=0.9216,Pmin⁡P0=0.92163/2=0.884736{a_{\min}\over a_0}=\left({3\times2\over2.5\times2.5}\right)^2=0.9216,\qquad {P_{\min}\over P_0}=0.9216^{3/2}=0.884736
周期关系来自开普勒定律,并利用了总质量不变。这个解析例子移走了初始供体约六分之一的质量,超出交互场景所展示的局部8%区间;它不是完整恒星演化轨迹。
自算最初供体与受体质量比0.4、1.5和3时,保守间距对移除比例的变化
保守圆轨道,质量与轨道角动量留在双星内。每条曲线直接计算质量账本;横轴是指定的移除量,不是恒星演化时间。

图中横轴一律是移除量相对初始供体质量的比例。不同最初质量比给出不同趋势,即使每一条都没有物质逃出体系。轻星送给重星时,质量乘积变小,维持相同轨道角动量就要求间距增大。

对应阅读:Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1

给流失物质也留一个位置

若受体只留住比例bb,剩下的物质离开双星,质量变化应写成三笔:供体减少xx,受体增加bxbx,体系外增加(1−b)x(1-b)x。最后一笔虽然不再属于双星,却不能从收支中消失。

dJo=−jesc(1−b) dx,jesc=γJoMdJ_o=-j_{\rm esc}(1-b)\,dx,\qquad j_{\rm esc}=\gamma{J_o\over M}
规定逃逸物质的比角动量。只有这个特定定义下,常数gamma才表示总轨道角动量除以总质量的倍数。
JoJo0=(MM0)γ,aa0=(MM0)2γ+1(Md0Ma0MdMa)2{J_o\over J_{o0}}=\left({M\over M_0}\right)^\gamma,\qquad {a\over a_0}=\left({M\over M_0}\right)^{2\gamma+1}\left({M_{d0}M_{a0}\over M_dM_a}\right)^2
假设bb和这个定义下的gammagamma都不随移除量改变。轨道与外逃物质共同携带最初的角动量。

相同的失质量比例,也可能给出不同间距:带走较多比角动量的流失,会更强地拉低剩余轨道角动量。它与质量重新分配造成的张开或收缩竞争。决定轨道还需要这部分物质携带的比角动量。

对应阅读:Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1

从哪里逃走,会留下不同的杠杆臂

各向同性流失仍会携带发射星本来具有的轨道运动。供体风与“气体先到受体附近、再各向同性发射”因此不同。用到质心的杠杆臂计算,二者在同一比角动量定义下分别给出:

γdonor=MaMd=1q,γreemission=MdMa=q\gamma_{\rm donor}={M_a\over M_d}={1\over q},\qquad\gamma_{\rm reemission}={M_d\over M_a}=q
这些gamma会随质量比改变,不能拿最初数值当成整个过程的常数。

完全不捕获的缓慢供体风有b=0b=0。把供体杠杆臂代回微分关系,得到dlna=−dlnMdln a=-dln M,也就是下面的极限。

aM=a0M0aM=a_0M_0
这是缓慢、各向同性、可维持近圆轨道的风极限。瞬间爆发的失质量可能改变偏心率,不能沿用这条缓慢轨道响应。

MESA的固定jdot_ml_check原生算例还把供体风、受体再发射和环流失混合起来。阅读它时,要连同各自的比角动量定义一起抄写;环半径或环的杠杆臂参数,并不自动等于这里的gamma。这里用独立ODE比较该原生算例的数学关系。恒星结构求解与整套MESA演化属于另外的计算层次。

对应阅读:Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1 · MESA 25.12.1, jdot_ml_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899

质量和角动量守恒,不等于轨道能守恒

Eorb=−GMdMa2aE_{\rm orb}=-{GM_dM_a\over2a}
把前面得到的质量与间距代入,可以记录轨道能变化。它是诊断量,不是包含恒星、气体、热和辐射的总能量。

即使采用保守转移,质量与轨道角动量都没有流到外面,轨道能也一般会变。真实体系中的自转、吸积、热和辐射提供了需要进一步处理的能量交换。本章只把指定的质量与角动量收支变成圆轨道响应,不借此声称气体能量已经解决。

同样,横轴“移除了多少质量”还没有时钟。要把它变成每年的转移率,需要额外观测或一个转移模型。后面的渐台二教程会用历史周期增长来做一个有条件的速率推算,而不是从流线回放快慢读出答案。

对应阅读:Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1 · Ivanova et al. (2013), Common envelope evolution; arXiv:1209.4302v2

练习与解释性答案

  1. 最初为3Modot3M_odot和2Modot2M_odot的保守圆轨道,在移走多少质量时最紧?最小间距比是多少?

    展开解释

    质量和固定时,质量乘积在相等处最大。令3−x=2+x3-x=2+x,得x=0.5Modotx=0.5M_odot;代入间距式得amin/a0=0.9216a_{min}/a_0=0.9216。这不决定转移是否能实际持续到那里。

  2. 受体留住一半,供体失去0.02Modot0.02M_odot。三笔质量应怎样记?

    展开解释

    供体变化−0.02Modot-0.02M_odot,受体变化+0.01Modot+0.01M_odot,外逃物质增加0.01Modot0.01M_odot。三者相加为零;轨道角动量还要另记外逃物质带走的部分。

  3. 为什么不能把缓慢风的间距关系直接用于瞬间失去同样质量的事件?

    展开解释

    缓慢关系假设轨道能持续调整、接近圆形。瞬时失质量时,位置和速度来不及按新圆轨道关系调整,结果可能具有偏心率,甚至不再束缚。相同的总失质量不保证相同的轨道过程。

动手检验预测

让相同的流失量走不同通道

  1. 把最初质量比设为三,移除比例设为8%,在展开区把受体保留比例设为零。
  2. 先选供体风,再选给定比角动量,并把其倍数设为二。比较间距曲线:为什么同样的流失质量可以对应张开与收缩?
  3. 导出两份记录,分别检查质量变化和外逃角动量,再解释轨道能如何进入更大的恒星与气体能量问题。

来源与计算范围

  1. Soberman, Phinney & van den Heuvel (1997), astro-ph/9703016v1

    第2、3节、角动量流失通道、响应指数及附录A。使用质量比前,先对照论文中供体与受体的命名。

  2. MESA 25.12.1, jdot_ml_check; commit 4a68bfbce0f6930c3f861d0611e914f02569a899

    原生inlist与run_binary_extras.f90提供混合流失的解析参照。这里用独立ODE检查数学关系,没有编译或通过MESA演化测试。原生代码为LGPL3或更新版本,本Lab未复制它。

  3. Ivanova et al. (2013), Common envelope evolution; arXiv:1209.4302v2

    能量近似、能量来源与去向、时标和核/遗骸边界讨论。标量终态不预测流动过程或抛射形态。

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

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