Bối cảnh và vấn đề nghiên cứu
Trong các nghiên cứu vật lý hạt nhân và ứng dụng phân tích phóng xạ, detector bán dẫn germanium có độ tinh khiết cao (HPGe - High Pure Germanium) đóng vai trò then chốt nhờ khả năng ghi nhận bức xạ tia gamma và tia X với độ phân giải năng lượng vượt trội (lên đến 0,1%). Tuy nhiên, việc khảo sát thực nghiệm trực tiếp hàm đáp ứng của detector và tối ưu hóa hệ phổ kế thường gặp khó khăn do hạn chế về trang thiết bị, điều kiện đo đạc và phông môi trường. Sự phát triển của các phương pháp mô phỏng máy tính, đặc biệt là phương pháp Monte Carlo thông qua chương trình MCNP, cung cấp công cụ tính toán vận chuyển bức xạ, hỗ trợ kiểm chứng và phân tích số liệu thực nghiệm.
Khóa luận "Nghiên cứu ảnh hưởng của các thông số mô phỏng đến hàm đáp ứng của detector HPGe" của tác giả Trương Thị Phương Thi (Khoa Vật lý, Trường Đại học Sư phạm TP.HCM, năm 2012, người hướng dẫn: ThS. Trần Thiện Thanh) tập trung vào việc mô hình hóa hệ phổ kế gamma sử dụng chương trình MCNP4C2 nhằm giải quyết các nhiệm vụ nghiên cứu cụ thể:
- Xây dựng tệp đầu vào (input file) mô tả cấu trúc hình học, thành phần vật liệu của hệ nguồn - đầu dò LEGe (Low Energy Germanium) và buồng chì che chắn trong chương trình MCNP4C2.
- Khảo sát mức độ ảnh hưởng của các thông số mô phỏng vật lý photon (
PHYS:P, chế độ vận chuyển electron MODE P E so với MODE P, chế độ tán xạ kết hợp NOCOH, cơ chế bức xạ hãm/tạo electron IDES) đến hiệu suất đỉnh của phổ gamma.
- So sánh, đối chiếu kết quả mô phỏng với số liệu chuẩn của Quados (2003) để kiểm tra độ tin cậy.
- Đánh giá sự phụ thuộc của hiệu suất đỉnh vào mật độ vật liệu tinh thể detector.
- Đánh giá mức độ suy giảm phông phóng xạ môi trường (độ lệch phổ) của ba cấu hình buồng chì chuẩn (cấu hình 737, 747 và 777).
- Đối tượng nghiên cứu: Đầu dò bán dẫn germanium năng lượng thấp (LEGe), nguồn phát bức xạ photon dạng đĩa, các lớp vật liệu che chắn và ba cấu hình buồng chì (737, 747, 777).
- Phạm vi nghiên cứu: Mô phỏng số trên phần mềm MCNP4C2; dải năng lượng khảo sát từ 15 keV đến 1 MeV; khoảng cách từ nguồn dạng đĩa đến cửa sổ Beryllium là 2 cm trong môi trường không khí.
Cơ sở lý thuyết và phương pháp
Khung lý thuyết và mô hình tương tác bức xạ
Khóa luận vận dụng lý thuyết vật lý hạt nhân về tương tác của photon với vật chất và nguyên lý hoạt động của detector bán dẫn chuyển tiếp P-I-N:
- Nguyên lý detector HPGe: Vùng I (vùng nghèo) nhạy với bức xạ ion hóa. Khi phân cực ngược, photon tương tác tạo cặp electron - lỗ trống; điện trường quét hạt mang điện về hai điện cực P và N tạo xung điện tích. Tinh thể Ge có vùng cấm hẹp nên được làm lạnh bằng nitơ lỏng ($LN_2$) ở nhiệt độ 77 K trong buồng chân không nhằm giảm dòng rò nhiệt.
- Đại lượng đặc trưng:
- Độ phân giải năng lượng: $R = \frac{\text{FWHM}}{H_0}$ (tính theo %).
- Hiệu suất toàn phần ($\varepsilon_t$): xác suất photon để lại năng lượng khác không trong vùng hoạt động.
- Hiệu suất đỉnh ($\varepsilon_p$): xác suất photon truyền toàn bộ năng lượng vào vùng hoạt động.
- Tỉ số đỉnh trên tổng: $P/T = \frac{\varepsilon_p}{\varepsilon_t}$.
- Cơ chế tương tác photon trong MCNP:
- Tán xạ Compton: Mô tả theo tiết diện vi phân Klein-Nishina hiệu chỉnh bởi thừa số dạng $I(Z,v)$.
- Tán xạ Thomson (tán xạ kết hợp): Không làm thay đổi năng lượng photon, sử dụng thừa số dạng $C^*(Z,v)$.
- Hấp thụ quang điện: Làm mất toàn bộ năng lượng photon, giải phóng electron quỹ đạo và phát photon huỳnh quang (tia X đặc trưng) hoặc điện tử Auger.
- Hiệu ứng tạo cặp: Xảy ra trong trường hạt nhân khi năng lượng photon $h\nu > 1{,}022\text{ MeV}$, tạo ra cặp $e^- - e^+$ và hai photon hủy hạt $0{,}511\text{ MeV}$.
Phương pháp nghiên cứu
Nghiên cứu sử dụng phương pháp thử thống kê Monte Carlo mô phỏng quá trình vận chuyển hạt từng bước:
- Phần mềm: Chương trình MCNP4C2 (Monte Carlo N-Particle Transport Code System).
- Kỹ thuật tính toán:
- Thiết lập hình học hệ đo gồm 9 ô mạng (cells) và 19 mặt biên (surfaces).
- Tally F8 (Pulse Height Tally): phân bố năng lượng theo độ cao xung tạo ra trong detector trên 8192 kênh phổ.
- Tùy chọn mở rộng năng lượng Gauss (GEB - Gaussian Energy Broadening) để làm khớp phổ mô phỏng với phổ thực nghiệm theo hàm thực nghiệm:
$$\text{FWHM} = a + b\sqrt{E} + cE$$
với các hệ số thực nghiệm: $a = 3{,}592 \times 10^{-4}$, $b = 1{,}832 \times 10^{-2}$, $c = 0{,}01918$.
- Nguồn dữ liệu: Thư viện tiết diện tương tác hạt nhân của MCNP4C2 và tập dữ liệu đối sánh quốc tế từ tài liệu Quados (2003).
Nội dung chính theo từng chương
Chương 1: Detector bán dẫn siêu tinh khiết HPGe
Chương này tổng quan về cấu tạo và tính chất vật lý của detector bán dẫn germanium:
- Phân tích ưu thế của detector HPGe so với detector Ge(Li) (vốn không ổn định ở nhiệt độ phòng do hiện tượng trôi dạt lithium ra khỏi vùng nội).
- Phân loại 7 nhóm detector germanium: Ultra-LEGe (năng lượng cực thấp), LEGe (năng lượng thấp), Coaxial Ge (đồng trục), REGe (điện cực ngược), Well (hình giếng), XtRa (dải rộng) và BEGe (năng lượng mở rộng). Đề tài xác định sử dụng loại LEGe cho dải năng lượng thấp và trung bình.
- Định nghĩa các tham số kỹ thuật: độ phân giải năng lượng (FWHM), hàm đáp ứng phân bố xung, hiệu suất ghi toàn phần, hiệu suất đỉnh và tỉ số $P/T$.
Chương 2: MCNP- Tập hợp các công cụ mô phỏng vận chuyển bức xạ bằng phương pháp Monte - Carlo
Nội dung chương trình bày cơ sở thuật toán và các khối lệnh điều khiển của MCNP:
- Bản chất thuật toán Monte Carlo: giải bài toán vận chuyển hạt ngẫu nhiên, sai số thống kê tỉ lệ nghịch với $\sqrt{N}$ (với $N$ là số sự kiện mô phỏng).
- Giới thiệu chương trình MCNP do Thomas N. Godfrey và cộng sự phát triển, xử lý dải năng lượng từ $10^{-11}\text{ MeV} - 20\text{ MeV}$ đối với neutron và $10^{-3}\text{ MeV} - 10^3\text{ MeV}$ đối với photon/electron.
- Hai chế độ xử lý vật lý photon:
- Xử lý đơn giản: Bỏ qua tán xạ kết hợp và huỳnh quang.
- Xử lý chi tiết: Tính đến tán xạ Thomson, huỳnh quang tia X sau hấp thụ quang điện.
- Cấu trúc tệp đầu vào MCNP: khai báo ô mạng (
cell cards), khai báo mặt (surface cards), thẻ nguồn (SDEF), thẻ loại hạt (MODE P, MODE P E), thẻ ghi nhận kết quả (F8:P, F8:P,E), thẻ vật liệu (M), thẻ điều kiện dừng (NPS, CTME) và xử lý phổ thực nghiệm (GEB).
Chương 3: Mô tả hệ phổ kế gamma HPGe dùng trong mô phỏng và kết quả mô phỏng
Chương này mô tả chi tiết các tham số hình học mô phỏng và phân tích kết quả tính toán:
- Thông số hình học hệ đo:
- Tinh thể Ge: bán kính vùng hoạt động $19{,}9\text{ mm}$, chiều cao $14{,}8\text{ mm}$, lớp chết bề mặt dày $100\ \mu\text{m}$.
- Vỏ kẹp nhôm (Al): phần kẹp phía sau bán kính $30\text{ mm}$, chiều cao $5\text{ mm}$; phần kẹp bên bán kính ngoài $30\text{ mm}$, bán kính trong $20\text{ mm}$, chiều cao $15\text{ mm}$.
- Cửa sổ Beryllium (Be): bán kính $26\text{ mm}$, chiều dày $1{,}5\text{ mm}$, đặt cách mặt đầu dò $5\text{ mm}$.
- Vỏ ngoài nhôm: dày $1{,}5\text{ mm}$, bán kính ngoài $41{,}5\text{ mm}$, chiều cao $58\text{ mm}$.
- Nguồn gamma: dạng đĩa bán kính $2{,}5\text{ mm}$, phát đẳng hướng trong không khí, cách cửa sổ Be $2\text{ cm}$.
- Thông số cấu hình buồng chì:
- Cấu hình 737: Lớp vỏ thép carbon dày $0{,}95\text{ cm}$, lớp chì dày $10\text{ cm}$, lớp lót trong gồm $1{,}5\text{ mm}\text{ Cu}$ và $1{,}0\text{ mm}\text{ Sn}$; có lỗ tròn ngang để kết nối bình nitơ lỏng.
- Cấu hình 747: Lớp vỏ thép carbon dày $0{,}95\text{ cm}$, lớp chì dày $10\text{ cm}$, lớp lót trong gồm $1{,}6\text{ mm}\text{ Cu}$ và $1{,}0\text{ mm}\text{ Sn}$; diện tích sàn $0{,}4\text{ m}^2$.
- Cấu hình 777: Cấu hình nền siêu thấp, lớp vỏ thép carbon dày $0{,}95\text{ cm}$, lớp chì dày $15\text{ cm}$ (tăng 50% bề dày so với 737/747), lớp lót trong gồm $1{,}5\text{ mm}\text{ Cu}$ siêu tinh khiết và $1{,}0\text{ mm}\text{ Sn}$.
- Phân chia ô mạng mô phỏng (9 ô): Ô 1 (nguồn), Ô 2 (tinh thể Ge), Ô 3 (lớp chết), Ô 4 (nhôm bên trong), Ô 5 (buồng chân không), Ô 6 (nhôm bên ngoài), Ô 7 (cửa sổ Be), Ô 8 (không gian xung quanh nguồn và detector, độ quan trọng = 1), Ô 9 (không gian ngoài, độ quan trọng = 0).
- Các thông số vật lý khảo sát trên thẻ
PHYS:P: EMCPE (ngưỡng năng lượng xử lý chi tiết, chọn $100\text{ MeV}$), IDES ($0$: sinh electron/TTB; $1$: không sinh electron), NOCOH ($0$: có tán xạ kết hợp; $1$: không có tán xạ kết hợp). Bỏ qua hiệu ứng Doppler (NODOP=1) và va chạm hạt nhân (ISPN=0) trong phiên bản MCNP4C2.
Kết quả và đóng góp
1. Ảnh hưởng của các thông số mô phỏng đến hiệu suất đỉnh
Kết quả khảo sát trên tinh thể detector có mật độ chuẩn $5{,}323\text{ g/cm}^3$ khi không đặt trong buồng chì:
| Trường hợp khảo sát |
Năng lượng |
Độ sai biệt hiệu suất đỉnh trong khóa luận (%) |
Số liệu tham chiếu Quados (2003) (%) |
Mức độ sai lệch so với Quados (%) |
Không có tán xạ kết hợp (NOCOH=1) so với mặc định |
15 keV |
4,69 |
3,40 |
37,94 |
|
30 keV |
0,61 |
0,70 |
12,86 |
Chế độ MODE P E so với MODE P |
15 keV |
+7,30 |
+5,50 |
32,73 |
|
30 keV |
+1,52 |
+1,40 |
8,57 |
|
750 keV |
-2,52 |
-1,50 |
68,00 |
|
1 MeV |
-3,02 |
-2,20 |
37,27 |
Nhận xét quy luật:
- Ở vùng năng lượng thấp (15 keV, 30 keV), thông số
NOCOH ảnh hưởng rõ nét vì tán xạ kết hợp làm lệch hướng photon năng lượng thấp mà không làm tiêu hao năng lượng. Chế độ MODE P E làm tăng hiệu suất đỉnh so với MODE P.
- Ở vùng năng lượng cao (750 keV, 1 MeV), thông số
IDES và việc vận chuyển electron làm giảm hiệu suất đỉnh do hiện tượng electron động năng cao thoát ra khỏi thể tích hoạt động của detector trước khi truyền hết năng lượng.
- Ảnh hưởng tổng thể của
NOCOH và IDES lên hiệu suất đỉnh nhỏ hơn 5%.
2. Sự phụ thuộc của hiệu suất đỉnh vào mật độ detector
- Khi mật độ detector tăng từ $5{,}0\text{ g/cm}^3$ đến các giá trị cao hơn, hiệu suất đỉnh ở dải năng lượng trung bình và cao tăng theo tỉ lệ thuận do xác suất tương tác tăng.
- Riêng tại các đỉnh năng lượng thấp (15 keV và 30 keV), khi mật độ tinh thể tăng thì hiệu suất đỉnh lại giảm. Nguyên nhân là mật độ của lớp chết bề mặt ($100\ \mu\text{m}$) tăng theo, dẫn đến phần lớn photon năng lượng thấp bị hấp thụ hoặc suy giảm năng lượng ngay tại lớp chết trước khi đi vào vùng hoạt động.
3. Đánh giá mức độ che chắn của các cấu hình buồng chì
Độ suy giảm phông môi trường được đánh giá qua độ lệch phổ (phần trăm sai khác giữa trường hợp có buồng chì và không có buồng chì). Số liệu so sánh tại năng lượng 1 MeV với mật độ detector $5{,}0\text{ g/cm}^3$:
| Cấu hình buồng chì |
Bề dày lớp chì (cm) |
Chi tiết lớp lót hấp thụ tia X của chì |
Độ lệch phổ tại 1 MeV (%) |
| Cấu hình 747 |
10,0 |
1,6 mm Cu + 1,0 mm Sn |
0,7799 |
| Cấu hình 737 |
10,0 |
1,5 mm Cu + 1,0 mm Sn |
0,8230 |
| Cấu hình 777 |
15,0 |
1,5 mm Cu (độ tinh khiết cao) + 1,0 mm Sn |
1,0814 |
Kết luận thực nghiệm mô phỏng:
- Cấu hình 777 cho khả năng giảm phông môi trường cao nhất nhờ bề dày lớp chì lớn hơn 50% ($15\text{ cm}$).
- Cấu hình 737 cho hiệu quả giảm phông tốt hơn cấu hình 747 dù cùng độ dày lớp chì $10\text{ cm}$, do khác biệt về kích thước hình học không gian bên trong.
- Lớp lót Cu và Sn trong cả 3 cấu hình triệt tiêu hiệu quả bức xạ tia X đặc trưng của chì phát sinh trong khoảng $75 - 85\text{ keV}$.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận
- Sai lệch số liệu: Kết quả mô phỏng độ sai biệt hiệu suất đỉnh của tác giả còn có độ chênh lệch đáng kể khi đối chiếu với số liệu chuẩn của Quados (2003), cụ thể sai lệch $37{,}94%$ ở $15\text{ keV}$ đối với trường hợp
NOCOH=1 và sai lệch lên đến $68{,}00%$ ở $750\text{ keV}$ đối với trường hợp MODE P E.
- Giới hạn phần mềm: Phiên bản MCNP4C2 sử dụng trong nghiên cứu chưa hỗ trợ xử lý mở rộng Doppler (hiệu ứng dãn nở Doppler do động lượng của electron liên kết trước va chạm).
Hướng nghiên cứu tiếp theo
- Khảo sát ảnh hưởng của các thông số mô phỏng đến hàm đáp ứng detector với các dạng hình học nguồn khác nhau (nguồn thể tích, nguồn điểm) và các dòng detector khác (Coaxial, Well, BEGe).
- Ứng dụng các phiên bản MCNP nâng cao (như MCNP5) để đánh giá đầy đủ ảnh hưởng của hiệu ứng dãn nở Doppler.
- Nghiên cứu sự phụ thuộc của hiệu suất ghi vào các thông số mở rộng: mật độ nguồn, khoảng cách nguồn - detector và kích thước hình học chi tiết của đầu dò.
Giá trị tham khảo
Tài liệu khóa luận cung cấp tài liệu tham khảo chuyên ngành cho sinh viên, học viên cao học và nghiên cứu viên thuộc các ngành:
- Vật lý hạt nhân, Kỹ thuật hạt nhân và Vật lý y khoa: Cung cấp thông số hình học chi tiết của đầu dò LEGe, các công thức thực nghiệm xác định FWHM và hàm làm khớp mở rộng đỉnh Gauss (
GEB).
- Mô phỏng bức xạ và an toàn bức xạ: Cung cấp phương pháp thiết lập tệp đầu vào (9 ô mạng, 19 mặt biên), cách sử dụng thẻ ghi nhận
F8 kết hợp xử lý thẻ vật lý PHYS:P (IDES, NOCOH, EMCPE), cùng số liệu kích thước chuẩn của 3 buồng chì che chắn thương mại (737, 747, 777).
Câu hỏi thường gặp
1. Khóa luận sử dụng loại đầu dò HPGe nào và kích thước hình học của tinh thể Ge là bao nhiêu?
Tác giả sử dụng đầu dò bán dẫn germanium năng lượng thấp LEGe. Tinh thể Ge siêu tinh khiết có bán kính vùng hoạt động $19{,}9\text{ mm}$, chiều cao $14{,}8\text{ mm}$, được bao bọc bởi lớp chết bề mặt dày $100\ \mu\text{m}$. Mặt trước đầu dò có cửa sổ mỏng Beryllium bán kính $26\text{ mm}$, dày $1{,}5\text{ mm}$, cách đầu dò $5\text{ mm}$.
2. Tại sao ở vùng năng lượng thấp (15 keV và 30 keV), hiệu suất đỉnh lại giảm khi mật độ tinh thể detector tăng?
Ở dải năng lượng thấp, khi mật độ khối của detector tăng thì mật độ của lớp chết bề mặt ($100\ \mu\text{m}$) cũng tăng tương ứng. Các photon năng lượng thấp (15 keV, 30 keV) có quãng chạy ngắn, dễ bị hấp thụ hoàn toàn hoặc tiêu hao một phần năng lượng ngay trong lớp chết trước khi đi vào thể tích hoạt động của đầu dò, dẫn đến số đếm tại đỉnh năng lượng toàn phần bị suy giảm.
3. Ý nghĩa của các thông số NOCOH và IDES trên thẻ PHYS:P trong MCNP là gì?
NOCOH: Điều khiển tán xạ kết hợp (Thomson). Giá trị NOCOH=0 cho phép xảy ra tán xạ kết hợp; NOCOH=1 tắt tán xạ kết hợp. Thông số này ảnh hưởng chủ yếu ở vùng năng lượng thấp (15 keV, 30 keV).IDES: Điều khiển cơ chế tạo electron và bức xạ hãm. Giá trị IDES=0 cho phép photon tạo ra electron trong MODE E hoặc xử lý bức xạ hãm bia dày (TTB); IDES=1 ngăn không cho tạo ra electron. Thông số này ảnh hưởng rõ nét ở vùng năng lượng cao (750 keV, 1 MeV).
4. Điểm khác biệt cơ bản giữa ba cấu hình buồng chì 737, 747 và 777 là gì?
- Cấu hình 737 và 747 có lớp chì dày $10\text{ cm}$ bọc ngoài bằng thép carbon $0{,}95\text{ cm}$. Cấu hình 737 có lỗ tròn ngang kết nối bình nitơ lỏng; cấu hình 747 thiết kế nhỏ gọn với diện tích sàn $0{,}4\text{ m}^2$.
- Cấu hình 777 là buồng chì phông siêu thấp với lớp chì dày $15\text{ cm}$ (dày hơn 50% so với 737 và 747) kết hợp lớp đồng tinh khiết cao $1{,}5\text{ mm}$ và thiếc $1{,}0\text{ mm}$, mang lại độ lệch phổ cao nhất ($1{,}0814%$ tại 1 MeV so với $0{,}8230%$ của 737 và $0{,}7799%$ của 747).
5. Vì sao nghiên cứu chưa khảo sát được ảnh hưởng của hiệu ứng dãn nở Doppler?
Chương trình mô phỏng được sử dụng trong khóa luận là MCNP4C2. Phiên bản này chưa tích hợp thư viện và mô hình tính toán hiệu ứng dãn nở Doppler (Doppler broadening) do động lượng của electron liên kết trước va chạm gây ra, buộc tác giả phải thiết lập bỏ qua thông số này (NODOP=1).
Kết luận
Khóa luận đã xây dựng mô hình mô phỏng hệ phổ kế gamma đầu dò HPGe (LEGe) kết hợp nguồn đĩa và các cấu hình buồng chì bằng chương trình MCNP4C2. Nghiên cứu đã định lượng mức độ ảnh hưởng của các thông số vật lý photon (NOCOH, IDES, MODE P E), chứng minh ảnh hưởng của chúng đến hiệu suất đỉnh dưới 5% và làm rõ quy luật biến thiên của hiệu suất theo mật độ detector cùng lớp chết bề mặt. Đồng thời, kết quả tính toán đã xác nhận cấu hình buồng chì 777 với bề dày chì 15 cm có hiệu quả che chắn giảm phông phóng xạ môi trường tốt nhất trong ba cấu hình được khảo sát.
