Tổng quan nghiên cứu
Nghiên cứu phản ứng quang hạt nhân trong vùng năng lượng cao đóng vai trò then chốt trong việc khám phá cấu trúc hạt nhân và cơ chế tương tác vi mô của vật chất. Trong các hệ thống máy gia tốc hiện đại và công nghệ lò phản ứng tiên tiến, vật liệu Molipđen (Mo) tự nhiên thường xuyên chịu sự bắn phá của các chùm bức xạ năng lượng cao, dẫn đến các biến đổi cấu trúc hạt nhân phức tạp. Tuy nhiên, cơ sở dữ liệu thực nghiệm về phản ứng quang hạt nhân do photon hãm gây ra trong dải năng lượng hàng gigaelectronvolt (GeV) trên thế giới vẫn còn tương đối hạn chế, chủ yếu tập trung vào chùm hạt proton hoặc neutron.
Luận văn thạc sĩ khoa học chuyên ngành Vật lý nguyên tử của tác giả Phạm Công Thuyên, được thực hiện dưới sự hướng dẫn khoa học của Tiến sĩ Kim Tiến Thành tại Trung tâm Vật lý hạt nhân thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam cùng Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, đã giải quyết khoảng trống học thuật này. Nghiên cứu tập trung xác định chính xác suất lượng và quy luật phân bố suất lượng của các phản ứng quang hạt nhân trên bia Mo tự nhiên dưới tác động của chùm photon hãm có năng lượng cực đại 2,5 GeV.
Thí nghiệm được triển khai thực địa tại Trung tâm Gia tốc Pohang (Hàn Quốc) trong khuôn khổ đề tài nghiên cứu cơ bản do Quỹ Phát triển Khoa học và Công nghệ Quốc gia (NAFOSTED, mã số 103.21) tài trợ. Kết quả nghiên cứu cung cấp một tập số liệu hạt nhân có độ tin cậy cao, hỗ trợ đắc lực cho việc kiểm chuẩn các mô hình lý thuyết vi mô, tối ưu hóa các quy trình tính toán an toàn bức xạ và định lượng độ hao mòn vật liệu trong các công trình kỹ thuật hạt nhân đương đại.
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 vận dụng hệ thống các mô hình lý thuyết hạt nhân kinh điển để phân tích toàn diện cơ chế phản ứng theo từng thang năng lượng:
- Mô hình hạt nhân hợp phần và cộng hưởng lưỡng cực khổng lồ (GDR): Ở vùng năng lượng photon thấp dưới 30 MeV, quá trình hấp thụ bức xạ điện từ diễn ra chủ yếu thông qua tương tác lưỡng cực điện, hình thành trạng thái hợp phần kích thích. Năng lượng kích thích cộng hưởng được mô tả qua sự kết hợp của hai mô hình Goldhaber-Teller và Steinwedel-Jensen theo biểu thức thực nghiệm W = 31,2 A^(-1/3) + 20,6 A^(-1/6) MeV. Tiết diện kích hoạt trong vùng này phản ánh xác suất phát ra một vài nucleon riêng lẻ.
- Cơ chế giả đơtron Levinger: Trong khoảng năng lượng từ 40 MeV đến 150 MeV, bước sóng photon ngắn hơn kích thước hạt nhân, photon có xu hướng tương tác trực tiếp với các cặp neutron - proton liên kết cục bộ (cụm giả đơtron) thay vì toàn bộ thể tích hạt nhân, với thông số Levinger L dao động từ 6,8 đến 10,3.
- Mô hình thác lũ nội hạt nhân Serber và phản ứng photospallation: Tại vùng năng lượng GeV, cơ chế phản ứng sinh nhiều hạt chiếm ưu thế. Theo mô hình của Serber, quá trình gồm hai giai đoạn nối tiếp: giai đoạn thác lũ nội hạt nhân diễn ra cực nhanh (khoảng 10^(-22) giây) phát ra các hạt sơ cấp, tiếp theo là giai đoạn bay hơi thống kê của hạt nhân kích thích dư cho đến khi giải phóng hết năng lượng bằng tia gamma.
- Công thức bán thực nghiệm 5 thông số Rudstam: Phân bố suất lượng các mảnh sản phẩm được mô hình hóa qua hàm phân bố điện tích và phân bố khối lượng, liên kết chặt chẽ với điện tích khả dĩ nhất và độ dốc theo số khối A.
Phương pháp nghiên cứu
Phương pháp thực nghiệm chủ đạo là kích hoạt bức xạ photon không phá hủy kết hợp với kỹ thuật đo phổ gamma độ phân giải cao:
- Mẫu nghiên cứu và kỹ thuật chọn mẫu: Mẫu Mo tự nhiên gồm 7 đồng vị bền được chế tạo bởi Reactor Experiments Inc dưới dạng đĩa tròn chuẩn hóa có đường kính 0,5 inch, độ dày 0,00394 inch và độ tinh khiết đạt 99,95%. Hai mẫu thí nghiệm độc lập có khối lượng đo lường chính xác là 0,1175 gam (ký hiệu Mo2) và 0,1191 gam (ký hiệu Mo3).
- Hệ thống gia tốc và chiếu xạ: Chùm electron năng lượng 2,5 GeV với độ rộng xung 1 ns, tần số 10 Hz từ máy gia tốc tuyến tính dài 160 mét tại Trung tâm Gia tốc Pohang được bắn phá vào bia chuyển đổi Vonfram (W) dày 0,2 mm để tạo ra phổ bức xạ hãm liên tục có thông lượng chùm hạt tích lũy lên tới 21,9 x 10^13 electron trong thời gian chiếu 240 phút.
- Hệ phổ kế gamma và chuẩn hóa: Sử dụng detector bán dẫn Gecmani siêu tinh khiết (HPGe) thể tích 100 cm³ có độ phân giải năng lượng 1,8 keV tại đỉnh 1332,2 keV của Co-60. Hệ đo được chuẩn hiệu suất ghi tuyệt đối bằng hai bộ nguồn chuẩn QCY và QCYK gồm 12 đồng vị bức xạ tại 5 khoảng cách hình học đo từ 0,5 cm đến 10 cm.
- Hệ thống hiệu chỉnh sai số đo: Nhằm đảm bảo tính chính xác tuyệt đối của suất lượng phản ứng, toàn bộ số liệu phổ được xử lý qua các thuật toán làm khớp đỉnh Gauss trên phần mềm chuyên dụng Gamma Vision và Genie 2000. Đồng thời, nghiên cứu đã tiến hành 5 phép hiệu chỉnh định lượng bắt buộc: hiệu chỉnh thời gian chết và chồng chập xung bằng máy phát xung chuẩn, hiệu chỉnh độ tự hấp thụ tia gamma trong bề dày mẫu kim loại, hiệu chỉnh hiệu ứng cộng đỉnh đối với các vạch gamma nối tầng, và hiệu chỉnh can nhiễu phóng xạ do sự phân rã của các đồng vị mẹ sang đồng vị con theo chuỗi phân rã phóng xạ.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình phân tích phổ gamma kích hoạt ở nhiều mốc thời gian làm nguội khác nhau (từ 45 phút, 1,86 ngày đến 35 ngày) và thời gian đo kéo dài từ 10 phút đến 2400 phút đã định danh và xác định chính xác suất lượng của 68 đồng vị phóng xạ cùng các trạng thái đồng phân hạt nhân (isomer):
- Kênh phản ứng phát một neutron đơn lẻ chiếm suất lượng cao nhất: Đồng vị Mo-99 tạo thành từ phản ứng quang hạt nhân trên Mo-100 có ngưỡng năng lượng 8,29 MeV đạt giá trị suất lượng lớn nhất là 143,4 x 10^(-10) phản ứng/giây, nhờ đóng góp áp đảo từ đỉnh cộng hưởng lưỡng cực khổng lồ nằm ở vùng năng lượng thấp của phổ bức xạ hãm.
- Các sản phẩm tách ít nucleon có hiệu suất hình thành vượt trội: Đồng vị đồng phân Nb-91m (chu kỳ bán rã 60,86 ngày) đạt suất lượng 13,815 x 10^(-10) phản ứng/giây; Zr-89 (chu kỳ bán rã 78,41 giờ) đạt 9,3277 x 10^(-10) phản ứng/giây; Nb-90 (chu kỳ bán rã 14,60 giờ) đạt 7,826 x 10^(-10) phản ứng/giây; và Sr-85 (chu kỳ bán rã 64,84 ngày) đạt 7,634 x 10^(-10) phản ứng/giây.
- Quy luật suy giảm suất lượng theo độ hụt khối lượng: Các đồng vị sản phẩm có số khối A càng xa hạt nhân bia ban đầu thì suất lượng giảm dần theo quy luật hàm mũ rõ rệt. Cụ thể, đồng vị Sr-80 sinh ra từ phản ứng phát ra 4 proton và từ 8 đến 16 neutron có ngưỡng năng lượng phản ứng lên tới trên 110 MeV chỉ đạt suất lượng 0,163 x 10^(-10) phản ứng/giây, tức thấp hơn gần 880 lần so với suất lượng của kênh tạo thành Mo-99.
Thảo luận kết quả
Dữ liệu thực nghiệm được so sánh trực tiếp với các tính toán mô phỏng theo công thức bán thực nghiệm 5 thông số của Rudstam:
- Tương quan giữa số liệu đo đạc và lý thuyết bán thực nghiệm: Đồ thị phân bố suất lượng theo số khối A trên hệ tọa độ bán logarit biểu diễn một đường cong dạng chuông đối xứng. Đỉnh của đường phân bố điện tích tập trung tại giá trị điện tích khả dĩ nhất, phản ánh trạng thái cân bằng tỷ lệ neutron/proton sau quá trình bay hơi hạt. Các kênh phản ứng sinh nhiều hạt (photospallation) phù hợp rất tốt với dự đoán lý thuyết của mô hình thác lũ nội hạt nhân kết hợp bay hơi thống kê.
- Trình bày dữ liệu qua đồ thị và bảng biểu: Sự tương quan giữa năng lượng ngưỡng phản ứng và suất lượng tạo thành được hệ thống hóa qua bảng ma trận 68 đồng vị. Phổ phân rã suy giảm hoạt độ theo thời gian phơi mẫu được minh họa bằng các biểu đồ đường cong phân rã gamma đa kênh, giúp tách biệt các đỉnh phổ chập đôi phức tạp như trường hợp vạch 231,67 keV của Y-85m và Sr-85m.
- Ý nghĩa vật lý: Kết quả khẳng định ở năng lượng 2,5 GeV, tương tác quang hạt nhân chuyển hóa hoàn toàn sang cơ chế thác lũ sinh hạt phức tạp, nhưng phần đuôi năng lượng thấp của phổ photon hãm vẫn đóng góp một tỷ trọng đáng kể vào các phản ứng cộng hưởng đơn giản.
Đề xuất và khuyến nghị
Dựa trên các phát hiện thực nghiệm và cơ sở dữ liệu thu được, nghiên cứu đưa ra 4 khuyến nghị hành động cụ thể:
- Cập nhật dữ liệu hạt nhân vào các thư viện quốc tế: Kiến nghị Cục Năng lượng nguyên tử và Viện Vật lý tiến hành tích hợp toàn bộ 68 số liệu suất lượng thực nghiệm đo đạc được trên bia Mo tự nhiên vào hệ cơ sở dữ liệu hạt nhân quốc tế (EXFOR/ENDF) trong vòng 12 tháng tới, nhằm cung cấp dữ liệu đối chuẩn chính xác cho cộng đồng nghiên cứu vật lý năng lượng cao toàn cầu.
- Tối ưu hóa thiết kế hệ thống che chắn an toàn bức xạ cho máy gia tốc: Đề xuất các trung tâm nghiên cứu bức xạ và cơ sở vận hành máy gia tốc electron áp dụng các hệ số phát xạ thực nghiệm này để hiệu chỉnh các mô hình mô phỏng Monte Carlo (như MCNP, FLUKA, GEANT4). Mục tiêu là giảm sai số tính toán liều bức xạ phát tán xuống dưới 5%, hoàn thành trong giai đoạn 2025-2026.
- Thúc đẩy lộ trình sản xuất đồng vị phóng xạ y tế Mo-99/Tc-99m không dùng Urani: Khuyến nghị các viện nghiên cứu hạt nhân phối hợp cùng các cơ sở y tế triển khai dự án thử nghiệm sản xuất đồng vị y học hạt nhân Mo-99 thông qua phản ứng quang hạt nhân trên máy gia tốc electron công suất lớn, thay thế dần phương pháp phân hạch Urani truyền thống nhằm nâng cao tính an toàn và giảm thiểu chất thải phóng xạ nguy hại.
- Chuẩn hóa quy trình xử lý can nhiễu và hiệu ứng cộng đỉnh: Đề xuất ban hành tài liệu hướng dẫn kỹ thuật chuẩn quốc gia về quy trình phân tích phổ gamma kích hoạt năng lượng cao, tập trung vào phương pháp hiệu chỉnh chuỗi phân rã đồng phân và hiệu ứng chồng chập đỉnh cho các phòng thí nghiệm đo lường bức xạ trước năm 2026.
Đối tượng nên tham khảo luận văn
Tài liệu luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:
- Kỹ sư và chuyên gia an toàn bức xạ tại các trung tâm gia tốc và lò phản ứng: Cung cấp thông số định lượng chính xác về suất lượng phản ứng để tính toán độ hoạt hóa của vật liệu Mo, dự báo tốc độ hao mòn cấu kiện và thiết kế hệ thống tường chắn bảo vệ tối ưu.
- Nhà nghiên cứu vật lý hạt nhân thực nghiệm và lý thuyết: Tài liệu là nguồn số liệu thực chứng quý giá để kiểm chứng độ chính xác của các mô hình vi mô, đặc biệt là cơ chế giả đơtron Levinger và mô hình thác lũ nội hạt nhân Serber ở thang năng lượng GeV.
- Chuyên gia y học hạt nhân và kỹ thuật dược chất phóng xạ: Giúp nắm bắt các thông số vật lý và ngưỡng năng lượng kích hoạt phục vụ định hướng nghiên cứu công nghệ điều chế đồng vị y tế Mo-99 độ tinh khiết cao trên máy gia tốc tuyến tính.
- Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Vật lý: Sử dụng làm tài liệu tham khảo chuẩn mực về phương pháp kích hoạt hạt nhân, kỹ thuật chuẩn phổ kế gamma bán dẫn HPGe và thuật toán xử lý sai số đo lường học thuật phức tạp.
Câu hỏi thường gặp
- Tại sao nghiên cứu lại sử dụng chùm photon hãm năng lượng cực đại 2,5 GeV thay vì các nguồn photon năng lượng thấp?
Photon năng lượng thấp chỉ kích thích các phản ứng đơn giản như phát một hoặc hai nucleon trong vùng cộng hưởng khổng lồ. Việc sử dụng photon hãm 2,5 GeV từ máy gia tốc tuyến tính Pohang cho phép mở rộng nghiên cứu sang vùng năng lượng GeV, kích hoạt các phản ứng sinh nhiều hạt (photospallation) phức tạp và kiểm chứng toàn diện mô hình thác lũ Serber. - Nguyên tố Molipđen (Mo) có vai trò và ý nghĩa thực tiễn như thế nào trong kỹ thuật hạt nhân?
Molipđen là kim loại chịu nhiệt, chống ăn mòn cao, được ứng dụng rộng rãi làm vật liệu cấu trúc trong lò phản ứng và bia máy gia tốc. Ngoài ra, đồng vị Mo-100 trong Mo tự nhiên là tiền chất quan trọng để điều chế đồng vị phóng xạ Mo-99 (phân rã ra Tc-99m) dùng trong hơn 80% các ca chẩn đoán hình ảnh y học hạt nhân. - Phương pháp nào được sử dụng để hiệu chỉnh các can nhiễu phóng xạ do phân rã dây chuyền?
Tác giả đã áp dụng hệ phương trình phân rã phóng xạ liên hoàn kết hợp tỷ số diện tích đỉnh phổ gamma và quy luật hàm mũ theo chu kỳ bán rã của từng cặp đồng vị mẹ - con (ví dụ cặp Nb-89m phân rã sang Zr-89 hoặc Y-87m phân rã sang Y-87), đảm bảo suất lượng tính ra phản ánh đúng phản ứng trực tiếp. - Công thức bán thực nghiệm Rudstam đóng vai trò gì trong việc đánh giá số liệu thực nghiệm?
Công thức Rudstam gồm 5 thông số cho phép mô hình hóa quy luật phân bố điện tích và khối lượng của các mảnh sản phẩm sau phản ứng spallation. Việc làm khớp số liệu thực nghiệm với công thức này giúp đánh giá độ tin cậy của phép đo và ngoại suy tiết diện phản ứng cho các đồng vị khó đo đạc trực tiếp. - Các yếu tố chính nào quyết định độ chính xác của việc xác định suất lượng phản ứng quang hạt nhân?
Độ chính xác phụ thuộc chặt chẽ vào độ phân giải của detector bán dẫn HPGe (1,8 keV), độ chính xác của đường cong chuẩn hiệu suất ghi tuyệt đối từ 12 nguồn chuẩn, việc kiểm soát thời gian chết của hệ đếm và các phép hiệu chỉnh định lượng về tự hấp thụ và cộng đỉnh.
Kết luận
- Luận văn đã xác định thực nghiệm đầy đủ suất lượng và phân bố suất lượng của 68 đồng vị phóng xạ tạo thành trên bia Mo tự nhiên dưới tác xạ của chùm photon hãm năng lượng cực đại 2,5 GeV.
- Khẳng định sự áp đảo của cơ chế phản ứng quang hạt nhân sinh nhiều hạt (photospallation) ở vùng năng lượng cao, đồng thời chỉ rõ sự đóng góp của vùng cộng hưởng GDR đối với kênh tạo Mo-99 có suất lượng cao nhất đạt 143,4 x 10^(-10) phản ứng/giây.
- Kiểm chứng thành công sự phù hợp giữa số liệu thực nghiệm và mô hình bán thực nghiệm Rudstam, làm sáng tỏ cơ chế phát xạ hạt qua mô hình thác lũ nội hạt nhân Serber.
- Đóng góp bộ số liệu hạt nhân thực nghiệm chất lượng cao cho cơ sở dữ liệu khoa học Việt Nam và quốc tế, phục vụ trực tiếp cho tính toán an toàn máy gia tốc và công nghệ đồng vị y tế.
- Trong giai đoạn tiếp theo, các nhà nghiên cứu và kỹ sư hạt nhân được khuyến khích khai thác bộ dữ liệu này để tích hợp vào các thuật toán mô phỏng hạt nhân chuyên sâu và phát triển các ứng dụng y học phóng xạ thực tiễn.