Tổng quan nghiên cứu

Chiến lược phát triển và ứng dụng năng lượng nguyên tử vì mục đích hòa bình theo Quyết định số 01/2006/QĐ-TTg của Thủ tướng Chính phủ đã đặt mục tiêu bảo đảm tự chủ 50% nhu cầu đồng vị và dược chất phóng xạ trong nước. Đến nay, cả nước đã hình thành mạng lưới gồm 35 cơ sở y học hạt nhân với hơn 45 thiết bị xạ hình hiện đại, 6 hệ máy gia tốc Cyclotron và Lò phản ứng hạt nhân Đà Lạt công suất 500 kW cung cấp từ 350 đến 400 Ci mỗi năm cho 25 trung tâm ung bướu trên toàn quốc.

Mặc dù vậy, đồng vị phóng xạ Iodine-125 với chu kỳ bán rã 59,4 ngày, phát bức xạ gamma mềm 35,5 keV và tia X đặc trưng 28 keV – vốn là tiêu chuẩn vàng 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 thể sản xuất trong nước mà hoàn toàn phụ thuộc vào nguồn nhập khẩu. Tại Bệnh viện Bạch Mai, phương pháp cấy hạt Iodine-125 kích thước 4,5 x 0,8 mm triển khai từ năm 2015 đã mang lại kết quả khả quan khi giúp hơn 10 bệnh nhân ung thư tuyến tiền liệt giảm 100% triệu chứng bệnh và bảo tồn trọn vẹn mô lành xung quanh khối u.

Xuất phát từ yêu cầu thực tiễn cấp bách đó, luận văn thạc sĩ chuyên ngành Vật lý nguyên tử tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội đã tập trung giải quyết bài toán: "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 hướng đến mục tiêu xây dựng mô hình toán học giải tích, lập trình thuật toán mô phỏng động học tạo thành và phân rã các đồng vị phóng xạ, từ đó xác định bộ thông số công nghệ tối ưu cho Lò phản ứng hạt nhân Đà Lạt ở các dải thông lượng từ 9.10^12 n/cm^2/s đến 2.10^13 n/cm^2/s. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật to lớn, tạo tiền đề khoa học vững chắc giúp giảm thiểu rủi ro vận hành mẫu nóng và tiết kiệm hàng tỷ đồng chi phí thử nghiệm thực tế.

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

Khung lý thuyết áp dụng

Quá trình sản xuất Iodine-125 dựa trên phản ứng kích hoạt nơtron nhiệt hạt nhân Xenon-124: 124Xe(n,γ)125Xe. Đồng vị trung gian Xenon-125 có chu kỳ bán rã 18 giờ, phân rã thông qua cơ chế bắt electron (EC) để chuyển hóa thành Iodine-125 có chu kỳ bán rã 59,4 ngày. Tiết diện kích hoạt nơtron nhiệt của Xenon-124 đạt 128 barn cùng tích phân cộng hưởng nơtron trên nhiệt đạt 3600 barn.

Song song với quá trình tạo thành hạt nhân mục tiêu, phản ứng phụ thứ cấp xảy ra do sự bắt nơtron của chính sản phẩm: 125I(n,γ)126I với tiết diện rất lớn lên tới 894 barn và tích phân cộng hưởng 13730 barn. Đồng vị tạp chất Iodine-126 có chu kỳ bán rã 13 ngày, phát tia gamma năng lượng cao gây độc tính phóng xạ nghiêm trọng đối với mô lành của bệnh nhân. Tiêu chuẩn quốc tế của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) đòi hỏi tỷ lệ tạp chất Iodine-126/Iodine-125 trong dược chất thành phẩm bắt buộc phải dưới mức 1% (tối ưu là dưới 0,9%) và hoạt độ riêng phải đạt trên 600 GBq/mg (tương đương trên 16000 Ci/g).

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

Nghiên cứu sử dụng nguồn dữ liệu vật lý hạt nhân chuẩn quốc tế từ IAEA, bao gồm hệ tiết diện phản ứng nơtron nhiệt, tích phân cộng hưởng và hằng số phân rã phóng xạ của 9 hạt nhân đồng vị trong chuỗi phản ứng Xenon và Iodine.

Về cỡ mẫu và phương pháp chọn mẫu mô phỏng, đề tài tiến hành khảo sát trên hai nhóm vật liệu bia chính: nhóm bia khí Xenon tự nhiên (hàm lượng Xenon-124 chiếm 0,095%) với khối lượng bia từ 15 gram đến 50 gram, và nhóm bia khí Xenon làm giàu (độ giàu Xenon-124 đạt từ 1% đến trên 99,9%) với khối lượng bia tiêu chuẩn từ 0,4 gram đến 0,5 gram. Cỡ mẫu và dải khối lượng này được lựa chọn dựa trên dung tích thể tích thực của ống chiếu nhôm chuẩn và giới hạn an toàn áp suất khí nén trong kênh chiếu lò phản ứng hạt nhân.

Phương pháp phân tích cốt lõi là thiết lập hệ phương trình vi phân động học phi tuyến mô tả tốc độ biến thiên mật độ hạt nhân theo thời gian. Thuật toán tích phân số Runge-Kutta bậc 4 (RK4) được lập trình hướng đối tượng trên ngôn ngữ Visual Basic 6.0 chạy trên hệ điều hành Windows. Lý do lựa chọn thuật toán RK4 là nhờ tính ổn định số vượt trội, sai số cục bộ ở cấp O(h^5), giúp tính toán chính xác đồng thời cả 9 chuỗi tích lũy hạt nhân trong khoảng thời gian chiếu xạ từ 10 đến 300 giờ và thời gian làm nguội từ 20 đến 90 ngày suốt giai đoạn nghiên cứu 2017 - 2019.

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

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

Kết quả kiểm chứng mô hình tính toán với số liệu công bố chuẩn của IAEA cho thấy độ tương thích xuất sắc với sai số tương đối dưới 0,3%. Cụ thể, khi chiếu 15 gram Xenon tự nhiên ở thông lượng 5.10^13 n/cm^2/s trong 200 giờ và làm nguội 40 ngày, sản lượng Iodine-125 tính toán đạt 29,074 GBq so với 29,15 GBq của IAEA; hàm lượng tạp chất Iodine-126 đạt 0,996% và hoạt độ riêng đạt 608,70 GBq/mg. Đối với bia 0,4 gram Xenon giàu trên 99% chiếu 10 giờ, sản lượng đạt 53,98 GBq so với số liệu 50,71 GBq của IAEA.

Đối với phương pháp sản xuất theo mẻ gián đoạn (Batch Process), khi chiếu 50 gram Xenon tự nhiên trong 20 ngày, việc tăng thông lượng nơtron từ 5.10^12 n/cm^2/s lên 10^14 n/cm^2/s giúp tổng hoạt độ Iodine-125 tăng vọt từ 1 Ci lên hơn 14 Ci. Tuy nhiên, thông lượng cao làm tăng mạnh tỷ lệ cháy hạt nhân, khiến hàm lượng tạp chất Iodine-126 tăng từ 1% lên tới 20%, đòi hỏi thời gian làm nguội kéo dài hơn 60 ngày. Trong khi đó, với bia 0,5 gram Xenon làm giàu 99,9% chiếu 17 giờ ở thông lượng 4.10^14 n/cm^2/s, hoạt độ thu được đạt tới 26111,59 mCi với tỷ lệ tạp chất Iodine-126 chỉ 0,82% và hoạt độ riêng đạt mức cực cao 17326,58 Ci/g.

Đối với phương pháp sản xuất theo chu trình lặp tuần hoàn (Loop Process), việc chiếu 15 gram Xenon tự nhiên ở thông lượng 10^14 n/cm^2/s theo chu kỳ 17 giờ chiếu xạ và 113 giờ phân rã qua 10 vòng lặp (tổng thời gian 47 ngày) đã tạo ra bước đột phá công nghệ: tổng hoạt độ Iodine-125 đạt trên 5,8 Ci (5809,23 mCi), tỷ lệ tạp chất Iodine-126 trong sản phẩm cuối cùng chỉ chiếm 0,413% (thấp hơn nhiều so với ngưỡng 1%) và hoạt độ riêng đạt 16786,08 Ci/g.

Thảo luận kết quả

Các kết quả phân tích động học được thể hiện trực quan thông qua các biểu đồ đường tích lũy hoạt độ và bảng đối sánh đa thông số. Đồ thị biểu diễn hoạt độ Iodine-125 trong buồng phân rã chỉ ra rằng thời gian làm nguội 113 giờ tương đương hơn 6 chu kỳ bán rã của Xenon-125, bảo đảm hơn 90% lượng Xenon-125 chuyển hóa hoàn toàn thành Iodine-125 trước khi thu hồi khí Xenon để tái chiếu xạ.

Hiện tượng suy giảm hoạt độ riêng từ 17218 Ci/g ở vòng lặp thứ nhất xuống còn 16786 Ci/g ở vòng lặp thứ 10 được giải thích là do sự tích tụ không thể tránh khỏi của đồng vị bền Iodine-127 với hàm lượng 11,1 phần triệu gram (ppm). So với các công bố quốc tế từ các trung tâm nghiên cứu tại Hoa Kỳ, Nga và Ấn Độ, phương pháp tuần hoàn khí Xenon liên tục giúp giảm thiểu tối đa hiện tượng tự hấp thụ và cháy nơtron của Iodine-125 trong vùng hoạt, cho phép tận dụng nguồn khí Xenon tự nhiên có chi phí thấp mà vẫn thu được dược chất đạt tiêu chuẩn y tế khắt khe.

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

Thứ nhất, Viện Nghiên cứu Hạt nhân Đà Lạt cần ưu tiên thiết kế và lắp đặt hệ thống chiếu xạ vòng lặp tuần hoàn kín chuyên dụng nối từ vùng hoạt lò phản ứng đến buồng làm lạnh thu hồi khí trong giai đoạn 2026 - 2028, hướng tới mục tiêu sản xuất ổn định từ 100 đến 150 Ci Iodine-125 mỗi năm đáp ứng nhu cầu nội địa.

Thứ hai, tối ưu hóa quy trình vận hành chiếu xạ thực tế tại Lò phản ứng hạt nhân Đà Lạt theo hai phương án: bố trí chiếu bia Xenon giàu (khối lượng 0,5 gram, độ giàu trên 99%) tại hốc bẫy nơtron trung tâm với thông lượng 2.10^13 n/cm^2/s trong thời gian từ 15 đến 17 giờ; hoặc bố trí chiếu bia Xenon tự nhiên (khối lượng 50 gram) tại kênh chiếu ướt với thông lượng 9.10^12 n/cm^2/s, bảo đảm thời gian làm nguội phân rã tối thiểu 113 giờ để giữ tỷ lệ tạp chất Iodine-126 dưới 0,5%.

Thứ ba, Bộ Khoa học và Công nghệ cùng Viện Năng lượng nguyên tử Việt Nam cần đầu tư nguồn kinh phí xây dựng hệ thống bẫy lạnh chân không để thu hồi khí Xenon-124 với hiệu suất thu hồi đạt trên 95%, giúp tiết giảm từ 30% đến 40% chi phí nguyên liệu đắt đỏ trong các đợt chiếu xạ dài ngày.

Thứ tư, Cục An toàn Bức xạ và Hạt nhân phối hợp với các bệnh viện chuyên khoa ung bướu như Bệnh viện Bạch Mai và Bệnh viện Trung ương Quân đội 108 ban hành quy chuẩn kỹ thuật quốc gia về đóng gói hạt titan chứa Iodine-125 chuẩn GMP-WHO trước năm 2027 nhằm hoàn thiện chuỗi giá trị từ sản xuất lò hạt nhân đến ứng dụng lâm sàng.

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

Nhóm kỹ sư vật lý lò và chuyên viên vận hành hạt nhân: Nắm vững bộ công cụ giải thuật và mô hình phân bố dòng nơtron để tính toán hiệu suất kích hoạt, thời gian cháy bia và tối ưu hóa thời gian dừng lò cho các kênh chiếu mẫu nóng.

Nhóm chuyên gia hóa - dược phóng xạ: Ứng dụng các quy trình công nghệ chưng cất, bẫy khí hóa lỏng và kỹ thuật hấp thụ Iodine-125 trong dung dịch kiềm NaOH 0,05M nhằm điều chế dung dịch NaI đạt hoạt độ riêng trên 16000 Ci/g.

Nhóm bác sĩ xạ trị và kỹ sư y vật lý tại các bệnh viện ung bướu: Nắm bắt chi tiết các thông số bức xạ học, phổ năng lượng tia gamma 35,5 keV và mức độ tinh khiết hạt nhân của hạt phóng xạ Iodine-125 để lập kế hoạch phân liều xạ trị áp sát chính xác cho khối u tuyến tiền liệt.

Nhóm học viên cao học và nghiên cứu sinh ngành Vật lý hạt nhân: Kế thừa thuật toán tích phân số Runge-Kutta bậc 4 và cấu trúc mã nguồn Visual Basic để mở rộng mô hình tính toán cho các đồng vị phóng xạ tiềm năng khác như Samarium-153, Strontium-89, Rhenium-186 và Rhenium-188.

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

Iodine-125 có ưu điểm vượt trội gì trong điều trị ung thư tuyến tiền liệt so với các đồng vị phóng xạ khác?
Iodine-125 phát tia gamma năng lượng thấp 35,5 keV với quãng chạy trong mô mềm chỉ vài milimét và chu kỳ bán rã 59,4 ngày. Đặc tính này cho phép tiêu diệt triệt để tế bào ung thư tại chỗ với liều hấp thụ tập trung trong vòng 6 tháng mà hoàn toàn không gây tổn thương các mô lành xung quanh như trực tràng và bàng quang.

Tạp chất Iodine-126 xuất hiện từ đâu và vì sao phải khống chế dưới 1%?
Iodine-126 sinh ra từ phản ứng phụ do Iodine-125 vừa tạo thành tiếp tục bắt nơtron trong lò. Do Iodine-126 phát bức xạ gamma năng lượng cao tới 1,4 MeV gây tổn thương mô lành nghiêm trọng, tiêu chuẩn y tế quốc tế bắt buộc phải kiểm soát tỷ lệ Iodine-126/Iodine-125 dưới 1% thông qua việc tối ưu thời gian chiếu và làm nguội.

Sản xuất Iodine-125 theo vòng lặp tuần hoàn vượt trội hơn sản xuất theo mẻ như thế nào?
Phương pháp vòng lặp tuần hoàn đưa khí Xenon ra khỏi vùng chiếu sau mỗi 17 giờ để phân rã trong buồng riêng, giúp sản phẩm Iodine-125 không bị nơtron tiếp tục bắn phá. Kỹ thuật này nâng hoạt độ riêng lên trên 16700 Ci/g, kiểm soát tạp chất Iodine-126 ở mức 0,413% ngay cả khi sử dụng nguồn bia Xenon tự nhiên.

Thời gian làm nguội 113 giờ trong buồng phân rã có ý nghĩa then chốt ra sao?
Chu kỳ bán rã của Xenon-125 là 18 giờ. Khoảng thời gian làm nguội 113 giờ tương đương hơn 6 chu kỳ bán hủy, giúp hơn 90% hạt nhân Xenon-125 phân rã hoàn toàn thành Iodine-125 và bám dính vào thành buồng hấp thụ trước khi khí Xenon chưa phân rã được bơm tuần hoàn trở lại vùng hoạt.

Chương trình tính toán trong luận văn có thể mở rộng cho các đồng vị nào khác?
Mô hình toán học hướng đối tượng và thuật toán Runge-Kutta bậc 4 trong luận văn được thiết kế dạng module mở rộng. Chương trình hoàn toàn có thể áp dụng trực tiếp để mô phỏng chính xác động học sản xuất các đồng vị y học quan trọng khác như Strontium-89, Samarium-153, Rhenium-186 và Rhenium-188 trong lò phản ứng hạt nhân.

Kết luận

  • Luận văn đã xây dựng thành công hệ phương trình vi phân mô tả chính xác quá trình biến đổi hạt nhân của 9 đồng vị trong chuỗi sản xuất Iodine-125 từ bia khí Xenon-124.
  • Ứng dụng xuất sắc thuật toán tích phân số Runge-Kutta bậc 4 trên phần mềm chuyên dụng, cho kết quả kiểm chứng hoàn toàn trùng khớp với số liệu chuẩn của IAEA với độ tin cậy cao.
  • Chứng minh tính ưu việt của phương pháp chiếu xạ vòng lặp tuần hoàn khi tạo ra hơn 5,8 Ci Iodine-125 với độ sạch tạp chất Iodine-126 đạt 0,413% và hoạt độ riêng vượt 16780 Ci/g.
  • Xác lập bộ thông số công nghệ chiếu xạ và phân rã tối ưu hóa hoàn toàn tương thích với cấu hình thực tế của Lò phản ứng hạt nhân Đà Lạt.
  • Cung cấp cơ sở lý thuyết và công cụ mô phỏng chuẩn xác phục vụ trực tiếp chiến lược tự chủ sản xuất dược chất phóng xạ phục vụ y tế tại Việt Nam giai đoạn 2026 - 2030.

Để thúc đẩy việc hiện thực hóa các kết quả tính toán vào thực tiễn, các nhà khoa học, kỹ sư hạt nhân và cơ sở y tế có thể kết nối, ứng dụng mô hình công nghệ này nhằm sớm đưa dây chuyền sản xuất hạt phóng xạ Iodine-125 nội địa vào vận hành thường quy.