Tổng quan nghiên cứu

Nhu cầu ứng dụng đồng vị phóng xạ trong y tế và công nghiệp tại Việt Nam ngày càng tăng cao theo định hướng của Quyết định số 01/2006/QĐ-TTg về Chiến lược phát triển ứng dụng năng lượng nguyên tử vì mục đích hòa bình. Đến năm 2019, toàn quốc đã vận hành 35 cơ sở y học hạt nhân cùng hơn 45 thiết bị xạ hình hiện đại như SPECT và PET-CT, song song với hệ thống 6 máy gia tốc Cyclotron tại các bệnh viện lớn. Mặc dù Lò phản ứng hạt nhân Đà Lạt (công suất danh định 500 kW) cung cấp đều đặn 350 đến 400 Ci mỗi năm cho 20 loại đồng vị và dược chất phóng xạ, việc sản xuất đồng vị Iodine-125 ($^{125}\text{I}$) vẫn chưa được thương mại hóa trong nước.

Đồng vị $^{125}\text{I}$ đóng vai trò cốt lõi trong phương pháp xạ trị áp sát cấy hạt phóng xạ (kích thước chuẩn 4,5 x 0,8 mm, phát tia gamma mềm 35 keV, chu kỳ bán rã 59,408 ngày). Điển hình tại Bệnh viện Bạch Mai, kỹ thuật này đã điều trị thành công cho hơn 10 bệnh nhân ung thư tuyến tiền liệt với tỷ lệ giảm 100% triệu chứng lâm sàng. Nhằm chấm dứt tình trạng phụ thuộc hoàn toàn vào nguồn nhập khẩu đắt đỏ, luận văn tập trung giải quyết bài toán mô hình hóa toán học và tính toán tối ưu các thông số công nghệ để sản xuất $^{125}\text{I}$ từ bia khí Xenon-124 ($^{124}\text{Xe}$) được chiếu xạ bởi dòng nơtron nhiệt trong lò phản ứng.

Mục tiêu cụ thể của công trình là xây dựng hệ thống phần mềm tính toán động học tích lũy ròng cho 9 hạt nhân phóng xạ, từ đó xác lập chế độ chiếu xạ và phân rã tối ưu tại hai vị trí trọng yếu của Lò Đà Lạt: hốc bẫy nơtron (thông lượng $\Phi = 2 \cdot 10^{13}\text{ n/cm}^2\text{/s}$) và kênh chiếu ướt ($\Phi = 9 \cdot 10^{12}\text{ n/cm}^2\text{/s}$). Kết quả nghiên cứu là luận cứ kỹ thuật trực tiếp để định hướng công nghệ sản xuất, giúp giảm thiểu rủi ro xử lý mẫu nóng và tối ưu hóa chi phí vận hành lò.

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, cơ chế phản ứng bắt nơtron nhiệt $(n,\gamma)$ và quy luật phân rã phóng xạ tự nhiên. Quá trình tạo thành $^{125}\text{I}$ diễn ra qua phản ứng hạt nhân cốt lõi:

$$^{124}\text{Xe}(n,\gamma)^{125}\text{Xe} \xrightarrow{EC, \beta^+, T_{1/2}=18\text{ h}} {}^{125}\text{I} \xrightarrow{EC, T_{1/2}=59,408\text{ d}} {}^{125}\text{Te (bền)}$$

Ba khái niệm then chốt chi phối toàn bộ quá trình tính toán bao gồm:

  • Tiết diện phản ứng vi mô ($\sigma$): Tiết diện bắt nơtron nhiệt của $^{124}\text{Xe}$ đạt $165 \pm 20\text{ barn}$, kèm theo tích phân cộng hưởng $I_0 = 3600\text{ barn}$.
  • Thông lượng nơtron ($\Phi$): Mật độ dòng hạt nơtron nhiệt tới bia trong trường lò phản ứng, dao động trong khoảng khảo sát từ $5 \cdot 10^{12}$ đến $10^{14}\text{ n/cm}^2\text{/s}$.
  • Hiện tượng cháy nơtron (Burn-up) và phản ứng cạnh tranh: Sự hấp thụ nơtron của chính sản phẩm sinh ra qua phản ứng $^{125}\text{I}(n,\gamma)^{126}\text{I}$ (tiết diện bắt $\sigma = 894\text{ barn}$, tích phân cộng hưởng $I_0 = 13730\text{ barn}$), tạo nên tạp chất không mong muốn là $^{126}\text{I}$ ($T_{1/2} = 13,0\text{ ngày}$).

Bên cạnh đó, khung lý thuyết phân rã bắt điện tử (Electron Capture - EC) giải thích đặc tính phát xạ năng lượng của $^{125}\text{I}$ gồm 6,68% tia gamma (E = 0,03549 MeV) và 93,32% tia X đặc trưng (E = 0,028 MeV), tạo nên hiệu quả sinh học cục bộ mà không gây tổn thương mô lành xung quanh.

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

Nghiên cứu thiết lập hệ 9 phương trình vi phân thường phi tuyến bậc nhất mô tả đồng thời tốc độ biến đổi và phân rã của các đồng vị trong chuỗi: $^{124}\text{Xe}$, $^{125}\text{Xe}$, $^{126}\text{Xe}$, $^{127}\text{Xe}$, $^{125}\text{I}$, $^{126}\text{I}$, $^{127}\text{I}$ và $^{128}\text{I}$.

  • Cỡ mẫu và thông số đầu vào: Khối lượng bia khảo sát gồm 15 đến 50 gram đối với khí Xenon tự nhiên (chứa 0,095% $^{124}\text{Xe}$) và 0,4 đến 0,5 gram đối với khí Xenon giàu ($^{124}\text{Xe} > 99%$).
  • Phương pháp phân tích số trị: Giải thuật tích phân số Runge-Kutta bậc 4 (RK4) được lựa chọn nhờ ưu thế vượt trội về độ ổn định, khả năng giải quyết các hệ phương trình vi phân cứng (stiff systems) đặc thù trong vật lý lò phản ứng và sai số cắt cụt cục bộ chỉ ở mức $O(h^5)$.
  • Công cụ thực thi và chuẩn hóa: Toàn bộ thuật toán được lập trình hướng đối tượng trên nền tảng Visual Basic 6.0 chạy trong môi trường Windows 7. Dữ liệu tính toán được đối chiếu, chuẩn hóa với các bộ số liệu thực nghiệm công bố bởi Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) theo các mốc thời gian chiếu 10 đến 300 giờ và thời gian làm nguội 20 đến 45 ngày.

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 thuật toán cho thấy độ tương thích vượt trội khi so sánh với số liệu chuẩn của IAEA:

  • Đối với bia 15 gram Xenon tự nhiên chiếu xạ ở thông lượng $\Phi = 5 \cdot 10^{13}\text{ n/cm}^2\text{/s}$ trong 200 giờ và phân rã 40 ngày, hoạt độ $^{125}\text{I}$ tính toán đạt 29,074 GBq (IAEA công bố 29,15 GBq, sai số dưới 0,26%), hàm lượng tạp chất $^{126}\text{I}$ đạt 0,996% và hoạt độ riêng đạt 608,70 GBq/mg. Khi tăng thời gian chiếu lên 300 giờ và phân rã 45 ngày, sản lượng đạt 39,81 GBq (so với 40,32 GBq của IAEA).
  • Đối với 0,4 gram bia Xenon giàu ($^{124}\text{Xe} > 99%$), sau 24 giờ chiếu và 20 ngày phân rã, sản lượng đạt 128,74 GBq với hoạt độ riêng 642,12 GBq/mg.

Khảo sát phương pháp sản xuất theo mẻ gián đoạn:

  • Sử dụng 50 gram Xenon tự nhiên: Khi tăng thông lượng $\Phi$ từ $5 \cdot 10^{12}\text{ n/cm}^2\text{/s}$ lên $10^{14}\text{ n/cm}^2\text{/s}$, hoạt độ $^{125}\text{I}$ tăng từ 1 Ci lên trên 14 Ci. Tuy nhiên, hàm lượng tạp chất $^{126}\text{I}$ tăng mạnh từ 1% lên mức 20% do hiệu ứng cháy nơtron của $^{125}\text{I}$.
  • Sử dụng 0,5 gram Xenon giàu 99,9%: Hoạt độ $^{125}\text{I}$ tăng tuyến tính theo độ giàu, đạt 26,11 Ci ở thời gian chiếu 17 giờ, với tỷ lệ tạp chất $^{126}\text{I}/^{125}\text{I}$ ở mức an toàn 0,82%. Nếu kéo dài thời gian chiếu lên 85 giờ, tỷ lệ tạp chất tăng vọt lên 10,19%.

Khảo sát chu trình chiếu xạ vòng lặp tuần hoàn (17 giờ chiếu, 113 giờ phân rã mỗi vòng):

  • Bia 15 gram Xenon tự nhiên: Sau 10 vòng lặp (tổng thời gian 47 ngày), tổng hoạt độ $^{125}\text{I}$ tích lũy đạt 5809,23 mCi (5,81 Ci), hàm lượng tạp chất $^{126}\text{I}$ giảm xuống chỉ còn 0,413% (< 1%), hoạt độ riêng đạt mức rất cao là 16786,08 Ci/g.
  • Bia 0,5 gram Xenon giàu 99,9%: Sau 10 vòng lặp, tổng hoạt độ $^{125}\text{I}$ trong buồng phân rã đạt mức kỷ lục 201507,7 mCi (201,5 Ci), trong khi hoạt độ tạp chất $^{126}\text{I}$ được giới hạn ở 833,32 mCi.

Thảo luận kết quả

Các phát hiện từ mô hình tính toán được minh họa rõ nét thông qua hệ thống đồ thị động học phân rã và bảng ma trận tối ưu hóa thời gian. Đồ thị tích lũy ròng chỉ ra rằng trong phương pháp mẻ gián đoạn, sự gia tăng thông lượng nơtron mang lại hiệu ứng kép: vừa tăng sản lượng $^{125}\text{I}$ vừa đẩy mạnh quá trình chuyển hóa $^{125}\text{I} \to {}^{126}\text{I}$. Tạp chất $^{126}\text{I}$ là đồng vị phát tia gamma cứng không mong muốn, đòi hỏi thời gian làm nguội kéo dài từ 30 đến 40 ngày để phân rã tự nhiên (do $T_{1/2} = 13,0\text{ ngày}$ ngắn hơn $59,4\text{ ngày}$ của $^{125}\text{I}$).

Ngược lại, biểu đồ phân tích chu trình vòng lặp tuần hoàn chứng minh tính ưu việt tuyệt đối về chất lượng hạt nhân. Bằng cách giới hạn thời gian chiếu ở mức 17 giờ (tương đương khoảng 90% chu kỳ bán rã 18 giờ của $^{125}\text{Xe}$), lượng khí $^{125}\text{Xe}$ tạo thành được nhanh chóng chuyển ra buồng làm lạnh độc lập để phân rã thành $^{125}\text{I}$. Nhờ việc tách sản phẩm ra khỏi trường nơtron, hiện tượng cháy hạt nhân bị triệt tiêu, giữ cho tỷ lệ tạp chất luôn dưới ngưỡng 0,5% ngay sau khi xử lý hóa học, loại bỏ hoàn toàn thời gian chờ đợi làm nguội kéo dài.

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

  1. Thiết kế và lắp đặt hệ thống chiếu xạ vòng lặp tuần hoàn tại Lò Đà Lạt: Viện Nghiên cứu Hạt nhân Đà Lạt cần chủ trì chế tạo mô hình vòng lặp khép kín kết nối kênh chiếu ướt ($\Phi = 9 \cdot 10^{12}\text{ n/cm}^2\text{/s}$) với buồng phân rã ngoài vùng hoạt. Mục tiêu đạt sản lượng từ 50 đến 100 Ci mỗi năm, hoàn thành lắp đặt thử nghiệm trong vòng 18 tháng.
  2. Ưu tiên đầu tư nguyên liệu bia khí Xenon giàu ($^{124}\text{Xe} > 99%$): Trung tâm Nghiên cứu và Điều chế Đồng vị phóng xạ cần triển khai công nghệ bẫy thu hồi khí Xenon tự động bằng hệ thống làm lạnh ngưng tụ và bơm nén tuần hoàn. Giải pháp này giúp tái sử dụng trên 95% lượng khí quý, giảm 40% chi phí nguyên liệu thô trong khung thời gian 12 tháng.
  3. Chuẩn hóa quy trình tinh chế hóa phóng xạ và kiểm tra chất lượng: Phòng Hóa phóng xạ cần ban hành quy chuẩn chiết tách $^{125}\text{I}$ bằng phương pháp hấp thụ trong dung dịch NaOH 0,05M, đảm bảo nồng độ tạp chất phóng xạ $^{126}\text{I} < 0,1%$ và hoạt độ riêng đạt trên 16000 Ci/g trước khi chia liều NaI thành phẩm, thực hiện định kỳ cho 100% các mẻ sản xuất.
  4. Hình thành chuỗi liên kết sản xuất và ứng dụng lâm sàng hạt xạ trị: Bộ Y tế cùng Bộ Khoa học và Công nghệ cần thiết lập đề án phối hợp giữa Lò phản ứng Đà Lạt và các bệnh viện hạt nhân tuyến đầu (Bệnh viện K, Bệnh viện Bạch Mai, Bệnh viện Trung ương Quân đội 108) nhằm chuyển giao công nghệ đóng gói hạt phóng xạ $4,5 \times 0,8\text{ mm}$ nội địa, đáp ứng 50% nhu cầu điều trị ung thư tuyến tiền liệt trong giai đoạn 2020 - 2025.

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

  • Kỹ sư vận lý lò và vận hành hạt nhân: Tiếp cận phương pháp tích phân số Runge-Kutta bậc 4 và bộ tham số tiết diện vi mô để tính toán thông lượng, lập kế hoạch nạp mẫu chiếu xạ tối ưu trong vùng hoạt lò phản ứng nghiên cứu.
  • Chuyên gia hóa dược phóng xạ và công nghệ bức xạ: Ứng dụng quy trình công nghệ 5 bước từ chiếu xạ, chờ phân rã 113 giờ, đến chưng cất hấp thụ trong dung dịch kiềm NaOH 0,05M để sản xuất dung dịch NaI thương phẩm.
  • Bác sĩ ung bướu và y học hạt nhân: Nắm vững đặc tính vật lý hạt nhân của $^{125}\text{I}$ (năng lượng phát xạ 35 keV, chu kỳ bán rã 59,4 ngày) phục vụ kỹ thuật cấy hạt phóng xạ điều trị ung thư tiền liệt tuyến và xạ trị áp sát u não.
  • Học viên cao học và nghiên cứu sinh ngành Vật lý hạt nhân: Sử dụng cấu trúc mã nguồn Visual Basic 6.0 và mô hình hệ phương trình vi phân làm tài liệu tham khảo để mở rộng tính toán cho các đồng vị y tế khác như $^{89}\text{Sr}$, $^{153}\text{Sm}$, $^{186}\text{Re}$ và $^{188}\text{Re}$.

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

Tại sao Iodine-125 lại được ưu tiên lựa chọn trong xạ trị cấy hạt ung thư tuyến tiền liệt?

$^{125}\text{I}$ phát tia gamma mềm 35 keV với quãng chạy trong mô chỉ vài milimét, kết hợp chu kỳ bán rã 59,408 ngày cho phép cung cấp liều hấp thụ tập trung cao để tiêu diệt tế bào ung thư mà không gây tổn thương mô lành lân cận. Hạt phóng xạ kích thước 4,5 x 0,8 mm được cấy vĩnh viễn an toàn, đạt hiệu quả điều trị 100% trên các ca lâm sàng.

Phản ứng hạt nhân cơ bản để sản xuất Iodine-125 trong lò phản ứng diễn ra như thế nào?

Quá trình bắt đầu bằng việc hạt nhân $^{124}\text{Xe}$ bắt nơtron nhiệt với tiết diện 165 barn để tạo ra $^{125}\text{Xe}$. Sau đó, $^{125}\text{Xe}$ phân rã bắt điện tử với chu kỳ bán rã 18 giờ để tạo thành sản phẩm mong muốn là đồng vị $^{125}\text{I}$.

Vì sao cần kiểm soát nghiêm ngặt hàm lượng tạp chất Iodine-126 dưới mức 1%?

$^{126}\text{I}$ phát tia gamma năng lượng cao (lên tới 0,666 MeV và 1,42 MeV), gây tăng liều chiếu xạ nguy hiểm cho các mô lành và nhân viên y tế. Việc duy trì tỷ lệ $^{126}\text{I}/^{125}\text{I} < 1%$ (và tối ưu $< 0,5%$) là tiêu chuẩn bắt buộc của Dược điển và IAEA đối với các chế phẩm dùng trong y tế.

Ưu điểm vượt trội của phương pháp chiếu xạ tuần hoàn so với phương pháp mẻ gián đoạn là gì?

Phương pháp tuần hoàn đưa khí Xenon ra buồng phân rã ngoài lò sau mỗi 17 giờ chiếu, giúp bảo vệ $^{125}\text{I}$ khỏi dòng nơtron, triệt tiêu phản ứng tạo tạp chất $^{126}\text{I}$. Nhờ đó, sản phẩm đạt hoạt độ riêng vượt trội (trên 16700 Ci/g) và có thể sử dụng ngay mà không cần chờ nguội 30 - 40 ngày.

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

Hoàn toàn đủ điều kiện. Hốc bẫy nơtron của Lò Đà Lạt đạt thông lượng $2 \cdot 10^{13}\text{ n/cm}^2\text{/s}$ và kênh chiếu ướt đạt $9 \cdot 10^{12}\text{ n/cm}^2\text{/s}$, hoàn toàn đáp ứng các thông số kỹ thuật tối ưu để sản xuất quy mô từ 5 đến 200 Ci cho mỗi đợt vận hành tuần hoàn.

Kết luận

  • Luận văn đã xây dựng thành công bộ công cụ phần mềm chuyên dụng trên nền tảng Visual Basic 6.0, ứng dụng giải thuật Runge-Kutta bậc 4 để mô phỏng chính xác động học 9 hạt nhân phóng xạ với sai số dưới 0,3% so với chuẩn IAEA.
  • Xác định phương pháp chiếu tuần hoàn (17 giờ chiếu, 113 giờ phân rã) là công nghệ tối ưu, cho phép sản xuất $^{125}\text{I}$ đạt hoạt độ riêng trên 16700 Ci/g và kiểm soát tạp chất $^{126}\text{I}$ ở mức 0,413%.
  • Chứng minh tính khả thi tuyệt đối của việc triển khai sản xuất $^{125}\text{I}$ tại Lò phản ứng hạt nhân Đà Lạt với thông lượng nơtron từ $9 \cdot 10^{12}$ đến $2 \cdot 10^{13}\text{ n/cm}^2\text{/s}$.
  • Thiết lập cơ sở dữ liệu khoa học định lượng giúp tối ưu hóa việc sử dụng bia Xenon giàu ($>99%$) và thiết kế buồng thu hồi khí trơ, giảm thiểu chi phí vận hành và rủi ro bức xạ.
  • Mở ra khả năng mở rộng thuật toán mô hình hóa cho hàng loạt đồng vị y học quan trọng khác như $^{89}\text{Sr}$, $^{153}\text{Sm}$ và $^{188}\text{Re}$.

Quý độc giả, nhà khoa học và các đơn vị y tế quan tâm có thể khai thác trực tiếp cơ sở dữ liệu và mô hình tính toán của luận văn để phục vụ công tác nghiên cứu, đào tạo hoặc ứng dụng triển khai sản xuất thực tế.