Bối cảnh và vấn đề nghiên cứu

Trong các ứng dụng kỹ thuật hạt nhân thuộc lĩnh vực y tế, công nghiệp, năng lượng và quan trắc phóng xạ môi trường, hệ phổ kế gamma sử dụng đầu dò nhấp nháy NaI(Tl) kết hợp với bộ phân tích đa kênh (MCA) là một trong những công cụ phát hiện bức xạ được sử dụng phổ biến nhờ hiệu suất ghi cao ở điều kiện nhiệt độ phòng. Tuy nhiên, trong các phép đo phổ phóng xạ định lượng, hiệu suất ghi thực tế của đầu dò không bao giờ đạt 100% do sự tương tác không hoàn toàn của photon trong thể tích nhạy và sự suy giảm bức xạ qua các lớp vỏ bảo vệ bên ngoài.

Để xác định đường cong hiệu suất đỉnh năng lượng toàn phần (FEPE - Full Energy Peak Efficiency) mà không cần tiến hành thực nghiệm phức tạp với nhiều nguồn phóng xạ nguy hiểm, các nhà nghiên cứu thường ứng dụng mô phỏng Monte Carlo thông qua chương trình MCNP5. Mặc dù vậy, kết quả tính toán hiệu suất từ mô phỏng thường tồn tại độ lệch đáng kể so với thực nghiệm nếu chỉ sử dụng các thông số kỹ thuật hình học và vật liệu danh định do nhà sản xuất cung cấp. Các nghiên cứu trước đây (như nghiên cứu của Hoàng Đức Tâm và cộng sự) chủ yếu khảo sát việc hiệu chỉnh bề dày của lớp phản xạ nhôm oxit ($Al_2O_3$) bao quanh tinh thể. Việc thiếu một quy trình chuẩn hóa nhằm tách rời và tối ưu hóa độc lập từng thông số cấu trúc của đầu dò đặt ra yêu cầu phải có phương pháp khảo sát mới.

Mục tiêu và nhiệm vụ nghiên cứu:

  1. Nghiên cứu 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 và nguyên lý vận hành của đầu dò nhấp nháy NaI(Tl).
  2. Xây dựng phương pháp mô phỏng vận chuyển hạt bằng chương trình MCNP5 nhằm mô hình hóa chính xác hệ phổ kế gamma.
  3. Đề xuất quy trình thực nghiệm và tính toán mới nhằm tách rời sự ảnh hưởng của 3 thông số kỹ thuật cốt lõi: mật độ lớp phản xạ $Al_2O_3$, bán kính tinh thể NaI(Tl) và chiều dài tinh thể NaI(Tl).
  4. Khớp hàm toán học, nội suy tìm các thông số kỹ thuật tối ưu từ dữ liệu thực nghiệm và đánh giá lại độ lệch hiệu suất FEPE giữa mô hình mô phỏng tối ưu hóa với số liệu thực nghiệm.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Đầu dò nhấp nháy tinh thể NaI(Tl) kích thước danh định $76{,}2 \times 76{,}2\text{ mm}$ ($3 \times 3\text{ inch}$) và các nguồn phóng xạ chuẩn điểm gồm $^{133}\text{Ba}$, $^{137}\text{Cs}$, $^{241}\text{Am}$, $^{152}\text{Eu}$, $^{22}\text{Na}$, $^{60}\text{Co}$.
  • Phạm vi không gian và thời gian: Nghiên cứu được thực hiện tại Phòng thí nghiệm Vật lý Hạt nhân, Khoa Vật lý, Trường Đại học Sư phạm TP. Hồ Chí Minh trong năm 2019; khảo sát dải năng lượng photon gamma trải dài từ $31\text{ keV}$ đến $1408\text{ keV}$.

Cơ sở lý thuyết và phương pháp

Khóa luận kết hợp chặt chẽ giữa cơ sở lý thuyết tương tác bức xạ, mô hình hóa toán học giải tích, phương pháp mô phỏng ngẫu nhiên Monte Carlo và phương pháp đo phổ gamma thực nghiệm.

1. Khung lý thuyết và mô hình tính toán

  • Cơ chế tương tác photon - vật chất: Khi chùm bức xạ gamma truyền qua vật chất có bề dày $d$, cường độ chùm tia suy giảm theo định luật hàm mũ: $$I = I_0 e^{-\mu d} = I_0 e^{-\mu_m \rho d}$$ Trong đó $\mu$ là hệ số suy giảm tuyến tính, $\mu_m$ là hệ số suy giảm khối lượng và $\rho$ là mật độ vật liệu bia. Quá trình mất năng lượng của photon trong tinh thể được chi phối bởi ba hiệu ứng cơ bản:
    • Hiệu ứng quang điện: Photon truyền toàn bộ năng lượng cho electron liên kết theo phương trình Einstein: $\frac{1}{2}mv_{\max}^2 = hf - \varepsilon_{lk}$.
    • Hiệu ứng tán xạ Compton: Va chạm giữa photon và electron tự do làm tán xạ photon và giải phóng electron với động năng $T_e = E - E'$.
    • Hiệu ứng tạo cặp: Xảy ra khi năng lượng photon $E \ge 2m_e c^2 \approx 1022\text{ keV}$, photon biến mất và tạo thành cặp electron - positron.
  • Cơ chế nhấp nháy trong tinh thể NaI(Tl): Tinh thể NaI tinh khiết được pha tạp Thallium ($Tl^+$) tạo ra các mức năng lượng hoạt hóa trung gian giữa vùng hóa trị và vùng dẫn. Cặp electron - lỗ trống sinh ra do ion hóa sẽ kích thích các ion $Tl^+ \to (Tl^+)^*$, sau đó tái hợp về trạng thái cơ bản và phát ra photon quang học nhìn thấy được để ống nhân quang điện (PMT) thu nhận.
  • Hiệu suất đỉnh năng lượng toàn phần (FEPE): Được tính toán bằng công thức: $$\varepsilon_{FEPE} = \frac{N}{A_0 e^{-\lambda t} I_\gamma t} = \varepsilon_{int} \cdot \varepsilon_g$$ Trong đó $N$ là diện tích đỉnh thuần (số đếm trừ phông), $A_0$ là hoạt độ ban đầu, $\lambda$ là hằng số phân rã, $I_\gamma$ là xác suất phát bức xạ, $t$ là thời gian đo sống, $\varepsilon_{int}$ là hiệu suất nội và $\varepsilon_g = \frac{\Omega}{4\pi}$ là hiệu suất hình học với góc khối $\Omega = 2\pi\left(1 - \frac{d}{\sqrt{d^2 + r^2}}\right)$.

2. Phương pháp nghiên cứu và công cụ

  • Mô phỏng Monte Carlo bằng chương trình MCNP5: Sử dụng mã MCNP5 (viết bằng Fortran 90) để mô phỏng sự truyền và lắng đọng năng lượng của photon. Tập tin đầu vào (file input) được xây dựng gồm các thẻ ô mạng (Cell Card), thẻ mặt (Surface Card), thẻ định nghĩa nguồn tổng quát (Data Card - SDEF) và thẻ tính toán xung phổ năng lượng Tally F8 ghi nhận photon trong ô mạng tinh thể (F8:p 1).
  • Phương pháp tách rời biến số và nội suy hàm khớp:
    • Mật độ lớp phản xạ $Al_2O_3$ ($\rho_{Al_2O_3}$): Sử dụng chùm photon năng lượng thấp ($31\text{ keV}$ từ $^{133}\text{Ba}$) có độ đâm xuyên rất ngắn và kết hợp ống chuẩn trực bằng đồng. Vì photon chỉ tương tác ở lớp ngoài cùng và không lan tới đáy hay biên tinh thể, ảnh hưởng của bán kính và chiều dài tinh thể bị triệt tiêu. Quan hệ giữa hiệu suất và mật độ được xấp xỉ tuyến tính qua khai triển MacLaurin: $\varepsilon_{FEPE} = \alpha - \beta \rho_{Al_2O_3}$.
    • Bán kính tinh thể ($R$): Sử dụng nguồn điểm năng lượng thấp ($31\text{ keV} - 121\text{ keV}$), không chuẩn trực, đặt cách đầu dò $40\text{ cm}$ ($d/R = 10{,}5$). Bức xạ chỉ phân bố tại mặt trước nên không chịu ảnh hưởng của chiều dài tinh thể; diện tích tiếp nhận phụ thuộc trực tiếp vào bán kính $R$ theo hàm tuyến tính $\varepsilon_{FEPE} = a + b R$. Bán kính tối ưu tổng hợp được tính theo phương pháp trung bình có trọng số sai số.
    • Chiều dài tinh thể ($L$): Sử dụng nguồn photon năng lượng cao ($662\text{ keV} - 1408\text{ keV}$), photon đâm xuyên sâu đến đáy tinh thể, từ đó khảo sát sự phụ thuộc tuyến tính $\varepsilon_{FEPE} = a + b L$.

Nội dung chính theo từng chương

Chương 1: Tương tác bức xạ gamma với vật chất

Chương 1 hệ thống hóa cơ sở lý thuyết vật lý hạt nhân về bức xạ điện từ bước sóng ngắn ($\lambda < 10^{-11}\text{ m}$). Tác giả phân tích chi tiết quy luật suy giảm cường độ chùm tia photon truyền qua môi trường vật chất theo hàm mũ của bề dày và mật độ. Ba cơ chế tương tác cơ bản (quang điện, tán xạ Compton và tạo cặp) được trình bày kèm theo các phương trình bảo toàn năng lượng, động lượng tương đối tính và các điều kiện biên của từng hiệu ứng. Đồng thời, chương này phân tích hình thái phổ bức xạ gamma thực nghiệm (điển hình ở năng lượng $1408\text{ keV}$), giải thích nguồn gốc hình thành đỉnh hấp thụ toàn phần, phổ liên tục Compton, mép Compton (Compton edge), đỉnh tán xạ ngược và các đỉnh thoát đơn/thoát đôi phát sinh từ quá trình hủy hạt positron ($511\text{ keV}$).

Chương 2: Phương pháp Monte Carlo và chương trình MCNP5

Chương 2 trình bày tổng quan lịch sử và nguyên lý toán học của phương pháp mô phỏng số ngẫu nhiên Monte Carlo, khởi nguồn từ bài toán chiếc kim của Georges-Louis Leclerc Comte de Buffon (1777), thí nghiệm lặp của R. Zilin’ski (5000 lần tung cho $\pi \approx 3{,}1236$) đến ứng dụng của John von Neumann và Stanislaw Ulam trong dự án Manhattan (1944). Chương này tập trung phân tích kiến trúc của chương trình MCNP5 do Phòng thí nghiệm Quốc gia Los Alamos phát triển từ các phiên bản tiền thân MCS (1963), MCN và MCNG (1973). Cấu trúc của một tập tin đầu vào MCNP5 được mô tả chi tiết:

  • Thẻ tiêu đề (Title Card).
  • Thẻ ô mạng (Cell Card): Định nghĩa số vật liệu, mật độ nguyên tử hoặc mật độ khối lượng ($g/\text{cm}^3$), toán tử giới hạn không gian hình học.
  • Thẻ mặt (Surface Card): Định nghĩa các mặt phẳng ($P, PX, PY, PZ$), mặt cầu ($S, SO, SX, SY, SZ$), mặt trụ ($C/X, C/Y, C/Z, CX, CY, CZ$) và mặt nón ($K/X, K/Y, K/Z, KX, KY, KZ$).
  • Thẻ dữ liệu (Data Card): Định nghĩa nguồn phát bức xạ tổng quát (SDEF), phân bố năng lượng rời rạc/liên tục (SI, SP, SB), vật liệu (Material Card - $M_m$ với mã đồng vị $ZZZ.AAA$) và bộ đếm ghi nhận năng lượng Tally F8.

Chương 3: Đầu dò NaI(Tl) và phương pháp xác định các thông số tối ưu của đầu dò NaI(Tl)

Chương 3 trình bày cấu tạo chi tiết, cơ chế tạo xung nhấp nháy của đầu dò NaI(Tl) và các công thức tính toán hiệu suất đỉnh năng lượng toàn phần $\varepsilon_{FEPE}$. Thông số hình học và vật liệu danh định của đầu dò do nhà sản xuất công bố được tổng hợp như sau:

Thành phần cấu trúc Bề dày / Bán kính (cm) Chiều dài (cm) Mật độ danh định ($g/\text{cm}^3$)
Lớp vỏ nhôm ngoài Bề dày: $0{,}10$ $8{,}14$ $2{,}699$
Lớp không khí đệm Bề dày: $0{,}10$ $7{,}79$ $0{,}001205$
Lớp phản xạ nhôm oxit ($Al_2O_3$) Bề dày: $0{,}16$ $7{,}78$ $0{,}55$
Tinh thể NaI(Tl) Bán kính: $3{,}81$ $7{,}62$ $3{,}67$

Tác giả thiết lập 3 mô hình thực nghiệm và mô phỏng tương ứng để tách rời từng thông số:

  • Thí nghiệm 1: Nguồn $^{133}\text{Ba}$ (hoạt độ $3{,}7 \cdot 10^9\text{ Bq}$ tại thời điểm $15/05/2013$, hằng số phân rã $\lambda = 2{,}09 \cdot 10^{-9}\text{ s}^{-1}$, thời gian đo $25500\text{ s}$ và $89600\text{ s}$) đặt trong ống chuẩn trực bằng đồng phía trước và bên cạnh đầu dò để xác định mật độ lớp phản xạ $Al_2O_3$ thông qua đỉnh $31\text{ keV}$. Dữ liệu từ NIST cho biết hệ số suy giảm khối của $Al_2O_3$ ở $30{,}973\text{ keV}$ là $\mu_m = 0{,}6408\text{ cm}^2/g$.
  • Thí nghiệm 2: Các nguồn photon năng lượng thấp ($^{133}\text{Ba}$, $^{137}\text{Cs}$, $^{241}\text{Am}$, $^{152}\text{Eu}$) phát các đỉnh $31\text{ keV}, 32\text{ keV}, 59\text{ keV}, 81\text{ keV}, 121\text{ keV}$ đặt cách đầu dò $40\text{ cm}$ (độ chính xác dịch chuyển $\pm 0{,}01\text{ mm}$), mô phỏng $10^7$ lịch sử hạt để khảo sát bán kính $R$ từ $3{,}72\text{ cm}$ đến $3{,}87\text{ cm}$.
  • Thí nghiệm 3: Các nguồn photon năng lượng cao ($^{137}\text{Cs}$, $^{22}\text{Na}$, $^{60}\text{Co}$, $^{152}\text{Eu}$) phát các đỉnh $662\text{ keV}, 964\text{ keV}, 1173\text{ keV}, 1274\text{ keV}, 1332\text{ keV}, 1408\text{ keV}$ đặt cách đầu dò $40\text{ cm}$ để khảo sát chiều dài tinh thể $L$ từ $7{,}48\text{ cm}$ đến $7{,}78\text{ cm}$.

Chương 4: Kết quả và thảo luận

Chương 4 công bố toàn bộ dữ liệu mô phỏng, phương trình khớp hàm, kết quả nội suy và so sánh đối chiếu giữa mô phỏng và thực nghiệm.

1. Kết quả xác định mật độ lớp phản xạ $Al_2O_3$

Mô phỏng thay đổi mật độ $\rho_{Al_2O_3}$ từ $0{,}4$ đến $3{,}6\text{ g/cm}^3$ cho đỉnh $31\text{ keV}$ thu được các hàm hồi quy tuyến tính:

  • Nguồn đặt phía trước đầu dò: $\varepsilon_{FEPE} \cdot 10^4 = (6{,}64 \pm 0{,}01) - (0{,}56 \pm 0{,}01)\rho$ ($R^2 = 0{,}9978$).
  • Nguồn đặt bên cạnh đầu dò: $\varepsilon_{FEPE} \cdot 10^4 = (6{,}56 \pm 0{,}02) - (0{,}56 \pm 0{,}01)\rho$ ($R^2 = 0{,}9977$).

Nội suy với hiệu suất thực nghiệm qua nhiều lần đo cho giá trị mật độ trung bình tối ưu là $\rho_{Al_2O_3} = 2{,}02 \pm 0{,}32\text{ g/cm}^3$, khác biệt rõ rệt so với giá trị danh định $0{,}55\text{ g/cm}^3$ của nhà sản xuất.

2. Kết quả xác định bán kính tối ưu của tinh thể NaI(Tl)

Khảo sát bán kính $R$ trong khoảng $3{,}72 - 3{,}87\text{ cm}$ cho 5 đỉnh năng lượng thấp cho thấy mối tương quan tuyến tính rất cao ($R^2 \ge 0{,}9990$):

  • Đỉnh $31\text{ keV}$: $\varepsilon \cdot 10^4 = -17{,}2 + 9{,}00 R$ ($R^2 = 0{,}9990$).
  • Đỉnh $32\text{ keV}$: $\varepsilon \cdot 10^4 = -18{,}30 + 9{,}40 R$ ($R^2 = 0{,}9990$).
  • Đỉnh $59\text{ keV}$: $\varepsilon \cdot 10^4 = -18{,}88 + 10{,}10 R$ ($R^2 = 0{,}9999$).
  • Đỉnh $81\text{ keV}$: $\varepsilon \cdot 10^4 = -17{,}7 + 9{,}41 R$ ($R^2 = 0{,}9999$).
  • Đỉnh $121\text{ keV}$: $\varepsilon \cdot 10^4 = -16{,}80 + 8{,}77 R$ ($R^2 = 0{,}9999$).

Bằng phương pháp tính trung bình có trọng số nghịch đảo phương sai sai số, giá trị bán kính tinh thể tối ưu được xác định là $R = 3{,}76 \pm 0{,}03\text{ cm}$ (so với giá trị danh định $3{,}81\text{ cm}$).

3. Kết quả xác định chiều dài tối ưu của tinh thể NaI(Tl)

Khảo sát chiều dài tinh thể từ $7{,}48\text{ cm}$ đến $7{,}78\text{ cm}$ với các đỉnh năng lượng cao cho thấy tương quan tuyến tính chặt chẽ:

  • Đỉnh $662\text{ keV}$: $\varepsilon \cdot 10^4 = 5{,}46 + 0{,}47 L$ ($R^2 = 0{,}9978$).
  • Đỉnh $964\text{ keV}$: $\varepsilon \cdot 10^4 = 3{,}065 + 0{,}478 L$ ($R^2 = 0{,}9999$).
  • Đỉnh $1173\text{ keV}$: $\varepsilon \cdot 10^4 = 2{,}271 + 0{,}455 L$ ($R^2 = 0{,}9999$).
  • Đỉnh $1274\text{ keV}$: $\varepsilon \cdot 10^4 = 2{,}064 + 0{,}432 L$ ($R^2 = 0{,}9995$).
  • Đỉnh $1332\text{ keV}$: $\varepsilon \cdot 10^4 = 1{,}869 + 0{,}434 L$ ($R^2 = 0{,}9999$).
  • Đỉnh $1408\text{ keV}$: $\varepsilon \cdot 10^4 = 1{,}723 + 0{,}422 L$ ($R^2 = 0{,}9999$).

Kết quả và đóng góp

1. Kết quả định lượng chính

Sau khi hiệu chỉnh đồng thời mật độ lớp phản xạ $\rho_{Al_2O_3} = 2{,}02\text{ g/cm}^3$ và bán kính tinh thể $R = 3{,}76\text{ cm}$, độ chênh lệch hiệu suất đỉnh toàn phần $\varepsilon_{FEPE}$ giữa mô phỏng MCNP5 và thực nghiệm giảm xuống dưới $2{,}16%$ trên toàn bộ dải năng lượng thấp.

Đỉnh năng lượng (keV) $\varepsilon_{FEPE}$ Thực nghiệm ($10^{-4}$) $\varepsilon_{FEPE}$ Mô phỏng ban đầu ($10^{-4}$) Độ lệch mô hình ban đầu (%) Bán kính tối ưu (cm) Độ lệch sau tối ưu (%)
31 $16{,}40$ $19{,}61$ $19{,}51%$ $3{,}73 \pm 0{,}07$ $1{,}39%$
32 $16{,}83$ $20{,}01$ $18{,}87%$ $3{,}74 \pm 0{,}07$ $1{,}05%$
59 $19{,}48$ $20{,}25$ $3{,}92%$ $3{,}77 \pm 0{,}06$ $2{,}16%$
81 $17{,}52$ $18{,}22$ $3{,}96%$ $3{,}74 \pm 0{,}06$ $0{,}65%$
121 $16{,}32$ $17{,}15$ $5{,}09%$ $3{,}78 \pm 0{,}06$ $1{,}10%$

2. Đóng góp khoa học của đề tài

  • Xây dựng quy trình phân rã tham số độc lập: Khác với các nghiên cứu trước đây vốn chỉ khảo sát ảnh hưởng tổng hợp hoặc chỉ điều chỉnh bề dày lớp phản xạ (như công bố của Hoàng Đức Tâm và cs chỉ thay đổi bề dày $Al_2O_3$), đề tài đã đưa ra phương pháp luận thực nghiệm kết hợp chuẩn trực và lựa chọn dải năng lượng thích hợp để cô lập sự ảnh hưởng của từng thông số: mật độ vật liệu phản xạ, bán kính mặt nhận bức xạ và chiều sâu tinh thể.
  • Hiệu chỉnh thông số thực tế cho đầu dò NaI(Tl): Chứng minh rằng mật độ thực tế của lớp $Al_2O_3$ trong đầu dò thương mại ($2{,}02\text{ g/cm}^3$) lớn hơn nhiều so với số liệu kỹ thuật danh định ($0{,}55\text{ g/cm}^3$), và bán kính hiệu dụng của tinh thể nhỏ hơn thông số công bố ($3{,}76\text{ cm}$ so với $3{,}81\text{ cm}$).
  • Nâng cao độ chính xác của mô phỏng: Đưa ra bộ thông số tối ưu giúp mô hình mô phỏng MCNP5 tiệm cận độ chính xác của các phép đo phổ gamma thực nghiệm, đặc biệt ở vùng năng lượng dưới $150\text{ keV}$ vốn có sai số mô phỏng ban đầu lên tới gần $20%$.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế trong văn bản: Nghiên cứu tập trung giải quyết chi tiết việc xác định mật độ lớp phản xạ và bán kính tinh thể ở vùng năng lượng thấp; việc tối ưu hóa đồng thời giá trị chiều dài tinh thể $L$ trên toàn bộ dải năng lượng cao ($662 - 1408\text{ keV}$) mới dừng lại ở bước khảo sát sự phụ thuộc tuyến tính và so sánh từng bước mà chưa tích hợp toàn diện thành một hàm mục tiêu đa biến duy nhất.
  • Hướng nghiên cứu tiếp theo: Tiếp tục hoàn thiện thuật toán tối ưu hóa đồng thời cả ba thông số (mật độ, bán kính, chiều dài) trên toàn bộ dải phổ liên tục từ $31\text{ keV}$ đến $1408\text{ keV}$, mở rộng phương pháp cho các loại đầu dò nhấp nháy khác hoặc đầu dò bán dẫn HPGe, và áp dụng mô hình tối ưu vào phân tích hoạt độ phóng xạ của các mẫu môi trường thực tế.

Giá trị tham khảo

  • Đối tượng thụ hưởng: Sinh viên, học viên cao học và nghiên cứu viên thuộc các chuyên ngành Vật lý Hạt nhân, Kỹ thuật Hạt nhân, Vật lý Y khoa và An toàn Bức xạ.
  • Nội dung giá trị:
    • Tài liệu tham khảo hữu ích về kỹ thuật xây dựng tập tin đầu vào (file input) trong chương trình MCNP5 (khai báo hình học Cell, Surface, nguồn SDEF và Tally F8).
    • Phương pháp luận tách biến thực nghiệm và thuật toán khớp hàm - nội suy sai số có thể ứng dụng trực tiếp để hiệu chuẩn các hệ phổ kế gamma trong các phòng thí nghiệm hạt nhân.

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

1. Tại sao nghiên cứu lại chọn đỉnh năng lượng 31 keV của nguồn 133Ba để xác định mật độ lớp phản xạ Al2O3?

Đỉnh năng lượng $31\text{ keV}$ là mức năng lượng thấp, có hệ số hấp thụ khối $\mu_m$ trong lớp $Al_2O_3$ lớn ($\approx 0{,}6408\text{ cm}^2/g$). Các photon này có quãng chạy ngắn, hầu như để lại toàn bộ năng lượng sát bề mặt ngoài tinh thể và không đâm xuyên tới đáy tinh thể (loại bỏ ảnh hưởng của chiều dài). Kết hợp với việc đặt nguồn trong ống chuẩn trực bằng đồng, góc khối tiếp xúc bị giới hạn nên hoàn toàn triệt tiêu ảnh hưởng của bán kính đầu dò.

2. Sự khác biệt giữa mật độ lớp phản xạ do nhà sản xuất cung cấp và mật độ tối ưu tìm được là bao nhiêu?

Nhà sản xuất cung cấp thông số mật độ danh định của lớp $Al_2O_3$ là $0{,}55\text{ g/cm}^3$. Kết quả mô phỏng kết hợp nội suy từ dữ liệu thực nghiệm cho thấy mật độ tối ưu thực tế là $2{,}02 \pm 0{,}32\text{ g/cm}^3$ (tính trung bình từ cấu hình đo mặt trước $1{,}99\text{ g/cm}^3$ và mặt bên $2{,}14\text{ g/cm}^3$).

3. Bán kính tối ưu của tinh thể NaI(Tl) được xác định bằng bao nhiêu và độ lệch hiệu suất giảm như thế nào?

Bán kính tối ưu được tính bằng phương pháp trung bình có trọng số sai số trên 5 đỉnh năng lượng thấp là $R = 3{,}76 \pm 0{,}03\text{ cm}$ (thay vì $3{,}81\text{ cm}$ danh định). Sau khi tối ưu bán kính và mật độ, độ lệch hiệu suất mô phỏng so với thực nghiệm giảm mạnh: đỉnh $31\text{ keV}$ giảm từ $19{,}51%$ xuống $1{,}39%$; đỉnh $32\text{ keV}$ giảm từ $18{,}87%$ xuống $1{,}05%$; đỉnh $59\text{ keV}$ từ $3{,}92%$ xuống $2{,}16%$; đỉnh $81\text{ keV}$ từ $3{,}96%$ xuống $0{,}65%$; đỉnh $121\text{ keV}$ từ $5{,}09%$ xuống $1{,}10%$.

4. Tally nào trong chương trình MCNP5 được sử dụng để mô phỏng đáp ứng của đầu dò NaI(Tl)?

Nghiên cứu sử dụng thẻ Tally F8 (F8:p), đóng vai trò như một detector ghi nhận xung và tính toán phổ phân bố năng lượng bị mất của photon (lắng đọng năng lượng toàn phần) trong ô mạng chứa vật liệu tinh thể NaI(Tl).

Kết luận

Khóa luận đã đề xuất và thực nghiệm thành công phương pháp mới nhằm xác định chính xác các thông số kỹ thuật cấu trúc của đầu dò nhấp nháy NaI(Tl) phục vụ mô phỏng Monte Carlo MCNP5. Bằng cách thiết lập các điều kiện đo chuẩn trực và phân vùng năng lượng thích hợp, tác giả đã tách rời và xác định được mật độ tối ưu của lớp phản xạ $Al_2O_3$ ($2{,}02 \pm 0{,}32\text{ g/cm}^3$) và bán kính tinh thể hiệu dụng ($3{,}76 \pm 0{,}03\text{ cm}$). Mô hình sau khi hiệu chỉnh đã khắc phục được sai số lớn ở vùng năng lượng thấp, đưa độ lệch giữa mô phỏng và thực nghiệm xuống dưới mức $2{,}16%$, cung cấp cơ sở dữ liệu tin cậy cho các bài toán phân tích phổ gamma định lượng.