Giới thiệu dự án

Bối cảnh và thực trạng đo lường bức xạ môi trường

Trong phân tích an toàn bức xạ và đánh giá phóng xạ môi trường (chiếm hơn 65% khối lượng công việc kiểm định hạt nhân theo thống kê của IAEA), việc định lượng chính xác hoạt độ phóng xạ của các đồng vị tự nhiên ($^{238}\text{U}$, $^{232}\text{Th}$, $^{40}\text{K}$) và nhân tạo ($^{137}\text{Cs}$) trong mẫu đất đá, trầm tích và vật liệu xây dựng là yêu cầu bắt buộc. Kỹ thuật phổ kế gamma sử dụng đầu dò bán dẫn Germanium siêu tinh khiết (High Purity Germanium - HPGe) là phương pháp tiêu chuẩn nhờ độ phân giải năng lượng vượt trội (FWHM $\le 2.0\text{ keV}$ tại $1332\text{ keV}$).

Tuy nhiên, nút thắt kỹ thuật lớn nhất trong phân tích phổ gamma mẫu khối thể tích lớn nằm ở việc chuẩn hóa Hiệu suất ghi đỉnh năng lượng toàn phần (Full Energy Peak Efficiency - FEPE). Hiệu suất đỉnh $\varepsilon(E)$ biến thiên phi tuyến theo dải năng lượng (từ $46.5\text{ keV}$ đến $2614.5\text{ keV}$) và phụ thuộc mật thiết vào:

  • Hiệu ứng tự hấp thụ (Self-absorption) phụ thuộc vào số nguyên tử hiệu dụng $Z_{\text{eff}}$ và mật độ khối $\rho$ của chất nền (matrix).
  • Hiệu ứng trùng phùng tổng (True Coincidence Summing - TCS) của các nhân phóng xạ phát nhiều tia gamma theo chuỗi tầng (cascade emission).
  • Sự suy giảm thể tích nhạy do biến đổi bề dày lớp chết (dead layer) của tinh thể Ge theo thời gian vận hành.
       [Photon Gamma (46.5 - 2614.5 keV)]
             [Tín hiệu xung FEPE]

Mục tiêu của dự án

  1. Xác định thực nghiệm FEPE: Thiết lập đường cong hiệu suất đỉnh năng lượng toàn phần trên hệ phổ kế HPGe GEM50P4-83 cho 4 mẫu chuẩn chuẩn hóa quốc tế: IAEA-RGU-1, IAEA-RGTh-1, IAEA-434, IAEA-447 trong dải năng lượng $46.5 - 2614.5\text{ keV}$.
  2. Xây dựng mô hình Monte Carlo MCNP6: Thiết kế tệp đầu vào 3D chuẩn xác từng milimet cho cấu trúc đầu dò HPGe, buồng chì che chắn và thành phần hóa học chi tiết của 4 mẫu chuẩn.
  3. Mô phỏng 3D tương tác photon: Chạy mô phỏng vận chuyển bức xạ với số lượng lịch sử hạt $NPS = 10^9$ hạt để thu thập phổ năng lượng mất mát qua Tally F8.
  4. Đánh giá và đối chuẩn sai lệch: Làm khớp giải tích đa thức bậc 5 $\ln(\varepsilon) = \sum_{i=0}^5 B_i (\ln E)^i$, phân tích độ lệch tương đối ($RD%$) giữa thực nghiệm và mô phỏng, định lượng ảnh hưởng của ma trận nguyên tố và lớp chết tinh thể.

Giải pháp và kết quả kỳ vọng

  • Phương pháp tiếp cận: Kết hợp đo thực nghiệm phông thấp với mô phỏng thống kê Monte Carlo bằng bộ mã MCNP6.
  • Chỉ số đo lường: Độ lệch khớp hàm thực nghiệm $RD \le 8.50%$ trên đa số các đỉnh; xác định nguyên nhân sai khác tại vùng năng lượng thấp ($<100\text{ keV}$) và vùng trùng phùng cascade ($>300\text{ keV}$).
  • Phạm vi nghiên cứu: Cấu hình hình học mẫu dạng hộp trụ ($\varnothing 7.3\text{ cm} \times 2.0\text{ cm}$), đặt áp sát bề mặt đầu dò (khoảng cách $0\text{ mm}$); dải mật độ mẫu từ $0.91\text{ g/cm}^3$ đến $1.55\text{ g/cm}^3$.

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

Phân tích hiện trạng các phương pháp chuẩn hiệu suất

Tiêu chí so sánh Chuẩn thực nghiệm nguồn chuẩn Phương pháp bán giải tích (Semi-empirical) Phương pháp Monte Carlo (MCNP6)
Chi phí vật tư Rất đắt đỏ (mua nhiều mẫu chuẩn IAEA) Thấp (sử dụng phần mềm thương mại) Rất thấp (mô phỏng số trên máy tính)
Tính linh hoạt hình học Rất kém (chỉ đo đúng hộp mẫu chuẩn) Trung bình (giới hạn một số hình dạng cơ bản) Không giới hạn (mọi dạng nguồn, Marinelli, trụ)
Bù trừ hiệu ứng ma trận Khó tùy biến thành phần nguyên tố Sai số lớn khi ma trận mẫu phức tạp Rất cao (khai báo chi tiết bảng thành phần XRF)
Thời gian thiết lập Tốn thời gian đo (86,400s - 172,800s/mẫu) Nhanh ($< 5$ phút tính toán) Vừa phải (phụ thuộc năng lực tính toán CPU)
Mức độ chính xác Chuẩn cơ sở thực tế Độ lệch $10 - 25%$ Độ khớp cao khi tinh chỉnh đúng lớp chết

Phân tích yêu cầu theo mô hình MoSCoW

  • Must Have: Mô hình hóa chính xác thông số hình học đầu dò GEM50P4-83; khai báo chi tiết thành phần nguyên tố 4 mẫu chuẩn (IAEA-RGU-1, IAEA-RGTh-1, IAEA-434, IAEA-447); tính toán diện tích đỉnh bằng phần mềm chuyên dụng; làm khớp hàm hiệu suất.
  • Should Have: Trừ phông phổ theo tỷ lệ thời gian đo thực; tính sai số truyền thống kê Poisson kết hợp sai số hoạt độ nguồn và xác suất phát bức xạ gamma ($I_\gamma$).
  • Could Have: Mô phỏng sự biến thiên của bề dày lớp chết mặt ngoài và mặt trong tinh thể Ge để tối ưu hóa mô hình.
  • Won't Have (Giai đoạn này): Tự động hóa tích hợp hiệu chỉnh trùng phùng tổng (TCS) bằng mã nguồn lai trực tiếp trong MCNP6.

Thiết kế kiến trúc hệ thống và công nghệ

Bảng thông số công nghệ và phiên bản (Technology Stack)

  • Hệ thu nhận phổ: Đầu dò HPGe GEM50P4-83 (ORTEC), hiệu suất tương đối 50%, làm lạnh bằng Nitơ lỏng ($77\text{ K}$).
  • Khối xử lý số tín hiệu: Digital Signal Analyzer Lynx DSA (Canberra/Mirion), tích hợp nguồn cao thế, tiền khuếch đại, bộ ADC 16,384 kênh.
  • Phần mềm xử lý phổ: Genie 2K (Canberra v3.4), Colegram (v1.2).
  • Mã nguồn mô phỏng hạt nhân: MCNP6 (Monte Carlo N-Particle v6.1 / v6.2, Los Alamos National Laboratory).
  • Phần mềm xử lý số liệu & Làm khớp: OriginPro (OriginLab v9.5/2018).

Thông số hình học chi tiết đầu dò HPGe GEM50P4-83

  • Đường kính tinh thể Ge: $65.9\text{ mm}$.
  • Chiều dài tinh thể Ge: $77.0\text{ mm}$.
  • Đường kính hốc tinh thể (Core hole): $11.5\text{ mm}$; Chiều sâu hốc: $64.9\text{ mm}$.
  • Bề dày lớp chết (Dead layer) mặt ngoài tinh thể: $0.7\text{ mm}$ (Ge).
  • Bề dày lớp chết mặt trong tinh thể: $0.0003\text{ mm}$ ($0.3\ \mu\text{m}$ Ge).
  • Bề dày vỏ bao bọc đầu dò (Endcap): $1.0\text{ mm}$ nhôm (Al); Khoảng cách đỉnh tinh thể đến vỏ: $5.0\text{ mm}$.

Implementation và kết quả

Cấu trúc dữ liệu đầu vào mô phỏng MCNP6

Mô hình mô phỏng được cấu trúc thành 3 khối thẻ (Cards) tiêu chuẩn trong tập tin đầu vào HPGe_GEM50P4_83.inp:

C ========================================================
C BLOCK 1: CELL CARDS - KHAI BAO O MANG HINH HOC
C ========================================================
1  1 -5.323  (-26 5 -27):(27 -24 5):(-5 6 29 -16):(-6 7 15 -16) IMP:P,E=1 $ Active Ge
2  2 -2.700  (1 -2 3 -4)                                        IMP:P,E=1 $ Al Cap
3  3 -1.550  (10 -11 12 -13)                                     IMP:P,E=1 $ Sample IAEA-RGU-1
4  0         #1 #2 #3 (-100)                                     IMP:P,E=1 $ Air Void
5  0         (100)                                               IMP:P,E=0 $ Graveyard

C ========================================================
C BLOCK 2: SURFACE CARDS - KHAI BAO MAT GIOI HAN
C ========================================================
1   PZ  0.0          $ Mat phang day vo Al
2   PZ  7.7          $ Mat phang dinh vo Al
3   CZ  3.5          $ Mat tru vo Al
5   PZ  0.5          $ Mat phang day tinh the Ge
6   PZ  8.2          $ Mat phang dinh tinh the Ge
7   CZ  3.295        $ Mat tru ngoai tinh the Ge (R=32.95mm)
15  CZ  0.575        $ Mat tru hoc tinh the Ge (R=5.75mm)
16  PZ  6.99         $ Mat phang day hoc tinh the Ge
100 SO  50.0         $ Mat cau the gioi ao

C ========================================================
C BLOCK 3: DATA CARDS - NGUON, VAT LIEU VA GHI NHAN TALLY
C ========================================================
MODE P E
SDEF  CEL=3 POS=0 0 8.5 RAD=D1 EXT=D2 AXS=0 0 1 ERG=D3 PAR=2
SI1   0.0 3.65       $ Ban kinh hop mau R=3.65cm
SP1   -21 1
SI2   0.0 2.00       $ Be day mau H=2.00cm
SP2   0   1
C Khai bao vat lieu mau IAEA-RGU-1 (Ti le phan tram khoi luong)
M3    6000  -0.0001  $ Carbon
      8016  -0.5316  $ Oxygen
      11023 -0.0002  $ Sodium
      12000 -0.0020  $ Magnesium
      13027 -0.0010  $ Aluminum
      14000 -0.4660  $ Silicon
F8:P  1              $ Tally ghi nhan nang luong mat mat trong cell Ge
E8    0.040 1000I 3.0
NPS   1000000000     $ 1 ty hat Monte Carlo

Thuật toán xác định hiệu suất và sai số lan truyền

Hiệu suất ghi nhận đỉnh năng lượng toàn phần thực nghiệm $\varepsilon_{\text{exp}}(E)$ được tính toán theo biểu thức giải tích:

$$\varepsilon_{\text{exp}}(E) = \frac{N_p(E)}{A \cdot t \cdot I_\gamma(E) \cdot m}$$

Trong đó:

  • $N_p(E)$: Diện tích đếm thực của đỉnh năng lượng sau khi trừ phông bằng Colegram.
  • $A$: Hoạt độ phóng xạ riêng của đồng vị chuẩn tại thời điểm đo ($\text{Bq/kg}$).
  • $t$: Thời gian đo sống (Live time) ($\text{s}$).
  • $I_\gamma(E)$: Xác suất phát bức xạ gamma của vạch năng lượng tương ứng.
  • $m$: Khối lượng mẫu đo ($\text{kg}$).

Đối với đỉnh chồng chập $186.2\text{ keV}$ giữa $^{226}\text{Ra}$ ($186.21\text{ keV}$, $I_\gamma = 3.555%$) và $^{235}\text{U}$ ($185.71\text{ keV}$, $I_\gamma = 57.2%$), hiệu suất được tách chiết theo tỷ lệ đóng góp thành phần:

$$\varepsilon(186.2\text{ keV}) = \frac{N_p(186.2)}{A(^{226}\text{Ra}) \cdot t \cdot I_\gamma(186.2) \cdot m} \times 0.572$$

Sai số tương đối của hiệu suất $\sigma_\varepsilon$ được xác định thông qua vi phân sai số toàn phần theo phân bố Poisson:

$$\left(\frac{\sigma_\varepsilon}{\varepsilon}\right)^2 = \left(\frac{\sigma_A}{A}\right)^2 + \left(\frac{\sigma_{N_p}}{N_p}\right)^2 + \left(\frac{\sigma_{I_\gamma}}{I_\gamma}\right)^2$$

Dữ liệu thực nghiệm và kết quả đối chuẩn (Benchmark Metrics)

1. Thông số hoạt độ và thời gian đo các mẫu chuẩn IAEA

Tên mẫu chuẩn Khối lượng $m$ (g) Mật độ $\rho$ ($\text{g/cm}^3$) Đồng vị chính khảo sát Hoạt độ $A$ ($\text{Bq/kg}$) Thời gian đo $t$ (s)
IAEA-RGU-1 $130.14$ $1.55$ $^{238}\text{U}$ chuỗi cân bằng $4940 \pm 30$ $86,400$
IAEA-RGTh-1 $118.60$ $1.42$ $^{232}\text{Th}$ chuỗi cân bằng $3250 \pm 90$ $86,400$
IAEA-434 $76.24$ $0.91$ $^{226}\text{Ra}$, $^{238}\text{U}$, $^{137}\text{Cs}$ $780 \pm 62$ ($^{226}\text{Ra}$) $172,800$
IAEA-447 $107.70$ $1.28$ $^{212}\text{Pb}$, $^{214}\text{Bi}$, $^{40}\text{K}$, $^{137}\text{Cs}$ $37.0 \pm 1.5$ ($^{212}\text{Pb}$) $172,800$

2. Kết quả đối sánh Hiệu suất đỉnh thực nghiệm, Làm khớp và Mô phỏng (Mẫu IAEA-RGU-1)

Độ lệch tương đối giữa giá trị thực nghiệm và mô phỏng được tính theo:

$$RD(%) = \frac{|\varepsilon_{\text{exp}} - \varepsilon_{\text{sim}}|}{\varepsilon_{\text{exp}}} \times 100$$

Chuỗi/Đồng vị Năng lượng $E$ (keV) Hiệu suất thực nghiệm $\varepsilon_{\text{exp}}$ Hiệu suất làm khớp $\varepsilon_{\text{fit}}$ Sai số làm khớp $RD_{\text{fit}}$ (%) Hiệu suất MCNP6 $\varepsilon_{\text{sim}}$ Sai số mô phỏng $RD_{\text{sim}}$ (%)
$^{210}\text{Pb}$ $46.54$ $0.0170 \pm 0.0002$ $0.0174$ $2.72$ $0.0136$ $27.99$
$^{234}\text{Th}$ $63.30$ $0.0584 \pm 0.0013$ $0.0538$ $8.50$ $0.0504$ $6.87$
$^{234}\text{Th}$ $92.38$ $0.0926 \pm 0.0081$ $0.0997$ $7.15$ $0.0946$ $5.38$
$^{226}\text{Ra}$ $186.21$ $0.0851 \pm 0.0007$ $0.0863$ $1.31$ $0.0996$ $13.36$
$^{214}\text{Pb}$ $242.00$ $0.0703 \pm 0.0005$ $0.0671$ $4.86$ $0.0856$ $21.67$
$^{214}\text{Pb}$ $295.22$ $0.0540 \pm 0.0004$ $0.0543$ $0.63$ $0.0745$ $27.10$
$^{214}\text{Pb}$ $351.93$ $0.0460 \pm 0.0003$ $0.0451$ $2.11$ $0.0655$ $31.16$
$^{214}\text{Bi}$ $609.31$ $0.0259 \pm 0.0002$ $0.0271$ $4.39$ $0.0441$ $38.59$
$^{214}\text{Bi}$ $1120.30$ $0.0182 \pm 0.0001$ $0.0183$ $0.56$ $0.0294$ $37.55$
$^{214}\text{Bi}$ $1238.11$ $0.0179 \pm 0.0001$ $0.0174$ $2.65$ $0.0275$ $36.72$
$^{214}\text{Bi}$ $1764.49$ $0.0142 \pm 0.0001$ $0.0142$ $0.17$ $0.0214$ $33.61$
$^{214}\text{Bi}$ $2204.21$ $0.0120 \pm 0.0001$ $0.0121$ $0.61$ $0.0180$ $33.08$

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

Cải tiến kỹ thuật so với các công trình công bố trước đây

Nghiên cứu mang lại những đóng góp khoa học rõ nét khi so sánh trực tiếp với các công bố của S. Baccouche (2012) và Farhad Masoudi (2018):

Độ sai biệt FEPE tại các đỉnh năng lượng đơn phát (63.3 keV, 92.4 keV):

Mô tả lớp vật liệu ma trận:
  1. Khảo sát toàn diện dải năng lượng cực rộng: Bao phủ từ vùng tia gamma mềm $46.54\text{ keV}$ ($^{210}\text{Pb}$) đến tia gamma cứng $2614.51\text{ keV}$ ($^{208}\text{Tl}$), chứng minh hiện tượng đảo chiều hiệu suất cực đại tại khoảng $90 - 150\text{ keV}$.
  2. Định lượng nguyên nhân sai lệch thực nghiệm - mô phỏng:
    • Tại vùng $E < 100\text{ keV}$: Độ sai biệt $RD$ giữa mô phỏng và thực nghiệm đạt mức lý tưởng ($\le 6.87%$ tại $63.3\text{ keV}$ và $92.4\text{ keV}$). Sự suy giảm hiệu suất thực tế ở $46.54\text{ keV}$ ($RD = 27.99%$) phản ánh chính xác hiện tượng tăng bề dày lớp chết mặt ngoài của tinh thể Ge sau thời gian dài sử dụng (vượt qua giá trị danh định $0.7\text{ mm}$ của ORTEC).
    • Tại vùng $E > 300\text{ keV}$: Giá trị mô phỏng MCNP6 luôn cao hơn thực nghiệm từ $20% - 44.41%$. Nguyên nhân cốt lõi được xác định là do hiệu ứng trùng phùng tổng (Coincidence Summing) trong chuỗi phân rã phức tạp của $^{214}\text{Bi}$, $^{214}\text{Pb}$, $^{228}\text{Ac}$, $^{208}\text{Tl}$ khi đặt mẫu ở hình học sát mặt đầu dò ($0\text{ cm}$).

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

Tình huống triển khai trong quan trắc phóng xạ thực tế

  • Phân tích mẫu đất nông nghiệp và môi trường: Giúp các phòng thí nghiệm đo bức xạ môi trường tính toán chính xác hàm lượng Bq/kg của $^{226}\text{Ra}$, $^{232}\text{Th}$, $^{40}\text{K}$ mà không cần phải mua sắm định kỳ các bộ mẫu chuẩn hóa học đắt tiền từ IAEA.
  • Giám sát an toàn phóng xạ quanh nhà máy điện hạt nhân và mỏ khai khoáng: Ứng dụng mô hình MCNP6 đã kiểm chứng để hiệu chỉnh hệ số tự hấp thụ cho các loại quặng kim loại nặng, đất hiếm có mật độ cao ($\rho > 2.0\text{ g/cm}^3$).

Yêu cầu triển khai hạ tầng kỹ thuật (Deployment Requirements)

  • Hệ thống phần cứng máy chủ mô phỏng: Workstation CPU Intel Xeon hoặc AMD EPYC $\ge 16$ cores, $32\text{ GB RAM}$, hỗ trợ chạy song song OpenMP/MPI để rút ngắn thời gian chạy $10^9$ hạt MCNP6 xuống dưới $45\text{ phút}$.
  • Hạ tầng đo lường: Phòng đo phông thấp được che chắn bằng buồng chì dày tối thiểu $10\text{ cm}$, lót mặt trong bằng cadmium ($1\text{ mm}$) và đồng ($1.5\text{ mm}$) để cắt tia X huỳnh quang của chì.

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

Hạn chế kỹ thuật hiện tại

  • Chưa tích hợp mã nguồn hiệu chỉnh True Coincidence Summing: Mã nguồn MCNP6 tiêu chuẩn chỉ tính toán xác suất vận chuyển đơn hạt photon độc lập, chưa tự động tái tạo xác suất chồng chập xung cascade theo thời gian thực của các đồng vị con phát nhiều tia gamma cùng lúc.
  • Giả định thông số lớp chết tĩnh: Chưa quét thực nghiệm (experimental scanning) bằng nguồn chuẩn điểm chùm hẹp để lập bản đồ biến thiên độ dày lớp chết thực tế theo toạ độ 3D của tinh thể Ge.

Kế hoạch nâng cấp và hướng nghiên cứu mở rộng

  1. Áp dụng phần mềm chuyên dụng EFFTRAN / GESPECOR: Kết hợp hệ số truyền hiệu suất (Efficiency Transfer) để tính toán ma trận hiệu chỉnh trùng phùng cho các mẫu hình học lớn.
  2. Khảo sát thực nghiệm hệ số tự hấp thụ $\mu/\rho$: Tiến hành đo xuyên qua mẫu với nguồn phát chùm chuẩn ($^{137}\text{Cs}$, $^{60}\text{Co}$) để đối chiếu trực tiếp với cơ sở dữ liệu tiết diện phản ứng photon trong thư viện MCPLIB của MCNP6.
  3. Tối ưu hóa thuật toán tinh chỉnh lớp chết: Viết script tự động hóa (Python + MCNP6 API) để tối ưu hóa đồng thời bề dày lớp chết mặt ngoài và bán kính tinh thể nhạy thông qua thuật toán di truyền (Genetic Algorithm).

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

                                [HỆ SINH THÁI HƯỞNG LỢI]
  • Sinh viên, học viên ngành Vật lý hạt nhân / Kỹ thuật hạt nhân: Tiếp cận quy trình phân tích phổ gamma thực tế từ bước thu nhận tín hiệu trên phần mềm chuyên dụng (Genie 2K, Colegram) đến mô phỏng vật lý thống kê.
  • Kỹ sư vận hành hệ phổ kế tại các trung tâm quan trắc: Nắm bắt phương pháp luận đánh giá sai số hệ thống do lớp chết đầu dò và ma trận mẫu, nâng cao độ tin cậy của chứng chỉ đo lường ISO/IEC 17025.
  • Các viện nghiên cứu và phòng thí nghiệm môi trường: Cắt giảm hơn $70%$ kinh phí duy trì mẫu chuẩn phóng xạ hàng năm nhờ làm chủ kỹ thuật mô phỏng truyền hiệu suất Monte Carlo.

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

1. Cần cấu hình máy tính như thế nào để chạy mô phỏng MCNP6 với 1 tỷ hạt ($10^9$ NPS)?

Để chạy bài toán photon với $NPS = 10^9$ hạt nhằm đảm bảo sai số thống kê của Tally F8 $< 1%$, hệ thống cần tối thiểu CPU 8 nhân / 16 luồng (xung nhịp $\ge 3.5\text{ GHz}$), RAM $16\text{ GB}$. Khi chạy đa luồng trên hệ điều hành Linux với lệnh mcnp6 tasks 16 inp=HPGe.inp, thời gian xử lý mất khoảng $1.5 - 2.5\text{ giờ}$.

2. Tại sao hiệu suất đỉnh tại năng lượng $46.54\text{ keV}$ của $^{210}\text{Pb}$ lại có độ sai biệt lớn ($27.99%$)?

Vùng năng lượng thấp ($< 60\text{ keV}$) cực kỳ nhạy cảm với sự suy giảm bức xạ qua lớp chết mặt ngoài tinh thể Ge ($0.7\text{ mm}$) và vỏ nhôm ($1.0\text{ mm}$). Trong thực tế, đầu dò sau nhiều năm hoạt động bị gia tăng lớp chết do sự khuếch tán tạp chất ở nhiệt độ phòng, làm suy giảm số đếm thực tế so với mô hình thiết kế ban đầu của nhà sản xuất.

3. Làm thế nào để xử lý sự chồng chập đỉnh tại vạch năng lượng $186.2\text{ keV}$?

Vạch $186.2\text{ keV}$ là sự chồng chập giữa $^{226}\text{Ra}$ ($186.21\text{ keV}$) và $^{235}\text{U}$ ($185.71\text{ keV}$). Cần sử dụng hệ phương pháp tách đóng góp bằng cách nhân diện tích đỉnh thực nghiệm với tỷ lệ phát nhánh chuẩn hóa $57.2%$ cho $^{226}\text{Ra}$ khi mẫu đạt trạng thái cân bằng thế kỷ (Secular Equilibrium sau 30 ngày lưu trữ kín).

4. Tại sao giá trị hiệu suất mô phỏng MCNP6 lại cao hơn thực nghiệm ở các vùng năng lượng $>300\text{ keV}$?

Ở hình học đo áp sát bề mặt đầu dò ($0\text{ cm}$), các đồng vị phát nhiều tia gamma đồng thời theo chuỗi (như $^{214}\text{Bi}$ phát vạch $609.3\text{ keV}$ cùng lúc với các tia khác) sẽ gây ra hiện tượng trùng phùng tổng (TCS). Hai photon đi vào đầu dò cùng lúc bị ghi nhận thành một xung có biên độ bằng tổng năng lượng, làm hụt diện tích đỉnh đơn lẻ trong thực nghiệm, trong khi mô phỏng độc lập từng photon không bị mất mát này nếu không cấu hình mô hình phân rã phức hợp.

5. Hàm khớp đa thức bậc 5 theo $\ln(E)$ có áp dụng được cho mọi dải năng lượng không?

Hàm đa thức $\ln(\varepsilon) = \sum_{i=0}^5 B_i (\ln E)^i$ đạt độ khớp xuất sắc ($R^2 > 0.99$, độ sai lệch $RD < 8.5%$) cho dải năng lượng từ $46.5\text{ keV}$ đến $2614.5\text{ keV}$ đối với đầu dò HPGe đồng trục. Đối với dải năng lượng cực thấp ($< 30\text{ keV}$), cần chuyển sang hàm khớp hữu tỷ hoặc chia dải năng lượng thành 2 nhánh trước và sau đỉnh suy giáp (K-edge absorption jump).

Kết luận

Khóa luận đã giải quyết trọn vẹn bài toán xác định và đối chuẩn hiệu suất ghi đỉnh năng lượng toàn phần (FEPE) của hệ phổ kế gamma HPGe GEM50P4-83 cho 4 hệ mẫu chuẩn môi trường đại diện IAEA. Thông qua việc kết hợp giữa phân tích thực nghiệm phông thấp và mô phỏng 3D Monte Carlo MCNP6 với $10^9$ lịch sử hạt, nghiên cứu đã:

  • Xây dựng thành công đường chuẩn hiệu suất giải tích chính xác cao trên toàn bộ dải năng lượng quan trắc $46.5 - 2614.5\text{ keV}$ với sai số làm khớp $\le 8.50%$.
  • Định lượng hóa ảnh hưởng của thành phần nguyên tố ma trận mẫu và làm sáng tỏ nguyên nhân sai khác giữa mô hình lý thuyết và thực nghiệm (chỉ ra vai trò chi phối của lớp chết tinh thể ở vùng năng lượng thấp và hiệu ứng trùng phùng tổng ở vùng năng lượng cao).
  • Cung cấp cơ sở dữ liệu vật liệu và bộ mã đầu vào MCNP6 chuẩn mực, tạo tiền đề ứng dụng trực tiếp cho các phòng thí nghiệm kiểm định an toàn bức xạ môi trường trên toàn quốc.