Chemistry Labs

核化学与放射化学

放射性同位素测年

放射性衰变常数与母体-子体同位素比如何构建地球年代学时钟,涵盖从考古碳-14 测年到等时线及地球年龄测定的广阔领域。

直觉直觉:封存在岩石中的放射性原子沙漏

沙漏之所以能够精确计时,是因为细沙以恒定速率由上部容器落入下部容器。在自然界中,放射性原子核构成了天然的微观原子沙漏:不稳定的母核依据严格的核量子跃迁概率衰变为放射成因的子核,这一过程完全不受外界温度、地层压力或化学键环境的影响。当岩浆冷却结晶为固体矿物或生物生命活动停止呼吸代谢的瞬间,沙漏被正式倒置启动:母体原子被牢固锁死在晶格节点中,而子体产物的持续累积便忠实记录下地质历史的浩瀚时间跨度。

选择地质年代计以观察其衰变轨迹:5 万年尺度内 14C 剩余份额的指数衰减(标有 5.73、11.5 和 50 千年标记)、直至 45.4 亿年地球年龄的放射成因 40Ar/40K 累积比值,以及追溯数十亿年行星演化历史的 207Pb/206Pb* 比值曲线。

中学阶段中学阶段:基本年代方程与碳-14 测年法

定义: 基本放射性年代方程与封闭体系假设

对于衰变常数为 λ\lambda 的放射性母核 PP,其衰变积累的放射成因子核数目为 D∗=P(eλt−1)D^* = P(e^{\lambda t} - 1)。解出结晶年龄 tt 得到 t=1λln⁡(1+D∗P)t = \frac{1}{\lambda}\ln\left(1 + \frac{D^*}{P}\right)。该方程所求得的年龄要代表真实的地质年代,前提是该矿物自冷却至其特有的封闭温度(TcT_c)以下以来,晶格系统必须保持严格的封闭体系(既无母体或子体元素的流失,也无外界原子的侵入)。

放射性碳测年依托于高层大气中的宇宙射线核反应:次级热中子轰击氮原子核产生碳-14:X714X27214N+X01X2021n→X614X26214C+X11X2121p\ce{^{14}_7N + ^1_0n -> ^{14}_6C + ^1_1p}。新生的碳-14 迅速被氧化为 X14X2214COX2\ce{^{14}CO2} 并通过植物光合作用进入全球食物网,使活体生物内的 X14X2214C/X12X2212C\ce{^{14}C/^{12}C} 原子比恒定维持在约 1.2×10−121.2 \times 10^{-12} 的稳态水平。生物死亡后新陈代谢立即终止,体内的 X14X2214C\ce{^{14}C} 便以 t1/2=5730±40 年t_{1/2} = 5730\pm 40\text{ 年}(剑桥半衰期)通过 β−\beta^- 衰变自发恢复为氮-14。按照国际地质年代学公约,未经校正的常规年龄采用利比半衰期(5568 年5568\text{ 年})并以“距今放射性碳年龄 BP”(以 1950 年为基准点)表示。鉴于古地磁场波动和太阳耀斑活动,必须应用基于树木年轮与湖泊纹泥建立的 IntCal20 树轮校准曲线转换为真实日历年。现代加速器质谱技术(AMS)能直接在毫克级微量样品上计数碳-14 粒子,使测年极限拓展至约 5 万年。

主要地质年代测定同位素体系与目标矿物
衰变体系衰变模式母核半衰期有效测年范围典型分析矿物与材料
X14X2214C→X14X2214N\ce{^{14}C -> ^{14}N}Beta- 负衰变5730 年300 至 50,000 年木炭、木材、古骨胶原、贝壳、泥炭
X40X2240K→X40X2240Ar\ce{^{40}K -> ^{40}Ar}轨道电子俘获(10.7%)12.48 亿年10 万年至 45 亿年透长石、黑云母、白云母、角闪石
X87X2287Rb→X87X2287Sr\ce{^{87}Rb->^{87}Sr}Beta- 负衰变488 亿年1000 万年至 46 亿年云母、钾长石、花岗岩类全岩等时线
X238X22238U→X206X22206Pb\ce{^{238}U->^{206}Pb}Alpha 与 Beta 级联衰变链44.68 亿年100 万年至 46 亿年锆石、斜锆石、独居石、榍石
X235X22235U→X207X22207Pb\ce{^{235}U->^{207}Pb}Alpha 与 Beta 级联衰变链7.038 亿年100 万年至 46 亿年与 238U 配对构建谐和曲线与铅铅等时线
X147X22147Sm→X143X22143Nd\ce{^{147}Sm->^{143}Nd}Alpha 衰变1060 亿年1 亿年至 46 亿年基性玄武岩、月球斜长岩、球粒陨石

例题: 考古木炭的放射性常规碳-14 年龄计算

某考古发掘遗址出土的一根古雪松木横梁,其实测比放射性活度恰好为工业革命前现代木材活度标准的 12.5%12.5\%(即 1/81/8)。试使用标准的剑桥物理半衰期 t1/2=5730 年t_{1/2} = 5730\text{ 年},计算该木制文物在进行树木年轮曲线校准前的常规同位素测年年龄。

解答

放射性剩余活度分数刚好为 N/N0=0.125=(1/2)3N/N_0 = 0.125 = (1/2)^3,表明母核经历了整整 3 个半衰期的衰变过程:t=3×t1/2=3×5730=17190 年t = 3 \times t_{1/2} = 3 \times 5730 = 17190\text{ 年}。代入指数衰变方程严格计算:t=−t1/2ln⁡2ln⁡(NN0)=−57300.69315×(−2.07944)≈17190 年t = -\frac{t_{1/2}}{\ln 2}\ln\left(\frac{N}{N_0}\right) = -\frac{5730}{0.69315}\times (-2.07944) \approx 17190\text{ 年}。在考古学实践中,这一未经校正的常规放射性碳年龄需进一步通过 IntCal20 曲线校准,以消除历史大气浓度波动,最终换算为日历公历年份(cal BP)。

大学大学阶段:K-Ar 分支衰变、40Ar/39Ar 阶梯加热与 Rb-Sr 等时线

t=1λln⁡(1+λλe40Ar∗40K)withλ=5.543×10−10 y−1, λe=0.580×10−10 y−1t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{^{40}\text{Ar}^*}{^{40}\text{K}}\right) \quad \text{with} \quad \lambda = 5.543 \times 10^{-10}\ \text{y}^{-1},\ \lambda_e = 0.580 \times 10^{-10}\ \text{y}^{-1}

钾-40 发生双通道分支衰变:轨道电子俘获衰变为惰性气体 X40X2240Ar\ce{^{40}Ar}(分支比 10.72%10.72\%,λe=0.580×10−10 年−1\lambda_e = 0.580 \times 10^{-10}\ \text{年}^{-1}),以及 β−\beta^- 负衰变为 X40X2240Ca\ce{^{40}Ca}(分支比 89.28%89.28\%,λβ=4.963×10−10 年−1\lambda_\beta = 4.963 \times 10^{-10}\ \text{年}^{-1}),总衰变常数 λ=5.543×10−10 年−1\lambda = 5.543 \times 10^{-10}\ \text{年}^{-1}(半衰期 12.48 亿年)。在地下炽热岩浆中,挥发性极强的氩气完全逸出脱气,使矿物初始氩含量清零(X40X2240Ar0≈0\ce{^{40}Ar}_0 \approx 0)。一旦岩浆喷发冷却为透长石或黑云母晶体,放射成因的 X40X2240Ar∗\ce{^{40}Ar}^* 便被封闭在矿物晶格中。现代 X40X2240Ar/X39X2239Ar\ce{^{40}Ar/^{39}Ar} 测年技术先将样品送入反应堆快中子轰击,通过 X39X2239K(n,p)X39X2239Ar\ce{^{39}K}(n,p)\ce{^{39}Ar} 反应建立同位素关联,再在惰性气体质谱计中进行激光阶梯升温释放分析,由平坦的年龄坪谱能够准确剔除后期变质加热导致的放射氩流失或过剩氩干扰。

87Sr86Sr=(87Sr86Sr)0+87Rb86Sr(eλt−1)\frac{^{87}\text{Sr}}{^{86}\text{Sr}} = \left(\frac{^{87}\text{Sr}}{^{86}\text{Sr}}\right)_0 + \frac{^{87}\text{Rb}}{^{86}\text{Sr}}\left(e^{\lambda t} - 1\right)

当矿物结晶时含有非零的初始子体同位素 D0D_0 时,单一样品的放射性测量在数学上无法求解年龄。等时线法(Isochron method)通过引入子体元素的非放射成因稳定同位素进行比值归一化,圆满解决了这一难题。例如在铷-锶测年体系(X87X2287Rb→X87X2287Sr\ce{^{87}Rb -> ^{87}Sr},t1/2=48.8 亿年t_{1/2} = 48.8\text{ 亿年},λ=1.42×10−11 年−1\lambda = 1.42 \times 10^{-11}\ \text{年}^{-1})中以稳定的 X86X2286Sr\ce{^{86}Sr} 作为分母。同一岩浆房同源结晶的多种矿物,由于结晶化学分异导致元素比 Rb/Sr\ce{Rb/Sr} 各不相同,但在固化瞬间共享完全均一的初始锶同位素比 (X87X2287Sr/X86X2286Sr)0(\ce{^{87}Sr/^{86}Sr})_0。以 X87X2287Rb/X86X2286Sr\ce{^{87}Rb/^{86}Sr} 为横坐标、X87X2287Sr/X86X2286Sr\ce{^{87}Sr/^{86}Sr} 为纵坐标作图,数据点连成一条笔直的等时线,其斜率 m=eλt−1≈λtm = e^{\lambda t} - 1 \approx \lambda t 直接决定了矿物的结晶年龄 tt,完全无需事先假设初始子体丰度。

例题: 火山凝灰岩透长石晶体的钾-氩同位素年代测定

对火山凝灰岩夹层中分选出的透长石单晶进行高精度同位素质谱分析,测得放射成因氩与母体钾的原子数比为 X40X2240Ar∗X40X2240K=0.0950\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}} = 0.0950。已知总衰变常数 λ=5.543×10−10 年−1\lambda = 5.543 \times 10^{-10}\ \text{年}^{-1},电子俘获分支常数 λe=0.580×10−10 年−1\lambda_e = 0.580 \times 10^{-10}\ \text{年}^{-1},试计算该火山岩层的结晶地质年龄(以百万年 Ma 为单位)。

解答

代入 K-Ar 分支衰变测年公式:t=1λln⁡(1+λλeX40X2240Ar∗X40X2240K)t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}}\right)。首先计算分支比修正系数:λλe=5.543×10−100.580×10−10≈9.5569\frac{\lambda}{\lambda_e} = \frac{5.543 \times 10^{-10}}{0.580 \times 10^{-10}} \approx 9.5569。代入对数项真数:1+9.5569×0.0950=1+0.9079=1.90791 + 9.5569 \times 0.0950 = 1 + 0.9079 = 1.9079。计算自然对数值:ln⁡(1.9079)≈0.6460\ln(1.9079) \approx 0.6460。除以总衰变常数 λ\lambda:t=0.64605.543×10−10 年−1≈1.1654×109 年=1165.4 Mat = \frac{0.6460}{5.543 \times 10^{-10}\ \text{年}^{-1}} \approx 1.1654 \times 10^9\text{ 年} = 1165.4\text{ Ma}(即 11.654 亿年)。

进阶高级阶段:U-Pb 谐和线、Jack Hills 锆石与地球的年龄

锆石(ZrSiOX4\ce{ZrSiO4})被公认为同位素地质年代学的金标准。在岩浆结晶过程中,锆石四方岛状硅酸盐晶格中的四价锆位点(r=84 pmr = 84\text{ pm})容许电荷与半径相近的 UX4+\ce{U^4+}(r=100 pmr = 100\text{ pm})以类质同象方式大量赋存,而以超过 10510^5 的排斥倍数彻底排斥半径较大且为二价的 PbX2+\ce{Pb^2+}(r=119 pmr = 119\text{ pm}),从而使新形成锆石中的初始普通铅几乎为零(Pb0≈0\ce{Pb}_0 \approx 0)。锆石内部同时运行着半衰期迥异的两个独立铀衰变时钟:X238X22238U→X206X22206Pb\ce{^{238}U -> ^{206}Pb}(t1/2=44.68 亿年t_{1/2} = 44.68\text{ 亿年},λ238=1.55125×10−10 年−1\lambda_{238} = 1.55125 \times 10^{-10}\ \text{年}^{-1})和 X235X22235U→X207X22207Pb\ce{^{235}U -> ^{207}Pb}(t1/2=7.038 亿年t_{1/2} = 7.038\text{ 亿年},λ235=9.8485×10−10 年−1\lambda_{235} = 9.8485 \times 10^{-10}\ \text{年}^{-1})。以 X207X22207Pb∗/X235X22235U\ce{^{207}Pb}^*/\ce{^{235}U} 为横轴、X206X22206Pb∗/X238X22238U\ce{^{206}Pb}^*/\ce{^{238}U} 为纵轴绘制的参数轨迹称为韦瑟里尔(Wetherill)谐和曲线。保持封闭体系的样品数据点严密落在谐和曲线上;若岩石在后期热事件中发生放射铅的扩散流失,测点将构成一条笔直的不谐和线(Discordia),其与谐和线的上交点代表岩浆初生结晶年龄,下交点则代表变质扰动年代。采自西澳大利亚杰克山(Jack Hills)的碎屑锆石测得了高达 44.04 亿年44.04\text{ 亿年} 的谐和年龄,成为地球冥古宙原始地壳存在液态水与长英质陆壳的直接见证。

对两个铀-铅衰变方程两端求比值,便巧妙消除了对母体元素铀绝对化学浓度的测量依赖:X207X22207Pb∗X206X22206Pb∗=X235X22235UX238X22238U(eλ235t−1eλ238t−1)=1137.88(eλ235t−1eλ238t−1)\frac{\ce{^{207}Pb}^*}{\ce{^{206}Pb}^*} = \frac{\ce{^{235}U}}{\ce{^{238}U}}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right) = \frac{1}{137.88}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right),其中 137.88137.88 为现代地球与太阳系的天然 X238X22238U/X235X22235U\ce{^{238}U/^{235}U} 同位素比值。1956 年,克莱尔·帕特森(Clair Patterson)测定了不含铀的峡谷暗黑(Canyon Diablo)铁陨石陨硫铁相(FeS\ce{FeS})中的铅同位素组成,确立了太阳系初生普通铅的原始基准。通过将多种石陨石与现代大洋深海沉积物铅同位素联合投在一条全球共生的铅铅等时线上,帕特森最终测定地球的年龄为 45.5±0.7 亿年45.5 \pm 0.7\text{ 亿年}(即 45.4 亿年45.4\text{ 亿年})。现代对原始球粒陨石中富钙铝包裹体(CAI)的超高精度铀-铅测年,将太阳星云最早固相物质凝聚时间精确锚定在 45.67 亿年45.67\text{ 亿年} 前。

参考文献

  • Atmospheric Helium Three and Radiocarbon from Cosmic Radiation · W. F. Libby, 1946
  • Precision Measurement of Half-Lives and Specific Activities of 235U and 238U · A. H. Jaffey, K. F. Flynn, L. E. Glendenin, W. C. Bentley, A. M. Essling, 1971
  • The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0-55 cal kBP) · P. J. Reimer et al., 2020