---
title: "周期轨道与恒星径向迁移"
description: "在固定棒势中实际射击周期母轨道，区分轨道形状与线性稳定，再理解短暂旋臂怎样留下径向迁移。"
canonical_url: "https://sandboxphysics.com/zh/stories/galactic-bar-orbits/periodic-orbits-and-radial-migration/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/galactic-bar-orbits/periodic-orbits-and-radial-migration/"
last_modified: "2026-11-05"
---

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

星系棒轨道 · L3

# 周期轨道与恒星径向迁移

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

Sandbox Physics · 2026-11-05

周期轨道要求恒星在指定参考系中，位置和速度一起回到起点。径向迁移关注的则是角动量或导引中心半径怎样变化。两者会在共振附近相遇，但仅凭一条看起来闭合的曲线，还不知道恒星是否被俘获，也不知道半径变化能否长久保留。

## 开始前需要什么

先读频率与力矩两章。现在会用到数值积分、求一个函数的零点，以及用微小位移研究轨道响应的线性近似。

## 释放一批星，与解出周期轨道是两种问题

交互实验从八个环释放示踪星，每环十六个方位，每处各放正负两种径向速度。每个样本权重都是二百五十六分之一。冷释放时，正负两种速度相同，因此实际有一百二十八个不同位置。这便于比较规定引力下的响应，起点没有被构造成平衡分布函数。

我们对四个惯性坐标 $(x,y,v_x,v_y)$ 积分 $8\pi$ 个时间单位。每一步先改半步速度，再推进位置，最后补上半步速度，这就是KDK跃蛙法。把步长减半后比较整段历史，可以发现只看末态能量遗漏的轨迹偏差。

细网格保存二百五十七帧。中细差除以三，是二阶方法进入收敛区以后的误差估计；另一种高阶自适应积分法DOP853和解析极限，提供另外的比较。每次真实漂移都检查半径，内侧穿越发生在两帧之间也会被发现。

对应阅读：Sandbox Physics independent galactic-bar scientific reader v1

## 借对称性，把返回问题缩成一个零点

固定 $A=0.01$、$\Omega_p=1$、$R_b=0.5$。在随棒坐标中，从 $(x,y,u,v)=(a,0,0,w)$ 出发；$u,v$ 是随棒坐标的变化率。这一组位于共转外侧的星，初始 $w$ 为负。沿轨道等到第一次从下向上穿过负 $x$ 侧的 $y=0$，检查此时 $u$ 是否为零。

返回速度随 $w$ 变化。先给出一段两端异号的区间，再求返回速度为零的初始值。利用反射对称，完整周期应是这段返回时间的两倍；随后仍要真正积完一整圈，检查四个状态量是否都回来。

以 $a=1.45$ 为例，新解得到 $w=-0.4483214884$，完整周期 $20.13125220$。初始惯性切向速度其实是 $w+\Omega_pa\simeq1.001678512$，不要误把负的随棒速度当成恒星在惯性系中反转。整周期参考返回误差小于 $10^{-9}$。

$$
\boldsymbol z(T)=\boldsymbol z(0),\qquad \boldsymbol z=(x,y,u,v),\qquad (v_x,v_y)_{t=0}=(u-\Omega_py,\ v+\Omega_px)
$$

闭合必须指定参考系，并包含速度。周期返回要求位置和速度在同一参考系中一起回来。

对应阅读：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

## 这些母轨道只比圆稍微扁一点

读者程序重新求出 $a=1.44,1.45,1.46$ 的三个邻近母轨道，以及 $a=1.89,1.90,1.91$ 的另外三个。内侧一组沿垂直棒轴方向略长，外侧一组沿平行棒轴方向略长；它们都在共转外侧，分别位于圆轨道外Lindblad参考位置的两边。

轨道离圆很近。另画半径变化，就能看见可量的响应，不必把真实轨道拉成夸张椭圆。六个解给出两处局部样本，没有穿过共振连续追踪整条分岔，也没有据图形给它们套上文献中的标准轨道族名称。

[![内外两个新求周期母轨道，相对棒转过一圈时的真实半径变化。](https://sandboxphysics.com/stories/galactic-bar-orbits/materials/periodic-parents-radial-response.svg)](https://sandboxphysics.com/stories/galactic-bar-orbits/materials/periodic-parents-radial-response.svg)

自有整周期积分。纵轴是当前半径减初始半径，横轴是相对棒顺时针走过的圈数。这里放大的是诊断曲线，没有拉伸空间轨道。

对应阅读：Sandbox Physics independent galactic-bar scientific reader v1

## 轻轻挪开一点，才能谈线性响应

先固定Jacobi值，扣除只让恒星沿同一轨道早到或晚到的相位差，再看剩下的微小偏移：它是一圈圈放大，还是在局部来回摆动？后者在线性近似下称为椭圆型。在母轨道的初态上加一个极小位移，走完一个母轨道周期。在线性近似下，矩阵 $M$ 把初始偏移映到最终偏移；它的特征值告诉我们，这些小偏移是逐周期放大，还是主要改变相位。

正则变量，是让运动满足标准哈密顿方程的一组成对坐标和动量。运动会保持这些成对变量之间的有向面积结构，称为辛结构。检查它是在检验线性映射应有的性质，特征值则另回答小扰动是否增长。 这是正则坐标下的技术几何检查，要求比单纯体积守恒更强。

把原运动方程对微小偏移展开，只留一阶项，就得到变分方程。读者程序分别用变分方程，以及正负微小位移后的独立积分来求这个映射。检查哈密顿辛结构前，要把随棒速度换成正则动量 $(p_x,p_y)=(u-\Omega_py,v+\Omega_px)$；随棒速度本身不是正则动量。

两个基准母轨道的非平凡特征值对都在单位圆上，并远离实轴端点，理论上为一的平凡特征值对也单独核对过。这支持固定Jacobi截面、去除相位方向后的线性椭圆型判断。有限扰动的非线性行为仍可能不同，另外四个邻轨没有继承这项稳定性标签。

对应阅读：Sandbox Physics independent galactic-bar scientific reader v1

## 向外走，与换一个导引半径不同

在轴对称势里，恒星可以不断向内、向外摆动，而 $L_z$ 不变，导引半径也不变。棒或旋臂施加力矩时，$L_z$ 才可能一同变化。

读者保存的示踪星群中，一个例子从参考导引半径 $R_g=1$ 出发，到末态有 $\Delta L_z\simeq-0.0596448$。取 $v_0=1$，末态参考导引半径就是 $0.9403552$。这是棒持续存在期间的一次首末比较，之后还可能继续交换。

Sellwood与Binney的短暂旋臂实验多走了一步。共转附近的星能相对波形走出马蹄状交换；若扰动一直存在，后续半程可能又把角动量换回去。旋臂若在交换过程中消退，恒星就可能把改变后的导引半径保留下来。一连串这样的事件，会把不同出生半径的恒星混在一起。

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

## 研究哪段历史，就解哪段历史

固定外棒势适合隔离引力、参考系和周期返回。增长或消退的旋臂需要完整时间历史；逐渐减速的活棒，还会与周围星盘和暗晕交换角动量。

2026年前沿研究中，向外移动的共转共振能带走部分被俘获恒星，并保留一些捕获顺序。它们的慢共振角在一个区域内来回摆动；这里六条周期母轨道的相对角度则绕棒走完一圈。两类问题都涉及共振，问的却是不同的相空间运动。

对应阅读：Sellwood & Binney (2002), Radial mixing in galactic discs, astro-ph/0203510v2 · Chiba, Fujii, Baba, Dubinski & Schönrich (2026), Tree-ring structure of Galactic bar resonance in N-body simulations, arXiv:2604.01294v1

## 练习与解释性答案

1. 一条轨迹回到原来位置，速度却不同，算周期轨道吗？

   展开解释

   不算。之后的运动取决于位置和速度。周期性要求在指定参考系内，完整状态一起返回。

2. 返回速度在 $w=-0.48$ 与 $w=-0.42$ 两端都为负，是否证明区间里没有周期根？

   展开解释

   不能。两端同号只说明这次括区不满足求根要求，非单调函数仍可能在里面有偶数个零点。原宽区间尝试保留为失败；另行预登记的窄区间才找出并验证母轨道。

3. 基准母轨道是线性椭圆型，旁边一个半径相近的母轨道也能直接称稳定吗？

   展开解释

   不能。基准检查针对自身轨道和Jacobi值附近的无穷小扰动。邻近母轨道要各自分析，有限振幅的非线性行为还需要更多证据。

## 动手检验预测

### [把母轨道与普通释放放在一起看](https://sandboxphysics.com/zh/labs/galactic-bar-orbits/)

1. 选垂直母轨道预设，固定其强度和模式速度，跟住附近一颗释放星。
2. 换到平行母轨道预设，用独立半径诊断比较真实径向幅度。
3. 改变一个物理输入，说明固定场母轨道参考的适用条件，再跟随新接入的示踪星群结果。

[上一篇：旋转棒怎样交换能量与角动量](https://sandboxphysics.com/zh/stories/galactic-bar-orbits/torque-energy-and-a-rotating-bar/) [下一篇：从保存状态重算一条棒势轨道](https://sandboxphysics.com/zh/stories/galactic-bar-orbits/replay-a-bar-orbit-from-saved-states/)

## 来源与计算范围

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. [Chiba, Fujii, Baba, Dubinski & Schönrich (2026), Tree-ring structure of Galactic bar resonance in N-body simulations, arXiv:2604.01294v1](https://arxiv.org/abs/2604.01294v1)

   固定2026年4月1日v1，PDF草稿标4月3日。完整二十一页含附录、参考文献已读，关键印刷页3、5、6、7、8、10、13、15、18已实际渲染查看。对应A/B活静态晕、C/D卫星对照、式5—11与附录A—D；不再发布原PDF或原图。

4. [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对照、圆轨道与本轮极限、规定增长功、六个新周期根和三张自有合成图。不含交互程序、完成轨道输入、原论文或观测数据。