Bối cảnh và vấn đề nghiên cứu
Trong kỹ thuật hạt nhân và khoa học vật liệu, kiểm tra không phá hủy (Non-Destructive Testing - NDT) là phương pháp quan trọng nhằm phát hiện các khuyết tật cấu trúc bên trong vật liệu mà không làm tổn hại đến hình dạng hay cơ tính của mẫu khảo sát. Một trong những kỹ thuật NDT cốt lõi là phương pháp chụp ảnh phóng xạ, hoạt động dựa trên nguyên lý suy giảm cường độ của chùm tia bức xạ khi xuyên qua các chiều dày vật liệu khác nhau. Để đánh giá định lượng khả năng làm suy giảm và che chắn bức xạ của các vật liệu dạng hợp chất, các thông số đặc trưng cho tương tác photon-vật chất cần được xác định chính xác, trong đó thông số nền tảng là nguyên tử số hiệu dụng ($Z_{eff}$).
Hiện nay, nhiều phương pháp thực nghiệm và lý thuyết đã được phát triển để xác định nguyên tử số hiệu dụng, tiêu biểu là phương pháp đo gamma tán xạ và phương pháp đo gamma truyền qua. Nghiên cứu của Kucuk và cộng sự (2014) đã thiết lập hệ đo gamma truyền qua sử dụng đầu dò nhấp nháy NaI(Tl) để xác định $Z_{eff}$ và mật độ electron hiệu dụng cho 5 loại polyme ở nhiều mức năng lượng, cho thấy sự phù hợp cao giữa số liệu thực nghiệm và tính toán lý thuyết. Polyme là nhóm vật liệu cao phân tử tổng hợp có ưu điểm độ bền cao, tính thẩm mỹ, dễ gia công chế tạo và đang được ứng dụng rộng rãi trong đời sống cũng như kỹ thuật (điển hình như việc phát triển bê tông polyme giảm phát thải khí nhà kính). Đồng thời, trong nghiên cứu trước đó của Huỳnh Đình Chương và cộng sự, việc dự đoán mật độ chất lỏng dựa trên tỷ lệ diện tích đỉnh tán xạ đơn so với nước đã giả định mật độ chỉ phụ thuộc vào năng lượng photon mà bỏ qua ảnh hưởng của thành phần nguyên tố. Do đó, việc xác định nguyên tử số hiệu dụng của các hợp chất polyme là yêu cầu cần thiết để hoàn thiện cơ sở dữ liệu và đánh giá mức độ ảnh hưởng của thành phần vật liệu lên phép đo mật độ.
Mục tiêu và nhiệm vụ nghiên cứu của khóa luận được xác lập theo các nội dung sau:
- Tổng quan cơ sở lý thuyết về các quá trình tương tác của bức xạ gamma với vật chất, định luật suy giảm bức xạ và các mô hình toán học xác định hệ số suy giảm khối cùng nguyên tử số hiệu dụng.
- Xây dựng mô hình mô phỏng hệ đo gamma truyền qua trên chương trình MCNP6, sử dụng nguồn phóng xạ $^{137}\text{Cs}$ (mức năng lượng 0,662 MeV) và đầu dò NaI(Tl) để ghi nhận xung bức xạ.
- Thu thập phổ năng lượng từ mô phỏng, xử lý phổ bằng phần mềm Colegram để xác định cường độ bức xạ truyền qua 14 loại vật liệu polyme và qua không khí.
- Tính toán hệ số suy giảm khối ($\mu_m$) và nguyên tử số hiệu dụng ($Z_{eff}$) bằng phương pháp mô phỏng Monte Carlo ($Z_{eff,\text{MC}}$) và phương pháp tính trực tiếp theo lý thuyết ($Z_{eff,\text{PI}}$).
- Đánh giá độ sai lệch tương đối giữa giá trị mô phỏng với dữ liệu chuẩn WinXCom của NIST, đồng thời so sánh kết quả tính toán với dữ liệu thực nghiệm đã công bố từ các công trình nghiên cứu quốc tế trước đây.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Nguyên tử số hiệu dụng ($Z_{eff}$) và hệ số suy giảm khối ($\mu_m$) của 14 loại vật liệu polyme thông dụng, bao gồm: ABS, HDPE, PMMA, PC, CR, PEEK, PET, PMP, POM, PP, PS, PVC, PE, PVDF.
- Phạm vi năng lượng và nguồn bức xạ: Nguồn bức xạ gamma đơn năng $^{137}\text{Cs}$ có mức năng lượng phát 661,657 keV (tương đương 0,662 MeV).
- Phạm vi không gian và thời gian: Khóa luận đượ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 Thành phố Hồ Chí Minh, hoàn thành vào năm 2020.
Cơ sở lý thuyết và phương pháp
Khung lý thuyết và các mô hình tính toán
Khóa luận xây dựng trên nền tảng vật lý bức xạ tương tác với vật chất, tập trung vào ba hiệu ứng cơ bản của photon gamma trong môi trường:
- Hiệu ứng quang điện: Photon gamma va chạm với electron liên kết trong nguyên tử và truyền toàn bộ năng lượng để giải phóng quang electron với động năng cực đại $E_e = E - E_{lk}$.
- Hiệu ứng Compton: Tán xạ của photon gamma lên electron tự do hoặc electron liên kết yếu ở lớp ngoài, photon bị lệch hướng một góc $\theta$ và năng lượng sau tán xạ tuân theo công thức:
$$E' = \frac{E}{1 + \frac{E}{m_e c^2}(1 - \cos\theta)}$$
- Hiệu ứng tạo cặp: Photon gamma có năng lượng lớn hơn hai lần năng lượng nghỉ của electron ($E \ge 2m_e c^2 = 1,022\text{ MeV}$) tương tác với điện trường hạt nhân tạo thành cặp electron-positron, sau đó positron hủy cặp với electron sinh ra hai photon $0,511\text{ MeV}$ bay ngược chiều nhau.
Định luật suy giảm chùm tia gamma hẹp đơn năng khi qua bề dày vật chất $x$ (cm) được mô tả theo hàm mũ:
$$I = I_0 e^{-\mu x} = I_0 e^{-\mu_m \rho x}$$
Trong đó $I_0$ là cường độ chùm tia tới, $I$ là cường độ chùm tia sau khi qua bia, $\mu$ ($\text{cm}^{-1}$) là hệ số suy giảm tuyến tính, $\rho$ ($\text{g/cm}^3$) là khối lượng riêng và $\mu_m = \mu/\rho$ ($\text{cm}^2/\text{g}$) là hệ số suy giảm khối. Từ đó, hệ số suy giảm khối thực nghiệm/mô phỏng được xác định qua công thức:
$$\mu_m = \frac{1}{\rho x} \ln\left(\frac{I_0}{I}\right)$$
Theo lý thuyết hỗn hợp đối với hợp chất, $\mu_m$ được tính bằng tổng đóng góp của từng nguyên tố thành phần:
$$\mu_m = \sum_i w_i (\mu_m)_i$$
với $w_i = \frac{n_i A_i}{\sum_j n_j A_j}$ là phần trăm khối lượng của nguyên tố thứ $i$.
Nguyên tử số hiệu dụng ($Z_{eff}$) của hợp chất được định nghĩa bằng tỷ số giữa tiết diện tương tác nguyên tử hiệu dụng ($\sigma_a$) và tiết diện tương tác electron hiệu dụng ($\sigma_{el}$):
$$Z_{eff} = \frac{\sigma_a}{\sigma_{el}}$$
Trong đó:
- Tiết diện tương tác phân tử: $\sigma_m = \frac{1}{N_A} \mu_m \sum_i n_i A_i$
- Tiết diện tương tác nguyên tử hiệu dụng: $\sigma_a = \frac{\sigma_m}{\sum_i n_i} = \frac{1}{N_A} \sum_i f_i A_i (\mu_m)_i$
- Tiết diện tương tác electron hiệu dụng: $\sigma_{el} = \frac{1}{N_A} \sum_i \frac{f_i A_i}{Z_i} (\mu_m)_i$
- $f_i = \frac{n_i}{\sum_j n_j}$ là tỷ lệ số nguyên tử của nguyên tố thứ $i$ trong hợp chất; $A_i$ và $Z_i$ lần lượt là khối lượng nguyên tử và nguyên tử số của nguyên tố thứ $i$.
Từ các biểu thức trên, phương pháp tính trực tiếp ($Z_{eff,\text{PI}}$) có dạng công thức toán học:
$$Z_{eff,\text{PI}} = \frac{\sum_i f_i A_i (\mu_m)_i}{\sum_j \frac{f_j A_j}{Z_j} (\mu_m)_j}$$
Bên cạnh đó, tác giả cũng hệ thống hóa phương pháp của chương trình XMuDat (Nowotny) áp dụng công thức $Z_{eff} = \left(\sum_i \alpha_i Z_i^{m-1}\right)^{1/(m-1)}$ với $3 \le m \le 5$ và phương pháp nội suy hàm logarit để xác định nguyên tử số tương đương ($Z_{eq}$).
Phương pháp nghiên cứu và công cụ kỹ thuật
Khóa luận áp dụng kết hợp hai phương pháp chính:
- Phương pháp tính trực tiếp theo lý thuyết: Tính toán $Z_{eff,\text{PI}}$ thông qua việc tra cứu hệ số suy giảm khối từng nguyên tố $(\mu_m)_i$ tại mức năng lượng $0,662\text{ MeV}$ trên cơ sở dữ liệu WinXCom của NIST.
- Phương pháp mô phỏng Monte Carlo: Sử dụng phần mềm MCNP6 (Monte Carlo N-Particle version 6) để mô phỏng hình học không gian 3 chiều của hệ đo gamma truyền qua. Nguồn phát $^{137}\text{Cs}$ được cấu hình phát $6 \times 10^9$ hạt photon đẳng hướng. Xung năng lượng tại tinh thể đầu dò $\text{NaI(Tl)}$ được thu nhận thông qua thẻ Tally F8.
- Phương pháp xử lý phổ: Dữ liệu phổ mô phỏng xuất ra từ MCNP6 được đưa vào phần mềm Colegram để bóc tách diện tích đỉnh năng lượng $662\text{ keV}$. Đỉnh phổ được khớp bằng hàm phân bố Gauss và nền phổ được làm khớp bằng hàm đa thức bậc 4:
$$P(x) = \frac{1}{\sigma \sqrt{2\pi}} e^{-\frac{(x - \bar{x})^2}{2\sigma^2}}, \quad y = a_0 + a_1 x + a_2 x^2 + a_3 x^3 + a_4 x^4$$
- Đánh giá sai số: Độ lệch tương đối (Relative Deviation - RD) giữa giá trị tính toán/mô phỏng so với dữ liệu chuẩn NIST được xác định theo biểu thức:
$$\text{RD } (%) = \frac{|\text{Giá trị chuẩn} - \text{Giá trị khảo sát}|}{\text{Giá trị chuẩn}} \times 100$$
Nguồn dữ liệu
- Dữ liệu tương tác photon: Thư viện WinXCom (NIST) và các thư viện hạt nhân liên tục ENDF, ENDL, ACTL tích hợp trong MCNP6.
- Thông số phân rã phóng xạ của $^{137}\text{Cs}$: Tài liệu kỹ thuật nguồn chuẩn từ hãng Eckert & Ziegler (Eckert & Ziegler Reference & Calibration Source Production Information).
- Dữ liệu mật độ và công thức hóa học của polyme: Dữ liệu WinXCom, tài liệu của Liên đoàn Nhựa Anh quốc (BPF) và dữ liệu kỹ thuật từ Công ty TNHH Thương mại Plastic IDO.
- Dữ liệu so sánh: Các công bố quốc tế của Kucuk và cộng sự (2014), Manohara và cộng sự, Singh và cộng sự (2014), Medhat (2012).
Nội dung chính theo từng chương
Chương 1. Tổng quan về nguyên tử số hiệu dụng
Chương này tập trung trình bày bản chất vật lý của bức xạ gamma, các cơ chế tương tác cơ bản khi photon đi qua môi trường vật chất và hệ thống các công thức xác định hệ số suy giảm khối lượng cùng nguyên tử số hiệu dụng. Tác giả phân tích chi tiết điều kiện xảy ra và quy luật năng lượng của hiệu ứng quang điện, hiệu ứng tán xạ Compton và hiệu ứng tạo cặp. Đồng thời, chương 1 thiết lập các biểu thức liên hệ từ định luật suy giảm hàm mũ đến các đại lượng tiết diện tương tác phân tử, tiết diện nguyên tử hiệu dụng và tiết diện electron hiệu dụng. Bốn phương pháp xác định $Z_{eff}$ được giới thiệu gồm: phương pháp tính trực tiếp, phương pháp Monte Carlo, phương pháp XMuDat và phương pháp nội suy tiết diện tương tác hàm logarit.
Chương 2. Phương pháp Monte Carlo và chương trình MCNP6
Chương 2 mô tả nguồn gốc, cơ sở toán xác suất của phương pháp Monte Carlo bắt đầu từ bài toán cây kim của Buffon (1777), bước đi ngẫu nhiên của Lord Rayleigh (1899) đến sự phát triển của nhóm nghiên cứu tại Phòng thí nghiệm Quốc gia Los Alamos. Tác giả giới thiệu cấu trúc và tính năng của chương trình MCNP6 (công bố năm 2013, hợp nhất từ MCNP và MCNPX), có khả năng mô phỏng vận chuyển 37 loại hạt cơ bản và hạt nhân trên nền các thư viện dữ liệu hạt nhân liên tục. Cấu trúc tập tin đầu vào (input file) của MCNP6 được phân tích chi tiết thành các khối thẻ bắt buộc:
- Thẻ định nghĩa ô mạng (Cell Cards): Khai báo số hiệu cell, số hiệu vật liệu, mật độ (mang dấu âm nếu là $\text{g/cm}^3$, dấu dương nếu là nguyên tử/$\text{barn-cm}$), các quan hệ hình học giao, hội, bù.
- Thẻ định nghĩa mặt (Surface Cards): Khai báo các mặt giới hạn không gian hình học, trong đó các mặt sử dụng trong đề tài gồm mặt phẳng vuông góc với trục tọa độ (PX, PY, PZ) và mặt trụ dọc theo trục Z (CZ).
- Thẻ định nghĩa dữ liệu (Data Cards): Cấu hình nguồn phát tổng quát SDEF (loại hạt photon PAR=2, năng lượng ERG, vị trí POS, hướng phát DIR, bán kính quét RAD, khoảng cách EXT), các thẻ phân bố SIn, SPn, SBn và thẻ ghi nhận xung Tally F8 đóng vai trò như đầu dò nhấp nháy ghi nhận năng lượng photon hấp thụ trong thể tích nhạy.
- Thẻ khai báo vật liệu (Material Cards): Khai báo thành phần hóa học thông qua ký hiệu ZAID (ZZZAAA) và tỷ lệ phần trăm số nguyên tử hoặc tỷ lệ khối lượng.
Chương 3. Xác định nguyên tử số hiệu dụng
Chương 3 là phần thực nghiệm mô phỏng và phân tích kết quả. Tác giả xây dựng cấu hình hình học của hệ đo gamma truyền qua với các kích thước cụ thể:
- Khoảng cách từ vị trí nguồn $^{137}\text{Cs}$ đến thanh vật liệu polyme là $28,43\text{ cm}$.
- Khoảng cách từ thanh vật liệu polyme đến bề mặt đầu dò $\text{NaI(Tl)}$ là $20,3\text{ cm}$.
- Nguồn bức xạ $^{137}\text{Cs}$: Năng lượng 661,657 keV, hiệu suất phát photon $84,99%$, số lượng lịch sử hạt chạy mô phỏng là 6 tỷ hạt ($6 \times 10^9$ particles).
- Mẫu polyme: Dạng thanh trụ đặt trong ống có đường kính $2\text{ cm}$.
- Cấu hình đầu dò nhấp nháy $\text{NaI(Tl)}$: Tinh thể nhạy có đường kính $76\text{ mm}$, chiều cao $76\text{ mm}$; lớp vỏ bảo vệ gồm lớp chì phía trên dày $1,5\text{ mm}$, lớp chì xung quanh dày $1,5\text{ mm}$, lớp silicon đệm dày $2\text{ mm}$ và lớp phản xạ phía trước đầu dò dày $3\text{ mm}$.
Bảng thông số 14 loại vật liệu polyme được khảo sát:
| Tên vật liệu |
Tên viết tắt |
Công thức hóa học |
Mật độ trung bình ($\text{g/cm}^3$) |
Nguồn trích dẫn mật độ |
| Acrylonitrin Butadien Styren |
ABS |
$\text{C}{15}\text{H}{17}\text{N}$ |
1,105 |
1,06 - 1,15 (Plastic IDO) |
| High-density polyethylene |
HDPE |
$\text{C}_2\text{H}_4$ |
0,970 |
0,96 - 0,98 (Plastic IDO) |
| Polymethylmethacrylate |
PMMA |
$\text{C}_5\text{H}_8\text{O}_2$ |
1,190 |
1,19 (WinXCom) |
| Polycarbonate |
PC |
$\text{C}{15}\text{H}{16}\text{O}_2$ |
1,240 |
1,20 - 1,28 (Plastic IDO) |
| Polychloroprene |
CR |
$\text{C}_4\text{H}_5\text{Cl}$ |
1,230 |
1,23 (Plastic IDO) |
| Polyetheretherketones |
PEEK |
$\text{C}{19}\text{H}{12}\text{O}_3$ |
1,350 |
1,3 - 1,4 (Plastic IDO) |
| Polyethyleneterephthalate |
PET |
$\text{C}_{10}\text{H}_8\text{O}_4$ |
1,415 |
1,33 - 1,5 (Plastic IDO) |
| Polymethylpentene |
PMP |
$\text{C}6\text{H}{12}$ |
0,830 |
0,83 (BPF) |
| Polyoxymethylene |
POM |
$\text{CH}_2\text{O}$ |
1,425 |
1,45 - 1,55 (WinXCom) |
| Polypropylene |
PP |
$\text{C}_3\text{H}_6$ |
0,930 |
0,93 (Plastic IDO) |
| Polystyrene |
PS |
$\text{C}_8\text{H}_8$ |
1,060 |
1,06 (WinXCom) |
| Polyvinylchloride |
PVC |
$\text{C}_2\text{H}_3\text{Cl}$ |
1,300 |
1,3 (Plastic IDO) |
| Polyethylene |
PE |
$\text{C}_2\text{H}_4$ |
0,940 |
0,94 (WinXCom) |
| Polyvinylidenefluoride |
PVDF |
$\text{C}_2\text{H}_2\text{F}_2$ |
1,820 |
1,79 - 1,85 (Plastic IDO) |
Sau khi chạy mô phỏng, tác giả ghi nhận cường độ photon truyền qua môi trường không khí (đóng vai trò cường độ ban đầu $I_0$) đạt $349.960.500$ xung tại mật độ không khí $0,001205\text{ g/cm}^3$. Tiếp theo, cường độ chùm photon truyền qua 14 mẫu polyme ($I$) được xác định thông qua việc khử phông và khớp đỉnh Gauss trên Colegram, từ đó tính ra hệ số suy giảm khối $\mu_m$ và nguyên tử số hiệu dụng $Z_{eff}$.
Kết quả và đóng góp
Số liệu kết quả tính toán và mô phỏng
Kết quả xác định hệ số suy giảm khối ($\mu_m$) và nguyên tử số hiệu dụng ($Z_{eff}$) của 14 loại vật liệu polyme ở mức năng lượng $0,662\text{ MeV}$ được tổng hợp chi tiết trong bảng sau:
| Vật liệu |
Cường độ truyền qua ($I$) |
$\mu_{m,\text{MC}}$ ($\text{cm}^2/\text{g}$) |
$\mu_{m,\text{NIST}}$ ($\text{cm}^2/\text{g}$) |
$\text{RD}_{\mu_m}$ (%) |
$Z_{eff,\text{PI}}$ |
$Z_{eff,\text{MC}}$ |
$\text{RD}{Z{eff}}$ (%) |
| ABS |
297.714.400 |
0,0732 |
0,0832 |
12,1 |
3,45 |
3,04 |
12,1 |
| HDPE |
301.193.400 |
0,0774 |
0,0880 |
12,1 |
2,67 |
2,34 |
12,1 |
| PMMA |
294.056.500 |
0,0731 |
0,0832 |
12,1 |
3,60 |
3,16 |
12,1 |
| PC |
293.131.700 |
0,0715 |
0,0825 |
13,4 |
3,70 |
3,20 |
13,4 |
| CR |
293.986.300 |
0,0708 |
0,0802 |
11,7 |
4,60 |
4,06 |
11,6 |
| PEEK |
289.268.800 |
0,0705 |
0,0803 |
12,2 |
4,41 |
3,88 |
12,1 |
| PET |
286.632.400 |
0,0705 |
0,0803 |
12,2 |
4,55 |
3,99 |
12,1 |
| PMP |
307.780.200 |
0,0774 |
0,0880 |
12,1 |
2,67 |
2,34 |
12,1 |
| POM |
284.862.300 |
0,0722 |
0,0822 |
12,1 |
4,00 |
3,51 |
12,1 |
| PP |
303.061.200 |
0,0774 |
0,0880 |
12,1 |
2,67 |
2,34 |
12,1 |
| PS |
299.829.100 |
0,0729 |
0,0830 |
12,1 |
3,50 |
3,08 |
12,1 |
| PVC |
291.739.100 |
0,0700 |
0,0790 |
11,4 |
5,33 |
4,72 |
11,4 |
| PE |
302.594.000 |
0,0774 |
0,0880 |
12,1 |
2,67 |
2,34 |
12,1 |
| PVDF |
273.554.800 |
0,0677 |
0,0771 |
12,2 |
5,33 |
4,68 |
12,3 |
So sánh đối chiếu giá trị $Z_{eff}$ với các công bố khoa học trước đây:
| Vật liệu |
$Z_{eff,\text{PI}}$ (Khóa luận) |
$Z_{eff,\text{MC}}$ (Khóa luận) |
Nghiên cứu khác |
Độ lệch $Z_{eff,\text{PI}}$ so với tác giả khác |
Độ lệch $Z_{eff,\text{MC}}$ so với tác giả khác |
| PMMA |
3,60 |
3,16 |
3,4 [22]; 3,52 [23] |
< 6% |
- |
| PC |
3,70 |
3,20 |
3,726 [24] |
< 6% |
Sai biệt lớn nhất: -16,3% |
| PP |
2,67 |
2,34 |
2,707 [23]; 2,66 [4] |
< 6% |
- |
| PS |
3,50 |
3,08 |
3,001 [24] |
14,3% |
- |
| PE |
2,67 |
2,34 |
2,710 [24] |
< 6% |
Sai biệt nhỏ nhất: 2,4% |
Nhận xét và đóng góp của nghiên cứu
- Phương pháp tính trực tiếp ($Z_{eff,\text{PI}}$): Cho kết quả rất sát với các công trình thực nghiệm quốc tế đối với 4/5 loại polyme khảo sát (PMMA, PC, PP, PE) với độ lệch đều dưới 6%. Riêng đối với polystyrene (PS), độ chênh lệch so với số liệu của Singh và cộng sự là 14,3%. Điều này khẳng định tính đúng đắn và độ tin cậy của thuật toán tính trực tiếp từ dữ liệu tiết diện nguyên tử WinXCom.
- Phương pháp mô phỏng Monte Carlo ($Z_{eff,\text{MC}}$): Toàn bộ hệ số suy giảm khối và nguyên tử số hiệu dụng thu được từ mô hình MCNP6 đều có độ lệch tương đối ổn định nằm trong khoảng từ 11,4% đến 13,4% (thấp nhất là PVC với 11,4% và cao nhất là PC với 13,4%). Độ lệch này mang tính hệ thống.
- Đóng góp khoa học: Khóa luận đã hệ thống hóa quy trình tính toán lý thuyết và quy trình xây dựng tập tin mô phỏng hệ đo gamma truyền qua trên MCNP6; cung cấp bộ số liệu tính toán chi tiết về $\mu_m$ và $Z_{eff}$ cho 14 loại hợp chất polyme thông dụng ở năng lượng $0,662\text{ MeV}$, đóng vai trò làm dữ liệu tham chiếu phục vụ cho hướng nghiên cứu đánh giá ảnh hưởng của thành phần vật liệu lên phương pháp đo mật độ chất lỏng trước đây của nhóm nghiên cứu.
Hạn chế và hướng nghiên cứu tiếp
Khóa luận thẳng thắn chỉ ra các điểm tồn tại và giải pháp hoàn thiện được tác giả đúc kết:
- Tác giả tự đánh giá khóa luận chưa đạt được hoàn toàn mục đích nghiên cứu đặt ra ban đầu do tồn tại độ sai lệch hệ thống từ 11,4% đến 13,4% giữa kết quả mô phỏng Monte Carlo và dữ liệu chuẩn WinXCom.
- Nguyên nhân cốt lõi của sai số: Trong quá trình xây dựng tập tin đầu vào MCNP6, việc khai báo thẻ vật liệu (Material Cards) chưa hoàn toàn tương thích và đồng nhất với cơ sở dữ liệu WinXCom; đồng thời các giá trị mật độ polyme được tổng hợp từ nhiều nguồn tài liệu kỹ thuật thương mại khác nhau chứ chưa được đo đạc thực nghiệm trực tiếp trên cùng một mẫu chuẩn.
- Các hướng nghiên cứu tiếp theo được đề xuất trong văn bản:
- Kiểm tra, hiệu chỉnh lại cấu trúc tập tin đầu vào của mô hình mô phỏng MCNP6 (đặc biệt là chuẩn hóa khai báo tỷ lệ đồng vị và mật độ trong ô mạng).
- Khai thác thêm các khía cạnh khác của các phương pháp xác định nguyên tử số hiệu dụng.
- Mở rộng áp dụng các phương pháp tính toán khác (như chương trình XMuDat, phương pháp nội suy hàm logarit tính $Z_{eq}$) để có thêm dữ liệu đa chiều phục vụ so sánh và đánh giá.
- Triển khai bố trí hệ thực nghiệm vật lý thực tế để đo đạc trực tiếp trên mẫu, từ đó kiểm chứng và đánh giá chính xác độ phù hợp của mô hình mô phỏng MCNP6.
Giá trị tham khảo
Báo cáo khóa luận có giá trị tham khảo rõ ràng cho sinh viên và học viên chuyên ngành Vật lý hạt nhân, Kỹ thuật hạt nhân, Vật lý y khoa và Khoa học vật liệu:
- Tài liệu hướng dẫn thiết lập MCNP6: Cung cấp chi tiết cú pháp viết tập tin đầu vào cho bài toán tương tác bức xạ gamma đơn năng, bao gồm cách khai báo các mặt chuẩn (PX, PY, PZ, CZ), cấu hình ô mạng, khai báo nguồn bức xạ dạng ống chuẩn trực (SDEF) và thiết lập đầu dò ghi nhận phổ nhấp nháy NaI(Tl) bằng thẻ Tally F8.
- Quy trình xử lý phổ bức xạ: Minh họa cụ thể phương pháp tách nền và khớp đỉnh năng lượng bằng cách kết hợp hàm Gauss và đa thức bậc 4 trên phần mềm Colegram.
- Cơ sở dữ liệu vật lý vật liệu: Đưa ra bảng tra cứu công thức hóa học, mật độ, hệ số suy giảm khối và nguyên tử số hiệu dụng của 14 hợp chất polyme tiêu biểu ở mức năng lượng $0,662\text{ MeV}$, hỗ trợ trực tiếp cho các nghiên cứu về vật liệu che chắn bức xạ nhẹ hoặc chế tạo phantom mô phỏng mô mềm sinh học.
Câu hỏi thường gặp
1. Tại sao nghiên cứu lại chọn nguồn phóng xạ $^{137}\text{Cs}$ và mức năng lượng 0,662 MeV?
Trong khóa luận, tác giả sử dụng nguồn $^{137}\text{Cs}$ vì đây là nguồn chuẩn phát photon gamma đơn năng với năng lượng đặc trưng $661,657\text{ keV}$ (xấp xỉ $0,662\text{ MeV}$) và hiệu suất phát cao ($84,99%$). Ở mức năng lượng này, tương tác chủ yếu của photon với các vật liệu có nguyên tử số thấp như polyme là hiệu ứng tán xạ Compton, rất thuận lợi cho việc kiểm tra định luật suy giảm bức xạ và xác định nguyên tử số hiệu dụng.
2. Có bao nhiêu hạt photon được phát ra trong mô phỏng MCNP6 và sử dụng thẻ tally nào để ghi phổ?
Trong chương trình MCNP6, tác giả thiết lập số lượng lịch sử hạt phát ra từ nguồn là 6 tỷ hạt ($6 \times 10^9$ hạt). Đầu dò nhấp nháy $\text{NaI(Tl)}$ được mô phỏng thu nhận xung năng lượng bằng thẻ Tally F8, cho phép ghi nhận sự phân bố phổ năng lượng bị mất trong thể tích nhạy cảm của tinh thể đầu dò.
3. Nguyên nhân chính dẫn đến độ lệch 11,4% – 13,4% giữa kết quả mô phỏng MCNP6 và chuẩn WinXCom là gì?
Tác giả giải trình nguyên nhân gây ra độ lệch này là do trong quá trình tạo tập tin đầu vào MCNP6, thành phần nguyên tố và tỷ lệ đồng vị của vật liệu polyme khai báo trong thẻ vật liệu chưa hoàn toàn tương thích với dữ liệu WinXCom, kết hợp với việc mật độ khối lượng của các loại polyme được thu thập từ nhiều nguồn tài liệu khác nhau thay vì đo trên một mẫu đồng nhất.
4. Trong 14 loại polyme, vật liệu nào có nguyên tử số hiệu dụng lý thuyết cao nhất và thấp nhất?
Theo kết quả tính toán trực tiếp ($Z_{eff,\text{PI}}$):
- Vật liệu có nguyên tử số hiệu dụng cao nhất là PVC ($\text{C}_2\text{H}_3\text{Cl}$) và PVDF ($\text{C}_2\text{H}_2\text{F}2$) với cùng giá trị $Z{eff,\text{PI}} = 5,33$ (do trong phân tử chứa nguyên tố Clo $Z=17$ và Flo $Z=9$).
- Vật liệu có nguyên tử số hiệu dụng thấp nhất là các polyme chỉ chứa Carbon và Hydro gồm HDPE, PMP, PP và PE với cùng giá trị $Z_{eff,\text{PI}} = 2,67$.
5. Phần mềm nào được sử dụng để bóc tách diện tích đỉnh phổ năng lượng từ dữ liệu MCNP6?
Dữ liệu xung năng lượng thu được từ thẻ Tally F8 của MCNP6 được xử lý bằng phần mềm Colegram. Trong đó, diện tích đỉnh năng lượng $662\text{ keV}$ được làm khớp bằng hàm phân bố Gauss và thành phần phông nền bức xạ được làm khớp bằng hàm đa thức bậc 4.
Kết luận
Khóa luận tốt nghiệp của sinh viên Trần Thị Mỹ Duyên đã thực hiện khảo sát toàn diện hệ số suy giảm khối và nguyên tử số hiệu dụng của 14 loại vật liệu polyme ở mức năng lượng $0,662\text{ MeV}$ thông qua phương pháp tính trực tiếp và mô phỏng Monte Carlo trên MCNP6. Kết quả tính lý thuyết trực tiếp cho thấy sự phù hợp cao với dữ liệu thực nghiệm đã công bố (độ lệch dưới 6% trên đa số mẫu), trong khi kết quả mô phỏng MCNP6 phản ánh sự suy giảm chùm tia mang tính hệ thống với độ lệch từ 11,4% đến 13,4% so với chuẩn WinXCom. Mặc dù cần tiếp tục hiệu chỉnh tham số khai báo vật liệu và tiến hành thực nghiệm đối chứng, công trình đã xây dựng thành công mô hình mô phỏng hệ truyền qua hoàn chỉnh và cung cấp bộ dữ liệu vật lý quan trọng phục vụ hướng nghiên cứu đánh giá mật độ vật chất trong kỹ thuật hạt nhân.
