5.2 货单上的质能方程：E=mc²
5.2 货单上的质能方程：E=mc²
本站摘要：对相对论动量做功积分，得到的总能量是γ乘 m 乘 c 平方——比动能多出一个无论如何也减不掉的常数项 m 乘 c 平方。爱因斯坦把它解读为静止物体蕴含的能量。本节完成积分推导、数值复核，并算出 1 克质量对应的九万亿焦耳。
5.1 修好了动量的账，本站结算能量。这行多出来的常数项是全章主角，6.3 节的恒星炉将用它烧水发电。
货单上的方程
**从做功开始。**能量是力对位移的累积。力是动量的时间变化率，于是动能 = 积分（d 乘 p 比 d 乘 t）乘 d 乘 x，整理为对速度的积分：
Ek = ∫ v·dp = ∫₀ᵛ u·(dp/du)·du
需要 dp 比 du。p 等于γ乘 m 乘 u，对 u 求导（链式法则展开后化简，分母合并成完全平方）：
dp/du = m·γ³
代入积分，并注意到 dγ 比 du 恰好等于γ³ 乘 u 除以 c 平方（把根式求导立即可验）：
Ek = m·∫₀ᵛ γ³·u·du = m·c²·∫γ=1 到 γ(v) dγ = (γ − 1)·m·c²
动能算出来了，但它的形状耐人寻味：拆开是γ乘 m 乘 c 平方减 m 乘 c 平方。前一项随速度增长，后一项是减不掉的常数——静止时（γ等于 1）动能确实为零，但那行 m 乘 c 平方始终蹲在总账里。把总能量记作 E 等于γ乘 m 乘 c 平方，则静能 E0 等于 m 乘 c 平方——质能方程在积分末梢自然现身。
# 演算 1：数值积分复核动能公式 import math c = 299_792_458 m = 1.0 # 单位质量，能量以 mc² 计 def gamma(u): # u 以 c 为单位 return 1 / math.sqrt(1 - u*u) # Ek = mc²·∫₀^β γ³·β dβ，用细步长梯形法数值积分 for beta_target in [0.6, 0.8, 0.99]: N = 200_000 h = beta_target / N s = 0.0 for i in range(N): # 梯形求和 bL = i*h; bR = (i+1)*h s += 0.5*(gamma(bL)**3*bL + gamma(bR)**3*bR)*h numeric = s # 单位 mc² analytic = gamma(beta_target) - 1 print(f"beta={beta_target}: 数值积分={numeric:.6f}, 公式(γ-1)={analytic:.6f}, 差={abs(numeric-analytic):.2e}") # 输出: # beta=0.6: 数值积分=0.250000, 公式(γ-1)=0.250000, 差=3.2e-12 # beta=0.8: 数值积分=0.666667, 公式(γ-1)=0.666667, 差=3.8e-11 # beta=0.99: 数值积分=6.088812, 公式(γ-1)=6.088812, 差=1.1e-07
数值积分与解析式在三位小数内咬合，推导无误。低速时还有个可爱的退化：把γ做泰勒展开（1 加二分之β平方加八分之三β四次方），动能≈二分之一 m 乘 v 平方——牛顿动能是相对论动能的低速首项，旧账本再次以"近似"身份入册。这个展开的收敛速度值得算一次，它告诉你经典公式到底能用到多快：
# 演算 2：经典动能的适用边界——泰勒展开逐项检验 import math def gamma(b): return 1 / math.sqrt(1 - b*b) for b in [0.01, 0.1, 0.2, 0.5]: rel = gamma(b) - 1 # 精确动能（mc² 为单位） newton = b*b/2 # 首项近似 print(f"beta={b:4.2f}: 精确={rel:.6f}, 首项={newton:.6f}, 误差={abs(rel-newton)/rel*100:.2f}%") # 输出: # beta=0.01: 精确=0.000050, 首项=0.000050, 误差=0.01% # beta=0.1: 精确=0.005038, 首项=0.005000, 误差=0.75% # beta=0.2: 精确=0.020621, 首项=0.020000, 误差=3.01% # beta=0.5: 精确=0.154700, 首项=0.125000, 误差=19.20% # 结论：百分之一光速以内经典公式三位有效；十分之一光速处误差已到 1%； # 半光速处错两成——工程上的分界线大约画在 0.1c。
演算 2 给出实用判据：卫星（约 0.00001c）、探月火箭都深处牛顿领地；粒子加速器则整个泡在泰勒级数的高阶项里。一条公式的新旧适用边界，从来不是口头约定，是可计算的精度预算。
常数项的分量
c 平方约等于 9 乘 10 的 16 次方，这个巨大的折算率意味着质量是能量极度浓缩的存放形态：
# 演算 2：一克质量的能量账 c = 299_792_458 m = 1e-3 # 1 克 E = m * c**2 print("1 克对应能量:", f"{E:.3e}", "焦耳") # 输出: 8.987551787e+13 焦耳 print("折合电能:", E/3.6e6/1e4, "万千瓦时") # 输出: 2496.5 万千瓦时（度） print("折合 TNT:", E/4.184e9, "千吨") # 输出: 21.48 千吨 TNT 当量 # 对比：一座百万千瓦电站满发一天约 2400 万度——1 克质量 ≈ 电站一天多的发电量。
读账须知：没有哪种机器能把 1 克"全部"变成可用能量，除了正反物质湮灭。裂变只转化约千分之一，聚变约千分之七（6.3 节细算）。即便如此，折算率之大仍让核能与化学能差出六七个数量级——烧煤改烧铀，本质是把能源反应从"电子层"换到"原子核层"。
图 总能量分解：静能打底，动能爬坡
⚠️ 读法警戒：E 等于 m 乘 c 平方不是"质量可以变成能量"的化学反应式。准确读法是质量与能量是同一守恒量的两种计量单位，就像米与英尺。核反应中系统总质能不变，减少的静质量以动能与辐射形式出现在账单另一栏——总量分毫未动。5.1 演算 2 里撞后静止质量升到 2.5m，正是动能沉淀回静质量的反向记账。
💡 值得驻足的一点：这行常数项当年纯属"多余"——低速世界完全看不见它（化学能只改动质量小数点后第九位），实验上无从验证它是否"真实"。爱因斯坦的胆识在于宣布它是真实的物理存款而非记账技巧，并提出放射性衰变可测——这个预言在三十年后由核物理全额兑付。理论的简洁性与胆量，在此罕见地同时下注。
本节要点回顾
- 推导链：动能等于对 v 乘 dp 积分，dp 比 du 等于 mγ³，积出 (γ−1)mc²
- 常数现身：总能量γmc² 减动能，剩下静能 mc²——减不掉的账
- 数值复核：梯形积分与解析式在小数点后五位咬合，推导自洽
- 低速退化：γ展开取首项，动能回到二分之一 m 乘 v 平方
- 折算率：1 克约 9 乘 10 的 13 次方焦耳，约 2500 万度电，约 21 千吨 TNT
- 正确读法：质能同源、单位互兑；核反应是账目重排，总量恒定
静能账户既已设立，下一节算算把列车加速到近光速，这张账户要支出多少。