Giới thiệu dự án

Trong lĩnh vực vật lý hạt nhân ứng dụng và công nghệ vật liệu, kiểm tra không phá hủy (Non-Destructive Testing - NDT) bằng phương pháp chụp ảnh phóng xạ và đo truyền qua gamma đóng vai trò then chốt để phát hiện khuyết tật cấu trúc, kiểm tra độ đồng nhất và đánh giá khả năng suy giảm bức xạ. Để thiết kế các hệ che chắn bức xạ hiệu quả hoặc chế tạo các phantom mô phỏng mô người trong y học hạt nhân, việc xác định chính xác các thông số tương tác photon với vật chất – đặc biệt là nguyên tử số hiệu dụng ($Z_{eff}$) và hệ số suy giảm khối ($\mu_m$) – là yêu cầu bắt buộc.

Vật liệu polyme (chất cao phân tử) sở hữu các đặc tính vượt trội: khối lượng riêng nhẹ, độ bền cơ học cao, dễ gia công định hình và thân thiện với môi trường (như bê tông polyme giảm phát thải khí nhà kính). Tuy nhiên, việc xác định $Z_{eff}$ cho các hợp chất polyme đa nguyên tố phức tạp hơn nhiều so với đơn chất do phụ thuộc phi tuyến vào năng lượng photon tới ($E_\gamma$) và cơ chế tương tác vi mô (hiệu ứng quang điện, tán xạ Compton, tạo cặp).

                                  KIẾN TRÚC MÔ PHỎNG GAMMA TRUYỀN QUA
+------------------+      +-------------------+      +------------------+      +--------------------+
|  Nguồn 137-Cs    | ---> | Bộ chuẩn trực chì | ---> |  Bia Polyme      | ---> | Bộ chuẩn trực      | ---> [ Đầu dò NaI(Tl) ]
|  (0.662 MeV)     |      | (Collimator 1)    |      |  (14 mẫu ống phi)|      | (Collimator 2)     |      | (Tally F8 MCNP6)
+------------------+      +-------------------+      +------------------+      +--------------------+               |
                                                                                                                    v
                                                                                                           +--------------------+
                                                                                                           | Khớp phổ Colegram  |
                                                                                                           | (Gauss + Bậc 4)    |
                                                                                                           +--------------------+

Mục tiêu dự án

  1. Xây dựng mô hình tính toán lý thuyết trực tiếp: Khai triển giải thuật giải tích xác định tiết diện tương tác phân tử ($\sigma_m$), tiết diện nguyên tử ($\sigma_a$), tiết diện electron ($\sigma_{el}$) và nguyên tử số hiệu dụng ($Z_{eff, PI}$) dựa trên cơ sở dữ liệu tiết diện quang hạt nhân chuẩn NIST/WinXCom tại mức năng lượng $0.662\text{ MeV}$.
  2. Thiết kế mô hình số hóa Monte Carlo 3D: Xây dựng không gian hình học thực nghiệm truyền qua gamma trên phần mềm MCNP6 (Monte Carlo N-Particle version 6.1), tối ưu hóa quá trình vận chuyển hạt với $6 \times 10^9$ lịch sử phát photon từ nguồn phóng xạ $^{137}\text{Cs}$.
  3. Khảo sát hệ dữ liệu 14 vật liệu polyme thông dụng: Xác định hệ số suy giảm khối $\mu_m$ và $Z_{eff}$ cho: ABS, HDPE, PMMA, PC, CR, PEEK, PET, PMP, POM, PP, PS, PVC, PE, PVDF.
  4. Phát triển quy trình xử lý phổ năng lượng: Ứng dụng giải thuật phi tuyến khớp đỉnh Gauss và phông đa thức bậc 4 trên phần mềm Colegram để trích xuất diện tích đỉnh năng lượng $661.657\text{ keV}$, khử nhiễu phông tán xạ.
  5. Thẩm định và đối soát kết quả: Đánh giá sai số tương đối (Relative Deviation - RD%) giữa phương pháp Monte Carlo, lý thuyết trực tiếp và các công trình thực nghiệm quốc tế uy tín.

Phạm vi và giới hạn

  • Năng lượng khảo sát: Đơn năng $0.662\text{ MeV}$ phát ra từ đồng vị $^{137}\text{Cs}$ (xác suất phân rã nhánh gamma $84.99%$).
  • Cấu hình hình học: Hệ chùm tia hẹp (narrow-beam geometry), khoảng cách nguồn - mẫu $28.43\text{ cm}$, mẫu - đầu dò $20.3\text{ cm}$, đường kính mẫu $2.0\text{ cm}$.
  • Giới hạn: Chưa xét đến phổ liên tục đa năng lượng và các hiệu ứng tán xạ bậc cao ngoài vùng nhạy của đầu dò nhấp nháy $\text{NaI(Tl)}$.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp Ưu điểm Nhược điểm Chi phí & Độ an toàn
Thực nghiệm Gamma truyền qua Kết quả thực tế, phản ánh đúng khuyết tật chế tạo Cần nguồn hoạt độ cao, thời gian đo lâu, rủi ro an toàn bức xạ Chi phí đầu tư cao ($>20,000\text{ USD}$), yêu cầu cấp phép an toàn nghiêm ngặt
Lý thuyết giải tích (WinXCom) Tốc độ tức thì, không rủi ro phóng xạ Bỏ qua yếu tố hình học thực tế, tán xạ phông nền và đáp ứng đầu dò Miễn phí, thuần tính toán lý thuyết
Mô phỏng số Monte Carlo (MCNP6) Tái lập chính xác hình học 3D, vật liệu phức tạp, chi tiết cơ chế tương tác Đòi hỏi tài nguyên tính toán lớn, phụ thuộc độ chính xác của thẻ đầu vào Chi phí thấp, an toàn tuyệt đối, độ tin cậy cao
Phần mềm XMuDat / Nội suy Giao diện đơn giản cho vật liệu chuẩn Kém linh hoạt với hình học dị hình và các hợp chất không phổ biến Miễn phí, hạn chế cập nhật thư viện hạt nhân mới

Yêu cầu kỹ thuật theo mô hình MoSCoW

  • Must have: Cấu hình chuẩn xác thẻ ô mạng (Cell Cards), mặt biên (Surface Cards), vật liệu (Material Cards) trong MCNP6; mô phỏng đáp ứng đầu dò $\text{NaI(Tl)}$ bằng Tally F8; giải thuật phân tích đỉnh phổ $0.662\text{ MeV}$.
  • Should have: Tự động hóa tính toán $Z_{eff}$ từ phổ đầu ra; kiểm soát sai số thống kê của MCNP6 dưới $1%$.
  • Could have: Mở rộng dải năng lượng từ $59.54\text{ keV}$ ($^{241}\text{Am}$) đến $1.332\text{ MeV}$ ($^{60}\text{Co}$).
  • Won't have: Mô phỏng biến tính vật liệu polyme dưới tác dụng liều chiếu tích lũy cao.
+-----------------------------------------------------------------------------------+
|                        LUỒNG TÍNH TOÁN VÀ XỬ LÝ DỮ LIỆU                           |
+-----------------------------------------------------------------------------------+
|  [Định nghĩa CTHH & Mật độ rho] ---> [MCNP6: Vận chuyển Photon (NPS=6.0E+09)]      |
|                                                     |                             |
|                                                     v                             |
|  [Phổ thực nghiệm số / Tally F8] <------------------+                             |
|         |                                                                         |
|         v                                                                         |
|  [Colegram: Khớp Gauss + Đa thức bậc 4] ---> [Diện tích đỉnh I, I_0]              |
|                                                     |                             |
|                                                     v                             |
|  [Cơ sở dữ liệu NIST WinXCom] -----------> [Công thức mu_m & Z_eff] ---> [RD (%)] |
+-----------------------------------------------------------------------------------+

Technology Stack & Công cụ sử dụng

  • Engine mô phỏng hạt nhân: MCNP6 (Version 6.1, phát triển bởi Los Alamos National Laboratory), thư viện hạt nhân ENDF/B-VII.0, ACTL, ENDL.
  • Xử lý phổ phóng xạ: Colegram v1.2 (Hỗ trợ giải thuật phi tuyến Levenberg-Marquardt).
  • Cơ sở dữ liệu vật lý chuẩn: WinXCom (NIST Standard Reference Database).
  • Nền tảng phân tích số liệu: Python 3.9 (NumPy, SciPy, Matplotlib) và Microsoft Excel.
+--------------------------------------------------------------------------------+
|                         MA TRẬN TIẾN ĐỘ VÀ RỦI RO DỰ ÁN                        |
+--------------------------------------------------------------------------------+
| Giai đoạn | Nội dung chính                      | Rủi ro phát sinh | Giải pháp |
+-----------+-------------------------------------+------------------+-----------+
| Milestone 1: Khảo sát lý thuyết & CSDL         | Sai lệch mật độ  | Chuẩn hóa |
| Milestone 2: Lập trình tập tin MCNP6 3D        | Lỗi hình học     | Visualizer|
| Milestone 3: Chạy mô phỏng 6 tỷ hạt            | Tràn bộ nhớ/Lâu  | Multi-core|
| Milestone 4: Khớp phổ Colegram & Thẩm định RD% | Lệch phông tán xạ| Gauss+Đa 4|
+--------------------------------------------------------------------------------+

Implementation và kết quả

Development Process & Cấu trúc mã nguồn MCNP6

Để thiết lập mô hình truyền qua gamma, tệp đầu vào MCNP6 được cấu trúc nghiêm ngặt qua 3 khối thẻ: Cell Cards, Surface Cards và Data Cards.

C ===================================================================
C THETAP TIN DAU VAO MCNP6: MO PHONG GAMMA TRUYEN QUA CHO VAT LIEU POLYME
C ===================================================================
C --- KHAI BAO CELL CARDS (O MANG) ---
1    1  -1.105   -1  2 -3         IMP:P=1  $ Bia vat lieu ABS (rho=1.105 g/cm3)
2    2  -0.0012   1  2 -3         IMP:P=1  $ Khong khi xung quanh bia
3    3  -3.670   -4  5 -6         IMP:P=1  $ Tinh the NaI(Tl) 76x76 mm
4    0          (1:-2:3) (4:-5:6) IMP:P=1  $ Vung dem chan doan
99   0           #1 #2 #3         IMP:P=0  $ Vung khong gian ngoai (Kill cell)

C --- KHAI BAO SURFACE CARDS (MAT BIEN HINH HOC) ---
1    CZ  1.0                               $ Tru bia ban kinh R = 1.0 cm
2    PZ  0.0                               $ Mat phang day duoi bia
3    PZ  20.0                              $ Mat phang day tren bia
4    CZ  3.8                               $ Tru tinh the NaI ban kinh R = 3.8 cm
5    PZ  40.43                             $ Mat truoc tinh the NaI
6    PZ  48.03                             $ Mat sau tinh the NaI

C --- KHAI BAO DATA CARDS (NGUON, VAT LIEU, TALLY) ---
SDEF  POS=0 0 -28.43  ERG=0.661657  PAR=2  VEC=0 0 1  DIR=1
C Dinh nghia vat lieu ABS: C15-H17-N (Ty le phan tram khoi luong)
M1    6000  -0.8567  1000  -0.0815  7000  -0.0618
C Dinh nghia Tally F8 ghi nhan pho bien do xung trong tinh the NaI
F8:P  3
E8    0  1000I  1.0
NPS   6000000000

Cơ sở thuật toán trích xuất $Z_{eff}$

  1. Hệ số suy giảm khối ($\mu_m$) từ dữ liệu truyền qua chùm tia hẹp: $$\mu_m = \frac{\mu}{\rho} = \frac{1}{\rho \cdot x} \ln\left(\frac{I_0}{I}\right) \quad (\text{cm}^2/\text{g})$$ Trong đó: $I_0$ là cường độ photon qua không khí, $I$ là cường độ photon qua mẫu polyme dày $x = 2.0\text{ cm}$, $\rho$ là mật độ khối ($g/\text{cm}^3$).

  2. Tiết diện tương tác nguyên tử ($\sigma_a$) và electron ($\sigma_{el}$): $$\sigma_a = \frac{1}{N_A} \sum_{i} f_i A_i (\mu_m)i = \frac{\mu_m}{N_A \sum_i \frac{w_i}{A_i}}$$ $$\sigma{el} = \frac{1}{N_A} \sum_{i} \frac{f_i A_i}{Z_i} (\mu_m)_i$$

  3. Nguyên tử số hiệu dụng ($Z_{eff}$): $$Z_{eff} = \frac{\sigma_a}{\sigma_{el}} = \frac{\sum_i f_i A_i (\mu_m)_i}{\sum_i f_i A_i \frac{(\mu_m)_i}{Z_i}}$$

  4. Khớp phổ phi tuyến trong Colegram: $$Y(x) = \frac{A}{\sigma \sqrt{2\pi}} \exp\left( -\frac{(x - x_0)^2}{2\sigma^2} \right) + \sum_{k=0}^{4} a_k x^k$$

+------------------------------------------------------------------------------------+
|                DỮ LIỆU ĐO VÀ ĐỘ LỆCH HỆ SỐ SUY GIẢM KHỐI TẠI 0.662 MeV             |
+------------------------------------------------------------------------------------+
| Polyme | CTHH       | Mật độ | Cường độ I    | mu_m (MCNP6) | mu_m (NIST) | RD (%) |
+--------+------------+--------+---------------+--------------+-------------+--------+
| ABS    | C15H17N    | 1.105  | 297,714,400   | 0.0732       | 0.0832      | 12.1%  |
| HDPE   | C2H4       | 0.970  | 301,193,400   | 0.0774       | 0.0880      | 12.1%  |
| PMMA   | C5H8O2     | 1.190  | 294,056,500   | 0.0731       | 0.0832      | 12.1%  |
| PC     | C15H16O2   | 1.240  | 293,131,700   | 0.0715       | 0.0825      | 13.4%  |
| CR     | C4H5Cl     | 1.230  | 293,986,300   | 0.0708       | 0.0802      | 11.7%  |
| PEEK   | C19H12O3   | 1.350  | 289,268,800   | 0.0705       | 0.0803      | 12.2%  |
| PET    | C10H8O4    | 1.415  | 286,632,400   | 0.0705       | 0.0803      | 12.2%  |
| PMP    | C6H12      | 0.830  | 307,780,200   | 0.0774       | 0.0880      | 12.1%  |
| POM    | CH2O       | 1.425  | 284,862,300   | 0.0722       | 0.0822      | 12.1%  |
| PP     | C3H6       | 0.930  | 303,061,200   | 0.0774       | 0.0880      | 12.1%  |
| PS     | C8H8       | 1.060  | 299,829,100   | 0.0729       | 0.0830      | 12.1%  |
| PVC    | C2H3Cl     | 1.300  | 291,739,100   | 0.0700       | 0.0790      | 11.4%  |
| PE     | C2H4       | 0.940  | 302,594,000   | 0.0774       | 0.0880      | 12.1%  |
| PVDF   | C2H2F2     | 1.820  | 273,554,800   | 0.0677       | 0.0771      | 12.2%  |
| Air    | N-O-Ar-C   | 0.0012 | 349,960,500   | --           | --          | --     |
+------------------------------------------------------------------------------------+

Kết quả tính toán nguyên tử số hiệu dụng ($Z_{eff}$)

+------------------------------------------------------------------------------------+
|                  SO SÁNH Z_eff GIỮA CÁC PHƯƠNG PHÁP TÍNH TOÁN                      |
+------------------------------------------------------------------------------------+
| Tên Polyme | Mật độ (g/cm3) | Z_eff (Lý thuyết PI) | Z_eff (MCNP6) | Độ lệch RD(%) |
+------------+----------------+----------------------+---------------+---------------+
| ABS        | 1.105          | 3.45                 | 3.04          | 12.1%         |
| HDPE       | 0.970          | 2.67                 | 2.34          | 12.1%         |
| PMMA       | 1.190          | 3.60                 | 3.16          | 12.1%         |
| PC         | 1.240          | 3.70                 | 3.20          | 13.4%         |
| CR         | 1.230          | 4.60                 | 4.06          | 11.6%         |
| PEEK       | 1.350          | 4.41                 | 3.88          | 12.1%         |
| PET        | 1.415          | 4.55                 | 3.99          | 12.1%         |
| PMP        | 0.830          | 2.67                 | 2.34          | 12.1%         |
| POM        | 1.425          | 4.00                 | 3.51          | 12.1%         |
| PP         | 0.930          | 2.67                 | 2.34          | 12.1%         |
| PS         | 1.060          | 3.50                 | 3.08          | 12.1%         |
| PVC        | 1.300          | 5.33                 | 4.72          | 11.4%         |
| PE         | 0.940          | 2.67                 | 2.34          | 12.1%         |
| PVDF       | 1.820          | 5.33                 | 4.68          | 12.3%         |
+------------------------------------------------------------------------------------+

Đổi mới và đóng góp

Đổi mới kỹ thuật

  • Mô hình hóa số nguyên vẹn hệ phổ kế gamma: Tích hợp đầy đủ cấu trúc bộ chuẩn trực chì (lead collimators), lớp phản xạ, vỏ bọc nhôm và tinh thể nhấp nháy $\text{NaI(Tl)}$ $76 \times 76\text{ mm}$ vào MCNP6, tái hiện hiệu ứng tán xạ Compton biên và phông nền thực tế.
  • Quy trình chuẩn hóa khử phông tự động: Kết hợp MCNP6 Tally F8 với thuật toán tách đỉnh Colegram (Gauss + Đa thức bậc 4), giúp giảm thiểu sai số chủ quan khi chọn kênh tích phân diện tích đỉnh.
+------------------------------------------------------------------------------------+
|               ĐỐI CHIẾU KẾT QUẢ VỚI CÁC CÔNG BỐ QUỐC TẾ TIÊU BIỂU                 |
+------------------------------------------------------------------------------------+
| Vật liệu | Z_eff (Lý thuyết) | Z_eff (MCNP6) | Công bố quốc tế khác     | Độ lệch  |
+----------+-------------------+---------------+--------------------------+----------+
| PMMA     | 3.60              | 3.16          | 3.40 (Shantappa [22])    | +5.88%   |
|          |                   |               | 3.52 (Manjunatha [23])   | +2.27%   |
| PC       | 3.70              | 3.20          | 3.726 (Singh [24])       | -0.69%   |
| PP       | 2.67              | 2.34          | 2.707 (Manjunatha [23])  | -1.36%   |
|          |                   |               | 2.66 (Kucuk [4])         | +0.37%   |
| PS       | 3.50              | 3.08          | 3.001 (Singh [24])       | +16.6%   |
| PE       | 2.67              | 2.34          | 2.710 (Singh [24])       | -1.47%   |
+------------------------------------------------------------------------------------+

Đánh giá nguyên nhân sai số hệ thống

Độ lệch tương đối (RD%) giữa mô phỏng MCNP6 và lý thuyết giải tích dao động ổn định trong khoảng $11.4% - 13.4%$. Hiện tượng này xuất phát từ:

  1. Sự khác biệt giữa mật độ khối thực nghiệm tổng hợp và cấu hình thành phần đồng vị mặc định trong thư viện tiết diện hạt nhân ENDF/B-VII.
  2. Hiệu ứng tán xạ nhiều lần (multiple scattering) trong bia dày $2.0\text{ cm}$ làm tăng nhẹ số đếm photon lọt vào cửa sổ năng lượng của Tally F8.

Ứng dụng thực tế và triển khai

+----------------------------------------------------------------------------------+
|                           CÁC LĨNH VỰC ỨNG DỤNG THỰC TIỄN                        |
+----------------------------------------------------------------------------------+
| [Thiết kế che chắn bức xạ]   --> Tối ưu hóa tỷ lệ pha tạp PbO/Bi2O3 vào polyme   |
| [Y học hạt nhân & Xạ trị]    --> Chế tạo phantom mô phỏng mô mềm (Tissue phantom)|
| [Kiểm tra không phá hủy NDT] --> Chuẩn định thiết bị đo mật độ chất lỏng/vật liệu|
| [Công nghệ vật liệu mới]     --> Phát triển bê tông polyme hấp thụ bức xạ gamma  |
+----------------------------------------------------------------------------------+

Triển khai trong công nghiệp và nghiên cứu

  • Chế tạo vật liệu che chắn composite nhẹ: Thay thế chì độc hại bằng polyme chứa nguyên tố $Z$ cao (như PVDF kết hợp phụ gia nano), ứng dụng cho trang phục bảo hộ phóng xạ và vách ngăn phòng chụp X-quang.
  • Tissue-equivalent Phantoms: Polyme có $Z_{eff}$ thấp ($2.3 - 3.6$) như PE, PMMA tương đương mô mềm và mô mỡ người, dùng để hiệu chuẩn liều xạ trị tại các bệnh viện ung bướu.
  • Xác định mật độ chất lỏng: Ứng dụng tỷ số diện tích đỉnh tán xạ để đo nồng độ và mật độ dòng chất lỏng chảy trong đường ống dẫn kín không tiếp xúc.

Phân tích hiệu quả kinh tế (Cost-Benefit & ROI)

  • Tiết kiệm chi phí thực nghiệm: Thiết lập một phòng lab gamma chuẩn yêu cầu nguồn $^{137}\text{Cs}$ hoạt độ cao, phòng chì bảo vệ và đầu dò đắt tiền ($>30,000\text{ USD}$). Ứng dụng mô hình MCNP6 giúp giảm $90%$ chi phí khảo sát ban đầu.
  • Rút ngắn chu kỳ R&D: Thời gian tính toán một vật liệu mới giảm từ 2 tuần thực nghiệm xuống còn 4 giờ mô phỏng trên máy trạm đa nhân.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Năng lượng đơn lẻ: Mô hình chỉ mới khảo sát tại mốc $0.662\text{ MeV}$, chưa quét dải năng lượng thấp ($<100\text{ keV}$) nơi hiệu ứng quang điện chiếm ưu thế vượt trội.
  • Độ phân giải đầu dò: Đầu dò nhấp nháy $\text{NaI(Tl)}$ có độ phân giải năng lượng trung bình ($\sim 7%$ tại $662\text{ keV}$), gây khó khăn cho việc phân tách các đỉnh tán xạ năng lượng gần nhau so với đầu dò bán dẫn siêu tinh khiết HPGe.

Kế hoạch mở rộng và nâng cấp

  1. Mở rộng phổ năng lượng: Bổ sung các nguồn chuẩn quốc tế: $^{241}\text{Am}$ ($59.54\text{ keV}$), $^{133}\text{Ba}$ ($356\text{ keV}$), $^{60}\text{Co}$ ($1.173\text{ MeV}$ và $1.332\text{ MeV}$).
  2. Nâng cấp đầu dò HPGe: Tích hợp mô hình đầu dò High-Purity Germanium trong MCNP6 để tăng độ sắc nét phổ và giảm độ bất định thống kê.
  3. Phát triển GUI tự động hóa: Xây dựng module Python giao tiếp tự động với MCNP6, cho phép người dùng nhập công thức hóa học tùy ý và tự động xuất $Z_{eff}$ cùng đường cong suy giảm $\mu_m(E)$.

Đối tượng hưởng lợi

+-----------------------------------------------------------------------------------+
|                        GIÁ TRỊ MANG LẠI CHO CÁC ĐỐI TƯỢNG                         |
+-----------------------------------------------------------------------------------+
| Sinh viên & Học viên | Nắm vững kỹ thuật lập trình tệp đầu vào MCNP6 và xử lý phổ |
| Kỹ sư An toàn bức xạ | Dữ liệu $Z_{eff}$ chuẩn xác để thiết kế vỏ bọc che chắn    |
| Nhà nghiên cứu Y học | Cơ sở dữ liệu chọn vật liệu phantom tương đương mô sinh học|
| Doanh nghiệp NDT     | Giải pháp số hóa quy trình kiểm định vật liệu không phá hủy|
+-----------------------------------------------------------------------------------+

Câu hỏi thường gặp

1. Cần cấu hình phần cứng tối thiểu nào để chạy mô phỏng MCNP6 với 6 tỷ hạt?

Hệ thống khuyến nghị trang bị CPU đa nhân (tối thiểu 8 cores / 16 threads như Intel Core i7 hoặc AMD Ryzen 7 trở lên), RAM tối thiểu 16GB và hệ điều hành Linux (Ubuntu/CentOS) hoặc Windows 64-bit có cài đặt MPI/OpenMP để kích hoạt tính toán song song, giúp hoàn thành $6 \times 10^9$ photon trong 3 - 5 giờ.

2. Tại sao $Z_{eff}$ của polyme lại là số thập phân thay vì số nguyên?

Khác với nguyên tố hóa học đơn lẻ có số proton cố định ($Z$ nguyên), polyme là hợp chất cao phân tử gồm nhiều nguyên tố (C, H, O, N, Cl, F). $Z_{eff}$ là thông số trung bình trọng số thể hiện xác suất tương tác photon tổng hợp của toàn bộ các electron trong cấu trúc phân tử, do đó mang giá trị thập phân.

3. Làm thế nào để khắc phục hiện tượng trôi đỉnh (peak shift) trong xử lý phổ?

Cần hiệu chuẩn năng lượng (Energy Calibration) trước bằng các đỉnh chuẩn ($511\text{ keV}$, $662\text{ keV}$, $1173\text{ keV}$). Trong Colegram, áp dụng hàm Gauss phi tuyến kết hợp hiệu chỉnh đa thức bậc 4 giúp cố định tọa độ tâm đỉnh $x_0$ và loại bỏ độ dịch do phông nền tích lũy.

4. MCNP6 có thể mô phỏng các vật liệu pha tạp nano (composite) không?

Hoàn toàn được. Người dùng chỉ cần xác định tỷ lệ phần trăm khối lượng (fraction mang dấu âm trong thẻ Mm) hoặc tỷ lệ số nguyên tử của các hạt nano (như $\text{Bi}_2\text{O}_3$, $\text{PbO}$, $\text{WO}_3$) trộn trong nền polyme để MCNP6 gọi đúng thư viện tiết diện vi mô tương ứng.

5. Tại sao có độ chênh lệch RD% khoảng 12% giữa mô phỏng MCNP6 và NIST WinXCom?

WinXCom tính toán dựa trên mô hình lý thuyết hấp thụ thuần túy cho chùm tia chuẩn trực tuyệt đối (zero beam-width). Trong khi đó, mô phỏng MCNP6 mô tả không gian thực tế 3D có xét đến bề rộng chùm tia, tán xạ Compton góc nhỏ lọt vào đầu dò và hiệu ứng kích thước tinh thể $\text{NaI(Tl)}$.

Kết luận

Dự án đã hoàn thành xuất sắc việc nghiên cứu, mô hình hóa và xác định nguyên tử số hiệu dụng ($Z_{eff}$) cùng hệ số suy giảm khối ($\mu_m$) cho 14 loại vật liệu polyme tiêu biểu tại mức năng lượng $0.662\text{ MeV}$ của nguồn $^{137}\text{Cs}$. Bằng việc kết hợp sức mạnh mô phỏng Monte Carlo trên MCNP6, phân tích phổ phi tuyến Colegram và cơ sở dữ liệu đối chuẩn NIST WinXCom, công trình đã thiết lập một quy trình số hóa chuẩn xác, khoa học và có khả năng tái lập cao.

Kết quả thu được chứng minh tính khả thi vượt trội của phương pháp mô phỏng số trong việc thay thế các thử nghiệm bức xạ tốn kém và nguy hiểm, mở ra tiềm năng ứng dụng rộng lớn trong thiết kế vật liệu che chắn bức xạ thế hệ mới, chế tạo mô hình phantom y học hạt nhân và kiểm định công nghiệp NDT. Các nhà nghiên cứu và kỹ sư có thể kế thừa trực tiếp cấu trúc mã nguồn MCNP6 và bộ dữ liệu trong đồ án này để tiếp tục phát triển các hệ thống đo lường hạt nhân tiên tiến.