Tổng quan nghiên cứu

Nhu cầu ứng dụng bức xạ và đồng vị phóng xạ trong y tế tại Việt Nam đang tăng trưởng mạnh mẽ với sự hình thành của mạng lưới 35 cơ sở y học hạt nhân cùng hơn 45 thiết bị xạ hình hiện đại trên cả nước. Căn cứ theo Quyết định số 01/2006/QĐ-TTg của Thủ tướng Chính phủ về chiến lược ứng dụng năng lượng nguyên tử vì mục đích hòa bình, mục tiêu trọng tâm đặt ra là phấn đấu tự chủ sản xuất 50% nhu cầu đồng vị phóng xạ và dược chất phóng xạ trong nước. Hiện tại, Lò nghiên cứu hạt nhân Đà Lạt với công suất định danh 500 kW đang vận hành cung cấp ổn định từ 350 đến 400 Ci mỗi năm cho 25 khoa và trung tâm ung bướu lớn trên toàn quốc. Các sản phẩm chủ lực bao gồm Iodine-131, Technetium-99m, Phosphor-32, Lutetium-177 và Samarium-153.

Mặc dù vậy, đồng vị phóng xạ Iodine-125—nguồn phát bức xạ gamma mềm đóng vai trò cốt lõi trong kỹ thuật cấy hạt phóng xạ điều trị ung thư tuyến tiền liệt và u não—vẫn chưa được sản xuất trong nước và hoàn toàn phụ thuộc vào nguồn nhập khẩu với chi phí đắt đỏ. Nhằm giải quyết triệt để rào cản công nghệ này, luận văn thạc sĩ chuyên ngành Vật lý nguyên tử của tác giả Đoàn Thị Thu Hiền tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội đã thực hiện đề tài tính toán quá trình sản xuất Iodine-125 từ khí Xenon-124 chiếu xạ bởi dòng nơtron nhiệt trong lò phản ứng. Nghiên cứu tập trung giải bài toán động học chuyển hóa hạt nhân tại Lò phản ứng hạt nhân Đà Lạt trong giai đoạn 2017–2019, thiết lập mô hình tính toán số lượng hạt nhân và tối ưu hóa thời gian chiếu nhằm giảm thiểu chi phí thử nghiệm thực tế và cung cấp cơ sở dữ liệu tin cậy cho việc triển khai sản xuất lâm sà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 trên nền tảng lý thuyết tương tác nơtron với vật chất và phản ứng hạt nhân bắt bức xạ nơtron nhiệt $(n,\gamma)$. Chu trình sản xuất Iodine-125 diễn ra qua phản ứng kích hoạt hạt nhân bia khí Xenon theo sơ đồ chuyển hóa: $^{124}\text{Xe}(n,\gamma)^{125}\text{Xe}$. Hạt nhân trung gian $^{125}\text{Xe}$ có chu kỳ bán rã $T_{1/2} = 18\text{ giờ}$, tự động phân rã theo kiểu bắt điện tử (EC) và phát bức xạ positron $\beta^+$ để tạo thành đồng vị mục tiêu $^{125}\text{I}$. Sau đó, $^{125}\text{I}$ tiếp tục phân rã với chu kỳ bán rã $T_{1/2} = 59.408\text{ ngày}$ về đồng vị bền $^{125}\text{Te}$, giải phóng tia gamma năng lượng thấp khoảng 35 keV. Tiết diện phản ứng bắt nơtron nhiệt của $^{124}\text{Xe}$ đạt giá trị rất cao là $\sigma = 165 \pm 20\text{ barn}$.

Khung lý thuyết thứ hai là mô hình tích lũy và phân rã phóng xạ vi phân liên tục dựa trên hệ phương trình vi phân Bateman. Mô hình này mô tả đồng thời tốc độ biến thiên số lượng hạt nhân bia, sự hình thành sản phẩm trung gian, sự tích lũy của $^{125}\text{I}$ và phản ứng phụ bắt nơtron thứ cấp tạo ra tạp chất phóng xạ nguy hại là $^{126}\text{I}$ ($T_{1/2} = 13\text{ ngày}$) theo kênh phản ứng $^{125}\text{I}(n,\gamma)^{126}\text{I}$. Việc kiểm soát chặt chẽ tỷ lệ hoạt độ giữa $^{126}\text{I}$ và $^{125}\text{I}$ là chỉ số chất lượng quyết định tính an toàn bức xạ y tế của sản phẩm cuối cùng.

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

Nguồn dữ liệu tính toán được thiết lập từ các thông số thực tế của cấu hình vùng hoạt Lò phản ứng hạt nhân Đà Lạt. Tác giả thực hiện khảo sát trên 2 nhóm kích thước mẫu bia tiêu biểu theo phương pháp chọn mẫu kỹ thuật định lượng có chủ đích: nhóm 1 sử dụng bia khí Xenon tự nhiên có khối lượng 50 g (tỷ lệ đồng vị $^{124}\text{Xe}$ chiếm 0.096%) và nhóm 2 sử dụng bia khí Xenon làm giàu khối lượng 0.5 g với độ giàu $^{124}\text{Xe}$ đạt trên 99%. Hai vị trí chiếu xạ thực nghiệm được khảo sát gồm kênh chiếu biên với thông lượng nơtron nhiệt $\Phi = 9 \times 10^{11}\text{ n/cm}^2/\text{s}$ và bẫy nơtron trung tâm với thông lượng nơtron nhiệt cực đại $\Phi = 2 \times 10^{13}\text{ n/cm}^2/\text{s}$.

Tác giả lựa chọn phương pháp tích phân số Runge-Kutta bậc 4 (RK4) để giải hệ phương trình vi phân động học phi tuyến mô tả tốc độ tích lũy ròng hạt nhân. Lý do lựa chọn thuật toán RK4 bắt nguồn từ yêu cầu đảm bảo sai số địa phương ở bậc $O(h^4)$, triệt tiêu hiện tượng mất ổn định số khi xử lý hệ phương trình chuỗi phân rã có độ chênh lệch thời gian bán rã lớn từ 18 giờ của $^{125}\text{Xe}$ đến gần 60 ngày của $^{125}\text{I}$. Toàn bộ thuật toán được đóng gói thành chương trình phần mềm với giao diện trực quan, tự động hóa các bước nhập dữ liệu đầu vào và kết xuất kết quả mô phỏng.

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

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

Thứ nhất, việc sử dụng bia khí Xenon làm giàu có độ tinh khiết $^{124}\text{Xe} > 99%$ tạo ra sự vượt trội hoàn toàn về chất lượng sản phẩm so với bia Xenon tự nhiên (chỉ chứa 0.096% $^{124}\text{Xe}$). Hoạt độ riêng của $^{125}\text{I}$ thu được từ bia làm giàu cao gấp hơn 1000 lần so với bia tự nhiên, đạt trạng thái không chất mang với hoạt độ riêng cực cao đáp ứng hoàn hảo tiêu chuẩn chế tạo dược chất phóng xạ.

Thứ hai, thời gian chiếu xạ tối ưu được xác định rõ ràng qua mô hình động học. Khi chiếu mẫu bia Xenon giàu 0.5 g tại bẫy nơtron ($\Phi = 2 \times 10^{13}\text{ n/cm}^2/\text{s}$), hoạt độ của $^{125}\text{I}$ tăng trưởng gần như tuyến tính trong 15 ngày đầu chiếu xạ trước khi tốc độ tích lũy chậm lại do sự cân bằng giữa tốc độ sinh và tốc độ bắt nơtron thứ cấp.

Thứ ba, quy trình làm nguội phân rã từ 3 đến 5 ngày sau khi kết thúc đợt chiếu 15 ngày giúp hàm lượng tạp chất phóng xạ nguy hại $^{126}\text{I}$ giảm xuống dưới ngưỡng an toàn $0.005%$ tổng hoạt độ.

Thứ tư, tính toán sản lượng quy mô vận hành lò 120 ngày mỗi năm cho thấy, chỉ với 1 mẻ chiếu 0.5 g khí Xenon giàu, Lò phản ứng hạt nhân Đà Lạt có thể tạo ra sản lượng hàng chục Curie $^{125}\text{I}$ sạch phóng xạ, đủ cung ứng cho nhu cầu xạ trị áp sát hàng năm của hàng trăm bệnh nhân trong nước.

Thảo luận kết quả

Các dữ liệu động học tích lũy hạt nhân được mô tả sinh động qua các đồ thị đường cong phân rã và bảng tổng hợp đối chiếu số liệu. Đồ thị hoạt độ theo thời gian chiếu xạ cho thấy khi kéo dài thời gian chiếu vượt quá điểm tối ưu 15 ngày, lượng $^{125}\text{I}$ tiếp tục bị kích hoạt chuyển hóa thành $^{126}\text{I}$. Đây là đồng vị không mong muốn vì nó phát ra các tia bức xạ gamma cứng năng lượng cao từ 388 keV đến 666 keV, gây tổn thương nghiêm trọng đến các mô lành xung quanh khối u khi cấy hạt phóng xạ vào cơ thể bệnh nhân. Do đó, việc chuyển mẫu khí ra buồng phân rã từ 3 đến 5 ngày là giải pháp then chốt để $^{125}\text{Xe}$ ($T_{1/2} = 18\text{ giờ}$) chuyển hóa hoàn toàn thành $^{125}\text{I}$ mà không sinh thêm tạp chất $^{126}\text{I}$.

Khi so sánh kết quả tính toán của phần mềm với bộ số liệu thực nghiệm chuẩn do Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) công bố đối với cả bia tự nhiên và bia giàu, mức độ sai lệch chỉ dao động dưới 3%. Sự tương thích cao này khẳng định tính chính xác tuyệt đối của mô hình toán học và thuật toán RK4 đã xây dựng, tạo cơ sở khoa học vững chắc để áp dụng vào thực tế vận hành lò phản ứng.

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

Thứ nhất, tiến hành thiết kế và chế tạo viên nang (capsule) bằng hợp kim nhôm chuyên dụng chịu áp suất cao nhằm đóng gói an toàn từ 0.5 g đến 1.0 g khí Xenon làm giàu $>99%$, đảm bảo độ kín tuyệt đối với tỷ lệ rò rỉ dưới 0.1% khi chiếu xạ trong môi trường bẫy nơtron. Chủ thể thực hiện là Viện Nghiên cứu Hạt nhân Đà Lạt với lộ trình triển khai trong 12 tháng.

Thứ hai, xây dựng và hoàn thiện dây chuyền chiết tách hóa học tự động trong buồng nóng (hotcell) chuyên dụng, đảm bảo hiệu suất thu hồi sản phẩm Iodine-125 đạt tối thiểu 95% và tỷ lệ thất thoát hoạt độ dưới 5%. Chủ thể thực hiện là Trung tâm Dược phóng xạ thuộc Viện Nghiên cứu Hạt nhân phối hợp với Bộ Khoa học và Công nghệ trong thời gian 18 tháng.

Thứ ba, triển khai đầu tư công nghệ chế tạo hạt vi xạ trị Iodine-125 bọc vỏ titan kích thước chuẩn $4.5 \times 0.8\text{ mm}$ phục vụ xạ trị áp sát u tuyến tiền liệt, hướng tới mục tiêu tự chủ 50% nhu cầu hạt phóng xạ tại các bệnh viện ung bướu vào năm 2028. Chủ thể thực hiện là liên danh giữa Bệnh viện Trung ương Quân đội 108, Bệnh viện Bạch Mai và Viện Năng lượng nguyên tử Việt Nam với thời gian thực hiện 24 tháng.

Thứ tư, chuẩn hóa và ban hành quy chuẩn kiểm định chất lượng dược chất phóng xạ, kiểm soát tỷ lệ tạp chất $^{126}\text{I} < 0.005%$ theo đúng tiêu chuẩn an toàn bức xạ quốc tế của IAEA. Chủ thể thực hiện là Cục An toàn bức xạ và hạt nhân phối hợp với Bộ Y tế trong thời gian 12 tháng.

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

Nhóm 1: Các nhà vật lý hạt nhân và kỹ sư vận hành lò phản ứng nghiên cứu. Luận văn cung cấp mô hình toán học số RK4 và thuật toán tính toán chuỗi phân rã, giúp tối ưu hóa chế độ chiếu xạ cho 278 lò nghiên cứu tương tự trên thế giới cũng như mở rộng cho các đồng vị y học khác như Samarium-153 và Lutetium-177.

Nhóm 2: Bác sĩ chuyên khoa ung bướu, bác sĩ y học hạt nhân và kỹ sư vật lý y khoa. Người đọc sẽ nắm bắt tường tận đặc tính vật lý của hạt phóng xạ $^{125}\text{I}$ (kích thước $4.5 \times 0.8\text{ mm}$, năng lượng gamma mềm 35 keV, chu kỳ bán rã 59.4 ngày) để xây dựng phác đồ tính liều xạ trị áp sát chính xác, tối ưu hóa hiệu quả điều trị u tuyến tiền liệt và u não.

Nhóm 3: Các nhà quản lý chính sách và chuyên gia hoạch định chiến lược năng lượng nguyên tử. Tài liệu cung cấp cơ sở kỹ thuật xác đáng để phê duyệt các dự án đầu tư sản xuất đồng vị tự chủ trong nước, phục vụ hiệu quả cho mạng lưới 35 cơ sở y học hạt nhân và 6 trung tâm máy gia tốc cyclotron trên toàn quốc.

Nhóm 4: Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Vật lý hạt nhân và Kỹ thuật hạt nhân. Luận văn là tài liệu tham khảo chuẩn mực về phương pháp mô phỏng số động học lò phản ứng và kỹ thuật lập trình giao diện tính toán tương tác bức xạ.

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

Tại sao Iodine-125 lại được ưu tiên hàng đầu trong kỹ thuật cấy hạt xạ trị ung thư tuyến tiền liệt?

Hạt phóng xạ Iodine-125 phát ra bức xạ gamma mềm năng lượng thấp 35 keV với tầm bắn trong mô bệnh chỉ vài milimét và có chu kỳ bán rã 59.4 ngày. Đặc tính này cho phép tập trung liều hấp thụ cực cao để tiêu diệt khối u tại chỗ mà không gây tổn thương mô lành xung quanh, đồng thời hạt có thể lưu lại an toàn vĩnh viễn trong cơ thể bệnh nhân mà không cần phẫu thuật lấy ra.

Phản ứng hạt nhân cơ bản nào được sử dụng để chế tạo Iodine-125 trong lò phản ứng?

Quá trình chế tạo dựa trên phản ứng bắt nơtron nhiệt của khí Xenon-124 có tiết diện phản ứng $\sigma = 165\text{ barn}$ theo sơ đồ $^{124}\text{Xe}(n,\gamma)^{125}\text{Xe}$. Đồng vị trung gian $^{125}\text{Xe}$ có chu kỳ bán rã 18 giờ sẽ tự động phân rã bắt điện tử để chuyển đổi thành sản phẩm mục tiêu $^{125}\text{I}$.

Vì sao cần phải kiểm soát nghiêm ngặt hàm lượng tạp chất Iodine-126 trong sản phẩm?

Iodine-126 sinh ra từ phản ứng bắt nơtron thứ cấp của $^{125}\text{I}$, phát ra bức xạ gamma cứng năng lượng cao từ 388 keV đến 666 keV. Nếu tỷ lệ $^{126}\text{I}$ vượt quá ngưỡng $0.005%$, bức xạ năng lượng cao sẽ đâm xuyên mạnh gây bỏng xạ và hoại tử các cơ quan lành lân cận vị trí cấy hạt của người bệnh.

Sử dụng bia khí Xenon làm giàu mang lại ưu thế gì so với Xenon tự nhiên?

Khí Xenon tự nhiên chỉ chứa 0.096% đồng vị $^{124}\text{Xe}$, trong khi bia làm giàu đạt độ tinh khiết trên 99%. Sử dụng 0.5 g Xenon giàu giúp nâng cao hoạt độ riêng của $^{125}\text{I}$ lên hơn 1000 lần, đạt trạng thái không chất mang và hạn chế tối đa sự hình thành các sản phẩm phụ không mong muốn.

Lò nghiên cứu hạt nhân Đà Lạt có đủ điều kiện thông lượng để sản xuất Iodine-125 không?

Hoàn toàn đáp ứng tốt. Với công suất định danh 500 kW, bẫy nơtron của Lò phản ứng hạt nhân Đà Lạt đạt thông lượng nơtron nhiệt lên tới $2 \times 10^{13}\text{ n/cm}^2/\text{s}$, đảm bảo khả năng chiếu xạ sản xuất hàng chục Curie $^{125}\text{I}$ tinh khiết mỗi năm phục vụ nhu cầu điều trị của hệ thống y tế trong nước.

Kết luận

  • Xây dựng thành công công cụ tính toán số học sử dụng thuật toán Runge-Kutta bậc 4 giải hệ phương trình động học tích lũy và phân rã với độ sai lệch dưới 3% so với chuẩn IAEA.
  • Xác định cấu hình chiếu xạ tối ưu với 0.5 g khí Xenon làm giàu $>99%$ tại bẫy nơtron nhiệt có thông lượng $2 \times 10^{13}\text{ n/cm}^2/\text{s}$.
  • Thiết lập chu trình sản xuất chuẩn gồm 15 ngày chiếu xạ kết hợp 3 đến 5 ngày lưu buồng phân rã để giảm tạp chất $^{126}\text{I}$ xuống dưới ngưỡng an toàn $0.005%$.
  • Chứng minh tính khả thi tuyệt đối của việc triển khai sản xuất thương mại Iodine-125 trên Lò phản ứng hạt nhân Đà Lạt công suất 500 kW.
  • Mở rộng thành công mô hình thuật toán cho việc nghiên cứu sản xuất các đồng vị y học quan trọng khác như Samarium-153 và Lutetium-177.

Đóng góp cốt lõi của luận văn là đã đặt nền móng lý thuyết và công nghệ hoàn chỉnh cho bài toán sản xuất Iodine-125 phục vụ điều trị ung thư tại Việt Nam. Trong lộ trình 2026–2028, việc nhanh chóng hoàn thiện buồng chiếu mẫu và dây chuyền đóng gói hạt phóng xạ chuẩn lâm sàng sẽ tạo bước đột phá trong tự chủ y học hạt nhân nước nhà. Hãy liên hệ và chia sẻ ngay công trình nghiên cứu giá trị này để cùng thúc đẩy các giải pháp công nghệ bức xạ tiên tiến cứu sống người bệnh!