双星演化与质量转移

概述#

银河系中约有 $40\%$–$70\%$ 的类太阳恒星处于双星或多星系统中（Raghavan et al. 2010，ApJS，190，1）；对于大质量 O 型星，这一比例超过 $70\%$（Sana et al. 2012，Science，337，444）。因此，大多数恒星并非孤立演化，而是与伴星共享引力束缚系统，它们的生命历程深刻地受到相互作用的影响。双星相互作用不仅创造了宇宙中一些最壮观的天象（Ia 型超新星、X 射线双星、毫秒脉冲星、引力波事件），也是宇宙化学演化的重要驱动力——Ia 型超新星是银河系铁元素的主要来源，而中子星并合提供了 r 过程元素。理解双星演化，是现代天体物理学的核心课题之一。

把恒星当成孤立球体来演化，只对大约三到六成的类太阳星近似成立，对 O 型星则几乎一开始就不成立。大质量恒星演化的外包层膨胀、质量损失和最终塌缩，都会被一颗足够近的伴星改写。恒星演化教科书里的单星轨迹，因此是极限情形，不是默认情形。Ia 型爆炸、X 射线双星和双致密星并合，都是这条相互作用链上的不同出口。

洛希瓣几何：双星系统的引力边界#

描述双星相互作用的核心概念是洛希瓣（Roche lobe）——每颗恒星的引力势占主导的空间区域。在同步旋转双星系统中，等势面形成两个相接的泪滴形区域，两颗星各占一个。两个泪滴的接触点称为拉格朗日点 L1，物质可以通过 L1 从一颗星流向另一颗星。洛希瓣半径的经典近似（Eggleton 1983）：

$$\frac{R_L}{a} = \frac{0.49 q^{2/3}}{0.6 q^{2/3} + \ln(1 + q^{1/3})}$$

其中 $a$ 是轨道分离，$q = M_1/M_2$ 是质量比。当一颗恒星演化膨胀（如进入巨星阶段）并充满洛希瓣，物质开始通过 L1 流向伴星，称为洛希瓣溢流（Roche lobe overflow, RLOF）。

洛希瓣不是一堵实体墙。它是共转参考系里有效势的分界面：引力加上离心势，在内拉格朗日点 L1 收成一个鞍点。物质一旦越过 L1，就进入伴星的势阱，通常先绕成盘，而不是沿径向砸下去。Eggleton 的公式把这滴泪的等效半径写成质量比的缓变函数，便于在演化程序里每一步判断“这颗星有没有满出来”。

溢流可以由两件事触发。施主星自己膨胀，主序结束后进入巨星或 AGB，半径涨过 $R_L$；或者轨道因磁制动、潮汐或引力波辐射而收缩，洛希瓣反过来去撞恒星表面。两种触发在观测上都可以表现为正在转移的密近双星，但后续稳不稳定完全不同。一旦自转与轨道不同步，或轨道明显偏心，这套等势面几何本身就要改写，不能再把瞬时的泪滴轮廓当成严格边界。

质量转移类型#

根据驱动质量转移的物理机制，RLOF 分为三种情况（Case A/B/C）：

质量转移的稳定性取决于两颗星的质量比和施主星的状态。如果质量转移速率过快，受主星无法处理流入的物质，系统可能进入共包层演化（Common Envelope Evolution, CEE）。

Case A 发生在两星还都很紧、施主仍在主序烧氢的时候，转移可以拖很久，轨道和两星结构一起被改。Case B 最常见：氢壳层点亮、恒星迅速胀成亚巨星或红巨星，半径在热时标上猛涨，洛希瓣来不及跟着变。Case C 更晚，氦核燃烧之后的 AGB 外包层又一次涨满，转移的是更疏松、更易被剥掉的包层。

稳定与否不看转移“看起来有多猛”，而看质量比和施主的半径—质量响应。施主更重时，丢掉质量往往让它相对洛希瓣更满，形成正反馈，轨道也因质量比反转前的收缩而更紧。受主若来不及热调整，物质在它周围堆积，系统就滑向共包层。同一质量比，在主序和红巨星阶段可以走出完全相反的结局。

共包层演化：轨道急剧收缩#

当质量转移变得不稳定（通常发生在大质量比 $q \gg 1$ 的系统或巨星的热时标上的快速 RLOF），施主星的整个外包层可能迅速膨胀并包裹伴星，形成共包层。在共包层中：

1. 伴星（和施主星的核心）在致密包层中旋转
2. 动力学摩擦将轨道动能转移给包层
3. 包层被逐渐加热并可能被驱逐出去
4. 如果包层成功驱逐，留下一个轨道大幅收缩的紧密双星（致密残骸+伴星）

能量守恒近似（$\alpha_{CE}$ 参数化）：

$$\alpha_{CE} \cdot \left(\frac{G M_{d,c} M_a}{2 a_f} - \frac{G M_d M_a}{2 a_i}\right) = E_{bind}$$

其中 $\alpha_{CE}$ 是包层驱逐效率（$0 < \alpha_{CE} \leq 1$），$a_i$ 和 $a_f$ 是初始和最终轨道分离。共包层演化是产生紧密双致密星系统（如双白矮星、白矮星+中子星）的关键机制，是引力波波源（如 LIGO 探测到的双黑洞、双中子星）形成的必要步骤（Ivanova et al. 2013，A&A Review，21，59）。共包层演化的物理机制至今尚未完全理解，$\alpha_{CE}$ 参数的不确定性是双致密星种群合成的最大误差来源。

共包层要解决的几何问题很具体。两颗大质量星在主序上可以相距很远，演化成巨星后半径会涨到把伴星吞进去；若不把轨道在包层里耗掉一大截，后来的两颗致密星根本进不了LIGO的频率窗口。能量形式把“轨道收缩释放的结合能”拿去付包层的束缚能，付得起就留下近距离双星，付不起就并合成一颗星。$\alpha_{CE}$ 把所有尚未算清的辐射、喷流和复合效率收成一个 0 到 1 之间的数。

这个数至今不是从三维流体力学唯一标定出来的。种群合成只能靠它扫描，再拿超新星率、X 射线双星和引力波并合率来反推。Ivanova 等人的综述把现状说成：方向清楚，微观机制仍缺。Webbink 较早把双白矮星并合同能量收支连在一起，既是 Ia 型前身星讨论的源头之一，也是共包层形式化的经典文献（Webbink 1984，ApJ，277，355）。

Ia 型超新星：白矮星的两种死法#

Ia 型超新星是宇宙距离测量的标准烛光，也是银河系中铁族元素的主要来源。它们被认为来自白矮星在双星系统中积累质量至钱德拉塞卡极限附近。但"给白矮星喂质量"的机制存在两种主要竞争理论：

单简并模型（Single Degenerate, SD）#

白矮星从正常伴星（主序星、巨星）通过 RLOF 积累氢/氦，在表面燃烧后逐渐增加碳氧核质量，最终在质量接近 $M_{Ch} \approx 1.4\,M_\odot$ 时发生热核爆炸。

• 优点：概念清晰，有直接物理机制
• 问题：观测上从未在 Ia 型超新星遗迹中找到预期的幸存伴星；SD 模型预测的新星活动也难以解释 Ia 型超新星的实际发生率

单简并把洛希瓣溢流写成一条可计算的喂食管道，却在遗迹里找不到那颗该留下的伴星。新星把吸积层反复烧掉，净增长效率可以低到涨不到 $M_{Ch} \approx 1.4\,M_\odot$。这不是说 SD 不可能，而是说它最干净的版本已经被观测逼到角落里。

双简并模型（Double Degenerate, DD）#

两颗白矮星在引力波辐射驱动下轨道衰减并最终并合，如果总质量超过钱德拉塞卡极限，则引发爆炸。

• 优点：双白矮星系统已有直接观测；与某些超亮和亚亮 Ia 型超新星的观测特征吻合
• 问题：质量比的并合可能产生轻度吸积诱导坍缩（AIC）形成中子星，而不是热核爆炸

双简并把共包层和引力波耗散串成一条不必留下普通伴星的路。它与某些超亮、亚亮事件对得上，却过不了“并合一定爆炸”这一关。质量比不合适时，并合产物可以坍缩成中子星，观测上就从 Ia 型样本里消失。因此 DD 也不是闭合答案。

目前学界共识是 Ia 型超新星可能来自多种前体系统（"多前体通道"，Maoz, Mannucci & Nelemans 2014，ARA&A，52，107），不同环境和时标对应不同渠道。

SD 要靠洛希瓣溢流把氢或氦稳定喂到白矮星上，速率太快会爆发新星把物质吐回去，太慢又涨不到钱德拉塞卡极限。DD 要先经过至少一次共包层，把两颗白矮星的轨道收到引力波能在宇宙年龄内并合的程度。

两条路都还没被观测关死，也没有哪一条能单独解释全部速率和化学。测距上的含义写在标准烛光与Ia超新星测距：通道混合若随红移变，标准化后的峰值光度就可能带上演化残差。爆炸物理本身见Ia 型超新星。

X 射线双星：活跃的质量转移系统#

X 射线双星（X-ray Binary, XRB）是质量转移正在进行的双星系统，其中致密天体（中子星或黑洞）从伴星吸积物质，物质落入的引力势能转化为 X 射线辐射：

$$L_X \sim \frac{G M \dot{M}}{R} \approx 10^{38} \left(\frac{\dot{M}}{10^{-8}\,M_\odot\,\text{yr}^{-1}}\right) \text{ erg/s}$$

X 射线双星分为两类：

• 高质量 X 射线双星（HMXB）：伴星为 OB 型大质量星，通过恒星风吸积（如 Cyg X-1）
• 低质量 X 射线双星（LMXB）：伴星为低质量星，通过 RLOF 吸积盘吸积（如半人马座 X-4）

分类的判据是伴星质量和转移几何，不是 X 射线有多亮。天鹅 X-1 可以在尚未充满洛希瓣时靠恒星风喂黑洞；半人马座 X-4 必须靠溢流维持盘。风吸积给出持续而不稳定的供应，溢流吸积给出盘和态变。把两类源的光谱拿来直接平均，会把两种几何揉成一种不存在的“典型双星”。

LMXB 中的中子星可以被吸积物质自旋加速到毫秒级旋转周期——这就是毫秒脉冲星（MSP）的形成机制，有时称为"再生脉冲星"（recycled pulsar）。毫秒脉冲星的自转极为稳定，可用作高精度时钟，在脉冲星计时阵列中探测引力波背景（PTA 计划）。

正在转移的系统把洛希瓣几何变成可拍摄的现场。高质量系统里，OB 伴星的风先于溢流把质量送走，致密星在风中开出弓形激波；低质量系统则几乎总是 RLOF，物质先经过吸积盘，引力势能在盘内边界转化成 X 射线。公式里的 $10^{38}\,\mathrm{erg/s}$ 对应约 $10^{-8}\,M_\odot\,\mathrm{yr}^{-1}$ 的吸积率，数量级与中子星爱丁顿光度相当，不是另一套任意标定。

HMXB 跟着年轻星协走，伴星寿命短，转移窗口也短；LMXB 可以在长得多的时标上慢慢吸，把中子星自旋到毫秒。再生之后的脉冲星离开吸积盘，成为计时阵列的钟。X 射线双星是这条链还在发光的阶段，脉冲星计时阵列是它熄灭之后留下的时钟网络。天鹅 X-1 与半人马座 X-4 分别站在风吸积和盘吸积两端，不能互相当成对方的近景。

双星对引力波天文学的贡献#

LIGO/Virgo/KAGRA 探测到的引力波事件几乎全部来自双致密星系统的并合：

• BBH（双黑洞并合）：来自大质量双星演化，经历共包层演化或类共包层演化（化学均一演化）
• BNS（双中子星并合）：GW170817 + AT2017gfo 千新星——中子星并合产生 r 过程元素的直接证据
• NSBH（中子星+黑洞并合）：2021 年首次探测（GW200105、GW200115）

这三类并合对应三条完全不同的双星生命史。双黑洞可以走共包层，也可以走化学均一演化，伴星在主序上就被潮汐锁成均匀自转，巨星膨胀被压住，轨道不必靠包层摩擦来收。双中子星必须连续两次超新星都不把系统踢散，成功率更低，所以事件更稀。中子星–黑洞要在质量比极端的情况下仍保住较轻的那颗，直到引力波把轨道收进地面干涉仪的频带。

截至 2023 年，LIGO/Virgo/KAGRA 第三次观测运行（O3）共检测到约 90 个引力波候选事件（GWTC-3），为双星演化理论提供了前所未有的约束（Abbott et al. 2023，Physical Review X，13，041039）。

地面干涉仪听见的，几乎都是双星演化已经走完之后的最后几圈。大质量双星要变成可并合的双黑洞，必须在巨星阶段用共包层或化学均一演化把轨道收紧，否则两颗星会在包层膨胀时提前并合，留不下两颗恒星级黑洞。双中子星更苛刻：两颗星都要先经历超新星而不把系统踢散，再在共包层后留下足够紧的轨道。

GW170817 把这条链从引力波波形接到光学千新星 AT2017gfo。双中子星并合与千新星给出 r 过程元素的直接现场，千新星与 r 过程把光谱里的镧系特征写进化学演化。O3 的约 90 个候选事件把并合率从数量级估计推进到可以反推 $\alpha_{CE}$ 和踢速度的程度，误差仍然被共包层参数主导。双中子星若测到可靠的电磁波对应体，还可以当作标准汽笛，与 Ia 型标准烛光走一条完全不同的距离路径。

跨域连接#

• 参考系与惯性：洛希瓣只在共转参考系里才是等势面——离心势必须与引力势合并，内拉格朗日点是这个有效势的鞍点。推论：溢流是一条穿过鞍点的窄流而非球对称的风；一旦自转与轨道不同步或轨道偏心，这套几何本身就失效。
• 控制论：转移是否稳定取决于扰动后的反馈符号：若失去质量反而使施主星相对洛希瓣更满，就是正反馈，转移失控并走向共包层。判据是两条半径—质量指数之差，而不是转移速率本身有多大——速率快只是失稳的结果，不是失稳的原因，这也是同样的质量比在不同演化阶段结局相反的缘故。
• 蒙特卡洛方法：共包层效率之类参数无法从第一原理算出，只能靠种群合成：抽样大量初始双星、按参数化规则演化、再与观测率比对。推论：结论对参数先验极敏感，这也是双致密星并合率预言里最大的误差来源。
• 光谱学：被剥离的施主星表面暴露的是碳氮氧循环加工过的物质，氮增丰而碳、氧亏损。这条丰度指纹能在完全看不到伴星的情况下证明发生过质量转移，因为单星演化不会把这层物质翻上表面。
• 标准烛光与 Ia 超新星测距：若 Ia 型超新星来自多种前体通道，其混合比例会随宿主星系的年龄与金属丰度变化，标准化后的峰值光度就可能带上随红移的系统漂移。双星演化的不确定性正是沿这条路径进入宇宙学参数的。

参考文献#

• Raghavan, D. et al. (2010). A Survey of Stellar Families: Multiplicity of Solar-type Stars. The Astrophysical Journal Supplement Series, 190, 1–42.
• Sana, H. et al. (2012). Binary Interaction Dominates the Evolution of Massive Stars. Science, 337, 444–446.
• Eggleton, P.P. (1983). Approximations to the radii of Roche lobes. The Astrophysical Journal, 268, 368.
• Ivanova, N. et al. (2013). Common Envelope Evolution: Where we stand and how we can move forward. Astronomy & Astrophysics Review, 21, 59.
• Maoz, D., Mannucci, F., & Nelemans, G. (2014). Observational Clues to the Progenitors of Type Ia Supernovae. Annual Review of Astronomy and Astrophysics, 52, 107–170.
• Abbott, R. et al. (LIGO/Virgo/KAGRA Collaboration) (2023). GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo during the Second Part of the Third Observing Run. Physical Review X, 13, 041039.
• Webbink, R. F. (1984). Double white dwarfs as progenitors of Type I supernovae and of R Coronae Borealis stars. The Astrophysical Journal, 277, 355–360.

延伸阅读#

• de Mink, S.E. et al. (2014). The Rotation Rates of Massive Stars: The Role of Binary Interaction. The Astrophysical Journal, 782, 7. — 双星相互作用对大质量恒星自转的影响
• Postnov, K.A. & Yungelson, L.R. (2014). The Evolution of Compact Binary Star Systems. Living Reviews in Relativity, 17, 3. — 双致密星演化的权威综述
• Tauris, T.M. & van den Heuvel, E.P.J. (2023). Physics of Binary Star Evolution. Princeton University Press. — 双星演化领域最新权威教材