恒星掠过旋转的引力图案

银河棒中的恒星轨道与共转

旋转的棒为什么能改变恒星的轨道？ 棒形引力场会沿恒星轨道的切向拉动它，改变角动量。恒星与棒相对转得多快，会影响这股拉力怎样累积。有些轨道能在随棒转动的画面中首尾闭合；规定引力场的计算可以展示这些运动。

拖动共转圈，改变棒的转速。选中一颗恒星，比较沿途的力矩与角动量变化。

所选恒星的角动量变化 $\Delta L_z/(R_0v_0)=-0.115$

所选恒星的径向速度 $v_R/v_0=0.0105$

近圆共转半径 $R_{\rm CR}/R_0=1$

带权方位响应 $Q_2=0.0855$

旋转的引力场

棒的转速

1.00

\Omega_{\rm p}R_0/v_0

势场强度

0.0400

A/v_0^2

回放时刻

12.6

t/(R_0/v_0)

拖青色共转柄或金色强度柄，选一颗恒星追踪它沿途受到的力矩。

正在计算完整的默认轨道。

参考系、单位与模型范围

显示切换只转动坐标。箭头沿当前显示参考系中的速度，长度经过压缩；颜色始终表示实际径向速度。

$R/R_0,\quad v/v_0,\quad t/(R_0/v_0)$

$R_{\rm b}=0.5R_0,\qquad t_{\rm end}=8\pi R_0/v_0$

所有粒子的权重相同。暖释放给成对粒子相反的径向速度，没有采用高斯分布。亮度只用于显示。

圆圈是背景场的近圆参考。在这组控制范围内，内Lindblad半径都落在排除的场区域中。场轴标记规定的引力图案，不表示测得的质量形状。

强度柄的长度是指示量，场半径保持固定。改变物理输入会从相同释放条件重算历史；部分暖释放与快速图案的组合会进入排除半径，因此被拒绝。

步长细化只是逐次计算的局部数值检查，没有证明连续输入范围内处处都有统一误差。最初的宽射击括区因两端同号而失败，记录仍保留；固定场参考采用另行登记的窄括区。

Dehnen 1999 v1：规定场与轨道族

Sellwood与Binney 2002 v2：瞬态波与迁移

模型假设

恒星是规定平面引力场中的无质量示踪粒子，引力场由对数背景与旋转四极项组成。恒星不反过来产生引力场，初始释放也没有拟合真实恒星分布或自洽星系。位置和速度采用参考单位，能量与角动量均按单位质量表示。冷释放中有位置、速度重合的粒子对；暖预设给成对粒子相反的初始径向速度。每段轨道中的场强与转速保持固定。改变物理参数会从相同初态重算整段历史，并比较两种积分步长。轨道进入排除的场半径区域时，请求会被拒绝，画面保留上一份完整结果；暖释放配上最快、最强的棒就是一个实际反例。圆圈标出背景场近圆轨道的共转与外Lindblad共振，偏心轨道不必以这些圆为精确边界。固定场中的闭合参考另用射击法求解，并向前积分完整周期。只有两条基准轨道有离线线性椭圆型证据，邻近解未判稳定，也没有非线性稳定结论。两组有限样本没有追踪整个分岔，不套用文献中的轨道族名称。粒子亮度和箭头长度用于显示，不代表光度或定量速度刻度。

先抓重点

核心问题: 恒星沿途受到的拉力，怎样累积成角动量变化？
一句话直觉: 力矩沿整条轨道积累。场以固定速度旋转时，能量与角动量各自可以变化，但它们的一个组合保持不变。

核心数学模型

规定的旋转引力场

\begin{gathered}\Phi=v_0^2\ln(R/R_0)+\Phi_{\rm b}\\\Phi_{\rm b}=A f(R)\cos[2(\phi-\Omega_{\rm p}t)]\\ f(R)=-\left(\frac{R_{\rm b}}R\right)^3,\ R>R_{\rm b}\end{gathered}

第一项给出平坦的圆周速度曲线，第二项有两个相反的较强拉力方向，并以固定速度转动。画面里的轨道始终留在固定场半径之外。

切向拉力改变角动量

\frac{dL_z}{dt}=-\frac{\partial\Phi}{\partial\phi},\qquad\Delta L_z=\int\tau_z\,dt

这里的角动量和下文的能量，都按单位质量计。指向中心的拉力对中心没有力矩。棒的那部分引力可以沿切向作用，正负号随着恒星相对棒的位置改变。时间图把实际角动量变化与沿途力矩积分放在一起比较。

共转比较的是相对转速

\begin{gathered}\Omega(R)=\frac{v_0}{R},\quad\kappa(R)=\frac{\sqrt2v_0}{R}\\ R_{\rm CR}=\frac{v_0}{\Omega_{\rm p}}\\ R_{\rm OLR}=\left(1+\frac1{\sqrt2}\right)R_{\rm CR}\end{gathered}

共转是近圆运动的恒星恰好跟上图案。外Lindblad参考半径把径向小振荡与相对转动联系起来。这些圆来自背景场的近圆近似；偏心轨道要看实际运动。

能量与角动量一起改变

E_J=E-\Omega_{\rm p}L_z,\qquad\Delta E=\Omega_{\rm p}\Delta L_z

只有同一个恒幅、匀速旋转场，才保持这个Jacobi组合。举例说，在参考单位中，图案转速为一，角动量增加0.05，能量就增加0.05。改变引力场后，要从原来的释放条件重新计算。

怎样量出方位响应

Q_2=\left|\sum_jw_j e^{2i(\phi_j-\Omega_{\rm p}t)}\right|,\qquad w_j=\frac1{256}

初始方位等间隔时，这个响应为零。粒子偏向两个相反方向，它就可能增大。每个粒子的权重相同；这个数描述方位分布，不是棒的质量或光度占比。

常见难点

看起来绕了一圈

容易误解图上的曲线回到原处，就是周期轨道。

正确理解周期父轨道可以作为附近运动的参照，它必须在棒参考系中同时回到原来的位置和速度。闭合参考预设先求这个条件，再向前积分整周期。内外两组是分别求得的局部样本。

在运动中改变棒

容易误解拖动棒的转速后，Jacobi量仍该接着原值保持不变。

正确理解这个守恒量只属于同一个稳态旋转场。改变物理参数后，程序从相同释放条件重算完整轨道。

动手实验

01
先去掉切向拉力
先暂停，把“势场强度”调到零，再回放保存的时刻。预想所选恒星的角动量变化会怎样。
你会看到： 冷释放中的近圆恒星几乎保持原有角动量。恢复棒以后，所选恒星受到时变的力矩，正负与净变化都取决于整条路径。
02
改变跟得上棒的位置
把青色共转柄向外拖，比较“共转移到外圈”和“更快的内圈图案”。
你会看到： 棒转得越慢，近圆共转半径越大。转速为0.7时，共转约在参考半径的1.43倍，外Lindblad参考约在2.44倍。恒星位置与方位响应曲线随新算出的运动一起改变。
03
比较两组真实周期解
先选“垂直棒的闭合参考”，再选“沿棒的闭合参考”，看舞台下方的径向起伏图。
你会看到： 同一个固定场中，内侧样本略向垂直棒的方向伸长，外侧样本略沿棒伸长。单独的起伏图能看清这点小差异。离线检验只给出了两条基准轨道的线性椭圆型证据，邻近解的稳定性仍未知。
04
读一颗恒星沿途受到的力矩
选一颗可见恒星，在某个保存时刻暂停并导出，比较角动量变化与力矩积分。
你会看到： 粒子、图和导出指向同一段完整轨道。力矩积分与角动量变化的紧密吻合，也来自离散踢移算法本身；轨道的物理准确性还要由独立积分参照检验。

银河棒中的恒星轨道与共转

旋转的拉力怎样改变角动量

核心数学模型

规定的旋转引力场

切向拉力改变角动量

共转比较的是相对转速

能量与角动量一起改变

怎样量出方位响应

常见难点

看起来绕了一圈

在运动中改变棒

动手实验

先去掉切向拉力

改变跟得上棒的位置

比较两组真实周期解

读一颗恒星沿途受到的力矩

相关实验与文章

星系棒、共振与恒星迁移