Tổng quan nghiên cứu

Nghiên cứu phản ứng hạt nhân ở dải năng lượng cực cao đóng vai trò then chốt trong việc giải mã cấu trúc vi mô của vật chất và cung cấp dữ liệu cơ bản cho ngành năng lượng nguyên tử. Trong khi các tương tác hạt nhân do nơtron và hạt tích điện gây ra đã được khảo sát tương đối chi tiết, các dữ liệu thực nghiệm về phản ứng quang hạt nhân trong vùng năng lượng gigaelectronvolt (GeV) hiện vẫn còn rất khan hiếm trên toàn cầu. Nhằm giải quyết khoảng trống học thuật này, luận văn thạc sĩ chuyên ngành Vật lý nguyên tử của tác giả Đồng Thị Kim Cúc, dưới sự hướng dẫn khoa học của Tiến sĩ Phạm Đức Khuê 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 phân tích thực nghiệm chuyên sâu về tương tác quang hạt nhân trên vật liệu Zirconi (Zr).

Mục tiêu cốt lõi của công trình là nhận diện toàn diện các sản phẩm phóng xạ tạo thành, xác định chính xác suất lượng phản ứng (reaction yield) và khảo sát quy luật phân bố suất lượng theo số khối hạt nhân. Đối tượng nghiên cứu là Zirconi tự nhiên – một kim loại có độ bền cơ học cao, khả năng chống ăn mòn hóa học xuất sắc và là vật liệu kết cấu vỏ bọc thanh nhiên liệu tối quan trọng trong các lò phản ứng hạt nhân. Nghiên cứu được triển khai thực nghiệm tại Trung tâm Gia tốc Pohang (POSTECH), Hàn Quốc, sử dụng chùm electron được gia tốc trên đường hầm dài 160 mét để tạo chùm bức xạ hãm có năng lượng cực đại 2,5 GeV với thông lượng tổng cộng đạt $2,19 \times 10^{14}$ electron trong 4 giờ chiếu xạ liên tục. Kết quả nghiên cứu không chỉ làm sáng tỏ cơ chế phản ứng photospallation năng lượng cao mà còn bổ sung hệ số liệu thực nghiệm quý giá, hỗ trợ trực tiếp cho công tác mô phỏng an toàn bức xạ và che chắn máy gia tốc thế hệ mới.

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

Khung lý thuyết áp dụng

Cơ chế tương tác giữa photon và hạt nhân nguyên tử biến đổi mạnh mẽ theo dải năng lượng của chùm bức xạ tới, được chia thành ba vùng lý thuyết cơ bản:

  • Cơ chế cộng hưởng lưỡng cực khổng lồ (GDR): Chiếm ưu thế ở dải năng lượng thấp từ 8 đến 40 MeV. Tại đây, bước sóng của photon tương đương kích thước hạt nhân, tạo ra dao động tập thể giữa khối proton và khối nơtron theo mô hình hai chất lỏng của Goldhaber-Teller và mô hình chồng chập trạng thái vi mô của Wilkinson.
  • Cơ chế giả đơtron (QDM): Giữ vai trò chủ đạo ở dải năng lượng trung bình từ 30 đến 140 MeV. Theo mô hình Levinger, bước sóng photon giảm ngắn lại và tương tác trực tiếp với các cụm cặp nơtron - proton (giả đơtron) bên trong hạt nhân thay vì toàn bộ hạt nhân hợp phần. Tiết diện phản ứng tỷ lệ thuận với số cặp n-p và thông số Levinger đặc trưng.
  • Cơ chế thác lũ nội hạt nhân và mô hình Rudstam: Ở vùng năng lượng cao trên 140 MeV (vượt ngưỡng sinh hạt pion), tương tác diễn ra theo mô hình hai giai đoạn của Serber. Giai đoạn đầu, photon va chạm trực tiếp với từng nucleon riêng lẻ tạo nên dòng thác lũ nội hạt nhân (intranuclear cascade). Giai đoạn tiếp theo, hạt nhân dư thừa năng lượng sẽ khử kích thích chủ yếu thông qua quá trình bay hơi (evaporation) các nucleon, tạo ra phản ứng quang sinh nhiều hạt (photospallation). Quy luật phân bố suất lượng sản phẩm theo điện tích Z và số khối A được mô tả chính xác qua công thức bán thực nghiệm 5 tham số của Rudstam.

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

Nghiên cứu sử dụng phương pháp kích hoạt hạt nhân (Nuclear Activation Analysis) kết hợp phép đo phổ gamma phân giải cao:

  • Mẫu nghiên cứu: Tấm kim loại Zirconi tự nhiên (natZr) có độ tinh khiết cao trên 99,8%, kích thước chuẩn hóa, bao gồm 5 đồng vị bền tự nhiên: 90Zr (51,45%), 91Zr (11,22%), 92Zr (17,15%), 94Zr (17,38%) và 96Zr (2,80%).
  • Cơ chế tạo chùm bức xạ: Bắn phá chùm electron 2,5 GeV từ máy gia tốc tuyến tính vào bia hãm Vonfram (W) dày 2 mm để sinh bức xạ hãm có phổ năng lượng liên tục với năng lượng đỉnh 2,5 GeV, độ rộng xung 1 nanogiây và tần số lặp 10 Hz.
  • Hệ thiết bị ghi nhận phổ: Sử dụng hai hệ phổ kế gamma đa kênh hiện đại của Trung tâm Gia tốc Pohang, gồm detector bán dẫn Gecmani siêu tinh khiết HPGe Canberra (thể tích 100 cm³, độ phân giải 1,8 keV tại đỉnh 1332 keV) và HPGe Ortec (thể tích 100 cm³, độ phân giải 1,9 keV). Khoảng cách đặt mẫu đo được chuẩn hóa ở hai vị trí 0,5 cm và 2,5 cm.
  • Quy trình đo và phân tích: Để tách biệt phổ của các đồng vị có chu kỳ bán rã khác nhau và giảm thiểu sai số hình học, các phép đo được thiết kế theo nhiều chu kỳ: thời gian làm nguội (phơi mẫu) biến thiên từ 180 phút đến 905 giờ, thời gian thu nhận phổ kéo dài từ 15 phút đến 48 giờ. Dữ liệu phổ được xử lý bằng các phần mềm chuyên dụng như Genie 2000, GammaVision và FitzPeaks Gamma Analysis kết hợp tra cứu thư viện số liệu hạt nhân quốc tế.

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 thực nghiệm phổ gamma đã ghi nhận nhiều kết quả khoa học nổi bật:

  • Độ bao phủ sản phẩm rộng lớn: Nhận diện chính xác 46 đồng vị phóng xạ khác nhau sinh ra từ các phản ứng quang hạt nhân trên bia natZr, với số khối trải dài từ A = 54 (đồng vị 54Mn) đến A = 95 (đồng vị 95Zr và 95Nb).
  • Phổ chu kỳ bán rã đa dạng: Đồng vị phóng xạ có thời gian sống ngắn nhất đo được là 79Rb với chu kỳ bán rã 22,9 phút, trong khi đồng vị có thời gian sống dài nhất là 54Mn với chu kỳ bán rã lên tới 312,3 ngày.
  • Kênh phản ứng photospallation cực mạnh: Phản ứng quang sinh nhiều hạt chiếm ưu thế tuyệt đối. Điển hình là phản ứng trên đồng vị 96Zr tạo thành hạt nhân 54Mn đã phát xạ đồng thời tới 42 nucleon (bao gồm 27 nơtron và 15 proton), với mức năng lượng ngưỡng tính toán lên đến 368,45 MeV.
  • Nhận diện các trạng thái đồng phân (Isomer) và quang sinh pion: Ghi nhận 10 sản phẩm ở trạng thái đồng phân kích thích bao gồm 90mY, 87mY, 86mY, 85mY, 84mY, 87mSr, 85mSr, 84mRb, 82mRb và 74mBr. Đặc biệt, luận văn đã phát hiện hai đồng vị 95Nb (ngưỡng 146,46 MeV) và 90Nb (ngưỡng từ 145,81 MeV đến 191,03 MeV) được sinh ra từ phản ứng quang sinh hạt pion mang điện tích âm $(\gamma, xn\pi^-)$.

Thảo luận kết quả

Dữ liệu thực nghiệm thu được phản ánh tính quy luật rõ rệt của phản ứng hạt nhân ở dải năng lượng GeV. Khi biểu diễn phân bố suất lượng theo số khối sản phẩm A trên đồ thị tọa độ bán logarit, đường cong thực nghiệm tạo thành dạng hình chuông đối xứng đặc trưng của cơ chế bay hơi nucleon sau quá trình thác lũ Serber. Dữ liệu này hoàn toàn tương thích và khớp nối chặt chẽ với đường cong lý thuyết xây dựng từ mô hình bán thực nghiệm Rudstam 5 tham số.

Trong quá trình phân tích số liệu, tác giả đã áp dụng các thuật toán hiệu chỉnh chuyên sâu để triệt tiêu sai số hệ thống:

  • Hiệu chỉnh hiệu ứng cộng đỉnh (Summing Effect): Được thực hiện nghiêm ngặt cho các vạch gamma nối tầng phát ra đồng thời, tiêu biểu như cặp vạch 898,042 keV (xác suất phát 93,7%) và 1836,063 keV (xác suất phát 99,2%) của 88Y, cũng như các vạch của 87Zr và 91Rb.
  • Xử lý can nhiễu đồng vị: Giải quyết triệt để sự trùng lặp vạch năng lượng giữa các đồng vị khác nhau, ví dụ như mức năng lượng 388,531 keV xuất hiện đồng thời ở cả 87Y và 87mSr, bằng cách phối hợp phân tích phân rã theo hàm số mũ thời gian và lựa chọn các vạch phổ đơn năng độc lập.

Đề 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, luận văn đưa ra 4 nhóm khuyến nghị có tính ứng dụng thực tiễn cao:

  • Tích hợp dữ liệu vào ngân hàng hạt nhân quốc tế: Viện Năng lượng Nguyên tử Việt Nam cùng các nhóm nghiên cứu cần chuẩn hóa và đăng ký tập số liệu suất lượng của 46 đồng vị phân tích từ bia Zr vào hệ thống cơ sở dữ liệu EXFOR của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) trong vòng 12 tháng tới.
  • Tối ưu hóa các chương trình mô phỏng bức xạ: Các kỹ sư thiết kế cần đưa các thông số bán thực nghiệm Rudstam đã tinh chỉnh từ dữ liệu bia Zr vào các bộ mã mô phỏng Monte Carlo (như Geant4, FLUKA, PHITS), nhằm nâng cao độ chính xác dự báo liều chiếu và an toàn bức xạ lên hơn 15% trong giai đoạn 2026–2028.
  • Nghiên cứu ứng dụng sản xuất đồng vị phóng xạ y học: Bộ Khoa học và Công nghệ cùng các viện nghiên cứu chuyên ngành cần xem xét khả năng tận dụng chùm bức xạ hãm từ máy gia tốc để sản xuất đồng vị 89Zr (chu kỳ bán rã 78,41 giờ) có độ tinh khiết cao phục vụ kỹ thuật ghi hình cắt lớp PET trong chẩn đoán ung thư trong lộ trình 2 năm tới.
  • Mở rộng dải thực nghiệm sang các bia vật liệu lân cận: Khuyến nghị các nhóm nghiên cứu vật lý hạt nhân tiếp tục thực hiện các thí nghiệm đo phổ kích hoạt trên các bia kim loại chuyển tiếp khác như Molypden (Mo) và Niobi (Nb) ở dải năng lượng từ 1,0 GeV đến 3,0 GeV trước năm 2030 nhằm hoàn thiện bức tranh toàn diện về tương tác photospallation.

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

Công trình học thuật này mang lại giá trị lý thuyết và ứng dụng rõ nét cho 4 nhóm đối tượng chính:

  • Các nhà nghiên cứu vật lý hạt nhân cơ bản: Khai thác tập dữ liệu 46 phản ứng quang hạt nhân với năng lượng ngưỡng từ 7,86 MeV đến 368,45 MeV để kiểm chứng và hiệu chỉnh các mô hình cấu trúc hạt nhân lý thuyết.
  • Kỹ sư an toàn bức xạ và công nghệ máy gia tốc: Ứng dụng số liệu suất lượng thực nghiệm vào việc tính toán chiều dày che chắn chì và bê tông, đảm bảo an toàn vận hành cho các máy gia tốc hạt năng lượng cao từ 1 GeV trở lên.
  • Chuyên gia vật liệu và công nghệ lò phản ứng: Đánh giá mức độ kích hoạt phóng xạ và độ bền bức xạ dài hạn của vỏ bọc hợp kim Zirconi khi làm việc trong môi trường bức xạ photon cường độ lớn.
  • Học viên cao học và sinh viên chuyên ngành Vật lý: Tiếp cận phương pháp luận mẫu mực về kỹ thuật đo phổ gamma đa kênh bằng detector bán dẫn HPGe có độ phân giải 1,8 keV, cùng quy trình xử lý phổ phức tạp qua các công cụ Genie 2000 và FitzPeaks.

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

1. Tại sao Zirconi (Zr) tự nhiên lại được lựa chọn làm bia nghiên cứu trong phản ứng quang hạt nhân?
Zirconi là vật liệu kết cấu then chốt trong ngành công nghiệp hạt nhân nhờ tiết diện bắt nơtron nhiệt nhỏ, độ bền cơ học cao và tính chống ăn mòn vượt trội. Việc nghiên cứu phản ứng trên 5 đồng vị bền của natZr giúp xác định chính xác mức độ kích hoạt phóng xạ của vật liệu vỏ thanh nhiên liệu dưới chùm bức xạ hãm 2,5 GeV.

2. Máy gia tốc tuyến tính 2,5 GeV tại Pohang đóng vai trò gì trong sự thành công của nghiên cứu?
Máy gia tốc Pohang dài 160 mét cung cấp chùm electron 2,5 GeV với dòng xung lớn hơn 2 A và tần số 10 Hz. Chùm electron này bắn phá bia Vonfram 2 mm tạo ra dòng photon hãm thông lượng cao $2,19 \times 10^{14}$ electron, cung cấp đủ năng lượng kích hoạt các phản ứng ngưỡng cao trên 300 MeV.

3. Làm thế nào để phân biệt chính xác 46 đồng vị phóng xạ trên cùng một mẫu đo?
Tác giả đã phối hợp đo phổ gamma ở nhiều thời điểm làm nguội khác nhau (từ 180 phút đến 905 giờ) trên detector bán dẫn HPGe có độ phân giải 1,8 keV. Dựa vào năng lượng vạch phổ đặc trưng và chu kỳ bán rã (từ 22,9 phút của 79Rb đến 312,3 ngày của 54Mn), phần mềm Genie 2000 và FitzPeaks đã định danh chính xác từng đồng vị.

4. Phản ứng photospallation ở năng lượng 2,5 GeV có đặc điểm gì khác so với năng lượng thấp?
Ở dải năng lượng thấp dưới 40 MeV, phản ứng chủ yếu phát xạ 1 đến 2 hạt theo cơ chế cộng hưởng khổng lồ GDR. Ngược lại, ở mức 2,5 GeV, phản ứng photospallation diễn ra theo cơ chế thác lũ Serber và bay hơi, giải phóng đồng thời tới 42 nucleon như phản ứng tạo ra 54Mn từ bia 96Zr.

5. Hiệu ứng cộng đỉnh trong phép đo phổ gamma được khắc phục như thế nào?
Hiệu ứng cộng đỉnh xuất hiện khi hai tia gamma nối tầng đến detector đồng thời (như cặp vạch 898,042 keV và 1836,063 keV của 88Y). Nhóm nghiên cứu đã xây dựng công thức hiệu chỉnh giải tích dựa trên hiệu suất ghi toàn phần và chuẩn hóa khoảng cách đo 2,5 cm để hạn chế sai số xuống dưới mức 5%.

Kết luận

Luận văn thạc sĩ của tác giả Đồng Thị Kim Cúc đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với những đóng góp khoa học nổi bật:

  • Nhận diện thành công 46 đồng vị phóng xạ tạo thành trên bia natZr dưới tác dụng của chùm bức xạ hãm năng lượng cực đại 2,5 GeV.
  • Khảo sát chi tiết cơ chế phản ứng photospallation năng lượng cao với kênh phản ứng phát xạ tối đa 27 nơtron và 15 proton ở năng lượng ngưỡng 368,45 MeV.
  • Đo đạc và ghi nhận chính xác 2 kênh phản ứng quang sinh hạt pion $(\gamma, xn\pi^-)$ tạo thành đồng vị 95Nb và 90Nb.
  • Chứng minh sự tương thích cao giữa số liệu thực nghiệm và hàm phân bố suất lượng bán thực nghiệm Rudstam 5 tham số trên dải khối lượng A = 54 đến 95.
  • Hoàn thiện quy trình thực nghiệm đo phổ kích hoạt gamma sử dụng detector bán dẫn HPGe phân giải cao 1,8 keV kết hợp các thuật toán hiệu chỉnh sai số tiên tiến.

Đây là tài liệu học thuật giá trị cao dành cho các nhà nghiên cứu, kỹ sư và học viên ngành Vật lý nguyên tử. Bạn đọc quan tâm có thể tham khảo toàn văn luận văn tại thư viện Đại học Quốc gia Hà Nội để khai thác trọn vẹn tập số liệu thực nghiệm quý giá này.