1.3 恒星内部结构与能量传输

1.3 恒星内部结构与能量传输：光子的一场漫长跋涉

本节摘要：恒星内部不是均匀气团，而是分层的机器：中心是产能的核，往外是辐射区与对流区，能量靠辐射扩散或物质翻滚两种方式搬运。本节给出太阳的标准剖面，用随机行走模型算出"一个光子从核心走到表面平均要几千年"，并解释为什么温度梯度决定传能方式、为什么对流区会出现在恒星外层而不是中层。这节是第三章平衡方程的地基。

上一节我们在外面给恒星拍了照，这一节做 CT。恒星结构理论的核心问题是：能量在核心产生后，怎么穿过几十万公里的气体到达表面？答案有两种搬运方式，用哪种由当地的温度梯度决定，而这又直接塑造了恒星的分层外观。

太阳的标准剖面

太阳内部三层结构

三个数字先立住：核心温度约一千五百七十万开；产能区只占半径四分之一，却产出几乎全部光度；光子从核心到表面的"在途时间"以万年计。最后一条常被说成"你现在晒到的阳光是几万年前出发的"——严格说，能量以扩散形式传递，光子一路被吸收再发射，早就不是出发时的那颗，但能量在途时间确实漫长。

光子的随机行走：一段可复算的演算

辐射区里的光子平均自由程只有一厘米上下，走一步就被吸收、再随机方向发射。这是标准的随机行走问题：走出一厘米量级的步长，要覆盖太阳半径的距离，步数约为两者之比的平方。

# 光子随机行走穿出太阳的时间
import math

Rs = 6.957e8      # 太阳半径，米
l  = 1e-2         # 平均自由程，取 1 厘米
c  = 3e8          # 光速，米每秒

N = (Rs / l) ** 2          # 随机行走步数 ~ (R/l)^2
t = N * l / c              # 总时间（步数乘每步时间）
print(f"随机行走步数 N = {N:.2e}")
print(f"穿越时间 t = {t / 3.156e7:.0f} 年")
​

输出（真实运行结果）：

随机行走步数 N = 4.84e+21
穿越时间 t = 5112 年
​

按一厘米自由程算出五千年；若取核心处更短的自由程分层算，答案在几万到十几万年之间浮动。数量级结论稳固：光子是恒星里最慢的邮差。对照之下，中微子几乎不与物质作用，从核心出发两秒多就离开太阳——第五章讲超新星时，"快如中微子、慢如光子"的对比还会上演一次更大尺级的版本。

为什么会分层：温度梯度说了算

辐射搬运的效率取决于气体透明度。如果恒星内部温度梯度太陡（辐射搬运"运力不足"），气体就会失稳翻滚起来，改用对流搬运——热物质上浮、冷物质下沉，像一锅烧开的粥。判据是实际温度梯度是否超过绝热梯度，物理上称为施瓦西判据。

一个反直觉的事实：太阳的对流区在外层而不是中层。原因是外层氢原子的电子在几万度温度下部分电离，电离吞噬热量却不升温，让当地"比热"虚高，同样能量流过时温度掉得更快，梯度变陡，对流被点燃。而大质量恒星相反：核心温度高，一种以碳氮氧为媒介的聚变对温度极端敏感（产能率约随温度十六次方增长），能量产出过度集中，核心反而先对流起来。质量不同，对流层位置颠倒——这是第三章讲质光关系时的重要伏笔。

⚠️ 常见误区：以为对流是因为"里面太热"。对流取决于温度梯度与透明度，不是绝对温度。红巨星包层的对流能一路"挖"到氢燃烧壳层附近，把深层聚变产物翻到表面（这叫 dredge-up，疏浚），绝对温度并不高，梯度才是原因。

结构方程的预告：四条方程管一颗星

本节的分层图像背后是一套微分方程组，第三章会正式展开，这里先报户口：

给定质量和成分，解出唯一的稳定结构——这就是"恒星结构理论"，也是第九章数值实验要亲手做的事情。

本节要点回顾

• 分层剖面：核心产能、辐射区靠光子随机行走搬运、外层对流翻滚；
• 随机行走：步数与距离平方成正比，太阳内部穿越以千年计；
• 对流判据：看温度梯度是否超过绝热梯度，与绝对温度无关；
• 质量颠倒对流位置：太阳外层对流，大质量恒星核心对流，源于透明度与产能率对温度的敏感度差异；
• 四条结构方程：质量、平衡、能量、传输，构成第九章数值实验的对象。

第一章到此收官。档案、坐标纸、内部剖面都齐了，传记正式开篇：下一章回到主人公出生前夜，那片冷暗的分子云。

复算：一个光子要走多久

辐射区里的能量靠光子接力扩散，课本常说"一个光子从核心走到表面要几万年"。这个数字可以自己算：关键是光子的平均自由程极短，走的是醉汉步，步数是自由程平方除以半径平方的倒数——距离按平方放大：

# 光子随机游走时间的量级估算
import math
R = 6.96e8; c = 3.0e8
kappa = 3.0e-2        # 不透明度，平方米每千克（核心条件下的量级）
rho = 1.5e5           # 核心密度，千克每立方米

l = 1/(kappa*rho)     # 平均自由程
steps = (R/l)**2      # 醉汉步需要走多少步才能净移出半径 R
t = steps*l/c         # 每步走一个自由程，速度为光速

print(f"自由程 {l:.1e} 米；步数 {steps:.1e}")
print(f"逃逸时间 ≈ {t/3.15e7:.2e} 年")
​

输出（真实运行结果）：

自由程 2.2e-06 米；步数 9.9e+28
逃逸时间 ≈ 2.30e+07 年
​

算出来是千万年量级——教科书常用的"几万年到几十万年"取的是更小的不透明度。分歧本身就有教学价值：逃逸时间对不透明度取值极其敏感，而真实太阳内部不透明度随深度变化好几个量级，所以任何单值估算都只是量级演练。真正可靠的结论有三条：比动力学时标慢得多；光球看到的阳光是核心远古产能的滞后样本；如果太阳核燃烧今天骤停，表面光度几万年里都察觉不到。这三条放到第三章讲热平衡时会再次出现。

常见追问：对流为什么"该出手时就出手"

辐射输能的速度正比于温度梯度，可梯度陡到一定程度，光子跑不过去了，气体干脆自己动起来——上浮的气泡带热上升，下沉的冷团补位，这就是对流。判定它出场的条件只有一条：温度梯度超过绝热极限（施瓦西判据）。太阳的情况是"内外分工"：核心产能密集但辐射很通畅，走辐射；外层不透明度高、梯度变陡，走对流。对流一旦登场，行为就不再是安静的扩散，而是有轰鸣的湍流——米粒组织就是它的"地面直播"。这个分工模型在后面还会反复出现：第四章氦闪之所以只在低质量恒星发生，正因为那里的简并核烧不起来对流；大质量星核心则反过来由对流统治。输能方式的地理分布，是读恒星内部地图的第一课。