---
title: "差动旋转与轨道共振"
description: "从平坦圆周速率曲线推导角速度和径向小振动，再计算共转与内外Lindblad共振的圆轨道参考半径。"
canonical_url: "https://sandboxphysics.com/zh/stories/galactic-bar-orbits/differential-rotation-and-orbital-resonances/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/galactic-bar-orbits/differential-rotation-and-orbital-resonances/"
last_modified: "2026-11-05"
---

[星系棒、共振与恒星迁移](https://sandboxphysics.com/zh/stories/galactic-bar-orbits/)

星系棒轨道 · L1

# 差动旋转与轨道共振

[打开旋转棒Lab](https://sandboxphysics.com/zh/labs/galactic-bar-orbits/)

Sandbox Physics · 2026-11-05

恒星沿轨道走多快，与旋臂图案每秒转过多少角度，是两个不同的问题。比较两者的角度变化，能知道恒星多久遇到一次图案；再把相遇频率与径向摆动频率相比，就能找到最简单的轨道共振。

## 开始前需要什么

先读概念篇。需要圆周运动和导数；导数可以先理解成局部斜率，二阶导数则描述曲线的弯曲程度。

## 先定好方向与时间单位

在右手平面上取坐标 $x,y$，方位角 $\phi$ 逆时针增加。半径 $R$ 是到中心的距离。恒星沿圆周的切向速度记作 $v_\phi$，每秒转过的角度则是角速度 $\Omega=d\phi/dt=v_\phi/R$。角速度的单位是时间的倒数。

选定长度 $R_0$ 和速度 $v_0$，时间单位便是 $R_0/v_0$。数值计算把两个尺度都取为一。此时模式速度为一，指的是 $v_0/R_0$，并非每个时间单位转一整圈；一整圈有 $2\pi$ 弧度。

$$
\Omega(R)=\frac{v_0}{R},\qquad T_{\rm orb}=\frac{2\pi R}{v_0}
$$

这里采用圆轨道速率不随半径变化的近似。它给出轴对称背景的圆轨道参考，棒加进来以后，真实恒星还会有径向运动。

对应阅读：Dehnen (1999v1/2000), The Effect of the Outer Lindblad Resonance of the Galactic Bar on the Local Stellar Velocity Distribution · Sandbox Physics independent galactic-bar scientific reader v1

## 一条由固定成员组成的线会绕紧

让两颗星从同一个方位角出发，分别在 $R_0$ 与 $2R_0$ 处走圆轨道。经过 $2\pi R_0/v_0$，内圈走完一圈，外圈只走完半圈。起初连着它们的径向线，已经被拉出 $\pi$ 的角度差。

刚性旋转图案的规则不同：相同时间里，整幅图案各处转过相同角度。集体引力起伏可以带着这样的模式速度传播，组成它的恒星却有各自的绕行速度。短暂旋臂也可能边绕边变。先固定一个模式，能把两种运动清楚地分开。

对应阅读：Sellwood & Binney (2002), Radial mixing in galactic discs, astro-ph/0203510v2

## 在随图案转动的视角数相遇次数

让坐标系以模式速度 $\Omega_p$ 转动，恒星的新方位角就是 $\phi'=\phi-\Omega_pt$。近圆轨道上，它每秒改变 $\Omega-\Omega_p$。若图案有 $m$ 条等价臂，恒星相对走过 $2\pi/m$ 就又遇到一条臂。

做一个完整算例。取双重对称图案，$\Omega_p=1$。在 $R=0.8$，圆轨道角速度是 $1.25$，相对角速度为 $0.25$，相遇间隔便是 $4\pi\simeq12.57$。换到 $R=1.2$，相对角速度为 $-1/6$，图案从后面赶过恒星，相遇间隔是 $6\pi\simeq18.85$。

在共转处，圆轨道参考的相遇时间趋于无穷。受扰动的真实恒星可能在图案附近来回摆动，称为天平动。它是否真的被俘获，还要解轨道；相遇很慢本身不能代替这一步。

$$
T_{\rm encounter}=\frac{2\pi}{m|\Omega-\Omega_p|},\qquad \Omega(R_{\rm CR})=\Omega_p
$$

绝对值让间隔为正；相对角速度的正负则说明谁在赶超谁。

对应阅读：Dehnen (1999v1/2000), The Effect of the Outer Lindblad Resonance of the Galactic Bar on the Local Stellar Velocity Distribution

## 轻轻向外推一下，恒星还会径向摆动

这一步先拿掉棒。平坦圆轨道速率对应引力势 $\Phi_0=v_0^2\ln(R/R_0)$。以下能量都按单位质量计。轴对称引力保留单位质量角动量 $L_z$，径向运动便可以用有效势 $U=\Phi_0+L_z^2/(2R^2)$ 描述；第二项来自切向运动。

圆轨道半径记为 $R_c$。在这里，有效势一阶导数为零，且 $L_z=R_cv_0$。把恒星稍微移开 $\delta R$，最低阶回复力等于有效势的曲率乘位移，再加负号。对有效势求两次导数，代入圆轨道条件，得到 $U''(R_c)=2v_0^2/R_c^2$。

这个小振动叫本轮运动，频率 $\kappa$ 是圆轨道角速度的 $\sqrt2$ 倍。独立读者程序在 $R_c=1.3R_0$ 给恒星一个 $10^{-4}v_0$ 的径向速度，再从相邻两次径向最大值（远心点）量周期，得到 $\kappa\simeq1.08785658\,v_0/R_0$。

一般的本轮频率式也来自同一步。圆轨道的径向引力条件为 $\Phi'_0=R\Omega^2$，再求一次导数得到 $\Phi''_0=\Omega^2+R\,d\Omega^2/dR$。有效势中的角动量项，在圆轨道处另贡献 $3\Omega^2$，相加就是下面的表达式。

这里给的是角频率，一次径向小振动需要 $2\pi/\kappa$，不是 $1/\kappa$。频率平方为正，代表这一局部径向近似存在回复振动；有限棒势会把径向和方位运动耦合起来，因此圆轨道频率只是研究共振的起点。

$$
\frac{d^2\delta R}{dt^2}+\kappa^2\delta R=0,\qquad \kappa^2=R\frac{d\Omega^2}{dR}+4\Omega^2=2\frac{v_0^2}{R^2}
$$

最右边用了平坦圆周速率曲线。振动方程只保留位移的一阶项，因此适合近圆轨道。

对应阅读：Sellwood & Binney (2002), Radial mixing in galactic discs, astro-ph/0203510v2 · Sandbox Physics independent galactic-bar scientific reader v1

## 把相遇频率与径向频率放在一起

双重对称扰动沿圆轨道重复的频率为 $2|\Omega-\Omega_p|$。它若等于 $\kappa$，每次径向拉扯就可能遇到本轮运动的相同阶段。内外Lindblad共振的圆轨道参考位置，由下面两个带正负号的方程决定。

取 $R_0=v_0=\Omega_p=1$，得到 $R_{\rm ILR}=1-1/\sqrt2\simeq0.2929$、$R_{\rm CR}=1$、$R_{\rm OLR}=1+1/\sqrt2\simeq1.7071$。交互实验只接受整条轨道始终在 $R=0.5$ 以外的结果，因此内Lindblad参考圆不在这条计算范围内。

这些位置来自线性、近圆近似。有限强度的棒会改变频率，并形成有宽度的共振区域。参考圆能帮助定位，真实俘获的边界仍要从非线性轨道判断。

$$
\Omega_p=\Omega-\frac{\kappa}{2}\quad({\rm ILR}),\qquad \Omega_p=\Omega\quad({\rm CR}),\qquad \Omega_p=\Omega+\frac{\kappa}{2}\quad({\rm OLR})
$$

ILR与OLR分别代表内、外Lindblad共振。式中的两个频率都来自未扰动的近圆轨道。

对应阅读：Dehnen (1999v1/2000), The Effect of the Outer Lindblad Resonance of the Galactic Bar on the Local Stellar Velocity Distribution · Sellwood & Binney (2002), Radial mixing in galactic discs, astro-ph/0203510v2

## 练习与解释性答案

1. 保持平坦圆轨道速率，把 $\Omega_p$ 加倍，三个圆轨道参考半径会怎样变？

   展开解释

   它们都与模式速度成反比，因此各自减半。对这一种平坦圆周速率曲线，三个参考半径之间的比例保持不变。

2. 某颗星遇到双臂的频率，恰好等于径向摆动频率，足以证明轨道闭合吗？

   展开解释

   这说明它满足一个线性共振条件。闭合还要求在指定参考系里，完整的位置和速度同时回到起点。有限振幅和其他运动频率也会影响返回，仍需积分检验。

## 动手检验预测

### [先预测共转圆会移到哪里](https://sandboxphysics.com/zh/labs/galactic-bar-orbits/)

1. 改变模式速度前，按反比关系预测共转半径怎样变。
2. 拖动共转柄，比较新结果中的半径和预测。
3. 选中参考圆附近的一颗星，跟住真实径向运动，不把参考圆当成俘获边界。

[上一篇：恒星为什么会穿过星系旋臂](https://sandboxphysics.com/zh/stories/galactic-bar-orbits/why-stars-cross-spiral-arms/) [下一篇：旋转棒怎样交换能量与角动量](https://sandboxphysics.com/zh/stories/galactic-bar-orbits/torque-energy-and-a-rotating-bar/)

## 来源与计算范围

1. [Dehnen (1999v1/2000), The Effect of the Outer Lindblad Resonance of the Galactic Bar on the Local Stellar Velocity Distribution](https://arxiv.org/abs/astro-ph/9911161v1)

   固定十四页astro-ph/9911161v1，参数预飞已完整核读；对应式3、第3节与原轨道图。自有模型采用固定棒半径0.5和逆时针方向，原文历史银河系推断不当作当前测量。

2. [Sellwood & Binney (2002), Radial mixing in galactic discs, astro-ph/0203510v2](https://arxiv.org/abs/astro-ph/0203510v2)

   固定十三页v2，已完整核读。第2节单位质量能量与作用量交换，第3—4节短暂自引力波；原模拟与这里规定外棒势的示踪释放分别使用。

3. [Sandbox Physics independent galactic-bar scientific reader v1](https://sandboxphysics.com/stories/galactic-bar-orbits/materials/galactic-bar-orbits-reader-science-v1.zip)

   自有零网络Python积分与保存状态分析：完整二百五十六示踪星群、十个独立DOP853对照、圆轨道与本轮极限、规定增长功、六个新周期根和三张自有合成图。不含交互程序、完成轨道输入、原论文或观测数据。