Chemistry Labs

Hoá học hạt nhân và phóng xạ

Xác định niên đại bằng đồng vị phóng xạ

Cách hằng số phân rã phóng xạ và tỉ số đồng vị mẹ-con tạo thành đồng hồ địa thời học, từ cacbon phóng xạ trong khảo cổ đến các đường đẳng thời và tuổi của Trái Đất.

Trực giácTrực giác: những đồng hồ cát phóng xạ đóng băng trong đá

Đồng hồ cát đo thời gian trôi qua nhờ những hạt cát chảy từ bầu trên xuống bầu dưới với tốc độ không đổi. Trong tự nhiên, các hạt nhân phóng xạ đóng vai trò như những đồng hồ cát nguyên tử: hạt nhân mẹ không bền chuyển hoá thành hạt nhân con theo xác suất lượng tử hạt nhân, hoàn toàn không bị ảnh hưởng bởi nhiệt độ môi trường, áp suất chôn vùi hay liên kết hoá học. Khi magma nguội thành tinh thể khoáng vật rắn hoặc sinh vật ngừng trao đổi chất, đồng hồ cát bắt đầu được lật: các nguyên tử mẹ bị khoá chặt trong mạng tinh thể và sự tích luỹ đều đặn của hạt nhân con ghi lại thời gian địa chất sâu thẳm.

Chọn một đồng hồ địa thời để theo dõi quỹ đạo phân rã: phần 14C còn lại theo hàm mũ qua 50 nghìn năm (với các mốc 5,73, 11,5 và 50 nghìn năm), sự tích luỹ 40Ar/40K cho đến tuổi 4,54 tỷ năm của Trái Đất, và tỉ số 207Pb/206Pb* ghi dấu tiến hoá hành tinh qua hàng tỷ năm.

Phổ thôngMức phổ thông: phương trình tuổi cơ bản và xác định niên đại bằng cacbon-14

Định nghĩa: Phương trình tuổi phóng xạ cơ bản và hệ kín

Với đồng vị mẹ phóng xạ PP phân rã thành đồng vị con D∗D^* với hằng số phân rã λ\lambda, số nguyên tử con tích luỹ là D∗=P(eλt−1)D^* = P(e^{\lambda t} - 1). Giải tìm tuổi tinh thể tt ta được t=1λln⁡(1+D∗P)t = \frac{1}{\lambda}\ln\left(1 + \frac{D^*}{P}\right). Để phương trình này phản ánh đúng tuổi địa chất thực tế, mẫu khoáng vật phải duy trì một hệ kín tuyệt đối (không mất mát hay nhận thêm nguyên tử mẹ hoặc con) kể từ khi nguội xuống dưới nhiệt độ đóng (TcT_c) đặc trưng của nó.

Định tuổi bằng cacbon phóng xạ dựa vào phản ứng va chạm tia vũ trụ ở thượng tầng khí quyển, nơi các nơtron nhiệt bắn phá nitơ: X714X27214N+X01X2021n→X614X26214C+X11X2121p\ce{^{14}_7N + ^1_0n -> ^{14}_6C + ^1_1p}. Lượng cacbon-14 sinh ra bị oxy hoá thành X14X2214COX2\ce{^{14}CO2} và đi vào sinh quyển qua quang hợp, duy trì tỉ số cân bằng động X14X2214C/X12X2212C≈1,2×10−12\ce{^{14}C/^{12}C} \approx 1{,}2 \times 10^{-12} trong mô sống. Khi sinh vật chết đi, sự trao đổi chất dừng lại và X14X2214C\ce{^{14}C} phân rã β−\beta^- trở lại thành X14X2214N\ce{^{14}N} với chu kỳ bán rã t1/2=5730±40 na˘mt_{1/2} = 5730\pm 40\text{ năm} (chu kỳ Cambridge). Theo quy ước quốc tế, tuổi thô được báo cáo bằng 'năm cacbon phóng xạ BP' (tính từ mốc 1950 CN) sử dụng chu kỳ Libby (5568 na˘m5568\text{ năm}). Những dao động trong hoạt động mặt trời và từ trường Trái Đất đòi hỏi đường cong hiệu chuẩn IntCal20 dựa trên vòng cây và trầm tích varve. Máy phổ khối gia tốc (AMS) đếm trực tiếp từng nguyên tử từ mẫu dưới miligam, mở rộng giới hạn đo đến ∼50000 na˘m\sim 50000\text{ năm}.

Các hệ đồng vị địa thời học chính và khoáng vật mục tiêu
Hệ phân rãKiểu phân rãChu kỳ bán rãDải đo hiệu dụngKhoáng vật và vật liệu tiêu biểu
X14X2214C→X14X2214N\ce{^{14}C -> ^{14}N}Phân rã beta-5730 năm300 đến 50.000 nămThan bùn, gỗ, xương, vỏ sò, than củi
X40X2240K→X40X2240Ar\ce{^{40}K -> ^{40}Ar}Bắt electron (10,7%)1,248 tỷ năm100 nghìn năm đến 4,5 tỷ nămSanidin, biotit, muscovit, hornblend
X87X2287Rb→X87X2287Sr\ce{^{87}Rb->^{87}Sr}Phân rã beta-48,8 tỷ năm10 triệu năm đến 4,6 tỷ nămMica, fenspat kali, đá magma toàn phần
X238X22238U→X206X22206Pb\ce{^{238}U->^{206}Pb}Chuỗi alpha và beta4,468 tỷ năm1 triệu năm đến 4,6 tỷ nămZircon, baddeleyit, monazit, titanit
X235X22235U→X207X22207Pb\ce{^{235}U->^{207}Pb}Chuỗi alpha và beta0,7038 tỷ năm1 triệu năm đến 4,6 tỷ nămGhép cặp với 238U trên đồ thị hoà hợp
X147X22147Sm→X143X22143Nd\ce{^{147}Sm->^{143}Nd}Phân rã alpha106 tỷ năm100 triệu năm đến 4,6 tỷ nămBasalt, anorthosit mặt trăng, thiên thạch

Ví dụ: Tính toán tuổi cacbon phóng xạ thô

Một dầm gỗ tuyết tùng cổ khai quật từ một di chỉ khảo cổ có hoạt độ riêng X14X2214C\ce{^{14}C} đo được bằng đúng 12,5%12{,}5\% (1/81/8) hoạt độ của gỗ thời kỳ tiền công nghiệp. Dùng chu kỳ bán rã vật lý chuẩn Cambridge t1/2=5730 na˘mt_{1/2} = 5730\text{ năm}, hãy tính tuổi phóng xạ thô của cổ vật gỗ này trước khi hiệu chuẩn vòng cây.

Lời giải

Tỉ số hoạt độ còn lại là N/N0=0,125=(1/2)3N/N_0 = 0{,}125 = (1/2)^3. Giá trị này ứng với đúng 3 chu kỳ bán rã: t=3×t1/2=3×5730=17190 na˘mt = 3 \times t_{1/2} = 3 \times 5730 = 17190\text{ năm}. Viết theo công thức phân rã luỹ thừa: t=−t1/2ln⁡2ln⁡(NN0)=−57300,69315×ln⁡(0,125)=−57300,69315×(−2,07944)≈17190 na˘mt = -\frac{t_{1/2}}{\ln 2}\ln\left(\frac{N}{N_0}\right) = -\frac{5730}{0{,}69315}\times \ln(0{,}125) = -\frac{5730}{0{,}69315}\times (-2{,}07944) \approx 17190\text{ năm}. Các nhà khảo cổ sau đó sẽ hiệu chuẩn tuổi thô này qua đường cong IntCal20 để ra tuổi lịch cal BP.

Đại họcĐại học: phân rã phân nhánh K-Ar, 40Ar/39Ar và đường đẳng thời 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}

Kali-40 phân rã phân nhánh: bắt electron thành khí hiếm X40X2240Ar\ce{^{40}Ar} (10,72%10{,}72\%, λe=0,580×10−10 na˘m−1\lambda_e = 0{,}580 \times 10^{-10}\ \text{năm}^{-1}) và phát xạ β−\beta^- thành X40X2240Ca\ce{^{40}Ca} (89,28%89{,}28\%, λβ=4,963×10−10 na˘m−1\lambda_\beta = 4{,}963 \times 10^{-10}\ \text{năm}^{-1}), cho tổng hằng số phân rã λ=5,543×10−10 na˘m−1\lambda = 5{,}543 \times 10^{-10}\ \text{năm}^{-1} (t1/2=1,248 tỷ na˘mt_{1/2} = 1{,}248\text{ tỷ năm}). Magma nóng chảy thoát hết khí argon ra ngoài, đưa tỉ số ban đầu về X40X2240Ar0≈0\ce{^{40}Ar}_0 \approx 0. Khi kết tinh, X40X2240Ar∗\ce{^{40}Ar}^* phóng xạ sinh ra bị giam chặt trong mạng tinh thể chứa kali như sanidin hay biotit. Kỹ thuật hiện đại X40X2240Ar/X39X2239Ar\ce{^{40}Ar/^{39}Ar} chiếu xạ mẫu trong lò phản ứng hạt nhân để chuyển hoá X39X2239K(n,p)X39X2239Ar\ce{^{39}K}(n,p)\ce{^{39}Ar}. Nung nóng từng bước bằng laser trong khối phổ kế khí hiếm giải phóng argon theo từng mức nhiệt, vẽ nên đồ thị cao nguyên tuổi giúp phát hiện mất mát argon hoặc argon dư bẫy trong mẫu.

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)

Khi khoáng vật chứa lượng đồng vị con ban đầu D0D_0 khác không tại thời điểm kết tinh, việc đo một mẫu đơn lẻ không thể tính được tuổi. Phương pháp đẳng thời (isochron) giải quyết bài toán này bằng cách chuẩn hoá theo một đồng vị bền không phóng xạ của nguyên tố con, chẳng hạn X86X2286Sr\ce{^{86}Sr} trong hệ rubidi-stronti (X87X2287Rb→X87X2287Sr\ce{^{87}Rb -> ^{87}Sr}, t1/2=48,8 tỷ na˘mt_{1/2} = 48{,}8\text{ tỷ năm}, λ=1,42×10−11 na˘m−1\lambda = 1{,}42 \times 10^{-11}\ \text{năm}^{-1}). Các khoáng vật khác nhau cùng kết tinh từ một lò magma có tỉ số nguyên tố Rb/Sr\ce{Rb/Sr} khác xa nhau do phân đoạn hoá học tinh thể, nhưng lại có cùng tỉ số đồng vị ban đầu (X87X2287Sr/X86X2286Sr)0(\ce{^{87}Sr/^{86}Sr})_0. Đồ thị biểu diễn X87X2287Sr/X86X2286Sr\ce{^{87}Sr/^{86}Sr} theo X87X2287Rb/X86X2286Sr\ce{^{87}Rb/^{86}Sr} tạo thành một đường thẳng có hệ số góc m=eλt−1≈λtm = e^{\lambda t} - 1 \approx \lambda t, cho phép xác định trực tiếp tuổi kết tinh tt mà không cần biết trước lượng con ban đầu.

Ví dụ: Xác định niên đại tinh thể sanidin núi lửa bằng phương pháp K-Ar

Phép đo khối phổ tỉ số đồng vị độ chính xác cao trên các tinh thể sanidin tách từ một lớp túp núi lửa cho tỉ số nguyên tử giữa argon phóng xạ và kali là X40X2240Ar∗X40X2240K=0,0950\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}} = 0{,}0950. Với các hằng số phân rã λ=5,543×10−10 na˘m−1\lambda = 5{,}543 \times 10^{-10}\ \text{năm}^{-1} và λe=0,580×10−10 na˘m−1\lambda_e = 0{,}580 \times 10^{-10}\ \text{năm}^{-1}, hãy tính tuổi kết tinh của lớp trầm tích núi lửa này theo triệu năm (Ma).

Lời giải

Thay số vào phương trình tuổi phân nhánh 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). Hệ số phân nhánh là λλ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. Tính biểu thức trong logarit: 1+9,5569×0,0950=1+0,9079=1,90791 + 9{,}5569 \times 0{,}0950 = 1 + 0{,}9079 = 1{,}9079. Giá trị logarit tự nhiên ln⁡(1,9079)≈0,6460\ln(1{,}9079) \approx 0{,}6460. Chia cho hằng số phân rã λ\lambda: t=0,64605,543×10−10 na˘m−1≈1,1654×109 na˘m=1165,4 Mat = \frac{0{,}6460}{5{,}543 \times 10^{-10}\ \text{năm}^{-1}} \approx 1{,}1654 \times 10^9\text{ năm} = 1165{,}4\text{ Ma} (1,165 tỷ na˘m1{,}165\text{ tỷ năm}).

Nâng caoNâng cao: đồ thị hoà hợp U-Pb, zircon Jack Hills và tuổi của Trái Đất

Zircon (ZrSiOX4\ce{ZrSiO4}) là tiêu chuẩn vàng của địa thời học. Trong quá trình magma kết tinh, vị trí zirconi (r=84 pmr = 84\text{ pm}) dung nạp UX4+\ce{U^4+} (r=100 pmr = 100\text{ pm}) vào mạng trực thoi nhưng loại trừ triệt để ion PbX2+\ce{Pb^2+} (r=119 pmr = 119\text{ pm}, điện tích +2+2) với hệ số phân bố vượt quá 10510^5, đưa lượng chì ban đầu về Pb0≈0\ce{Pb}_0 \approx 0. Vì zircon chứa hai đồng hồ phân rã độc lập với chu kỳ bán rã rất khác nhau—X238X22238U→X206X22206Pb\ce{^{238}U -> ^{206}Pb} (t1/2=4,468 tỷ na˘mt_{1/2} = 4{,}468\text{ tỷ năm}, λ238=1,55125×10−10 na˘m−1\lambda_{238} = 1{,}55125 \times 10^{-10}\ \text{năm}^{-1}) và X235X22235U→X207X22207Pb\ce{^{235}U -> ^{207}Pb} (t1/2=0,7038 tỷ na˘mt_{1/2} = 0{,}7038\text{ tỷ năm}, λ235=9,8485×10−10 na˘m−1\lambda_{235} = 9{,}8485 \times 10^{-10}\ \text{năm}^{-1})—đồ thị X206X22206Pb∗/X238X22238U\ce{^{206}Pb}^*/\ce{^{238}U} theo X207X22207Pb∗/X235X22235U\ce{^{207}Pb}^*/\ce{^{235}U} vạch nên đường cong hoà hợp Wetherill. Các tinh thể không bị xáo trộn sẽ nằm hoà hợp trên đường cong; mẫu bị mất chì do biến chất sẽ tạo thành đường bất hoà hợp tuyến tính mà giao điểm trên cho tuổi kết tinh nguyên sinh và giao điểm dưới cho tuổi biến chất. Những hạt zircon trầm tích tại Jack Hills ở Tây Úc cho tuổi hoà hợp 4,404 tỷ na˘m4{,}404\text{ tỷ năm}, bảo tồn dấu tích vỏ Trái Đất kỷ Hoả minh.

Lấy tỉ số hai phương trình phân rã của chì phóng xạ sẽ triệt tiêu nhu cầu đo nồng độ urani tuyệt đối: 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), với 137,88137{,}88 là tỉ số X238X22238U/X235X22235U\ce{^{238}U/^{235}U} tự nhiên ngày nay. Năm 1956, Clair Patterson đã xác định thành phần chì nguyên thuỷ của Hệ Mặt Trời bằng cách phân tích pha troilit (FeS\ce{FeS}) không chứa urani trong thiên thạch sắt Canyon Diablo. Đặt số liệu chì thiên thạch cùng chì trầm tích đáy đại dương lên một đường đẳng thời Pb-Pb duy nhất, Patterson đã tính ra tuổi của Trái Đất là 4,55±0,07 tỷ na˘m4{,}55 \pm 0{,}07\text{ tỷ năm} (4,54 tỷ na˘m4{,}54\text{ tỷ năm}). Các phân tích U-Pb độ chính xác cao ngày nay trên các thể vùi giàu canxi-nhôm (CAI) trong thiên thạch chondrit định tuổi thời điểm ngưng tụ của đĩa tiền hành tinh Hệ Mặt Trời là 4,567 tỷ na˘m4{,}567\text{ tỷ năm}.

Tài liệu tham khảo

  • 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