Bối cảnh và vấn đề nghiên cứu
Kỹ thuật kiểm tra không phá hủy (Non-Destructive Testing - NDT) đóng vai trò quan trọng trong việc đánh giá khuyết tật cấu trúc, tính chất vật liệu và kiểm soát chất lượng công nghiệp mà không làm tổn hại đến đối tượng kiểm tra. Trong số các phương pháp NDT, phương pháp bức xạ gamma truyền qua được ứng dụng rộng rãi để xác định các tham số tương tác bức xạ như hệ số suy giảm khối ($\mu_m$), tiết diện tương tác, mật độ electron ($N_{eff}$) và nguyên tử số hiệu dụng ($Z_{eff}$). Nguyên tử số $Z$ là tham số cơ bản đối với các đơn nguyên tố, nhưng đối với các hợp chất và dung dịch phức hợp, một giá trị $Z$ đơn lẻ không thể đại diện cho toàn bộ hệ thống; thay vào đó, tham số $Z_{eff}$ phụ thuộc vào năng lượng photon được sử dụng để mô tả khả năng tương tác và che chắn bức xạ.
Mặc dù $Z_{eff}$ đã được khảo sát nhiều trên các vật liệu như polymer, hợp kim, thủy tinh, vật liệu xây dựng và mẫu sinh học, các nghiên cứu tập trung vào dung dịch chất lỏng — đặc biệt là các axit vô cơ và hữu cơ — vẫn còn hạn chế. Bên cạnh đó, phương pháp đo mật độ khối ($\rho$) truyền thống bằng cân thể tích thường tiềm ẩn nguy cơ mất an toàn hóa học đối với các dung dịch axit mạnh, độc hại hoặc dễ bay hơi. Việc ứng dụng kỹ thuật gamma truyền qua cho phép xác định mật độ chất lỏng trực tiếp trong ống chứa kín mà không cần tiếp xúc hóa chất.
Đề tài đặt ra 3 mục tiêu nghiên cứu cụ thể:
- Tính toán giá trị nguyên tử số hiệu dụng $Z_{eff}$ cho các loại axit bằng phương pháp tính trực tiếp và phương pháp nội suy dựa trên cơ sở dữ liệu của Viện Tiêu chuẩn và Kỹ thuật Quốc gia Hoa Kỳ (NIST).
- Đo đạc thực nghiệm và mô phỏng Monte Carlo bằng chương trình MCNP6 để xác định hệ số suy giảm khối lượng $\mu_m$, từ đó suy ra $Z_{eff}$ và đánh giá độ sai lệch giữa mô phỏng, thực nghiệm và dữ liệu chuẩn.
- Xây dựng đường chuẩn mật độ cho dung dịch axit thông qua phương pháp mô phỏng và thực nghiệm nhằm xác định mật độ chất lỏng bằng kỹ thuật bức xạ truyền qua.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: 6 loại axit thực nghiệm (axit acetic, formic, lactic, nitric, phosphoric, sulfuric) cùng nước cất ($H_2O$ làm mẫu đối chứng), và 17 loại axit bổ sung trong mô phỏng (tổng cộng 24 chất khảo sát, mật độ từ 0,430 g/cm³ đến 1,839 g/cm³).
- Nguồn phóng xạ: Nguồn đồng vị $^{137}\text{Cs}$ phát bức xạ gamma đơn năng lượng 661,7 keV.
- Phạm vi thời gian và không gian: Nghiên cứu thực nghiệm được thực hiện tại Bộ môn Vật lý Hạt nhân, Khoa Vật lý, Trường Đại học Sư phạm TP.HCM (hoàn thành tháng 7 năm 2020), duy trì ở nhiệt độ môi trường 24°C, thời gian đo mỗi mẫu là 28.800 giây.
Cơ sở lý thuyết và phương pháp
Khóa luận thiết lập hệ thống cơ sở lý thuyết dựa trên các định luật tương tác bức xạ gamma với vật chất:
- Cơ chế tương tác photon: Quá trình hấp thụ quang điện ($\sigma_{photo} \sim Z^5/E^{7/2}$), tán xạ Compton ($\sigma_{Compton} \sim Z/E$), và hiệu ứng tạo cặp ($\sigma_{pair} \sim Z^2\ln(E)$). Ở mức năng lượng 661,7 keV của nguồn $^{137}\text{Cs}$, hiệu ứng tán xạ Compton chiếm ưu thế chủ đạo.
- Định luật suy giảm Beer-Lambert: Cường độ chùm tia gamma hẹp truyền qua lớp vật chất có bề dày $x$ và mật độ $\rho$ tuân theo phương trình $I = I_0 e^{-\mu x} = I_0 e^{-\mu_m \rho x}$, trong đó $\mu$ là hệ số suy giảm tuyến tính (cm⁻¹) và $\mu_m = \mu/\rho$ là hệ số suy giảm khối (cm²/g). Với chùm tia rộng, phương trình được hiệu chỉnh bằng hệ số tích lũy $B(E, Z)$; tuy nhiên, trong bố trí chùm tia chuẩn trực hẹp của nghiên cứu, hệ số $B$ được xem bằng 1.
- Quy tắc pha trộn (Mixture Rule): Hệ số suy giảm khối lý thuyết của hợp chất được tính theo tổng trọng số nồng độ khối lượng của các nguyên tố thành phần: $\mu_m = \sum w_i (\mu_m)_i$.
- Các mô hình tính $Z_{eff}$:
- Phương pháp tính trực tiếp: Tính từ thành phần phần trăm nguyên tử $f_i$, khối lượng nguyên tử $A_i$, nguyên tử số $Z_i$ và hệ số suy giảm khối lượng nguyên tố.
- Phương pháp tính theo tiết diện nguyên tử: Xác định tiết diện nguyên tử toàn phần $\sigma_a = (\mu/\rho)_{Compound} / (N_A \sum w_i/A_i)$ và tiết diện electron hiệu dụng $\sigma_e = (1/N_A) \sum [f_i A_i (\mu/\rho)i / Z_i]$, từ đó $Z{eff} = \sigma_a / \sigma_e$.
- Các công thức bán nghiệm khác: Công thức chương trình máy tính XMuDat, công thức Spiers ($m=3,94$), Tsai - Cho ($m=3,4$), Manninen - Koikkalainen, và Glasser.
Phương pháp nghiên cứu kết hợp giữa mô phỏng tính toán và đo đạc thực nghiệm:
- Phương pháp mô phỏng Monte Carlo: Sử dụng mã MCNP6 phiên bản 1.0 của Phòng thí nghiệm Quốc gia Los Alamos (LANL), vận chuyển $6 \times 10^9$ lịch sử hạt với phân bố phổ độ cao xung (Pulse Height Spectrum - PHS) 8192 kênh. Độ phân giải năng lượng được định hình qua hàm thực nghiệm $FWHM(MeV) = a + b\sqrt{E + cE^2}$ với các tham số khớp $a = -0,00526\text{ MeV}$, $b = 0,06096\text{ MeV}^{1/2}$, $c = 0,01401\text{ MeV}^{-1}$.
- Hệ đo thực nghiệm: Sử dụng đầu dò nhấp nháy $\text{NaI(Tl)}$ kích thước $76,2\text{ mm} \times 76,2\text{ mm}$ do Amptek sản xuất, kết hợp bộ phân tích đa kênh MCA 8192 kênh, chuẩn trực bằng khối chì (bán kính chuẩn trực nguồn $0,251\text{ cm}$, bán kính chuẩn trực đầu dò $0,50\text{ cm}$), căn chỉnh đồng trục bằng tia laser đỏ.
- Nguồn dữ liệu và xử lý số liệu: Dữ liệu chuẩn đối chứng trích xuất từ cơ sở dữ liệu XCOM/XMuDat của NIST. Xử lý phổ gamma bằng phần mềm Colegram (khớp đỉnh phân bố Gauss, phông nền dạng bậc One Step bằng thuật toán Marquardt-Levenberg) và Origin. Sai số thực nghiệm được đánh giá theo độ lệch chuẩn trung bình và phương pháp truyền sai số với độ tin cậy 99,99%.
Nội dung chính theo từng chương
Chương 1: Tổng quan về gamma truyền qua
Chương 1 hệ thống hóa lịch sử phát hiện tia gamma từ Paul Villard (1900) đến các nghiên cứu định danh của Ernest Rutherford (1903). Tác giả trình bày bản chất sóng điện từ năng lượng cao của gamma và ba cơ chế tương tác cơ bản với vật chất (quang điện, Compton, tạo cặp), nhấn mạnh sự phụ thuộc của tiết diện tương tác vào năng lượng photon $E$ và nguyên tử số $Z$. Cấu trúc phổ năng lượng gamma thực tế bao gồm đỉnh hấp thụ toàn phần (FEA), vùng tán xạ đơn (SS), tán xạ bội (MS), tán xạ ngược (BS) và các đỉnh thoát (SPE, DPE).
Chương này phân tích toán học định luật Beer-Lambert cho chùm tia chuẩn trực hẹp và chùm tia rộng kèm hệ số tích lũy $B(E,Z)$. Cơ chế xác định mật độ chất lỏng được thiết lập qua mối liên hệ tuyến tính $\ln(I) = \ln(I_0) - \mu_m \rho x$, tạo cơ sở lý thuyết cho kỹ thuật đo mật độ không tiếp xúc.
Chương 2: Phương pháp xác định nguyên tử số hiệu dụng
Chương 2 trình bày chi tiết các phương pháp xác định $Z_{eff}$ của hợp chất. Tác giả phân chia thành hai nhóm tiếp cận chính:
- Phương pháp trực tiếp và phương pháp tiết diện nguyên tử: Dựa trên hàm lượng nguyên tố và hệ số suy giảm khối lượng chuẩn từ NIST kết hợp quy tắc Mixture Rule để tính tiết diện phân tử $\sigma_m$, tiết diện nguyên tử $\sigma_a$, tiết diện electron $\sigma_e$, từ đó suy ra $Z_{eff} = \sigma_a/\sigma_e$.
- Phương pháp nội suy và các công thức bán nghiệm: Phân tích cơ sở của chương trình XMuDat và các phương trình do Spiers ($m=3,94$), Tsai - Cho ($m=3,4$), Manninen - Koikkalainen, và Glasser đề xuất.
Khóa luận xác định sử dụng phương pháp tính trực tiếp và phương pháp tính theo tiết diện nguyên tử làm hai trục tính toán chính để so sánh kết quả mô phỏng và thực nghiệm.
Chương 3: Khảo sát nguyên tử số hiệu dụng và mật độ của axit bằng mô phỏng Monte Carlo và thực nghiệm
Chương 3 mô tả chi tiết quy trình bố trí hệ đo thực nghiệm và xây dựng mô hình tính toán trên MCNP6:
- Bố trí thực nghiệm: Nguồn $^{137}\text{Cs}$ (hoạt độ thời điểm đo $0,76 \pm 0,03\ \mu\text{Ci}$) đặt trong khối chì $24\text{ cm} \times 22\text{ cm} \times 5\text{ cm}$; khoảng cách từ nguồn và mặt đầu dò đến tâm ống mẫu là $12,4\text{ cm}$. Mẫu chất lỏng đựng trong ống đong thủy tinh 100 ml có bán kính trong đo bằng panme là $r = 1,340 \pm 0,00898\text{ cm}$.
Thông số vật liệu thực nghiệm chính:
- Axit Acetic (C2H4O2, 1,045 g/cm3) - Xilong
- Axit Formic (CH2O2, 1,220 g/cm3) - Xilong
- Axit Lactic (C3H6O3, 1,206 g/cm3) - Xilong
- Axit Nitric (HNO3, 1,513 g/cm3) - Merck
- Axit Phosphoric (H3PO4, 1,390 g/cm3) - VWR Chemicals BDH
- Axit Sulfuric (H2SO4, 1,839 g/cm3) - Merck
- Nước cất (H2O, 1,000 g/cm3) - Merck
- Mô phỏng MCNP6: Thiết lập hình học 3D đồng nhất với thực nghiệm, mô phỏng 24 loại axit và nước cất với cấu trúc đầu vào khai báo đầy đủ thành phần hóa học, khối lượng riêng, cấu tạo đầu dò $\text{NaI(Tl)}$ có vỏ bọc $\text{Al}$, cửa sổ $\text{SiO}_2$, lớp phản xạ $\text{Al}_2\text{O}_3$ và lớp silicon.
- Xử lý phổ và tính sai số: Phổ PHS được trích xuất ROI, trừ phông môi trường, làm khớp bằng phần mềm Colegram kết hợp Origin. Sai số của bán kính ống đo và tỉ số đếm $R = N_{Acid}/N_{H2O}$ được tính toán qua công thức truyền sai số ở độ tin cậy 99,99%.
Chương 4: Xác định nguyên tử số hiệu dụng và mật độ của một số axit
Chương 4 tổng hợp toàn bộ số liệu thu được từ thực nghiệm và mô phỏng MCNP6:
- Tính toán hệ số suy giảm khối $\mu_m$ và $Z_{eff}$ cho 6 axit thực nghiệm và nước cất, so sánh với giá trị chuẩn NIST.
- Tính toán $\mu_m$ và $Z_{eff}$ cho 24 chất lỏng mô phỏng MCNP6, đối chiếu với số liệu của các nghiên cứu trước đây tại mức năng lượng 1 MeV.
- Xây dựng đường chuẩn mật độ tương quan giữa tỉ số đếm $R_{MP} = N_{Acid}/N_{H2O}$ và mật độ $\rho$, kiểm chứng độ chính xác khi nội suy mật độ thực nghiệm.
Kết quả và đóng góp
1. Kết quả đo đạc thực nghiệm và đối chiếu chuẩn NIST
Hệ số suy giảm khối lượng $\mu_m$ và nguyên tử số hiệu dụng $Z_{eff}$ đo bằng thực nghiệm trên hệ phổ kế gamma $\text{NaI(Tl)}$ đạt độ tương thích cao với cơ sở dữ liệu NIST:
| STT |
Vật liệu |
Mật độ $\rho$ (g/cm³) |
$\mu_m$ NIST (cm²/g) |
$\mu_m$ Thực nghiệm (cm²/g) |
$Z_{eff}$ NIST |
$Z_{eff}$ Thực nghiệm |
Độ lệch $Z_{eff}$ (%) |
| 1 |
Nước cất ($H_2O$) |
1,000 |
0,086 |
$0,085 \pm 0,006$ |
3,33 |
$3,31 \pm 0,09$ |
0,57 |
| 2 |
Axit Acetic ($C_2H_4O_2$) |
1,045 |
0,082 |
$0,085 \pm 0,006$ |
4,00 |
$3,83 \pm 0,07$ |
4,28 |
| 3 |
Axit Formic ($CH_2O_2$) |
1,220 |
0,080 |
$0,079 \pm 0,005$ |
4,80 |
$4,70 \pm 0,07$ |
2,02 |
| 4 |
Axit Lactic ($C_3H_6O_3$) |
1,206 |
0,082 |
$0,079 \pm 0,005$ |
4,00 |
$4,13 \pm 0,09$ |
3,24 |
| 5 |
Axit Nitric ($HNO_3$) |
1,513 |
0,078 |
$0,080 \pm 0,005$ |
6,43 |
$6,42 \pm 0,07$ |
0,07 |
| 6 |
Axit Phosphoric ($H_3PO_4$) |
1,390 |
0,079 |
$0,087 \pm 0,004$ |
6,19 |
$6,04 \pm 0,05$ |
2,44 |
| 7 |
Axit Sulfuric ($H_2SO_4$) |
1,839 |
0,079 |
$0,077 \pm 0,004$ |
7,14 |
$7,87 \pm 0,06$ |
10,17 |
Tất cả các giá trị $Z_{eff}$ thực nghiệm đều thỏa mãn điều kiện giới hạn vật lý $Z_{min} < Z_{eff} < Z_{max}$ của các nguyên tố cấu thành. Axit Nitric có độ sai biệt nhỏ nhất so với NIST (0,07%), trong khi Axit Sulfuric có độ lệch lớn nhất (10,17%).
2. Kết quả mô phỏng MCNP6
Đối với 24 chất lỏng mô phỏng, độ lệch tương đối của $\mu_m$ so với NIST phân bố đồng đều dưới 8% (dao động từ 3,6% ở Histine đến 7,4% ở Perchloric). Giá trị $Z_{eff}$ tính từ mô phỏng MCNP6 có độ sai lệch từ 3,43% đến 6,87% so với chuẩn NIST.
| Loại axit |
$Z_{eff}$ NIST |
$Z_{eff}$ Thực nghiệm |
$Z_{eff}$ Mô phỏng |
Sai lệch TN - MP (%) |
| Axit Acetic |
4,00 |
3,83 |
3,75 |
2,09 |
| Axit Formic |
4,80 |
4,70 |
4,49 |
4,47 |
| Axit Lactic |
4,00 |
4,13 |
3,74 |
9,44 |
| Axit Nitric |
6,43 |
6,42 |
5,87 |
8,57 |
| Axit Phosphoric |
6,19 |
6,04 |
5,76 |
4,64 |
| Axit Sulfuric |
7,14 |
7,87 |
6,65 |
15,50 |
3. Xây dựng đường chuẩn xác định mật độ chất lỏng
Đề tài đã thiết lập thành công đường chuẩn mật độ bán thực nghiệm biểu diễn logarit tự nhiên của tỉ số đếm $\ln(R)$ theo mật độ $\rho$ ($R = N_{Acid}/N_{H2O}$). Đường chuẩn cho phép xác định mật độ của các dung dịch axit kín với độ chính xác cao, chứng minh tính khả thi của kỹ thuật gamma truyền qua trong việc đo đạc an toàn các chất lỏng độc hại.
4. Đóng góp của đề tài
- Bổ sung cơ sở dữ liệu thực nghiệm và mô phỏng về nguyên tử số hiệu dụng $Z_{eff}$ và hệ số suy giảm khối $\mu_m$ cho nhóm dung dịch axit vô cơ và hữu cơ — một mảng dữ liệu còn thiếu hụt so với các nghiên cứu trên vật liệu rắn.
- Chuẩn hóa quy trình xử lý phổ gamma thực nghiệm bằng Colegram kết hợp mã mô phỏng MCNP6 trong nghiên cứu tương tác bức xạ với chất lỏng.
Hạn chế và hướng nghiên cứu tiếp
Giá trị tham khảo
Khóa luận là tài liệu tham khảo chuyên khảo hữu ích cho:
- Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Vật lý hạt nhân, Kỹ thuật hạt nhân, và An toàn bức xạ trong việc nắm bắt phương pháp đo phổ gamma truyền qua và kỹ thuật xử lý phổ đỉnh thực nghiệm bằng Colegram.
- Các kỹ sư NDT và cán bộ nghiên cứu công nghệ bức xạ quan tâm đến giải pháp đo lường mật độ không tiếp xúc đối với chất lỏng ăn mòn, độc hại hoặc các hợp chất trong đường ống kín.
- Người nghiên cứu mô phỏng MCNP6 có thể tham khảo tệp dữ liệu đầu vào mẫu (input file) trong phần phụ lục cùng phương pháp khớp hàm bề rộng năng lượng FWHM.
Câu hỏi thường gặp
1. Nguồn bức xạ và đầu dò nào được sử dụng trong nghiên cứu thực nghiệm?
Tác giả sử dụng nguồn phóng xạ đồng vị $^{137}\text{Cs}$ phát tia gamma đơn năng 661,7 keV với xác suất phát 84,99% (hoạt độ đo đạc thực tế $0,76 \pm 0,03\ \mu\text{Ci}$) và đầu dò nhấp nháy $\text{NaI(Tl)}$ kích thước $76,2\text{ mm} \times 76,2\text{ mm}$ do hãng Amptek cung cấp, kết hợp bộ phân tích đa kênh MCA 8192 kênh.
2. Các hàm toán học nào được dùng để xử lý phổ gamma trong phần mềm Colegram?
Trong phần mềm Colegram, tác giả sử dụng hàm phân bố Gaussian $G(x) = \frac{A}{\sigma \sqrt{2\pi}} \exp\left[-\frac{(x-x_0)^2}{2\sigma^2}\right]$ để làm khớp đỉnh hấp thụ toàn phần và hàm bậc One Step $STEP(X) = B + \frac{S}{2} \text{erfc}\left[\frac{x-x_0}{\sigma\sqrt{2}}\right]$ để mô tả hình dạng phông nền dưới đỉnh.
3. Quy chuẩn nào được sử dụng để xác định giá trị $Z_{eff}$ tính toán là hợp lệ?
Một giá trị nguyên tử số hiệu dụng $Z_{eff}$ của vật liệu hợp chất được xem là hợp lệ khi nằm trong khoảng giới hạn $Z_{min} < Z_{eff} < Z_{max}$, trong đó $Z_{min}$ và $Z_{max}$ lần lượt là nguyên tử số nhỏ nhất và lớn nhất của các nguyên tố cấu thành hợp chất (ví dụ với $H_2SO_4$ gồm H ($Z=1$), O ($Z=8$), S ($Z=16$), giá trị $Z_{eff} = 7,87$ nằm trong khoảng 1 đến 16).
4. Tại sao hệ số tích lũy $B(E, Z)$ lại được lấy bằng 1 trong các tính toán của khóa luận?
Tác giả sử dụng hai khối chuẩn trực bằng chì đặt trước nguồn (bán kính $0,251\text{ cm}$) và trước đầu dò (bán kính $0,50\text{ cm}$) nhằm tạo ra chùm tia hẹp và song song, hạn chế tối đa các photon tán xạ nhiều lần đi vào thể tích nhạy của đầu dò. Nhờ đó, đóng góp của các thành phần gamma thứ cấp là không đáng kể và hệ số tích lũy $B(E,Z)$ được xấp xỉ bằng 1.
5. Số lượng lịch sử hạt và hàm FWHM được thiết lập như thế nào trong mô phỏng MCNP6?
Trong chương trình MCNP6, tác giả thực hiện mô phỏng với tổng số lịch sử hạt là $6 \times 10^9$ hạt. Hàm bề rộng nửa độ cao FWHM được khai báo theo công thức $FWHM(MeV) = a + b\sqrt{E + cE^2}$ với các hệ số khớp thực nghiệm gồm $a = -0,00526\text{ MeV}$, $b = 0,06096\text{ MeV}^{1/2}$, và $c = 0,01401\text{ MeV}^{-1}$.
Kết luận
Khóa luận của tác giả Nguyễn Thu Hằng đã thực hiện thành công việc xác định nguyên tử số hiệu dụng $Z_{eff}$ và hệ số suy giảm khối $\mu_m$ cho 6 loại axit thực nghiệm và 24 loại axit mô phỏng bằng phương pháp gamma truyền qua với nguồn $^{137}\text{Cs}$. Kết quả thực nghiệm và mô phỏng MCNP6 cho thấy sự phù hợp chặt chẽ với cơ sở dữ liệu chuẩn của NIST, đồng thời chứng minh mối tương quan tuyến tính giữa mật độ chất lỏng và logarit tỉ số đếm bức xạ. Công trình cung cấp giải pháp đo mật độ dung dịch axit an toàn trong ống kín và bổ sung dữ liệu thực nghiệm có giá trị cho lĩnh vực vật lý bức xạ và kiểm tra không phá hủy.
