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ạ P phân rã thành đồng vị con D∗ với hằng số phân rã λ, số nguyên tử con tích luỹ là D∗=P(eλt−1). Giải tìm tuổi tinh thể t ta được t=λ1ln(1+PD∗). Để 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 (Tc) đặ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+X01X2021nX614X26214C+X11X2121p. Lượng cacbon-14 sinh ra bị oxy hoá thành X14X2214COX2 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 trong mô sống. Khi sinh vật chết đi, sự trao đổi chất dừng lại và X14X2214C phân rã β− trở lại thành X14X2214N với chu kỳ bán rã t1/2=5730±40 na˘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˘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.
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ụng
Khoáng vật và vật liệu tiêu biểu
X14X2214CX14X2214N
Phân rã beta-
5730 năm
300 đến 50.000 năm
Than bùn, gỗ, xương, vỏ sò, than củi
X40X2240KX40X2240Ar
Bắt electron (10,7%)
1,248 tỷ năm
100 nghìn năm đến 4,5 tỷ năm
Sanidin, biotit, muscovit, hornblend
X87X2287RbX87X2287Sr
Phân rã beta-
48,8 tỷ năm
10 triệu năm đến 4,6 tỷ năm
Mica, fenspat kali, đá magma toàn phần
X238X22238UX206X22206Pb
Chuỗi alpha và beta
4,468 tỷ năm
1 triệu năm đến 4,6 tỷ năm
Zircon, baddeleyit, monazit, titanit
X235X22235UX207X22207Pb
Chuỗi alpha và beta
0,7038 tỷ năm
1 triệu năm đến 4,6 tỷ năm
Ghép cặp với 238U trên đồ thị hoà hợp
X147X22147SmX143X22143Nd
Phân rã alpha
106 tỷ năm
100 triệu năm đến 4,6 tỷ năm
Basalt, 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 đo được bằng đúng 12,5% (1/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˘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)3. Giá trị này ứng với đúng 3 chu kỳ bán rã: t=3×t1/2=3×5730=17190 na˘m. Viết theo công thức phân rã luỹ thừa: t=−ln2t1/2ln(N0N)=−0,693155730×ln(0,125)=−0,693155730×(−2,07944)≈17190 na˘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
Kali-40 phân rã phân nhánh: bắt electron thành khí hiếm X40X2240Ar (10,72%, λe=0,580×10−10na˘m−1) và phát xạ β− thành X40X2240Ca (89,28%, λβ=4,963×10−10na˘m−1), cho tổng hằng số phân rã λ=5,543×10−10na˘m−1 (t1/2=1,248 tỷ na˘m). Magma nóng chảy thoát hết khí argon ra ngoài, đưa tỉ số ban đầu về X40X2240Ar0≈0. Khi kết tinh, X40X2240Ar∗ 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 chiếu xạ mẫu trong lò phản ứng hạt nhân để chuyển hoá X39X2239K(n,p)X39X2239Ar. 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.
86Sr87Sr=(86Sr87Sr)0+86Sr87Rb(eλt−1)
Khi khoáng vật chứa lượng đồng vị con ban đầu D0 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 trong hệ rubidi-stronti (X87X2287RbX87X2287Sr, t1/2=48,8 tỷ na˘m, λ=1,42×10−11na˘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 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. Đồ thị biểu diễn X87X2287Sr/X86X2286Sr theo X87X2287Rb/X86X2286Sr tạo thành một đường thẳng có hệ số góc m=eλt−1≈λt, cho phép xác định trực tiếp tuổi kết tinh t 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à X40X2240KX40X2240Ar∗=0,0950. Với các hằng số phân rã λ=5,543×10−10na˘m−1 và λe=0,580×10−10na˘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=λ1ln(1+λeλX40X2240KX40X2240Ar∗). Hệ số phân nhánh là λeλ=0,580×10−105,543×10−10≈9,5569. Tính biểu thức trong logarit: 1+9,5569×0,0950=1+0,9079=1,9079. Giá trị logarit tự nhiên ln(1,9079)≈0,6460. Chia cho hằng số phân rã λ: t=5,543×10−10na˘m−10,6460≈1,1654×109 na˘m=1165,4 Ma (1,165 tỷ na˘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) 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 pm) dung nạp UX4+ (r=100 pm) vào mạng trực thoi nhưng loại trừ triệt để ion PbX2+ (r=119 pm, điện tích +2) với hệ số phân bố vượt quá 105, đưa lượng chì ban đầu về Pb0≈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—X238X22238UX206X22206Pb (t1/2=4,468 tỷ na˘m, λ238=1,55125×10−10na˘m−1) và X235X22235UX207X22207Pb (t1/2=0,7038 tỷ na˘m, λ235=9,8485×10−10na˘m−1)—đồ thị X206X22206Pb∗/X238X22238U theo X207X22207Pb∗/X235X22235U 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˘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: X206X22206Pb∗X207X22207Pb∗=X238X22238UX235X22235U(eλ238t−1eλ235t−1)=137,881(eλ238t−1eλ235t−1), với 137,88 là tỉ số X238X22238U/X235X22235U 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) 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˘m (4,54 tỷ na˘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˘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