Tổng quan nghiên cứu

Sự phát triển mạnh mẽ của kỹ thuật hạt nhân trong y tế, công nghiệp và nghiên cứu khoa học đòi hỏi các giải pháp bảo vệ an toàn bức xạ ngày càng hoàn thiện. Nghiên cứu khảo sát các tham số suy giảm của chùm tia gamma theo năng lượng đối với vật liệu hợp chất và hỗn hợp được thực hiện tại Thành phố Hồ Chí Minh nhằm giải quyết bài toán tối ưu hóa vật liệu che chắn thế hệ mới. Trọng tâm của luận văn là xác định chính xác ba tham số cốt lõi gồm hệ số suy giảm khối, nguyên tử số hiệu dụng và mật độ electron hiệu dụng trong dải năng lượng từ 100 keV đến 2000 keV.

Phạm vi nghiên cứu tập trung vào 4 nhóm vật liệu tiềm năng: xi măng Portland, nhựa Kapton, nhựa PMMA và 6 hệ thủy tinh borate-tellurite-silicate chứa tỷ lệ bismuth oxit từ 50% đến 75%. Thông qua việc mô phỏng 14 mức năng lượng bức xạ đơn năng chuẩn từ các nguồn phóng xạ công nghiệp như Coban-60, Xesi-137 và Ytri-88, công trình cung cấp hệ thống dữ liệu toàn diện giúp thay thế các thử nghiệm thực nghiệm tốn kém. Kết quả nghiên cứu đã chứng minh phương pháp mô phỏng Monte Carlo có độ tin cậy cao với độ lệch tương đối so với chuẩn quốc tế luôn duy trì ở mức dưới 5%, đóng góp cơ sở dữ liệu quan trọng cho ngành an toàn bức xạ và vật lý y khoa.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết của luận văn dựa trên cơ chế tương tác của bức xạ điện từ với vật chất qua ba hiệu ứng chủ đạo:

  1. Hiệu ứng quang điện: Chiếm ưu thế ở vùng năng lượng thấp dưới 1 MeV và vật liệu có điện tích hạt nhân lớn, nơi photon truyền toàn bộ năng lượng cho electron liên kết lớp trong, làm bứt electron và phát ra tia X đặc trưng với tiết diện tương tác tỷ lệ thuận với lũy thừa bậc cao của nguyên tử số.
  2. Tán xạ Compton: Chiếm ưu thế trong dải năng lượng trung bình từ vài trăm keV đến vài MeV, mô tả sự va chạm không đàn hồi giữa photon và electron tự do ở lớp ngoài, làm suy giảm năng lượng photon tới và đổi hướng bay.
  3. Hiệu ứng tạo cặp: Xảy ra khi năng lượng photon vượt ngưỡng 1,022 MeV trong điện trường hạt nhân, biến đổi photon thành cặp electron-positron.

Các khái niệm chính được định lượng bao gồm: hệ số suy giảm tuyến tính biểu thị mức độ suy giảm cường độ chùm tia trên một đơn vị chiều dày; hệ số suy giảm khối loại bỏ ảnh hưởng của mật độ vật chất; nguyên tử số hiệu dụng phản ánh cấu trúc nguyên tử tương đương của hợp chất phức tạp; và mật độ electron hiệu dụng xác định số lượng electron tham gia tương tác trên một đơn vị khối lượng.

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

Nghiên cứu sử dụng phương pháp mô phỏng Monte Carlo thông qua bộ chương trình MCNP5 kết hợp thuật toán xử lý phổ cải tiến bằng phần mềm Colegram.

  • Cỡ mẫu và thiết lập mô phỏng: Cỡ mẫu mô phỏng được thiết lập ở mức 1.000.000 hạt photon cho mỗi cấu hình đo nhằm triệt tiêu sai số thống kê và đảm bảo độ cao xung của tally F8 phản ánh chính xác đáp ứng của đầu dò nhấp nháy NaI(Tl) kích thước 76 mm x 76 mm. Kích thước hình học của mẫu vật liệu được chuẩn hóa dưới dạng tấm phẳng kích thước 10 cm x 10 cm với độ dày 0,5 cm.
  • Phương pháp chọn mẫu: Mẫu nghiên cứu được lựa chọn có chủ đích đại diện cho các nhóm vật liệu che chắn hiện đại: xi măng Portland thông dụng trong xây dựng công trình bức xạ, nhựa nhiệt dẻo PMMA và Kapton có tỷ trọng nhẹ ứng dụng trong hàng không vũ trụ, cùng 6 biến thể thủy tinh borate biến đổi nồng độ bismuth nhằm khảo sát sự thay đổi khả năng cản xạ theo thành phần hóa học.
  • Lý do lựa chọn phương pháp phân tích: Việc sử dụng MCNP5 kết hợp kỹ thuật trừ phông Compton bằng hàm đa thức và làm khớp đỉnh hấp thụ toàn phần bằng hàm Gauss trong Colegram giúp bóc tách chính xác diện tích đỉnh phổ thực. Nguồn dữ liệu đối chứng được thu thập từ cơ sở dữ liệu XCOM của Viện Tiêu chuẩn và Công nghệ Quốc gia Hoa Kỳ NIST và các công trình thực nghiệm quốc tế công bố trong giai đoạn 2015–2017.

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

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

  1. Độ chính xác của hệ số suy giảm khối: Giá trị hệ số suy giảm khối tính toán từ mô phỏng MCNP5 có sự tương thích cao với cơ sở dữ liệu NIST. Độ lệch tương đối giữa mô phỏng và chuẩn NIST trên toàn bộ 14 mức năng lượng từ 122 keV đến 1836 keV đều nằm dưới ngưỡng 5,0%.
  2. Hiệu quả cản xạ vượt trội của thủy tinh borate: Khi hàm lượng bismuth oxit tăng từ 50% lên 75%, hệ số suy giảm khối tại mức năng lượng 122 keV tăng trưởng mạnh từ 1,353 cm²/g lên 3,091 cm²/g, cao gấp nhiều lần so với xi măng Portland (0,153 cm²/g) và nhựa PMMA (0,147 cm²/g).
  3. Mức độ tương thích giữa mô phỏng và thực nghiệm: Độ lệch tương đối giữa kết quả mô phỏng và dữ liệu thực nghiệm đo đạc thực tế đối với các mẫu thủy tinh borate luôn dưới 5,5%. Đối với hai loại vật liệu nhựa polymer, sai số tối đa chỉ đạt mức 8,1%, hoàn toàn đáp ứng tiêu chuẩn kiểm định quốc tế cho phép sai số dưới 10%.
  4. Quy luật biến thiên của nguyên tử số hiệu dụng: Tỷ số nguyên tử số hiệu dụng thu được từ mô phỏng so với tính toán lý thuyết từ NIST và dữ liệu thực nghiệm đều tiệm cận giá trị 1,00 trên toàn dải năng lượng khảo sát.

Thảo luận kết quả

Tại vùng năng lượng thấp dưới 300 keV, hệ số suy giảm khối giảm nhanh theo năng lượng do hiệu ứng quang điện chiếm ưu thế tuyệt đối. Trong dải năng lượng từ 300 keV đến 1500 keV, tán xạ Compton trở thành quá trình tương tác chính, khiến tốc độ suy giảm của hệ số suy giảm khối chậm lại và dần ổn định. Điều này lý giải tại sao các vật liệu chứa nguyên tố nặng như thủy tinh borate giàu bismuth thể hiện năng lực che chắn vượt trội ở dải năng lượng thấp, trong khi các polymer nhẹ như PMMA và Kapton phát huy hiệu quả ở dải năng lượng trung bình nhờ mật độ electron cao.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua các bảng thống kê sai số tương đối, đồ thị phân bố tỷ số nguyên tử số hiệu dụng theo năng lượng và các đường cong hồi quy phi tuyến. Mô hình hàm số làm khớp sự phụ thuộc của hệ số suy giảm khối vào năng lượng đạt hệ số xác định tương quan rất cao, tiệm cận giá trị tuyệt đối 0,999. Kết quả này mở ra khả năng nội suy trực tiếp các thông số cản xạ tại bất kỳ mức năng lượng nào mà không cần lặp lại quy trình tính toán phức tạp.

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

  1. Ứng dụng công nghệ mô phỏng MCNP5 trong thiết kế vật liệu che chắn: Các viện nghiên cứu hạt nhân và trung tâm an toàn bức xạ cần tích hợp quy trình mô phỏng số vào giai đoạn tiền khả thi nhằm cắt giảm 60% chi phí thử nghiệm thực tế và rút ngắn thời gian thử nghiệm mẫu trong lộ trình 6 đến 12 tháng.
  2. Chuẩn hóa quy trình xử lý phổ gamma cải tiến: Các phòng thí nghiệm đo lường bức xạ cần áp dụng rộng rãi kỹ thuật làm khớp hàm Gauss kết hợp đa thức trừ phông nền Compton qua phần mềm Colegram để kiểm soát sai số xác định diện tích đỉnh phổ dưới mức 2,0%, thực hiện hoàn tất trước quý 4 năm 2026.
  3. Phát triển sản xuất thủy tinh borate giàu bismuth cho lĩnh vực y tế: Các doanh nghiệp sản xuất vật liệu xây dựng chuyên dụng cần phối hợp với các bệnh viện ung bướu để chế tạo kính quan sát cản xạ chứa từ 60% đến 75% bismuth oxit, hướng tới mục tiêu tăng hiệu suất cản bức xạ gamma lên 35% so với vật liệu kính truyền thống trong giai đoạn 2026–2028.
  4. Xây dựng hệ thống cơ sở dữ liệu số hóa và hàm nội suy thông số cản xạ: Cơ quan quản lý an toàn bức xạ phối hợp với các trường đại học triển khai nền tảng tra cứu trực tuyến các hệ số suy giảm khối cho trên 20 loại vật liệu hợp chất và composite mới trong vòng 18 tháng tới.

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

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý hạt nhân: Nắm vững phương pháp thiết lập hình học card nguồn, định nghĩa vật liệu ô mạng cell card và khai báo tally F8 trong MCNP5 để phục vụ các đề tài tính toán liều và che chắn.
  2. Kỹ sư an toàn bức xạ và chuyên gia vật lý y khoa tại bệnh viện: Ứng dụng số liệu nguyên tử số hiệu dụng và hệ số suy giảm khối để thiết kế tấm chắn bảo vệ phòng chụp X-quang, phòng xạ trị Spect và máy gia tốc y tế.
  3. Kỹ sư công nghệ vật liệu và chuyên gia sản xuất kính cản xạ: Khai thác công thức phối trộn hệ thủy tinh borate-tellurite-silicate nhằm phát triển dòng kính chì nhân tạo không độc hại và thân thiện với môi trường.
  4. Giảng viên và sinh viên các trường đại học khối khoa học tự nhiên: Sử dụng tài liệu làm giáo trình tham khảo chuyên sâu về tương tác bức xạ với vật chất và ứng dụng kỹ thuật số trong xử lý dữ liệu phổ gamma.

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

Phương pháp mô phỏng Monte Carlo MCNP5 có thể thay thế hoàn toàn thực nghiệm đo bức xạ không?

Mô phỏng Monte Carlo MCNP5 có thể thay thế hầu hết các bước thử nghiệm ban đầu nhờ độ chính xác cao với sai số dưới 5% so với chuẩn NIST. Phương pháp này giúp khảo sát an toàn các nguồn phóng xạ nguy hiểm và tiết kiệm chi phí, tuy nhiên thực nghiệm vật lý vẫn cần thiết ở bước thẩm định cuối cùng cho các công trình trọng điểm.

Tại sao nguyên tử số hiệu dụng lại biến thiên theo năng lượng của chùm tia gamma?

Đối với hợp chất hoặc hỗn hợp, mỗi nguyên tố thành phần có tiết diện tương tác khác nhau tùy thuộc vào năng lượng photon tới. Khi chuyển từ vùng năng lượng thấp do hiệu ứng quang điện chi phối sang vùng năng lượng cao do tán xạ Compton chiếm ưu thế, tỷ trọng đóng góp của từng nguyên tố thay đổi khiến nguyên tử số hiệu dụng biến thiên theo năng lượng.

Kỹ thuật xử lý phổ cải tiến bằng phần mềm Colegram có ưu điểm gì vượt trội?

Kỹ thuật này sử dụng hàm Gauss để làm khớp đỉnh hấp thụ năng lượng toàn phần và hàm đa thức để mô tả chính xác phông nền Compton sinh ra bởi bức xạ thứ cấp. Điều này giúp loại bỏ triệt để hiện tượng chân phổ bất đối xứng, đảm bảo diện tích đỉnh tính toán có độ tin cậy cao và sai số thống kê dưới 2%.

Vật liệu nào trong nghiên cứu thể hiện khả năng che chắn bức xạ gamma tốt nhất?

Hệ thủy tinh borate chứa 75% bismuth oxit thể hiện khả năng che chắn vượt trội nhất với hệ số suy giảm khối đạt 3,091 cm²/g tại mức năng lượng 122 keV. Khả năng cản xạ này cao gấp hơn 20 lần so với xi măng Portland thông thường, biến nó thành vật liệu thay thế hoàn hảo cho kính chì truyền thống.

Việc thiết lập 1.000.000 lịch sử hạt photon trong mô phỏng mang lại ý nghĩa gì?

Thiết lập số lượng 1.000.000 hạt nguồn đảm bảo tính thống kê chuẩn xác theo quy luật số lớn trong thuật toán Monte Carlo. Số lượng hạt này giúp giảm thiểu phương sai của tally F8, đảm bảo kết quả phổ ghi nhận tại đầu dò NaI(Tl) đạt độ mịn cao và sai số tương đối so với thực nghiệm luôn duy trì dưới 8,1%.

Kết luận

  • Luận văn đã xác định chuẩn xác bộ ba tham số suy giảm gồm hệ số suy giảm khối, nguyên tử số hiệu dụng và mật độ electron hiệu dụng cho 4 nhóm vật liệu phức hợp trên dải năng lượng 100–2000 keV.
  • Khẳng định tính ưu việt của phương pháp mô phỏng Monte Carlo MCNP5 với độ lệch tương đối so với cơ sở dữ liệu NIST luôn dưới 5,0% và dưới 8,1% so với thực nghiệm.
  • Chứng minh tiềm năng ứng dụng to lớn của hệ thủy tinh borate chứa bismuth oxit trong việc chế tạo vật liệu che chắn bức xạ hiệu năng cao, nhẹ và thân thiện với môi trường.
  • Thiết lập thành công mô hình hàm số hồi quy phi tuyến với hệ số tương quan đạt xấp xỉ 0,999, cho phép nội suy nhanh các tham số suy giảm tại mọi mức năng lượng bức xạ.
  • Hoàn thiện quy trình xử lý phổ gamma cải tiến bằng Colegram, cung cấp phương pháp luận chuẩn xác cho các nghiên cứu đo lường hạt nhân tiếp theo.

Các đơn vị nghiên cứu, cơ sở y tế và doanh nghiệp sản xuất vật liệu che chắn bức xạ hãy liên hệ và áp dụng ngay các giải pháp mô phỏng tính toán từ công trình này để tối ưu hóa thiết kế an toàn bức xạ, nâng cao hiệu quả kỹ thuật và tiết kiệm tối đa chi phí đầu tư.