Tổng quan nghiên cứu

Trong lĩnh vực đo lường bức xạ hạt nhân và chẩn đoán hình ảnh y khoa hiện đại, các thiết bị cảm biến bán dẫn đang khẳng định vị thế ưu việt nhờ kích thước nhỏ gọn và độ phân giải năng lượng cao. Nhu cầu ghi nhận chính xác bức xạ photon đòi hỏi vật liệu đầu dò phải có khả năng tương tác mạnh, trong đó chất bán dẫn hợp chất Cadmium Telluride (CdTe) nổi lên như một giải pháp đột phá với số nguyên tử hiệu dụng trung bình xấp xỉ 50 và độ rộng vùng cấm từ 1,44 eV đến 1,52 eV tại nhiệt độ phòng 300 K. Tuy nhiên, việc khảo sát thực nghiệm các đặc tính ghi nhận bức xạ thường tiêu tốn kinh phí lớn cho thiết bị chuẩn, đồng thời tiềm ẩn nguy cơ phơi nhiễm phóng xạ đối với người thao tác.

Vấn đề cốt lõi của nghiên cứu là xác định chính xác đường cong hiệu suất ghi nhận bức xạ gamma và đánh giá toàn diện các yếu tố ảnh hưởng như năng lượng photon tới, kích thước hình học tinh thể và khoảng cách đo. Mục tiêu cụ thể của luận văn là xây dựng hệ mô phỏng hoàn chỉnh cấu trúc đầu dò bán dẫn XR-100CdTe của hãng Amptek bằng bộ công cụ Geant4, ứng dụng phương pháp Monte Carlo để khảo sát đường cong hiệu suất trong dải năng lượng khảo sát từ 7,8 keV đến 2 MeV.

Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh trong thời gian từ tháng 8 đến tháng 12 năm 2018. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp một quy trình tính toán số chuẩn xác, giúp giảm thiểu 100% rủi ro bức xạ trong quá trình thiết kế thử nghiệm và rút ngắn khoảng 70% thời gian hiệu chuẩn phổ kế thực nghiệm, đồng thời mở ra hướng ứng dụng hiệu quả trong y học hạt nhân và quan trắc môi trường.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên hệ thống cơ sở lý thuyết vững chắc về vật lý tương tác bức xạ và vật lý chất bán dẫn:

  • Lý thuyết tương tác giữa photon gamma và vật chất: Khi photon gamma đi qua môi trường vật chất, 4 quá trình tương tác chính xảy ra tùy thuộc vào năng lượng photon và số nguyên tử Z của môi trường:
    1. Hấp thụ quang điện: Photon truyền toàn bộ năng lượng cho electron liên kết, giải phóng electron quang điện với động năng xác định theo công thức năng lượng photon trừ năng lượng liên kết. Tiết diện tương tác quang điện tỷ lệ thuận với hàm bậc 5 của số nguyên tử ($Z^5$) và giảm mạnh theo năng lượng photon, chiếm vai trò chủ đạo ở dải năng lượng dưới 100 keV.
    2. Tán xạ Compton: Photon va chạm không xuyên tâm với electron, truyền một phần động năng cho electron giật lùi và tán xạ ra photon thứ cấp có bước sóng dài hơn. Quá trình này chiếm ưu thế trong khoảng năng lượng từ 100 keV đến 1 MeV.
    3. Quá trình tạo cặp: Xảy ra khi photon có năng lượng lớn hơn ngưỡng 1,022 MeV (tương ứng hai lần khối lượng nghỉ của electron) tương tác với trường Coulomb hạt nhân để tạo ra cặp electron - positron.
    4. Tán xạ Rayleigh: Quá trình tán xạ đàn hồi không làm mất năng lượng của photon, chủ yếu xuất hiện ở vùng năng lượng rất thấp.
  • Lý thuyết đầu dò bán dẫn CdTe: Tinh thể CdTe có cấu trúc mạng zincblende đối xứng lập phương, khối lượng phân tử 240,01 g/mol và điểm nóng chảy 1092 °C. Bức xạ ion hóa đi vào thể tích hoạt động tạo ra các cặp electron - lỗ trống với năng lượng ion hóa trung bình chỉ 4,43 eV trên mỗi cặp hạt dẫn (thấp hơn nhiều so với mức khoảng 30 eV trong buồng ion hóa khí và 100 eV trong đầu dò nhấp nháy NaI(Tl)). Dưới tác dụng của điện áp phân cực cao +500 V, các hạt mang điện được quét về hai điện cực để tạo thành xung tín hiệu có biên độ tỷ lệ với năng lượng hấp thụ.
  • Hệ thống khái niệm hiệu suất: Bao gồm hiệu suất tuyệt đối (tỷ số giữa số xung ghi nhận và tổng số photon phát ra từ nguồn), hiệu suất nội (tỷ số xung ghi nhận trên số photon thực tế tới bề mặt đầu dò), hiệu suất toàn phần và hiệu suất đỉnh năng lượng toàn phần (xác suất photon giải phóng 100% năng lượng trong vùng hoạt).

Phương pháp nghiên cứu

  • Nguồn dữ liệu và cấu hình hệ thống: Mô hình hóa chính xác đầu dò thương mại XR-100CdTe của Amptek với kích thước tinh thể CdTe hoạt động là $5 \times 5 \times 1\text{ mm}$ (diện tích bề mặt $25\text{ mm}^2$), độ phân giải năng lượng danh định dưới 1,5 keV FWHM tại đỉnh 122 keV của Co-57, được bao bọc bởi lớp cửa sổ Beryllium dày 100 µm và màng Graphite dày 25 µm.
  • Cỡ mẫu và phương pháp chọn mẫu: Sử dụng phương pháp lấy mẫu ngẫu nhiên Monte Carlo để mô phỏng tương tác của chùm hạt bức xạ. Cỡ mẫu mô phỏng được thiết lập ở mức từ $10^6$ đến $10^7$ lịch sử hạt phát ra từ nguồn hình học dạng đĩa nhằm đảm bảo tính hội tụ theo Định lý giới hạn trung tâm và Luật số lớn, giúp sai số thống kê duy trì ở mức dưới 1%.
  • Phương pháp phân tích và lý do lựa chọn: Nghiên cứu sử dụng bộ công cụ Geant4 (Geometry ANd Tracking) phiên bản chuẩn của CERN trên nền tảng ngôn ngữ C++ kết hợp khung phân tích dữ liệu ROOT. Lý do lựa chọn Geant4 là khả năng mô tả hình học ba chiều phức tạp thông qua các lớp chuyên biệt (G4UserDetectorConstruction, G4UserPhysicsList, G4UserPrimaryGeneratorAction), hỗ trợ đầy đủ các gói vật lý tương tác điện từ tiêu chuẩn và năng lượng thấp, đồng thời cho phép thiết lập ngưỡng cắt vi phân 1 mm giúp tối ưu hóa thời gian tính toán mà vẫn đảm bảo độ chuẩn xác tuyệt đối của phổ năng lượng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  1. Đặc tính đáp ứng của đường cong hiệu suất theo năng lượng: Trong dải năng lượng thấp từ 10 keV đến 40 keV, hiệu suất đỉnh của đầu dò CdTe đạt giá trị cực đại xấp xỉ 100% nhờ ưu thế hấp thụ quang điện tuyệt đối của vật liệu có số nguyên tử cao ($Z \approx 50$). Tuy nhiên, khi năng lượng photon tăng từ 50 keV lên 100 keV, hiệu suất ghi nhận giảm xuống khoảng 65%, và tiếp tục suy giảm mạnh xuống dưới mức 15% ở ngưỡng năng lượng 500 keV do tiết diện quang điện giảm nhanh và tán xạ Compton trở nên chiếm ưu thế.
  2. Ảnh hưởng rõ rệt của bề dày tinh thể bán dẫn: Việc khảo sát đường cong hiệu suất với các độ dày khác nhau cho thấy khi gia tăng bề dày tinh thể CdTe từ 1 mm lên 2,25 mm, hiệu suất ghi nhận tại dải năng lượng trung bình (từ 100 keV đến 300 keV) được cải thiện đáng kể, tăng thêm từ 35% đến 55%. Bề dày lớn hơn giúp gia tăng quãng đường tương tác của photon, từ đó nâng cao xác suất hấp thụ toàn phần.
  3. Sự vượt trội tuyệt đối của CdTe so với đầu dò Silic (Si): Kết quả so sánh trực tiếp ở cùng điều kiện hình học cho thấy tại mức năng lượng 60 keV, đầu dò CdTe dày 1 mm đạt hiệu suất ghi nhận cao gấp hơn 4 lần so với đầu dò Silic dày 500 µm (vốn chỉ đạt hiệu suất dưới 15%). Điều này khẳng định ưu thế của vật liệu có mật độ cao và số nguyên tử Z lớn trong việc giữ lại năng lượng photon.
  4. Quy luật suy giảm hiệu suất theo khoảng cách và lớp hấp thụ bề mặt: Hiệu suất tuyệt đối của hệ đo suy giảm tỷ lệ nghịch với bình phương khoảng cách nguồn - đầu dò đúng theo quy luật biến thiên của góc khối hình học. Đồng thời, lớp tiếp xúc cửa sổ Beryllium (100 µm) hoặc Nhôm chỉ làm suy giảm nhẹ khoảng 2% đến 5% số đếm ở vùng năng lượng siêu thấp dưới 10 keV và gần như không gây ảnh hưởng đến dải năng lượng cao hơn.

Thảo luận kết quả

Các kết quả mô phỏng phản ánh chính xác bản chất vật lý của quá trình bẫy hạt mang điện (hole trapping) trong chất bán dẫn hợp chất CdTe. Do độ linh động và thời gian sống của lỗ trống trong CdTe thấp hơn đáng kể so với electron, các tương tác diễn ra gần cực dương sẽ cho hiệu suất thu gom điện tích cao hơn các tương tác gần cực âm.

Dữ liệu mô phỏng được biểu diễn trực quan thông qua các bảng số liệu chi tiết và đồ thị hàm log-log thể hiện xác suất tương tác trong dải năng lượng từ 1 keV đến 1 MeV, cùng biểu đồ tuyến tính phân tích vùng năng lượng 10 keV đến 250 keV. Khi đối chiếu với thông số kỹ thuật xuất xưởng của nhà sản xuất Amptek cho thiết bị XR-100CdTe, đường cong hiệu suất mô phỏng bằng Geant4 có độ trùng khớp cao với độ lệch tương đối chỉ dao động trong khoảng 2% đến 4%. Điều này chứng minh độ tin cậy tuyệt đối của hệ mô phỏng được xây dựng.

Đề xuất và khuyến nghị

  1. Chuẩn hóa quy trình mô phỏng Geant4 trong thiết kế phổ kế hạt nhân: Các viện nghiên cứu và trường đại học chuyên ngành cần ban hành hướng dẫn thực hành chuẩn về mô phỏng Geant4 trong vòng 6 tháng tới. Việc tối ưu hóa các thông số ngưỡng cắt vi phân ở mức 1 mm và khai báo chính xác vật liệu NIST sẽ giúp giảm sai số tính toán xuống dưới 2%, tạo tiền đề vững chắc cho việc thiết kế các hệ đo mới.
  2. Tối ưu hóa kích thước hình học tinh thể CdTe cho ứng dụng y tế: Các nhóm nghiên cứu phát triển thiết bị chẩn đoán hình ảnh y sinh cần nâng độ dày tinh thể CdTe lên mức từ 2 mm đến 3 mm trong lộ trình 12 tháng, nhằm tăng hiệu suất phát hiện bức xạ gamma 140 keV của đồng vị phóng xạ Technetium-99m thêm tối thiểu 40% mà không làm suy giảm đáng kể độ phân giải năng lượng.
  3. Tích hợp vi mạch xử lý bù bẫy điện tích: Các kỹ sư điện tử hạt nhân cần nghiên cứu phát triển mạch tiền khuếch đại thông minh ứng dụng bộ lọc số FPGA trong giai đoạn 18 tháng tới, hướng đến mục tiêu hiệu chỉnh hiện tượng mất mát điện tích do bẫy lỗ trống, giữ độ phân giải năng lượng luôn đạt mức FWHM dưới 1,5 keV tại đỉnh 122 keV.
  4. Duy trì chế độ làm mát nhiệt điện ổn định cho hệ đầu dò: Các cơ sở vận hành thiết bị phân tích huỳnh quang tia X (XRF) và phổ kế gamma cần trang bị bộ làm mát Peltier hai giai đoạn đảm bảo nhiệt độ đầu dò luôn duy trì dưới -10 °C với công suất tiêu thụ dưới 1 W, nhằm triệt tiêu dòng rò nhiệt xuống dưới mức $5 \times 10^{-3}\text{ counts/s}$ và bảo toàn độ nhạy tín hiệu 0,82 mV/keV.

Đối tượng nên tham khảo luận văn

  • Học viên cao học và nghiên cứu sinh ngành Vật lý Kỹ thuật, Vật lý Hạt nhân: Cung cấp tài liệu tham khảo giá trị về phương pháp lập trình C++ trên nền tảng Geant4, kỹ thuật khai báo các lớp đối tượng hình học không gian ba chiều và thuật toán xử lý dữ liệu phổ bằng phần mềm ROOT.
  • Kỹ sư công nghệ y sinh và chuyên gia chẩn đoán hình ảnh: Nắm vững các đặc tính vật lý của đầu dò bán dẫn CdTe để ứng dụng vào việc tối ưu hóa thiết kế đầu thu cho hệ thống chụp cắt lớp vi tính bức xạ đơn photon (SPECT) và máy chụp X-quang kỹ thuật số liều thấp.
  • Cán bộ quản lý an toàn bức xạ và quan trắc môi trường phóng xạ: Ứng dụng quy trình mô phỏng để xây dựng đường cong hiệu chuẩn phổ kế gamma hiện trường nhanh chóng mà không cần phụ thuộc vào nguồn chuẩn phóng xạ hoạt độ cao.
  • Các nhà sản xuất và phát triển cảm biến đo lường bức xạ: Khai thác các dữ liệu so sánh định lượng về độ dày tinh thể, vật liệu cửa sổ bảo vệ và điện áp phân cực để cải tiến cấu trúc cơ - điện tử cho các thiết bị đo phóng xạ cầm tay thế hệ mới.

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

  • Tại sao đầu dò bán dẫn CdTe lại vượt trội hơn đầu dò Silic (Si) trong đo bức xạ gamma? Đầu dò CdTe vượt trội nhờ số nguyên tử hiệu dụng trung bình $Z \approx 50$ (cao hơn nhiều so với $Z=14$ của Silic), giúp tiết diện hấp thụ quang điện (vốn tỷ lệ thuận với $Z^5$) tăng gấp hàng trăm lần. Do đó, tại cùng độ dày 1 mm và mức năng lượng 60 keV, hiệu suất ghi nhận của CdTe đạt trên 80% trong khi Silic chỉ đạt dưới 15%.
  • Phương pháp Monte Carlo trong chương trình Geant4 hoạt động theo nguyên lý nào? Phương pháp Monte Carlo sử dụng các chuỗi số giả ngẫu nhiên để mô phỏng lại từng bước tương tác xác suất của từng hạt photon riêng lẻ theo hàm mật độ phân bố xác suất vật lý. Khi số lượng lịch sử hạt mô phỏng đạt từ $10^6$ đến $10^7$ sự kiện, giá trị trung bình thống kê sẽ tiệm cận chính xác với quy luật thực nghiệm thực tế.
  • Vì sao đầu dò XR-100CdTe cần được cấp điện áp +500 V và làm lạnh nhiệt điện? Điện trường mạnh từ mức điện áp +500 V giúp thu gom nhanh chóng các hạt mang điện về điện cực trước khi chúng bị tái hợp hoặc bẫy lại. Do điện áp cao có thể gây dòng rò lớn ở nhiệt độ phòng, hệ thống bắt buộc phải được làm lạnh xuống dưới -10 °C qua bộ làm mát Peltier để giảm dòng rò xuống mức tối thiểu dưới $1\ \mu\text{A}$.
  • Sự khác biệt căn bản giữa hiệu suất đỉnh và hiệu suất toàn phần là gì? Hiệu suất đỉnh là xác suất photon giải phóng 100% năng lượng ban đầu trong vùng hoạt của đầu dò (tạo nên đỉnh năng lượng đơn vạch), trong khi hiệu suất toàn phần tính đến mọi tương tác để lại bất kỳ mức năng lượng nào khác 0. Trong phân tích định lượng đồng vị phóng xạ, hiệu suất đỉnh là thông số quan trọng nhất.
  • Độ dày của tinh thể CdTe ảnh hưởng như thế nào đến khả năng ghi nhận bức xạ? Độ dày tinh thể quyết định quãng đường tương tác của photon. Đối với photon năng lượng cao trên 100 keV, việc tăng độ dày từ 1 mm lên 2,25 mm giúp xác suất va chạm và hấp thụ quang điện tăng lên rõ rệt, từ đó nâng cao hiệu suất ghi nhận của đầu dò thêm từ 35% đến 55%.

Kết luận

  • Xây dựng thành công hệ mô phỏng hình học và vật lý hoàn chỉnh cho đầu dò bán dẫn XR-100CdTe kích thước $5 \times 5 \times 1\text{ mm}$ trên nền tảng Geant4.
  • Xác lập chính xác đường cong hiệu suất bức xạ gamma trong dải năng lượng mở rộng từ 7,8 keV đến 2 MeV, chứng minh hiệu suất đạt cực đại xấp xỉ 100% ở vùng năng lượng dưới 40 keV.
  • Đánh giá định lượng tính ưu việt của vật liệu CdTe với hiệu suất ghi nhận vượt trội gấp hơn 4 lần so với đầu dò Silic ở mức năng lượng 60 keV.
  • Khảo sát chi tiết mối quan hệ giữa hiệu suất với khoảng cách đo và độ dày tinh thể, cho thấy độ dày 2,25 mm giúp tăng hiệu suất dải năng lượng trung bình thêm tới 55%.
  • Kiểm chứng kết quả mô phỏng với thông số kỹ thuật thực nghiệm của nhà sản xuất Amptek với sai số tương đối cực nhỏ chỉ từ 2% đến 4%.

Đóng góp lớn nhất của công trình là cung cấp một giải pháp công nghệ số chuẩn xác, an toàn và tiết kiệm kinh phí trong việc nghiên cứu, đánh giá đặc tính của các hệ phổ kế hạt nhân bán dẫn tiên tiến.

Trong giai đoạn tiếp theo từ năm 2024 đến 2026, hướng nghiên cứu sẽ tiếp tục mở rộng mô phỏng cho các hệ đầu dò mảng đa kênh CdZnTe dạng pixelated và tích hợp các thuật toán xử lý tín hiệu xung số thời gian thực. Hãy kết nối và tham khảo toàn diện cấu hình mô phỏng C++ trong luận văn để áp dụng hiệu quả vào các dự án nghiên cứu công nghệ bức xạ và vật lý kỹ thuật của bạn.