Chemistry Labs

12年级

半衰期与放射性衰变定律

原子核存续的数学规律:指数衰变定律、贝克勒尔活度表征、母子体衰变贝特曼方程组以及锝医用发生器的暂时平衡动力学。

直觉直觉:由纯粹概率构建的时间之钟

想象每分钟同时抛掷一万枚硬币,并将正面朝上的硬币全部移走。在第一分钟内,大约一半的硬币消失,剩下约五千枚;在第二分钟内,所剩硬币又减少一半,剩下约两千五百枚。没有任何一枚硬币能够预知自己何时掷出正面,一枚已经存在了十分钟的硬币在下一分钟翻面的概率与一枚新铸造的硬币毫无二致。放射性原子核的行为完全遵从这种无记忆性的统计法则:每个不稳定原子核在任意一秒内发生衰变的概率恒定不变,与其已经存续了多久完全无关。

由于宏观物质包含极其庞大的原子核系统(仅一微克同位素就包含数万亿个原子),这种微观的纯随机性经统计平均后,汇聚成一条极其平滑、完全可预测的指数衰变曲线。存活的原子核数目在每一个被称为半衰期的固定时间间隔内均减半。无论是追踪几小时内衰减的医用显像示踪剂,还是测定历经数十亿年演化的铀矿地质年龄,指数衰变定律都为科学探索提供了恒定精准的原子级计时器。

交互式衰变谱图。可分别观察通用的指数半衰减曲线(数据集 0)、在约 23 小时出现子体峰值的 ⁹⁹Mo/⁹⁹ᵐTc 医用发生器母子体动力学曲线(数据集 1),以及标注有数千年衰变时间节点的放射性碳衰变曲线(数据集 2)。

中学阶段基础核心:半衰期、衰变常数与放射性活度

定义: 放射性衰变定律与放射性活度

放射性核素的衰变率 −dNdt-\frac{\mathrm{d}N}{\mathrm{d}t} 正比于该时刻存在的放射性核数 N(t)N(t):−dNdt=λN-\frac{\mathrm{d}N}{\mathrm{d}t} = \lambda N,其中 λ\lambda 为特征衰变常数(s−1\text{s}^{-1})。积分后可得 N(t)=N0e−λt=N02−t/t1/2N(t) = N_0 \mathrm{e}^{-\lambda t} = N_0 2^{-t / t_{1/2}}。半衰期 t1/2=ln⁡2λ≈0.69315λt_{1/2} = \frac{\ln 2}{\lambda} \approx \frac{0.69315}{\lambda} 为初始核数目减半所需的时间。放射性活度 A(t)=λN(t)=A0e−λtA(t) = \lambda N(t) = A_0 \mathrm{e}^{-\lambda t} 度量单位时间内的衰变次数,单位为贝克勒尔(1 Bq=1 次衰变/秒1\text{ Bq} = 1\text{ 次衰变/秒})或居里(1 Ci=3.7×1010 Bq1\text{ Ci} = 3.7 \times 10^{10}\text{ Bq})。

N(t)=N0e−λt=N0 2−t/t1/2,A(t)=λN(t),t1/2=ln⁡2λ,τ=1λ≈1.4427 t1/2N(t) = N_0 \mathrm{e}^{-\lambda t} = N_0 \, 2^{-t / t_{1/2}}, \qquad A(t) = \lambda N(t), \qquad t_{1/2} = \frac{\ln 2}{\lambda}, \qquad \tau = \frac{1}{\lambda} \approx 1.4427 \, t_{1/2}
代表性放射性核素及其半衰期
核素半衰期主要衰变方式典型应用领域
X614X26214C\ce{^{14}_6C}5730 年β−\beta^-有机考古文物的碳-14 年代测定
X53131X2532131I\ce{^{131}_{53}I}8.02 天β−,γ\beta^-, \gamma甲状腺疾病的靶向放射治疗与显像诊断
X4399mX243299mTc\ce{^{99m}_{43}Tc}6.01 小时同质异能跃迁 (γ,140.5 keV\gamma, 140.5\text{ keV})临床 SPECT 诊断医学显像的标准核素
X2760X227260Co\ce{^{60}_{27}Co}5.27 年β−,γ\beta^-, \gamma工业伽马射线辐照灭菌与外照射放射治疗
X55137X2552137Cs\ce{^{137}_{55}Cs}30.1 年β−,γ\beta^-, \gamma裂变产物环境示踪与辐射测量仪器校准源
X88226X2882226Ra\ce{^{226}_{88}Ra}约 1600 年α,γ\alpha, \gamma历史居里活度单位基准源及早期近距离治疗
X92238X2922238U\ce{^{238}_{92}U}44.7 亿年α\alpha原始地球年代学测定基准与核反应堆燃料

例题: 放射性碳定年与残存份额计算

从某考古遗址出土的木质文物测得其 14C^{14}\text{C} 活度仅为现代活木的 25%25\%(14C^{14}\text{C} 的半衰期 t1/2=5730 年t_{1/2} = 5730\text{ 年})。求已过去的半衰期次数并推算文物的年代。此外,计算活度衰减至初始值的 0.1%0.1\%(1/10001/1000)以下所需经历的半衰期次数。

解答

残存份额:N/N0=0.25=(1/2)2N/N_0 = 0.25 = (1/2)^2,故历经的半衰期次数为 n=2n = 2 次。 文物年代:t=2×5730=11460 年t = 2 \times 5730 = 11460\text{ 年}(约 1.15 万年)。 活度降至 1/10001/1000 所需次数:(1/2)n=10−3  ⟹  n=ln⁡1000ln⁡2≈6.90780.69315≈9.97 个半衰期(1/2)^n = 10^{-3} \implies n = \frac{\ln 1000}{\ln 2} \approx \frac{6.9078}{0.69315} \approx 9.97\text{ 个半衰期}(约为 10 个半衰期,约合 5.71 万年)。

例题: 微克级碘-131 样品的放射性活度计算

某无载体放射性药物安瓿中含有 1.0 μg1.0\ \mu\text{g} 纯碘-131(X53131X2532131I\ce{^{131}_{53}I},摩尔质量 M≈130.9 g/molM \approx 130.9\text{ g/mol},半衰期 t1/2=8.02 天t_{1/2} = 8.02\text{ 天})。求其衰变常数 λ\lambda(s−1\text{s}^{-1})、放射性核总数 NN 以及初始放射性活度(以 Bq 和 GBq 表示)。

解答

以秒为单位的半衰期:t1/2=8.02×86400 s≈6.929×105 st_{1/2} = 8.02 \times 86400\text{ s} \approx 6.929 \times 10^5\text{ s}。 衰变常数:λ=ln⁡2t1/2=0.693156.929×105≈1.000×10−6 s−1\lambda = \frac{\ln 2}{t_{1/2}} = \frac{0.69315}{6.929 \times 10^5} \approx 1.000 \times 10^{-6}\text{ s}^{-1}。 放射性核总数:N=1.0×10−6 g130.9 g/mol×6.022×1023 mol−1≈4.60×1015 个N = \frac{1.0 \times 10^{-6}\text{ g}}{130.9\text{ g/mol}} \times 6.022 \times 10^{23}\text{ mol}^{-1} \approx 4.60 \times 10^{15}\text{ 个}。 放射性活度:A=λN=1.000×10−6 s−1×4.60×1015≈4.60×109 Bq=4.60 GBqA = \lambda N = 1.000 \times 10^{-6}\text{ s}^{-1} \times 4.60 \times 10^{15} \approx 4.60 \times 10^9\text{ Bq} = 4.60\text{ GBq}(约合 124 mCi124\text{ mCi})。

大学大学阶段:贝特曼方程组与放射平衡

dN2dt=λ1N1−λ2N2  ⟹  N2(t)=λ1λ2−λ1N1,0(e−λ1t−e−λ2t),tmax⁡=ln⁡(λ2/λ1)λ2−λ1\frac{\mathrm{d}N_2}{\mathrm{d}t} = \lambda_1 N_1 - \lambda_2 N_2 \implies N_2(t) = \frac{\lambda_1}{\lambda_2 - \lambda_1} N_{1,0} \left( \mathrm{e}^{-\lambda_1 t} - \mathrm{e}^{-\lambda_2 t} \right), \qquad t_{\max} = \frac{\ln(\lambda_2 / \lambda_1)}{\lambda_2 - \lambda_1}

当放射性母核 1 衰变生成放射性子核 2 时,子核的瞬时数量遵循哈里·贝特曼联立微分方程组。根据两者的半衰期比例,主要形成两种平衡状态: 1. 长期平衡(永久平衡):母核半衰期远大于子核半衰期(λ1≪λ2\lambda_1 \ll \lambda_2,例如 X226X22226Ra\ce{^{226}Ra} 的 t1/2=1600 年t_{1/2} = 1600\text{ 年},而子核 X222X22222Rn\ce{^{222}Rn} 仅为 3.82 天3.82\text{ 天})。经子核的数个半衰期后,e−λ2t→0\mathrm{e}^{-\lambda_2 t} \to 0,子核放射性活度与母核完全相等:A2(t)→A1(t)A_2(t) \to A_1(t)。 2. 暂时平衡:母核寿命大于子核寿命但处于可比时间尺度(λ1<λ2\lambda_1 < \lambda_2,典型代表为医用 X99X2299Mo→X99mX2299mTc\ce{^{99}Mo} \to \ce{^{99m}Tc} 发生器,t1/2t_{1/2} 分别为 65.94 小时65.94\text{ 小时} 与 6.01 小时6.01\text{ 小时})。子核活度会逐渐超过母核,并在 tmax⁡=ln⁡(λ2/λ1)λ2−λ1≈22.8 小时t_{\max} = \frac{\ln(\lambda_2 / \lambda_1)}{\lambda_2 - \lambda_1} \approx 22.8\text{ 小时}(约 23 小时)达到峰值,此后两者以母核的半衰期同步衰减:A2(t)→λ2λ2−λ1A1(t)≈1.10 A1(t)A_2(t) \to \frac{\lambda_2}{\lambda_2 - \lambda_1} A_1(t) \approx 1.10 \, A_1(t)。

在放射性药物学与环境生物物理学中,引入生物体内的放射性核素同时经历物理衰变与生物代谢清除两种独立途径。有效衰变常数等于物理衰变率与生物清除率的直接代数和(λeff=λphys+λbio\lambda_{\text{eff}} = \lambda_{\text{phys}} + \lambda_{\text{bio}}),由此导出有效半衰期的倒数调和关系:1Teff=1Tphys+1Tbio\frac{1}{T_{\text{eff}}} = \frac{1}{T_{\text{phys}}} + \frac{1}{T_{\text{bio}}}。若某个核素存在多种相互竞争的衰变途径(分枝衰变,分枝常数分别为 λa\lambda_a 和 λb\lambda_b),则总衰变常数为 λ=λa+λb\lambda = \lambda_a + \lambda_b,分支 aa 的分支比即为 fa=λa/λf_a = \lambda_a / \lambda。

进阶高阶物理深度:泊松统计与量子隧穿

由于单个放射性衰变是在庞大核子总体中独立发生的泊松事件,测得的 NN 次衰变计数必然带有不可消除的统计标准差 σ=N\sigma = \sqrt{N} 以及相对不确定度 σ/N=1/N\sigma / N = 1/\sqrt{N}。在低本底辐射测量或旧石器时代碳-14 年代测定中,累积足够多的有效计数是压低统计噪声的唯一物理途径。此外,核素的衰变常数在同位素图谱上跨越了二十个数量级以上。1928 年,乔治·伽莫夫(以及格尼与康登)揭示了 α 衰变半衰期源于量子力学势垒隧穿效应:核内预先形成的 α 粒子以约每秒 102110^{21} 次的高频撞击库仑势垒,指数隧穿几率 P∝exp⁡[−2πη]P \propto \exp[-2\pi\eta] 完美解释了盖革-努塔尔经验定律,即 α 粒子动能的微小增加会引发半衰期跨越十几个数量级的剧烈骤降。

参考文献

  • Introductory Nuclear Physics · K. S. Krane, 1987
  • Radiochemistry and Nuclear Chemistry, 4th Edition · G. Choppin, J.-O. Liljenzin, J. Rydberg, C. Ekberg, 2013
  • Atmospheric Helium Three and Radiocarbon from Cosmic Radiation · W. F. Libby, 1946