--- title: "两颗恒星之间的窄口在哪里" description: "从两颗星的引力和旋转出发，找到内侧鞍点、闭合Roche瓣与体积半径，并弄清弹道能回答什么。" canonical_url: "https://sandboxphysics.com/zh/stories/binary-mass-transfer/the-passage-at-l1/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/binary-mass-transfer/the-passage-at-l1/" last_modified: "2026-11-03" --- [双星物质转移与公共包层](https://sandboxphysics.com/zh/stories/binary-mass-transfer/) AC29 · L1 # 两颗恒星之间的窄口在哪里 Sandbox Physics · 2026-11-03 窄口的位置由两颗星的质量和旋转共同决定；它到供体中心的距离，却不等于供体Roche瓣的半径。先把力与坐标摆清楚，就能看出这两个看似相近的长度为什么不同。 ## 开始前需要什么 会使用万有引力和圆周运动；理解向量分量。后半段用到等势面与体积积分，相关定义会在文中给出。 ## 跟着双星一起转 在惯性坐标系里，两颗星一直移动。为了研究某一时刻气体能否仍单独围住供体，换到与圆轨道一起转的坐标系：两颗点质量固定在横轴上，绕它们转的气体仍然可以运动。这个换法只是让背景简洁，并没有把气体固定住。 取供体质量为$M_d$、受体质量为$M_a$，两星中心间距为$a$。横轴从供体指向受体，纵轴沿轨道平面，垂直轴沿轨道角动量，构成右手系。原点放在共同质心，不能随手放在两星几何中点；只有等质量时两者才重合。 $$ q={M_d\over M_a},\quad \mu={M_a\over M_d+M_a}={1\over1+q},\quad x_d=-\mu,\quad x_a=1-\mu $$ 这里的坐标已经除以间距。供体在左，受体在右；质量比始终是供体除以受体。 再把总质量和引力常数也归一化。于是$G=M_d+M_a=a=1$，圆轨道角频率同样为一。若恢复物理量，要使用$Omega^2=G(M_d+M_a)/a^3$。这些单位描述当前的冻结轨道。质量移除量只有在另给转移率之后，才对应经过时间。 对应阅读：Eggleton (1983), Approximations to the radii of Roche lobes ## 先看受力，再写势 静止在共转画面里的小气团，受到两颗星的引力；旋转坐标还带来离心项。两颗质量越近，引力项变化越陡。离心项则取决于到旋转轴的距离：沿垂直轴走远，没有增加到旋转轴的距离，因此不能凭对称感再添一个垂直方向的离心项。 $$ \mathcal U={1-\mu\over r_d}+{\mu\over r_a}+{x^2+y^2\over2} $$ 定义$r_d^2=(x+mu)^2+y^2+z^2$、$r_a^2=(x-1+mu)^2+y^2+z^2$。这里采用正的引力项，$mathcal U$是通常Roche势的负值。 $$ \ddot x=\partial_x\mathcal U+2\dot y,\qquad\ddot y=\partial_y\mathcal U-2\dot x,\qquad\ddot z=\partial_z\mathcal U $$ 运动气团还受到科里奥利项。找到静止点时速度为零；真正积分流线时，就不能把这些速度项删掉。 符号看似只是约定，受力却不能跟着含混。这里势的梯度直接进入加速度；若改用负的物理有效势，加速度就要取它的负梯度。只把势改号、仍沿原来的梯度推粒子，会得到完全不同的路线。 对应阅读：Eggleton (1983), Approximations to the radii of Roche lobes ## 窄口是鞍点，瓣半径来自整个体积 $$ 0=x-{(1-\mu)(x+\mu)\over|x+\mu|^3}-{\mu(x-1+\mu)\over|x-1+\mu|^3} $$ 在$-mu