---
title: "怎样从混合模读出核心转动"
description: "逐步读Li2024的模态拟合和近简并检验，正确换算目录转速，再追踪混合模身份与表面修正。"
canonical_url: "https://sandboxphysics.com/zh/stories/stellar-seismology/mixed-modes-core-rotation/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/stellar-seismology/mixed-modes-core-rotation/"
last_modified: "2026-11-03"
---

[读懂恒星内部](https://sandboxphysics.com/zh/stories/stellar-seismology/)

AC26 · L4

# 怎样从混合模读出核心转动

Sandbox Physics · 2026-11-03

混合模对核心和包层的反应权重不同，频率分裂便能用来分开估计两处平均转速。打开Li、Deheuvels和Ballot2024年的arXiv2405.12116v1，第2.2至2.5节、图6和附录E。图6检验了一个关键问题：自转本身会改变邻近模态的耦合，简单对称分裂不一定够用。

## 开始前需要什么

先读信号分析章和真实数据教程。下面在实际用到时说明角频率、惯量权重、方位阶数和近简并。

## 先看论文究竟拟合了哪些峰

在读转速之前，先打开第2.2节。作者首先用洛伦兹形拟合径向和四极、以声学运动为主的模态。式1的eta定义为半高半宽，半高全宽是它的两倍，混用这两个量会改变对阻尼的解释。

拟合后的径向、四极频率建立了声学排列，其中曲率和小频移都重要，真实演化恒星不能总用完美等间隔阶梯代替。偶极纯声学参照还多一个频移参数，因为它在相邻径向模态之间的位置会随演化改变。

Zenodo文件保留了这些不同阶段：fitting\_detail给出拟合量和误差，偶极bin\_points则只是在功率谱里选出的粗略起点。这些保留许多小数位的值用来引导搜索，频率误差仍需通过拟合估计。

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants · Li, Deheuvels & Ballot, Zenodo13207485 v2

## 带着单位读两个传播腔的相位

第2.3节把混合模连到包层中的声学腔和核心中的浮力腔。每个腔的相位描述波怎样容纳其中，耦合因子则表示两处沟通得多强。解出相位关系，就能得到围绕纯声学排列的许多混合模。

当倒数频率要与以秒计的重力模周期间隔相减、相除时，频率必须先用赫兹。190微赫兹的倒数约5263秒；若周期间隔约86秒，对应的是几十个重力波阶数，不是远小于一个周期。

相位约定尤其要注意。式9后面明确提醒：目录采用的重力相位偏置，与正确纯重力弱耦合定义相差二分之一。读论文方程时保留它的整套目录约定，跨到另一种定义比较时，再作作者指出的半周期换算；不能把两套定义混进同一个零耦合检验。

数值求根也要分清根和正切极点。改用连续的正弦、余弦形式，或按相位分段，就能避免把极点认成模态。检验脱耦极限和预期分支数，比画出一片像样的密集梳状谱更能说明方法是否正确。

$$
\tan\theta_{\rm p}=q\tan\theta_{\rm g}
$$

$q$为耦合因子。这是实际浮力/声学混合模的渐近关系，与可操纵的径向声腔分开使用。

$$
\theta_{\rm p}=\pi\frac{\nu-\nu_{\rm p}}{\Delta\nu},\qquad\theta_{\rm g}=\frac{\pi}{\Delta\Pi_1}\left(\frac1\nu-\frac1{\nu_{\rm g}}\right)
$$

声学参照是纯p模频率，重力参照采用论文声明的周期/相位约定。只有频率、周期单位一致，相位才无量纲。

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants

## 让不同模态给核心、包层不同权重

自转会分开方位阶数$m$不同的模态，一个偶极多重峰最多可含三个分量。论文第一种方法，把分裂写成核心与包层平均转速的加权组合。

权重$\zeta$是重力腔中的模态惯量占总惯量的比例。惯量可以理解为：给定表面位移归一化后，有多少内部物质参与运动。重力主导模态的权重接近一，对核心更敏感；声学主导模态则更能反映包层。

读目录前，先检查二和圆周率。Omega\_g、Omega\_p给出角转速，以每秒百万分之一弧度计，来源说明要求除以$2\pi$，才得到循环频率微赫兹。纯重力模的分裂公式还多一个二分之一。

例如KIC3860179的目录角频率数值为3.717347和0.510525，换成循环转速约为0.592和0.0813微赫兹。按一阶近似，纯重力主导分裂约0.296微赫兹，不是3.717微赫兹。这里是在换算外部目录结果，没有用本教程275天子集重新测转速。

$$
\delta\nu_{\rm rot}=\zeta\frac{\Omega_{\rm core}}{4\pi}+(1-\zeta)\frac{\Omega_{\rm env}}{2\pi}
$$

这是相应模态与耦合条件下的一阶关系。$\Omega$是弧度每秒的角频率，分裂则是循环频率。

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants · Li, Deheuvels & Ballot, Zenodo13207485 v2

## 用图6看清对称分裂何时失效

近简并指未扰动频率彼此足够接近，自转扰动能够混合它们的本征函数。转速仍然可以比振动频率小得多，却已经足以影响相邻模态之间的小间隔。

第二种方法先按方位阶数移动纯p模和纯g模频率，再为每个分量解耦合关系，让耦合排列也随自转改变。先算好混合模、最后加上对称分裂，一般无法表示同样的变化。

图6采用KIC11245496作为参照。红点是另一套近简并计算产生的合成多重峰不对称，灰色结果来自第二种方法的拟合。背景显示惯量权重，较强的不对称集中在声学主导模态附近。先沿着这些轴、符号和背景读图，再判断一个负包层转速估计意味着什么。

附录E给出独立矩阵观点：把位移展开在一组正交归一的相邻未扰动模态上，对角矩阵放角频率平方，自转算子则同时贡献对角和非对角扰动。对整个矩阵求本征值，才保留模态之间的混合；得到的是角频率平方，还需要相应开方、换算。

$$
\nu_{{\rm g},m}=\nu_{\rm g}-m\frac{\Omega_{\rm core}}{4\pi},\qquad\nu_{{\rm p},m}=\nu_{\rm p}-m\frac{\Omega_{\rm env}}{2\pi}
$$

这里沿用论文的方位符号，移动后的参照频率进入耦合计算，不是给所有最终观测峰统一加一个频移。

$$
M=A+R,\qquad A_{ij}=\omega_{0,i}^{2}\delta_{ij},\qquad R_{ij}=\langle\delta\mathcal L_{\rm rot}\boldsymbol\xi_i,\boldsymbol\xi_j\rangle
$$

内积按质量加权。要复现论文那颗参照星，需要实际恒星本征函数和旋转核。

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants

## 独立算一个小型近简并例子

配套mixed-mode-limits程序取两个正交的p/g基模，建立正定频率平方耦合矩阵，比较厄米本征值求解与独立二次特征方程根，并检验零耦合能返回纯模态。

示意纯模态只相差0.03微赫兹时，旋转多重峰不对称约达0.170微赫兹；把两频率拉开到3微赫兹，最大不对称降到约0.000448。邻近模态分得更开，混合便弱得多。

这个两模例子检验机制和数值极限，没有校准到论文耦合因子的换算，也没有完整恒星本征函数旋转核，更不声称复现图6。它帮助我们先看清非对角混合为什么会改变答案，再进入更大的恒星计算。

```
python3 science/mixed-mode-limits.py --input input/helmi__HD128279__HD_128279_observed_modes.csv --out rerun/mixed-mode-limits.json
```

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants

## 读总体结论时保留选择条件

论文从2495颗具有清晰混合模排列的恒星开始，识别800组双峰和1206组三峰，得到2006颗核心转速。其余恒星没有提供同样已分辨的分裂信息，没测出分裂不能改写为测得零转速。

对243颗包层转速显著大于零的恒星，论文给出约20的核心/包层转速比，并且分布很宽。这句话属于探测到包层转动的子样本，不能直接给每颗红巨星分配一个普遍比例。模态可见性、倾角、观测分辨能力和拟合方法，都会影响谁能进入这组。

原目录还保留缺失哨兵，部分辅助误差项和光度说明需要额外辨认。练习原样保留这些值，未解决的辅助字段就不进入拟合，不暗自修正符号或把光度解释成另一种量。拿目录建立新的总体先验前，应先读来源说明。

对应阅读：Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants · Li, Deheuvels & Ballot, Zenodo13207485 v2

## 跟着表面项检查模态身份

再回到Lindsay2025附录A和图5。恒星模型的近表面层会移动声学频率，Ball–Gizon一类表面项，用逆频率和三次频率项，再除以声明的模态惯量表示这种变化。

对于混合模，移动声学参照也可能改变它与重力腔的耦合。给每个最终峰盲加平滑频移，未必能保持身份。HD128279表格就是一个具体提醒：未修正偶极列的行序，与修正偶极列不同，逐行相减不能识别物理表面频移。

两项研究的对象也要各自保留。Kepler目录研究红巨星转动总体，HD128279的TESS例子则比较贫金属模型与年龄估计，它们由模态物理相连。某颗星的条件模型年龄，不会替另一颗星提供转速误差。

下次读同类论文，可以先写下模态辨认、角频率/循环频率、相位约定、惯量归一化、分辨能力和近简并处理，再读醒目的转速结论。这些选择决定了数字究竟在回答什么。

$$
\delta\nu_{\rm surf}=\frac{a_{-1}(\nu/\nu_{\rm ac})^{-1}+a_3(\nu/\nu_{\rm ac})^3}{\mathcal I}
$$

按这个归一化，系数带频率单位，$\mathcal I$及其表面归一也须说明。原论文有勘误；包里采用发表模型表作外部参照，没有重新复现完整表面项计算。

对应阅读：Lindsay et al. 2025, Precise Asteroseismic Ages for the Helmi Streams · Ball & Gizon 2014, A new correction of stellar oscillation frequencies for near-surface effects · Lindsay et al., Zenodo15102877 v1 model/observed-mode archive

## 练习与解释性答案

1. 目录角转速为$4\times10^{-6}\,\mathrm{rad\,s^{-1}}$，循环频率与一阶偶极公式里的纯重力分裂各是多少？

   展开解释

   先除以$2\pi$，得到约0.637微赫兹的循环转速；纯重力分裂再取一半，约0.318微赫兹。角转速、循环转速与模态分裂是三个不同量。

2. 自转相对振动频率仍很小，为什么也可能造成强近简并混合？

   展开解释

   混合要比较的是邻近模态间很小的频率差。相对190微赫兹微不足道的移动，可以超过0.03微赫兹间隔，因此非对角耦合仍会显著改变本征函数。

3. 为什么一个来源偶极种子还不够测核心转速？

   展开解释

   它没有拟合频率误差，也未保证多重峰身份。转速推断还需要方位分量、模态惯量与耦合信息，并适当处理近简并和观测窗口。

## 动手检验预测

### [改变内部一层，读出模态变化](https://sandboxphysics.com/zh/labs/stellar-oscillation-modes/)

1. 查看真实Kepler谱，找出相邻径向参照之间有多少候选峰。
2. 运行独立两模计算，比较纯频率接近和分开时的变化。
3. 回到图6，比较不对称多重峰与对称一阶分裂，并指出测量转速还需要哪些真实恒星模态信息。

[上一篇：用TESS光变检验红巨星模型](https://sandboxphysics.com/zh/stories/stellar-seismology/tess-red-giant-model/)

## 来源与计算范围

1. [Lindsay et al. 2025, Precise Asteroseismic Ages for the Helmi Streams](https://arxiv.org/html/2507.01091v1)

   固定arXiv2507.01091v1；第II.2、III、V.3、V.4节，图3/5、表3/4/6和附录A；已发表于ApJ 989, 189，完整优化和年龄误差属于外部发表结果。

2. [Lindsay et al., Zenodo15102877 v1 model/observed-mode archive](https://zenodo.org/records/15102877)

   元数据明确CC-BY4；选取的HD128279观测、history/index、八个GYRE结构和发表模态表完整核验；只有一条初始物理轨迹，不含840轨迹优化集。

3. [Li, Deheuvels & Ballot 2024, Core and envelope rotation rates for 2006 red giants](https://arxiv.org/html/2405.12116v1)

   固定2405.12116v1；第2.2至2.5节、近简并模拟图6和附录E，已发表于A&A 688, A184；角频率/循环频率单位和g模相位约定须分清。

4. [Li, Deheuvels & Ballot, Zenodo13207485 v2](https://zenodo.org/records/13207485)

   CC-BY4旋转目录及KIC3860179选定模态文件；偶极bin\_points仅是粗略种子，l0/l2拟合表另含误差，含混的辅助列不用于拟合。

5. [Ball & Gizon 2014, A new correction of stellar oscillation frequencies for near-surface effects](https://arxiv.org/abs/1408.0986)

   按模态惯量加权的逆频率/三次频率表面项；具体应用读Lindsay2025附录A。原论文有勘误DOI10.1051/0004-6361/201424325e；课程不声称独立复现全文。