恒星结构与流体静力平衡

概述#

一颗恒星为什么不会在自身引力下坍缩，也不会因内部压力而爆炸？答案在于流体静力平衡（hydrostatic equilibrium）——一种引力向内拉与压力向外推之间精确而脆弱的平衡。理解这一平衡是现代恒星物理学的基石。它告诉我们：恒星不是一团燃烧的气体，而是一台精密调节的热核机器，其内部结构遵循可以用数学方程描述的严格规律。

常见误解是把恒星看成一团着了火的气。燃烧是表面看到的结果，支撑来自压力梯度。没有向外递减的压强，上层会被压垮；没有引力，气体会散掉。另一类误解是平衡一旦建立就永远不动。脉动、壳层闪燃、核心坍缩都是短暂偏离；主序上绝大部分时间，偏离会在动力学时标内被拉回来。恒星演化走的是热核时标和开尔文–亥姆霍兹时标，比流体静力恢复慢得多。

流体静力平衡方程#

考虑恒星内部距中心距离 $r$ 处的一薄层壳体。这层壳体受到两个力的作用：

• 向内的引力：来自半径 $r$ 以内所有质量 $M(r)$ 对该壳体的引力
• 向外的压力差：下方（更高压力）和上方（更低压力）之间的压力差

当两者平衡时，得到流体静力平衡方程：

$$\frac{dP}{dr} = -\frac{G M(r) \rho(r)}{r^2}$$

其中 $P$ 是压力，$\rho$ 是密度，$G$ 是万有引力常数。这个方程意味着：越靠近恒星中心，上方压着的质量越多，所需的支撑压力也越大。

对于太阳，中心压力约为 $2.5 \times 10^{16}$ Pa（约 $2.5 \times 10^{11}$ 个大气压），中心温度约为 $1.57 \times 10^7$ K（源自标准太阳模型，Bahcall et al. 2005，The Astrophysical Journal，621，L85）。这些极端条件正是维持核聚变反应所必需的。

方程右边的负号表示压强必须向外递减。某一体积元若被轻微压缩，密度升高、压强升高，净力把它推回；若被拉开，压强梯度不够托住重量，气体回落。恢复的时钟是自由下落时标 $\sqrt{R^3/GM}$，对太阳大约一小时量级。太阳看起来千年不变，是因为我们看到的是热核时标，不是流体静力时标。日震看到的声波，正是在这条平衡附近的小振动，详见太阳结构与活动。

把同一方程写到行星或气体巨星上，形式不变，只是 (M(r)) 与 ( ho) 的数值换了场景。恒星特有的地方在于：压强梯度必须维持到足以点燃或曾经点燃热核反应的中心温度，否则「发光的气体球」会沿开尔文–亥姆霍兹时标收缩，而不是停在主序上。

质量连续性方程#

与流体静力平衡方程配套的是质量连续性方程，描述质量随半径的分布：

$$\frac{dM(r)}{dr} = 4\pi r^2 \rho(r)$$

这两个方程构成了描述恒星结构的基础方程组的前两条。

质量连续说的是：半径 $r$ 处往外走一小步，包进去的质量等于当地密度乘球壳体积。中心边界是 $M(0)=0$，表面是 $M(R)=M$。给定密度分布，积分一次就得到 $M(r)$；反过来，测到 $M(r)$ 的斜率就等于当地密度。流体静力方程里的 $M(r)$ 必须和这条连续方程自洽，不能随手指定一个质量再另给一套密度。

两条方程还不够闭合。压强要靠状态方程写成密度和温度的函数，温度又要靠能量传输和产能。所以“结构”不是单独解一条静力平衡，而是四条方程加一张状态方程一起迭代。质量连续保证球对称质量不丢不造；静力平衡保证这团质量不会在下一步自由落体。缺了连续方程，静力平衡右边的 $M(r)$ 就没有定义。

数值上，质量连续与静力平衡常被合成「以质量为自变量」的形式：把 (r) 写成 (M) 的函数，便于从中心积分到表面。坐标怎么选不改物理，只改边界条件的写法——中心 (r=0)、(L=0)，表面光学深度趋于零——两点边值问题的本质仍在。

压力的物理来源#

恒星内部维持压力的机制随密度和温度不同而变化：

热压力（气体压力）#

大多数主序星的主要支撑来自理想气体压力：

$$P_{gas} = \frac{\rho k_B T}{\mu m_H}$$

其中 $k_B$ 是玻尔兹曼常数，$T$ 是温度，$\mu$ 是平均分子量，$m_H$ 是氢原子质量。恒星核心的高温（数百万至数十亿开尔文）产生了足以抵抗引力的热压力。

热压来自粒子热运动的动量交换，微观图像见气体动理论。氢聚成氦以后，同样质量里的粒子数变少，$\mu$ 上升，同温同密下压强下降。核心必须升温才能继续托住外层，主序星的光度因此会随年龄缓慢增加。这不是额外假设，是理想气体状态方程加质量连续和静力平衡的直接推论。

辐射压力#

在质量大于约 $10\,M_\odot$ 的大质量恒星中，光子的辐射压力变得重要：

$$P_{rad} = \frac{4\sigma T^4}{3c}$$

其中 $\sigma$ 是斯特藩-玻尔兹曼常数，$c$ 是光速。对于质量远大于 $100\,M_\odot$ 的极超质量恒星，辐射压力甚至超过气体压力，导致恒星接近爱丁顿极限——辐射力等于引力的临界状态。

辐射压按 $T^4$ 走，质量越大核心越热，辐射压占总压的份额就越高。爱丁顿光度假定球对称、电子散射不透明度。主序大质量星可以长期靠近这条线，但不等于已经炸开：真实包层有对流、有孔隙，光子可以从低密度通道漏走。一旦局部光度超过当地爱丁顿值，外层会膨胀、脉动或抛出物质，流体静力平衡在包层变成近似。大质量恒星演化里的强星风，有一部分就是辐射压把静力平衡顶出了稳定区。

辐射压份额升高时，绝热指数被拉低，恒星对扰动更「软」。这不是说大质量星随时解体，而是说包层更容易膨胀、更容易触发强星风与脉动；流体静力平衡从「几乎刚性」变成「需要持续质量损失来维持的准平衡」。

简并压力#

当恒星演化到晚期、核燃烧停止后，普通热压力无法维持支撑。此时，量子力学的泡利不相容原理提供了额外压力。电子简并压力（白矮星）或中子简并压力（中子星）与温度无关，是由量子态占据的极限决定的。

简并压不靠温度。费米子不能挤进同一个量子态，密度升高时动量必须升高，于是有压强。泡利不相容原理因此能在核聚变停掉之后继续托住白矮星。Chandrasekhar（1931，The Astrophysical Journal，74，81）证明，相对论电子简并对应一条最大质量；对每个核子对应一个电子的碳氧组成，这条上限约为 $1.4\,M_\odot$。超过它，电子简并不能再提供静力平衡，核心继续坍缩。

中子简并要把广义相对论写进静力平衡，得到托尔曼–奥本海默–沃尔科夫方程。Oppenheimer 与 Volkoff（1939，Physical Review，55，374）用理想中子费米气算出的上限只有约 $0.7\,M_\odot$。核力排斥会把上限抬高，现代观测到的中子星已超过两倍太阳质量，具体数字取决于致密物质状态方程。简并压与热压的关键差别是：冷却不会让白矮星瘪下去，主序星一旦失去热源却必须收缩。

能量传输方程#

仅有压力平衡还不够——恒星核心产生的能量必须传输到表面并辐射出去。能量传输遵循两种主要机制：

辐射传输#

能量以光子形式向外传播，传输效率由不透明度 $\kappa$（opacity）决定：

$$\frac{dT}{dr} = -\frac{3\kappa \rho L(r)}{64 \pi \sigma T^3 r^2}$$

在太阳内部约 $70\%$ 的半径范围内（辐射区），能量主要通过辐射传输。光子在辐射区内经历约 $10^5$ 年的随机游走才能到达表面——尽管光速极快，但物质的高不透明度使光子每次只能移动极短的自由程。

辐射传输方程把当地光度、不透明度和温度梯度锁在一起。不透明度高，同样的光度就要求更陡的温度梯度。梯度一旦陡过绝热值，辐射通道不够用，对流接手。辐射、对流与导热的分工见热传导对流与辐射，恒星内部只是把同一套选择规则放到球对称引力场里。

随机游走的「十万年」是太阳辐射区的数量级说法：自由程极短、散射次数极多，能量向外走的有效速度远低于光速。日震与中微子之所以重要，正是因为它们分别穿过声波通道与弱相互作用通道，绕开了这条光子迷宫。

对流传输#

当辐射传输无法有效输出核心产生的能量时（温度梯度超过绝热梯度），物质通过对流运动传输热量——热的气泡上升，冷的气泡下沉。太阳外部约 $30\%$ 的半径是对流区，这里的湍流对流产生了太阳表面的米粒组织（granulation）和超米粒组织（supergranulation）。

施瓦西判据（Schwarzschild 判据）：当实际温度梯度超过绝热温度梯度时，发生对流不稳定性：

$$\nabla_{rad} > \nabla_{ad}$$

其中 $\nabla_{rad}$ 是维持辐射传输所需的对数温度梯度，$\nabla_{ad}$ 是绝热对数温度梯度。

对流一旦打开，温度梯度被钉在绝热值附近，能量靠物质宏观运动搬走，而不是靠光子一步一步随机走。低质量主序星不透明度高、产能集中度低，可以全对流；太阳只在外包层对流；大质量星则往往是对流核加辐射包层。对流改变的是化学混合和温度剖面，流体静力平衡本身仍然成立，只是压强梯度要按新的密度、温度分布重算。

核能产生方程#

恒星结构的第四个基本方程描述能量产生：

$$\frac{dL(r)}{dr} = 4\pi r^2 \rho(r) \epsilon(r)$$

其中 $L(r)$ 是半径 $r$ 以内的光度，$\epsilon$ 是单位质量的核能产生率。$\epsilon$ 对温度极为敏感：质子-质子链的产能率约 $\propto T^4$，而 CNO 循环约 $\propto T^{17}$，硅燃烧约 $\propto T^{40}$。这种极高的温度敏感性使得恒星核心的热核反应具有天然的自我调节能力——温度稍有上升，产能率激增，压力增大，恒星膨胀降温，形成负反馈。

产能方程是质量连续在能量上的对应物：每往外走一层，光度增加等于该层的核产能。恒星核合成决定 $\epsilon$ 用哪条反应网，结构方程则决定温度够不够点燃它。负反馈能工作，前提仍是流体静力平衡：膨胀必须真的降低中心温度，而维里定理保证引力束缚的理想气体正是如此。

维里定理与恒星的自我调节#

维里定理（Virial theorem）是流体静力平衡的深层含义：对于引力束缚的系统，热能等于引力势能绝对值的一半：

$$E_{thermal} = -\frac{1}{2} E_{grav}$$

这导致了一个反直觉但极为重要的结论：当恒星失去能量时，它的温度反而升高。收缩 → 引力势能释放 → 一半转化为热能（温度升高）、一半辐射出去。这解释了为什么主序前恒星（原恒星）在缓慢收缩过程中会变热，最终点燃核聚变。

简并星没有这条路。电子或中子的费米能已经定了压强，再损失热量只是降低热运动的那一小部分，半径几乎不动。所以白矮星冷却变暗，主序前的太阳型星收缩变热。恒星形成末段的赫亚什轨迹，正是维里定理在还没点火时的应用。

维里定理还规定了「点火」的门槛：原恒星收缩释放的能量一半用来升温，一半辐射掉；当中心温度够到质子隧穿库仑势垒的速率时，核产能才接手开尔文–亥姆霍兹收缩。主序不是外加开关，是静力平衡加维里关系走到核时标的自然结果。

标准太阳模型#

恒星结构方程组（流体静力平衡 + 质量连续 + 能量传输 + 核能产生）加上合适的状态方程和不透明度表，可以用数值方法求解，得到标准太阳模型（Standard Solar Model, SSM）。

标准太阳模型的预言已通过多种途径检验：

• 日震学（helioseismology）：通过分析太阳表面的声学振荡（压力模式和重力模式），精确探测太阳内部的密度和声速分布。日震学测量与标准太阳模型的预言吻合到 $0.1\%$ 的精度（但在对流区底部附近存在约 $1\%$ 的偏差，称为"太阳丰度问题"）
• 太阳中微子：核心核聚变产生的中微子已被地面中微子探测器（Super-Kamiokande、SNO、Borexino）探测到。历史上的"太阳中微子问题"（测量值只有预言的三分之一）通过发现中微子振荡而解决，同时确认了标准太阳模型的正确性

太阳丰度问题提醒我们：结构方程本身可以很稳，输入的不透明度与重元素丰度仍能在对流区底留下约百分之一的声速偏差。日震学检验的是整条解，不是单独一条静力平衡公式；差在表上，往往差在原子物理与光谱定标。

恒星结构方程的更广泛应用#

这套方程不仅适用于太阳，而且适用于从棕矮星到超质量恒星的一切恒星：

表里换的是状态方程，不是另写一套静力平衡。从棕矮星到极超质量星，左边仍是 $\mathrm{d}P/\mathrm{d}r = -GM\rho/r^2$，变的是 $P(\rho,T)$ 由谁主导。质量连续同样不变：$M(r)$ 永远是密度的球积分。换天体种类，先问压强从哪来，再问产能开不开，不要先改引力。

为什么这很重要#

流体静力平衡方程看似简单，却是宇宙中所有恒星存在的根本原因。它决定了恒星能活多久（质量越大、平衡越脆弱、寿命越短）、核心能到达多高的温度（决定能合成哪些元素）、以及最终会以何种方式死亡。理解恒星结构，就是理解宇宙中绝大多数发光物质的本质。标准太阳模型的成功，是20世纪天体物理学与粒子物理学合作的最佳范例之一。

热压托住主序，辐射压给大质量星划出一条接近爱丁顿的边界，简并压在核燃料用尽后接手。质量连续则保证这些压强作用在正确的 $M(r)$ 上。四句话合在一起，才是“恒星为什么既不塌也不炸”的完整答案。

跨域连接#

• 大气层结构：地球大气的压强递减服从同一条方程——压强梯度等于密度乘重力，等温时给出以标高为尺度的指数衰减。恒星的唯一差别是温度随深度剧增，标高不再恒定，压强才能跨越十几个数量级。推论：测出某处的压强梯度就等于称出了该处以内的质量。
• 微分方程：四条结构方程构成的是两点边值问题（中心处质量与光度为零，表面处压强与温度趋零），不是初值问题，因此只能用松弛迭代求解。这也解释了为什么给定质量与成分后解几乎唯一——主序才会是一条窄带而不是一片云。
• 气体动理论：热压强公式里的平均分子量随氢聚变成氦而上升，同温同密下压强反而下降，核心必须升温补偿。推论：主序星的光度会随年龄缓慢上升，而这一步完全不需要引入任何新物理。
• 原子结构：能量传输方程里的不透明度主要来自重元素离子的束缚—自由跃迁，取决于电子壳层结构。于是氧、氖的丰度定标一变，不透明度表就变，结构解跟着变——这正是日震学在对流区底部约百分之一偏差的来源。
• 中微子探测器：光子在辐射区随机游走十万年量级才出来，携带的是路径信息而非核心的即时状态；只有中微子几乎无阻碍穿出，能直接检验产能率。历史上的太阳中微子亏缺最终由中微子振荡解释，同时反过来确认了标准太阳模型。

参考文献#

• Bahcall, J.N., Serenelli, A.M., & Basu, S. (2005). New Solar Opacities, Abundances, Helioseismology, and Neutrino Fluxes. The Astrophysical Journal, 621, L85–L88.
• Chandrasekhar, S. (1939). An Introduction to the Study of Stellar Structure. University of Chicago Press.
• Hansen, C.J., Kawaler, S.D., & Trimble, V. (2004). Stellar Interiors: Physical Principles, Structure, and Evolution (2nd ed.). Springer.
• Christensen-Dalsgaard, J. (2002). Helioseismology. Reviews of Modern Physics, 74, 1073–1129.
• Kippenhahn, R., Weigert, A., & Weiss, A. (2012). Stellar Structure and Evolution (2nd ed.). Springer.
• Chandrasekhar, S. (1931). The maximum mass of ideal white dwarfs. The Astrophysical Journal, 74, 81–82.
• Oppenheimer, J.R. & Volkoff, G.M. (1939). On massive neutron cores. Physical Review, 55, 374–381.

延伸阅读#

• Prialnik, D. (2009). An Introduction to the Theory of Stellar Structure and Evolution (2nd ed.). Cambridge University Press. — 适合入门的教科书，推导清晰
• Basu, S. & Antia, H.M. (2008). Helioseismology and Solar Abundances. Physics Reports, 457, 217–283. — 日震学约束标准太阳模型的综述
• Eddington, A.S. (1926). The Internal Constitution of the Stars. Cambridge University Press. — 历史性经典，理解流体静力平衡的思想来源