跳到主要内容
Sandbox Physics
读懂恒星内部

AC26 · L2

声音怎样穿过分层恒星

Sandbox Physics ·

只改变内部一层，不同模态的频率会移动不同幅度。把界面两侧的波接起来，就能算出这种变化，也能用精确拼接解检验数值求解器。频率排列之所以有用，正因为它比单个声程时间携带更多结构信息。

让旅程的一部分变慢

保留常密度球形声腔和外缘压力释放条件。核心从中心延伸到半径aa，声速为c1c_1；外层一直到RR，声速为c2c_2。两处声速都为正。在这个模型里，改变声速相当于改变局部可压缩性，密度仍然不变。

采用无量纲半径x=r/Rx=r/R、界面位置b=a/Rb=a/R和声速比v=c1/c2v=c_1/c_2。时间用R/c2R/c_2计，普通频率用c2/Rc_2/R计，这样可以把内部结构与整体尺寸、声速尺度分开。无量纲频率为一，不一定就是一赫兹。

从中心走到外缘的声学时间，等于两层耗时相加。它给出一个高阶频率间隔估计，但没有要求相邻低阶模态也必须严格保持同样的距离。

τ=ac1+R−ac2,τc2R=bv+1−b\tau=\frac{a}{c_1}+\frac{R-a}{c_2},\qquad \frac{\tau c_2}{R}=\frac{b}{v}+1-b
每一项都是路程除以当地波速，界面移动时，这个无量纲传播时间连续变化。
ΔνasymRc2≃12(b/v+1−b)\frac{\Delta\nu_{\rm asym}R}{c_2}\simeq\frac{1}{2(b/v+1-b)}
下标提醒我们这是高阶声学估计，内部跃变的反射会扰动单个频率间隔。

对应阅读：MESA School 2026, Lab 3: Asteroseismology

在界面把两侧的波接起来

如果界面没有额外的表面作用，压力必须连续，气体的法向速度也必须连续。由运动方程可知，压力梯度除以密度应在两侧连续。这里密度相同，压力梯度本身便连续。

在各自均匀的一层内，变换后的ψ=rP\psi=rP仍满足正弦方程。到了界面，ψ\psi和径向导数都连续。任选一个非零整体归一化，把中心规则解传播到界面，再接着穿过外层，就能找到整个球体的形状。

如果密度也跳变，压力梯度的连续条件就会改变，不能再直接沿用ψ\psi导数连续。此时应回到下面保留密度因子的压力方程。

[P]=0,[P′ρ]=0[P]=0,\qquad\left[\frac{P^\prime}{\rho}\right]=0
方括号表示界面两侧之差，两式分别来自压力和法向速度的连续性。
−ddr(r2ρP′)=ω2r2ρc2P-\frac{d}{dr}\left(\frac{r^2}{\rho}P^\prime\right)=\omega^2\frac{r^2}{\rho c^2}P
这条压力Sturm–Liouville方程保留了密度变化时需要的因子，浏览器可操纵部分仍采用常密度。

对应阅读：MESA School 2026, Lab 3: Asteroseismology

精确求解两层频率

在无量纲方程里取R=1R=1，记s=ωR/c2s=\omega R/c_2，核心波数便是s/vs/v，外层波数是ss。把中心规则解的变换导数归一为一，传播到界面时的函数值与导数都能直接写出。

在均匀层里，知道一端的值和导数，就能算出另一端。下方外层传播式因此给出外边界的变换压力，它的零点就是允许模态。表达式只含正弦和余弦，不会把正切函数的极点误当成根。

声程近似可以帮助估计根大致在哪里。真正找多个模态时，应扫描足够细的范围，给每个变号区间设括号，再逐个细化，并检验节点数。只给非线性求解器换几个初值，可能反复找到同一个根。

ψ(b)=sin⁡(sb/v)s/v,ψ′(b)=cos⁡(sb/v)\psi(b)=\frac{\sin(sb/v)}{s/v},\qquad\psi^\prime(b)=\cos(sb/v)
它们来自核心规则解，整体归一化不影响允许的频率。
F(s)=ψ(b)cos⁡[s(1−b)]+ψ′(b)ssin⁡[s(1−b)]=0F(s)=\psi(b)\cos[s(1-b)]+\frac{\psi^\prime(b)}{s}\sin[s(1-b)]=0
允许模态让外边界残差为零，把求得的根除以2π2\pi，才是无量纲普通频率。

对应阅读：MESA School 2026, Lab 3: Asteroseismology

完整算一组不等间隔

取b=0.37b=0.37、v=0.4v=0.4。无量纲声程时间为1.555，对应的高阶间隔估计约为0.322。但前三个精确普通频率分别为0.298308、0.686279和0.919733，前两个间隔约为0.388和0.233，彼此不同，也不等于声程近似。

再放回一个物理例子：球半径一米，外层声速每秒300米，核心声速就是每秒120米。把无量纲频率乘以每秒300，前三阶约为89.5、206和276赫兹。前两阶相差约116赫兹，接下来的差却只有约70.0赫兹。

波不仅受到外边界限制，还会在内部界面反射。有些模态在慢核心里留下较强的压力变化，有些则给这一区域不同权重。一个平均声程无法囊括低阶模态的全部相位变化，不齐的间隔正是结构留下的信息。

import numpy as np
from scipy.optimize import brentq
b, v = 0.37, 0.4
def boundary(s):
    value = np.sin(s*b/v)/(s/v)
    derivative = np.cos(s*b/v)
    return value*np.cos(s*(1-b)) + derivative*np.sin(s*(1-b))/s
first_root = brentq(boundary, 1.0, 2.0)
print(first_root/(2*np.pi))  # 0.298307644...

对应阅读：MESA School 2026, Lab 3: Asteroseismology

把微分方程变成能检验的矩阵

一般的正声速剖面未必有两层解析解。可以用网格上的分段线性“帽子函数”近似变换后的形状：每个函数只在附近几个网格点非零。把方程乘以各个帽子函数并积分，就得到刚度矩阵和带权质量矩阵。

刚度矩阵衡量空间梯度，质量矩阵按声速平方的倒数给各层加权。式中的v(x)v(x)表示当地声速除以外层声速，在核心等于vv，在外层等于一。广义本征值问题寻找能让两项平衡的形状，本征值就是无量纲角频率的平方。两个二次型均为正，让频率平方保持正值。

如果界面落在一个网格单元内部，积分就在那里分开。把界面凑到最近网格点，会让控制移动时出现阶梯跳变；连续物理问题没有这种“吸附”，数值方法也不该凭空加上。

对上面的例子，独立一致有限元参照与精确传播根比较，前八模最大频率误差约为：64格时1.10%，128格时0.272%，256格时0.0682%。浏览器采用另行检验的正质量集中离散；192格、289个结构点的测试里，前八模最坏误差约0.179%。两者是不同数值方案和不同证据，不能混为一次收敛计算。

Kij=∫01ϕi′ϕj′ dx,Mij=∫01ϕiϕjv(x)2 dx,Ka=λMaK_{ij}=\int_0^1\phi_i^\prime\phi_j^\prime\,dx,\qquad M_{ij}=\int_0^1\frac{\phi_i\phi_j}{v(x)^2}\,dx,\qquad K\mathbf a=\lambda M\mathbf a
ϕi\phi_i是局部帽子函数，a\mathbf a包含形状系数，λ=(ωR/c2)2\lambda=(\omega R/c_2)^2。

对应阅读：MESA School 2026, Lab 3: Asteroseismology

形状也要一起检查

频率误差小，不代表模态身份一定对。还要确认频率排序正确，第nn个径向压力模态有n−1n-1个内部节点，外缘压力为零，中心满足规则条件，离散方程残差也应很小。

不同模态还满足质量权重下的正交关系，即带权重叠积分为零，并不是两条曲线永不相交。推导办法是：把每个本征方程乘以另一个模态，相减后积分。边界项消失，留下频率平方之差乘以带权重叠；频率不同，重叠便必须为零。

频率、形状和边界通过这些独立检验后，图中的节点才有可靠的物理含义。显示的节点圆环来自提供模型频率的同一组本征函数。

∫01ψm(x)ψn(x)v(x)2 dx=0(m≠n)\int_0^1\frac{\psi_m(x)\psi_n(x)}{v(x)^2}\,dx=0\qquad(m\ne n)
这个权重来自常密度变换后的波动方程，积分使用求得的本征函数。

对应阅读：MESA School 2026, Lab 3: Asteroseismology

一个音调仍容许不同内部

比较两个精确模型：第一个界面为0.37、核心声速比0.65；第二个界面为0.60、核心声速比约0.797655。内部结构明显不同，基本无量纲频率却完全相同，都约为0.421445。

它们第三阶频率约为1.26389和1.29771，多一个分辨清楚的模态便有机会区分这两个例子。不过，这仍不能证明任意连续声速剖面只有一个解。每个测量模态给出有限约束，自由变化的剖面却有许多自由度。

真实恒星建模还加入质量、成分、混合和大气物理，所以结构比较必须说明允许哪些模型。一条固定演化轨迹只能回答给定轨迹下的问题，要推断更一般的年龄或成分，还需要相应模型集合及不确定性。

对应阅读：MESA School 2026, Lab 3: Asteroseismology · Lindsay et al. 2025, Precise Asteroseismic Ages for the Helmi Streams

练习与解释性答案

  1. 界面穿过一个网格单元时，为什么要拆开这个单元的积分？

    展开解释

    界面两侧的声速倒数平方权重不同。在真实位置拆分，才能保留正确权重并让本征值连续变化；把界面凑到网格点，会人为地把结构改成一格一格跳动。

  2. 慢核心例子里，为什么不能把所有频率间隔强行调成高阶声程近似？

    展开解释

    声程近似概括总传播时间，低阶模态还携带界面相位和不同的分层权重。强制等间隔会删掉真实结构信息，也会掩盖与精确传播解的差别。

  3. 两个结构的基本频率相同，单独把这个频率测得更准，是否一定能区分它们？

    展开解释

    不能。例子里它们匹配的是同一个精确基本频率，提高这个数的精度并不消除多解；需要另一个模态，或者独立的结构观测量。

动手检验预测

改变内部一层，读出模态变化

  1. 选择“强烈的内部反射”，比较测量间隔备选值与橙色模型线之间的实际距离。
  2. 小幅拖动界面，观察选中模态的节点怎样连续移动。
  3. 回到整齐阶梯，解释哪处界面反射消失了，再列出描述真实恒星还需要加入的物理。

来源与计算范围

  1. MESA School 2026, Lab 3: Asteroseismology

    官方课堂的振动频率、声程积分与模型演化；这里可操纵的常密度声腔另作明确的教学模型。

  2. Lindsay et al. 2025, Precise Asteroseismic Ages for the Helmi Streams

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

反馈这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策