第九章 亲手写传记：数值实验
第九章 · 亲手写传记：数值实验
章节摘要：前八章读的是别人写的恒星传记，本章换读者执笔：从四条结构方程出发跑一遍最简恒星模型，用最小二乘拟合复核质量-光度关系，认识专业演化代码的工具链形态，最后盘点核天体物理的未解之谜。目标不是替代专业软件，而是把"量级演算"升级成"可复现的数值实验"。
学习目标
读完本章，你应当能够：
- 写出恒星结构四方程各自的物理含义与积分顺序，并解释为什么中心边界条件最难给；
- 用桌面积分重现太阳中心压强与温度的量级，理解网格与步长的取舍；
- 拟合主序质量-光度关系并解释指数约三点五的物理来源；
- 说出专业演化工具链的输入-输出形态与常见使用误区；
- 列举核天体物理的几个前沿悬念，说明它们各自卡在哪一步。
核心概念速览
一句金句：恒星结构方程只有四条，难的不是方程，是每一层都换一个物态方程与不透明度表的耐心。
子章节导航
第一节"从四条方程开始"把第一章的物理语言落成可积分的代码。第二节"光度-质量关系拟合"教读者用观测数据亲手复核第三章的指数。第三节"演化模型工具链"面向想深入者，介绍专业代码的输入输出形态。第四节"核天体物理与未解之谜"收束全册，把还没讲完的故事列成研究清单。
子章节之间的逻辑关系
一句话论点：先复算、再拟合、后对照——三步走完，读者就从读传记的人变成能写传记的人。
四条方程（9.1 节） │ 提供骨架 ▼ 拟合质量光度关系（9.2 节） │ 观测数据上验货 ▼ 专业工具链（9.3 节） │ 从玩具到工业级 ▼ 未解之谜（9.4 节） │ 传记的下一章由读者续写
本章知识点清单
- 四条结构方程的物理对应：质量守恒、流体静力学平衡、能量守恒、能量传输；
- 中心与表面两类边界条件的" shooting "逻辑：猜中心参数，积分到表面，不匹配就回猜；
- 质量-光度关系的拟合方法与指数的物理解读；
- 演化代码的输入参数组（质量、丰度、自转、金属丰度）与输出轨迹；
- 太阳丰度问题、锂问题、第一代恒星三个前沿悬念的卡点。
易错点与常见误读
积分步长不均匀会把陡峭的核心梯度抹平，中心压强偏低数倍；把拟合关系当定律外推，会在极端质量处全军覆没；专业代码不是黑箱，输入文件每一段都对应一条方程，读不懂输入就别信输出。此外要区分"结构模型"（某一时刻的快照）与"演化模型"（时间序列的连拍），本章第一节算的是快照，第三节工具链拍的是连拍。
与前后章的衔接
本章把全册的量级演算集中返厂：第一章的时标、第三章的寿命、第五章的能量账，都在这里获得代码版本。演算全程只用桌面算力，但方法论与专业科研同构——先量级、后拟合、再全模型，这三步是任何计算物理课题的通用节拍。
三节作业的递进关系
本章三步走不是并列的三个专题，而是一条能力递进线：9.1 练的是"从方程到代码"（复现量级），9.2 练的是"从数据到结论"（拟合与闭环），9.3 练的是"从玩具到工业"（工具链与版本感）。三节作业建议按序完成，每节的花费都在一个下午以内，但衔接处的检查点要守住：9.1 的积分不收敛，别急着进 9.2；9.2 的残差没归因，9.3 的参数敏感性实验会做不明白。递进线走完，读者的收获不是三个知识点，而是一套可以迁移到任何计算课题的工作流。
演算环境的约定
为避免环境差异淹没物理结论，本章统一约定：代码用任何带浮点运算的桌面环境即可复跑，示例用常见的科学计算风格书写但只依赖标准库；输出块是真实运行结果，读者复跑的偏差应只在末位有效数字。单位纪律单列一条：全程国际单位制或在行内显式标注太阳单位，混用单位是数值实验第一大事故源。这些约定看起来琐碎，却是"可复现"三个字的地基——第九章的全部教学价值都压在这块地基上。