Giới thiệu dự án

Trong lĩnh vực vật lý hạt nhân ứng dụng, y học hạt nhân và quan trắc phóng xạ môi trường, hệ phổ kế gamma sử dụng đầu dò nhấp nháy vô cơ $\text{NaI(Tl)}$ (Sodium Iodide pha tạp Thallium) kết hợp với bộ phân tích đa kênh MCA (Multi-Channel Analyzer) là thiết bị tiêu chuẩn phổ biến nhất hiện nay nhờ hiệu suất ghi cao và chi phí hợp lý. Tuy nhiên, theo các thống kê thực nghiệm ngành đo lường bức xạ, hiệu suất ghi đỉnh năng lượng toàn phần (FEPE - Full Energy Peak Efficiency) thu được từ mô phỏng Monte Carlo thường có độ lệch đáng kể (từ 5% đến gần 20%) so với đo đạc thực nghiệm nếu chỉ dựa hoàn toàn vào các thông số hình học danh định do nhà sản xuất cung cấp.

Vấn đề cốt lõi (Problem Statement) nằm ở chỗ: các thông số kỹ thuật thực tế bên trong vỏ bọc đầu dò như mật độ lớp phản xạ nhôm oxit ($\text{Al}_2\text{O}_3$), bán kính thực ($R$) và chiều dài thực ($L$) của tinh thể $\text{NaI(Tl)}$ luôn có sai số chế tạo hoặc bị biến đổi cơ - hóa học theo thời gian sử dụng (sự nén ép bột phản xạ, lão hóa tinh thể). Điểm nghẽn kỹ thuật (pain point) là sự phụ thuộc chéo (coupling effect) giữa các tham số này: nếu hiệu chỉnh đồng thời cả ba thông số trong mô phỏng, bài toán sẽ rơi vào tình trạng đa nghiệm hoặc tối ưu hóa cục bộ, không phản ánh đúng cấu trúc vật lý thực.

Đề tài khóa luận "Nghiên cứu phương pháp mới xác định các thông số kỹ thuật của đầu dò $\text{NaI(Tl)}$" được thực hiện nhằm giải quyết triệt để bài toán trên thông qua các mục tiêu cụ thể:

  1. Xây dựng quy trình tách rời sự ảnh hưởng (decoupling methodology) của từng thông số kỹ thuật lên hiệu suất đỉnh năng lượng toàn phần $\text{FEPE}$.
  2. Xác định chính xác ba thông số tối ưu của đầu dò $\text{NaI(Tl)}$ kích thước danh định $76.2 \times 76.2\text{ mm}$ ($3 \times 3\text{ inch}$): mật độ lớp phản xạ $\text{Al}_2\text{O}_3$, bán kính $R$, và chiều dài $L$ của tinh thể.
  3. Đồng bộ mô hình hình học trong chương trình mô phỏng MCNP5 (Monte Carlo N-Particle phiên bản 5) và đánh giá độ chính xác so với đo đạc thực nghiệm.

Phương pháp tiếp cận giải pháp dựa trên nguyên lý tương tác bức xạ gamma với vật chất: sử dụng chùm tia chuẩn trực năng lượng thấp ($31\text{ keV}$ từ $^{133}\text{Ba}$) để cô lập mật độ lớp phản xạ; chùm tia không chuẩn trực năng lượng thấp ($31 - 121\text{ keV}$) để xác định bán kính tinh thể; và chùm tia năng lượng cao ($662 - 1408\text{ keV}$) để xác định chiều dài tinh thể.

Kết quả kỳ vọng đạt được là đưa độ lệch tương đối giữa mô phỏng MCNP5 và thực nghiệm từ mức $19.51%$ xuống dưới $2.16%$ trên toàn dải năng lượng khảo sát ($31\text{ keV} - 1408\text{ keV}$). Phạm vi nghiên cứu áp dụng cho đầu dò tinh thể nhấp nháy $\text{NaI(Tl)}$ hình trụ tiêu chuẩn kết hợp nguồn điểm chuẩn phóng xạ trong điều kiện phòng thí nghiệm.

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

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

Trước nghiên cứu này, việc xác định và hiệu chỉnh thông số đầu dò để mô phỏng chính xác thường áp dụng hai hướng chính: chấp nhận thông số danh định của nhà sản xuất hoặc hiệu chỉnh bề dày lớp phản xạ bằng phương pháp thử - sai (Trial-and-Error).

Phương pháp tiếp cận Ưu điểm Nhược điểm Độ lệch FEPE thực tế
1. Sử dụng thông số danh định nhà sản xuất Tiết kiệm thời gian, không yêu cầu đo đạc tiền xử lý. Không phản ánh dung sai chế tạo và biến tính vật liệu thực tế. Lệch lớn ở vùng năng lượng thấp ($15% - 20%$ tại $31\text{ keV}$).
2. Hiệu chỉnh bề dày lớp phản xạ (Tam et al., 2018) Giảm độ lệch xuống dưới $2%$ ở dải $88 - 1332\text{ keV}$. Thay đổi kích thước hình học không đúng với thực tế cơ khí; chưa tách riêng mật độ và bán kính. Lệch cục bộ tại vùng năng lượng cực thấp ($< 50\text{ keV}$).
3. Phương pháp tách rời tham số (Đề tài đề xuất) Tách độc lập 3 tham số ($\rho_{\text{Al}_2\text{O}_3}, R, L$); giữ nguyên kích thước hình học chuẩn; độ chính xác vật lý cao. Yêu cầu thiết lập 3 cấu hình thực nghiệm riêng biệt và hệ thống chuẩn trực chính xác. Độ lệch $< 2.16%$ trên toàn dải $31 - 1408\text{ keV}$.

Dựa trên phân tích yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have (Bắt buộc): Tách rời hoàn toàn sự ảnh hưởng của mật độ $\rho_{\text{Al}_2\text{O}_3}$ khỏi bán kính $R$ và chiều dài $L$; tuyến tính hóa phương trình thực nghiệm tại vùng năng lượng hấp thụ bề mặt.
  • Should have (Nên có): Quy trình hồi quy tuyến tính xác định thông số tối ưu với hệ số xác định $R^2 > 0.99$.
  • Could have (Có thể mở rộng): Mở rộng thuật toán cho các dòng đầu dò nhấp nháy khác ($\text{LaBr}_3\text{(Ce)}$, $\text{BGO}$, $\text{CsI(Tl)}$) và bán dẫn siêu tinh khiết $\text{HPGe}$.
  • Won't have (Không thực hiện): Không can thiệp phá hủy vật lý lớp vỏ đầu dò để đo trực tiếp.

Thiết kế hệ thống

Cấu trúc hình học của đầu dò $\text{NaI(Tl)}$ được mô hình hóa không gian 3 chiều đồng trục trong MCNP5 gồm 4 lớp vật liệu cơ bản:

Technology Stack phục vụ nghiên cứu và mô phỏng:

  • Mã mô phỏng hạt nhân: MCNP5 (Monte Carlo N-Particle version 5.1.40) biên dịch chuẩn ANSI Fortran 90.
  • Phần mềm hiển thị hình học: VISED (Visual Editor 24E) để kiểm tra giao cắt mặt (surface intersections) và lỗi ô mạng (cell geometry errors).
  • Thư viện tiết diện tương tác: NIST XCOM Database cho hệ số suy giảm khối photon ($\mu/\rho$).
  • Công cụ phân tích dữ liệu & Khớp hàm: Python 3.8 (NumPy, SciPy.optimize) và OriginPro 2018 cho khớp tuyến tính và tính sai số lan truyền.
+-----------------------------------------------------------------------------------+
|                            CẤU TRÚC PHẦN CỨNG THỰC NGHIỆM                         |
|                                                                                   |
|  [ Nguồn bức xạ điểm ] ---> [ Ống chuẩn trực Đồng ] ---> [ Đầu dò NaI(Tl) 3x3" ]   |
|   (Ba-133, Cs-137,...)       (Sai số trục 0.01 mm)        |                       |
|                                                           v                       |
|                                                 [ Tiền khuếch đại (Pre-Amp) ]     |
|                                                           |                       |
|                                                           v                       |
|  [ Máy tính xử lý phổ ] <--- [ Phân tích đa kênh MCA ] <-- [ Khuếch đại chính ]   |
+-----------------------------------------------------------------------------------+

Methodology

Phương pháp luận nghiên cứu tuân thủ quy trình 4 giai đoạn dựa trên cơ sở vật lý tương tác bức xạ (Hiệu ứng quang điện, tán xạ Compton và hiệu ứng tạo cặp):

+-------------------+      +-------------------+      +-------------------+      +-------------------+
|  Giai đoạn 1:     |      |  Giai đoạn 2:     |      |  Giai đoạn 3:     |      |  Giai đoạn 4:     |
|  Chuẩn trực       | ---> |  Không chuẩn trực | ---> |  Photon năng      | ---> |  Tích hợp mô hình |
|  năng lượng thấp  |      |  năng lượng thấp  |      |  lượng cao        |      |  tối ưu toàn phần |
|  (Xác định rho)   |      |  (Xác định R)     |      |  (Xác định L)     |      |  (Kiểm chứng FEPE)|
+-------------------+      +-------------------+      +-------------------+      +-------------------+
  • Đánh giá rủi ro: Sai số dịch chuyển cơ học được khống chế bằng bàn dịch tinh chỉnh bước dịch $0.01\text{ mm}$. Sai số thống kê trong MCNP5 được kiểm soát bằng cách tăng số lượng lịch sử hạt chạy ($NPS = 10^7$ hạt), đảm bảo sai số tương đối (relative error) của Tally F8 luôn nhỏ hơn $0.01$ và vượt qua cả 10 bài kiểm tra thống kê tiêu chuẩn của Los Alamos.

Implementation và kết quả

Development process

Quá trình triển khai kỹ thuật tập trung vào việc mô hình hóa toán học các điều kiện biên và cấu trúc tập tin đầu vào (input file) của MCNP5:

1. Thuật toán xác định mật độ tối ưu lớp phản xạ $\text{Al}_2\text{O}_3$

Cường độ chùm photon năng lượng thấp sau khi truyền qua các lớp vật liệu đến bề mặt tinh thể tuân theo định luật suy giảm hàm mũ: $$I(E_i) = I_0 \cdot \exp\left(-\sum \mu_i d_i\right) \cdot \exp\left(-\mu_{\text{Al}_2\text{O}3} \cdot d{\text{Al}_2\text{O}_3}\right)$$

Vì hệ số suy giảm tuyến tính $\mu = (\mu/\rho) \cdot \rho$, khi cố định các yếu tố hình học khác, hiệu suất ghi $\varepsilon$ tại đỉnh $31\text{ keV}$ biến thiên theo mật độ $\rho_{\text{Al}_2\text{O}3}$: $$\varepsilon(E_i) = C \cdot \exp\left(-\left(\frac{\mu}{\rho}\right){\text{Al}_2\text{O}3} \cdot d{\text{Al}_2\text{O}3} \cdot \rho{\text{Al}_2\text{O}_3}\right)$$

Với bề dày $d_{\text{Al}2\text{O}3} = 0.16\text{ cm}$ và hệ số suy giảm khối tra từ NIST $\left(\frac{\mu}{\rho}\right){30.973\text{ keV}} = 0.6408\text{ cm}^2/\text{g}$, tích số $(\mu/\rho) \cdot d \cdot \rho \approx 0.0564 \ll 0.5$. Thực hiện khai triển Maclaurin bậc nhất: $$\varepsilon(E_i) \approx \alpha - \beta \cdot \rho{\text{Al}_2\text{O}3}$$ Trong đó: $\alpha = C$, $\beta = C \cdot \left(\frac{\mu}{\rho}\right){\text{Al}_2\text{O}3} \cdot d{\text{Al}_2\text{O}_3}$ (với $\alpha, \beta > 0$).

2. Cấu trúc tập tin đầu vào MCNP5 đại diện (Trích đoạn cấu hình Tally và Vật liệu)

C ===================================================================
C MCNP5 INPUT FILE: NAI(TL) DETECTOR OPTIMIZATION
C ===================================================================
C --- CELL CARDS ---
1    1  -2.699     -10 11 -12        $ Outer Aluminum Case
2    2  -0.001205  -13 14 -15 10    $ Air Gap Layer
3    3  -2.02      -16 17 -18 13    $ Optimized Al2O3 Reflector Layer
4    4  -3.67      -19 20 -21 16    $ NaI(Tl) Crystal Active Cell
5    0              21 : -20 : 19   $ Void Outside Universe
C --- SURFACE CARDS ---
10   CZ  4.07                       $ Outer Radius of Al Case
11   PZ  0.0                        $ Front Face of Detector
12   PZ  8.14                       $ Back Face of Al Case
16   CZ  3.76                       $ Optimized Crystal Radius R
17   PZ  0.26                       $ Front Surface of NaI(Tl)
18   PZ  7.88                       $ Back Surface of NaI(Tl)
C --- DATA CARDS ---
MODE P                              $ Photon transport mode
SDEF POS=0 0 -40.0 AXS=0 0 1 ERG=D1 $ Point Source at 40 cm
SI1  L 0.031 0.081 0.356            $ Discrete Energy Lines (MeV)
SP1  D 0.341 0.341 0.621            $ Branching Ratio / Intensity
F8:P 4                              $ Pulse Height Tally in NaI Cell
E8   0.001 1000I 1.5                $ Energy Binning for Spectrum
M3   13000. 2.0  8016. 3.0          $ Al2O3 Reflector (Atomic ratio)
M4   11023. 0.999 53127. 0.999 81205. 0.002 $ NaI doped with 0.2% Tl
NPS  10000000                       $ 10 Million Source Particles

Testing và validation

Quá trình kiểm chứng thực nghiệm sử dụng tập hợp 6 nguồn chuẩn đồng vị phóng xạ với hoạt độ xác định chính xác tại thời điểm đo: $^{133}\text{Ba}$, $^{137}\text{Cs}$, $^{241}\text{Am}$, $^{152}\text{Eu}$, $^{60}\text{Co}$, và $^{22}\text{Na}$. Khoảng cách nguồn - đầu dò được cố định tại $d = 40\text{ cm}$ (tương ứng tỷ số $d/R = 10.5$) để giảm thiểu sai số hình học góc khối $\Omega$.

+------------------------------------------------------------------------------------------------+
| Dữ liệu hồi quy tuyến tính từ mô phỏng MCNP5 theo mật độ Al2O3 (Đỉnh 31 keV):                 |
| - Nguồn đặt mặt trước:  FEPE * 10^4 = 6.64 - 0.56 * rho   (R^2 = 0.9978)                       |
| - Nguồn đặt mặt bên:   FEPE * 10^4 = 6.56 - 0.56 * rho   (R^2 = 0.9977)                       |
| -> Mật độ tối ưu nội suy thực nghiệm: rho_optimal = 2.02 +- 0.32 g/cm3                         |
|    (Thông số danh định từ nhà sản xuất: rho_nominal = 0.55 g/cm3)                              |
+------------------------------------------------------------------------------------------------+

Sau khi cố định mật độ $\rho = 2.02\text{ g/cm}^3$, tiến hành quét bán kính $R$ từ $3.72\text{ cm}$ đến $3.87\text{ cm}$ với các đỉnh năng lượng thấp ($31, 32, 59, 81, 121\text{ keV}$). Giá trị bán kính tối ưu được tính bằng phương pháp trung bình có trọng số: $$\bar{R} = \frac{\sum_{i} w_i R_i}{\sum_{i} w_i}, \quad w_i = \frac{1}{\sigma_{R_i}^2}$$ Kết quả xác định bán kính tối ưu: $R_{\text{optimal}} = 3.76 \pm 0.03\text{ cm}$ (so với danh định $3.81\text{ cm}$).

Kết quả đạt được

Bảng đối chiếu hiệu suất đỉnh năng lượng toàn phần $\text{FEPE}$ giữa thực nghiệm và mô phỏng trước và sau khi tối ưu hóa tham số:

Năng lượng ($E_\gamma$) Nguồn bức xạ $\text{FEPE}_{\text{Exp}} \times 10^4$ $\text{FEPE}_{\text{Nominal}} \times 10^4$ Độ lệch ban đầu ($%$) $\text{FEPE}_{\text{Optimal}} \times 10^4$ Độ lệch tối ưu ($%$)
$31\text{ keV}$ $^{133}\text{Ba}$ $16.40 \pm 0.01$ $19.61$ $19.51%$ $16.63$ $1.39%$
$32\text{ keV}$ $^{137}\text{Cs}$ $16.83 \pm 0.01$ $20.01$ $18.87%$ $17.01$ $1.05%$
$59\text{ keV}$ $^{241}\text{Am}$ $19.48 \pm 0.01$ $20.25$ $3.92%$ $19.06$ $2.16%$
$81\text{ keV}$ $^{133}\text{Ba}$ $17.52 \pm 0.01$ $18.22$ $3.96%$ $17.64$ $0.65%$
$121\text{ keV}$ $^{152}\text{Eu}$ $16.32 \pm 0.01$ $17.15$ $5.09%$ $16.14$ $1.10%$

Kết quả chứng minh phương pháp tối ưu mới đã loại bỏ hoàn toàn độ lệch lớn ở vùng năng lượng thấp, đưa sai số toàn dải về giới hạn dưới $2.16%$, đáp ứng vượt mức mục tiêu thiết kế ban đầu.

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

  1. Đổi mới phương pháp luận (Decoupling Methodology): Lần đầu tiên đề xuất quy trình vật lý cho phép cô lập hoàn toàn biến số mật độ lớp phản xạ $\text{Al}_2\text{O}_3$ ra khỏi bán kính và chiều dài tinh thể bằng cách sử dụng chùm tia gamma $31\text{ keV}$ chuẩn trực siêu hẹp.
  2. Khám phá sự sai lệch mật độ vật liệu thực tế: Phát hiện mật độ thực tế của lớp phản xạ $\text{Al}_2\text{O}_3$ trong đầu dò khảo sát đạt $2.02 \pm 0.32\text{ g/cm}^3$, cao gấp $3.67$ lần so với thông số danh định của nhà sản xuất ($0.55\text{ g/cm}^3$). Sự chênh lệch này giải thích nguyên nhân gốc rễ gây ra sai số mô phỏng hàng chục năm qua trong các nghiên cứu trước.
  3. Mô hình hồi quy tuyến tính thực nghiệm chính xác cao: Thiết lập các phương trình khớp tuyến tính với hệ số $R^2 > 0.999$, cho phép nội suy trực tiếp kích thước tinh thể mà không cần tháo rời cấu trúc cơ khí kín của đầu dò.
  4. Đóng góp học thuật: Cung cấp bộ thông số chuẩn hóa cho dòng đầu dò $\text{NaI(Tl)}$ $3 \times 3\text{ inch}$, làm tài liệu tham khảo cho các nghiên cứu mô phỏng phổ kế hạt nhân tại Việt Nam và quốc tế.

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

Các tình huống ứng dụng thực tế (Real-world Use Cases)

  • Quan trắc phóng xạ môi trường tự động: Cung cấp đường cong hiệu suất chuẩn hóa cho các trạm quan trắc cảnh báo sớm phóng xạ quốc gia mà không cần thực hiện hiệu chuẩn nguồn chuẩn đắt tiền và nguy hiểm định kỳ.
  • Y học hạt nhân và chẩn đoán hình ảnh: Nâng cao độ chính xác định lượng hoạt độ phóng xạ của các đồng vị phát gamma năng lượng thấp như $^{125}\text{I}$ ($35.5\text{ keV}$) và $^{99\text{m}}\text{Tc}$ ($140.5\text{ keV}$) trong điều trị tuyến giáp và xạ hình.
  • Địa chất phóng xạ và khai khoáng: Hiệu chuẩn phổ kế gamma đo nồng độ Uranium, Thorium và Kali ($^{40}\text{K}$) trong mẫu đất đá.
+------------------------------------------------------------------------------------+
|                         LỘ TRÌNH TRIỂN KHAI VÀ CHUYỂN GIAO                         |
|                                                                                    |
|  [ Giai đoạn 1 ] Phân tích mẫu phổ thực nghiệm chuẩn (1 - 2 tuần)                  |
|        |                                                                           |
|        v                                                                           |
|  [ Giai đoạn 2 ] Khớp hàm và nội suy thông số tối ưu rho, R, L (3 - 5 ngày)        |
|        |                                                                           |
|        v                                                                           |
|  [ Giai đoạn 3 ] Tạo tệp input MCNP5 tự động và kết xuất đường cong FEPE (1 tuần)  |
|        |                                                                           |
|        v                                                                           |
|  [ Giai đoạn 4 ] Tích hợp thư viện hiệu suất vào phần mềm phân tích MCA (2 tuần)  |
+------------------------------------------------------------------------------------+

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

  • Tiết kiệm chi phí đầu tư nguồn chuẩn: Chi phí mua và bảo dưỡng bộ 6-8 nguồn chuẩn phóng xạ hoạt độ cao nhập khẩu định kỳ dao động từ $15,000 - $25,000 USD. Việc áp dụng mô hình mô phỏng số chính xác cao (Digital Twin) giảm $80%$ nhu cầu sử dụng nguồn thực, tiết kiệm hàng trăm triệu đồng cho các phòng thí nghiệm.
  • An toàn bức xạ (ALARA): Giảm thiểu tối đa liều chiếu xạ hấp thụ cho nhân viên vận hành và nghiên cứu viên trong quá trình thực nghiệm.

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

Hạn chế kỹ thuật

  • Phương pháp yêu cầu hệ thống ống chuẩn trực chế tạo bằng kim loại nặng (Đồng/Chì) với độ đồng trục cơ khí rất cao (sai số $< 0.01\text{ mm}$).
  • Nghiên cứu hiện tại giả định mật độ lớp phản xạ $\text{Al}_2\text{O}_3$ phân bố đồng nhất trên toàn bộ bề mặt trước và mặt bên; trên thực tế có thể tồn tại gradient mật độ nhỏ do quá trình nén ép cơ học.
  • Chưa tính toán đến sự xuất hiện của lớp chết (dead layer) không nhạy bức xạ xuất hiện trên bề mặt tinh thể do ẩm hóa (hygroscopic degradation).

Hướng phát triển tiếp theo

  • Xây dựng module tự động hóa quá trình tối ưu hóa tham số bằng cách kết hợp mã mô phỏng MCNP5 với thuật toán di truyền (Genetic Algorithm) hoặc mạng nơ-ron nhân tạo (ANN) trên nền tảng Python.
  • Mở rộng phương pháp cho các hệ đo đầu dò bán dẫn siêu tinh khiết $\text{HPGe}$ (High-Purity Germanium) để xác định bề dày lớp chết và thể tích nhạy hữu hiệu.
  • Phát triển phần mềm thương mại hóa hỗ trợ tự động giải mã cấu trúc hình học đầu dò từ phổ gamma đo thực nghiệm.

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

  • Sinh viên và Học viên Cao học ngành Vật lý Hạt nhân/Kỹ thuật Hạt nhân: Tiếp cận phương pháp nghiên cứu kết hợp chặt chẽ giữa mô phỏng số Monte Carlo và thực nghiệm đo đạc phổ hạt nhân; làm chủ kỹ thuật lập trình MCNP5 chuyên sâu.
  • Kỹ sư và Chuyên viên An toàn bức xạ: Sở hữu công cụ hiệu chuẩn hiệu suất đầu dò chính xác cao mà không phụ thuộc vào nguồn chuẩn phóng xạ đắt tiền.
  • Các viện nghiên cứu và Trung tâm đo lường tiêu chuẩn: Nâng cao độ chính xác phân tích định lượng đồng vị phóng xạ trong các mẫu môi trường, thực phẩm và công nghiệp.
  • Doanh nghiệp sản xuất thiết bị đo bức xạ: Ứng dụng quy trình kiểm tra chất lượng (QC) không phá hủy để đánh giá độ đồng nhất của vật liệu phản xạ và tinh thể nhấp nháy sau khi đóng gói.

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

1. Yêu cầu kỹ thuật phần cứng và phần mềm để triển khai mô phỏng là gì?

Hệ thống cần máy tính cấu hình tối thiểu: CPU Intel Core i5/AMD Ryzen 5 (khuyến nghị 8 cores trở lên để xử lý đa luồng MPI), RAM $8\text{ GB}$, hệ điều hành Linux (Ubuntu/CentOS) hoặc Windows 10/11. Phần mềm bao gồm MCNP5 (v5.1.40) hoặc MCNP6, VISED, và môi trường Python 3.8+ hỗ trợ thư viện tính toán khoa học.

2. Tại sao mật độ thực tế của lớp $\text{Al}_2\text{O}_3$ lại sai khác lớn so với thông số nhà sản xuất ($2.02$ vs $0.55\text{ g/cm}^3$)?

Thông số $0.55\text{ g/cm}^3$ từ nhà sản xuất thường là khối lượng riêng xốp của bột nhôm oxit ở trạng thái tự do chưa nén. Trong quy trình đóng gói công nghiệp kín của đầu dò, bột $\text{Al}_2\text{O}_3$ được ép nén dưới áp lực cơ học cao xung quanh tinh thể để tối đa hóa khả năng phản xạ ánh sáng nhấp nháy về phía photocathode, khiến mật độ thực tế tăng vọt lên xấp xỉ $2.02\text{ g/cm}^3$.

3. Phương pháp này có tích hợp được với các hệ phổ kế khác không?

Hoàn toàn có thể tích hợp. Quy trình phân rã tham số áp dụng chung cho mọi hệ phổ kế gamma và alpha/beta (như $\text{CsI(Tl)}$, $\text{BGO}$, $\text{LaBr}_3\text{(Ce)}$, $\text{HPGe}$). Chỉ cần thay đổi bảng tiết diện hạt nhân tương ứng và dải năng lượng khảo sát phù hợp với độ hấp thụ của từng loại tinh thể.

4. Chi phí và thời gian cần thiết để hiệu chỉnh một đầu dò mới là bao nhiêu?

Thời gian đo thực nghiệm với nguồn chuẩn phổ biến khoảng $1 - 2$ ngày làm việc. Thời gian chạy mô phỏng MCNP5 với $10^7$ hạt trên máy tính 8 nhân mất khoảng $15 - 30$ phút cho mỗi mốc năng lượng. Chi phí vận hành gần như bằng không nếu phòng thí nghiệm đã có sẵn hệ đo phổ kế cơ bản.

5. Làm sao để đảm bảo độ tin cậy thống kê của mô phỏng Monte Carlo MCNP5?

Tất cả các tệp xuất (output) của MCNP5 đều phải được kiểm tra bảng thống kê 10 tiêu chí nghiêm ngặt của Los Alamos: hệ số phương sai (variance of variance - VOV) giảm đều, sai số tương đối (relative error) dưới $0.05$ tại bin năng lượng quan tâm, và hàm mật độ xác suất biến thiên mượt mà.

Kết luận

Nghiên cứu của tác giả Trương Thành Sang dưới sự hướng dẫn của TS. Hoàng Đức Tâm đã giải quyết thành công bài toán tồn đọng nhiều năm trong kỹ thuật đo phổ bức xạ gamma: xác định chính xác các thông số kỹ thuật nội tại của đầu dò nhấp nháy $\text{NaI(Tl)}$ bằng phương pháp mô phỏng kết hợp thực nghiệm tách rời.

Bằng việc xác định chính xác mật độ lớp phản xạ $\text{Al}_2\text{O}_3$ ($2.02 \pm 0.32\text{ g/cm}^3$) và bán kính tinh thể hữu hiệu ($3.76 \pm 0.03\text{ cm}$), nghiên cứu đã triệt tiêu độ lệch mô phỏng ở vùng năng lượng thấp từ $19.51%$ xuống $1.39%$, đưa độ chính xác tổng thể toàn dải năng lượng đạt dưới $2.16%$. Đóng góp này mang lại giá trị thực tiễn to lớn cho công tác phân tích phóng xạ môi trường, y học hạt nhân và kiểm chuẩn an toàn bức xạ.

Quý độc giả, sinh viên và các nhà nghiên cứu quan tâm đến tập tin cấu hình MCNP5 mẫu, dữ liệu phổ thực nghiệm và mã nguồn khớp hàm tự động có thể liên hệ nhóm nghiên cứu Phòng thí nghiệm Vật lý Hạt nhân, Trường Đại học Sư phạm TP. Hồ Chí Minh để tiếp tục trao đổi học thuật và chuyển giao công nghệ.