Tổng quan về luận án

Nghiên cứu cấu trúc hạt nhân và cơ chế tương tác ở vùng năng lượng trung bình là một trong những hướng đi then chốt của vật lý hạt nhân hiện đại. Luận án tiến sĩ vật lý chuyên ngành Vật lý nguyên tử và hạt nhân (mã số: 9440106) của Nghiên cứu sinh Nguyễn Thị Xuân với đề tài "Nghiên cứu các đặc trưng của một số phản ứng hạt nhân với bức xạ hãm năng lượng cực đại sau vùng cộng hưởng lưỡng cực khổng lồ và proton năng lượng tới 45 MeV" đã giải quyết những khoảng trống tri thức quan trọng trong cơ sở dữ liệu phản ứng hạt nhân quốc tế.

Về bối cảnh khoa học, phần lớn các công trình nghiên cứu phản ứng quang hạt nhân trên thế giới trước đây chủ yếu tập trung trong vùng cộng hưởng lưỡng cực khổng lồ (Giant Dipole Resonance - GDR) với năng lượng photon dưới 30 MeV, nơi cơ chế hạt nhân hợp phần chiếm ưu thế tuyệt đối. Khi năng lượng kích thích vượt qua vùng GDR (từ 30 MeV đến 70 MeV) và với chùm proton năng lượng tới 45 MeV, các quá trình phản ứng trở nên phức tạp do sự xuất hiện đồng thời và cạnh tranh của nhiều kênh phản ứng đa hạt $(\gamma, xnyp)$, phát xạ hạt tiền cân bằng (pre-equilibrium) và cơ chế tương tác trực tiếp (direct reaction).

Khoảng trống nghiên cứu thực nghiệm (research gap) được chỉ ra rõ nét qua thống kê từ Cơ sở dữ liệu phản ứng hạt nhân thực nghiệm EXFOR (IAEA): "có tới 1350 mục cho phản ứng $(\gamma,n)$, 313 mục cho $(\gamma,2n)$, 64 mục cho $(\gamma,3n)$ và chỉ 14 cho $(\gamma,4n)$; đồng thời có 789 mục cho phản ứng $(\gamma,p)$ nhưng chỉ có 7 mục cho phản ứng $(\gamma,2p)$". Thư viện dữ liệu quang hạt nhân của IAEA cũng chỉ chứa 146 mục tiết diện được đánh giá, thiếu vắng trầm trọng dữ liệu cho các phản ứng sinh nhiều hạt ở năng lượng trung bình.

Để giải quyết khoảng trống này, luận án tập trung vào ba câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. RQ1 & H1: Tiết diện tích phân, suất lượng phản ứng và tỷ số suất lượng ứng phân thay đổi như thế nào theo năng lượng chùm bức xạ hãm từ 50 đến 70 MeV trên các bia nguyên tố tự nhiên? Giả thuyết: Khi năng lượng photon vượt trên đỉnh GDR, cơ chế giả đơtron (Quasi-Deuteron Model - QDM) và tương tác trực tiếp bắt đầu chi phối, làm gia tăng suất lượng các kênh sinh nhiều nucleon và tỷ số tạo thành trạng thái ứng phân spin cao.
  2. RQ2 & H2: Sự phụ thuộc spin của mật độ mức hạt nhân ảnh hưởng ra sao đến tỷ số ứng phân (Isomeric Ratio - IR) của các cặp đồng phân $^{137m,g}\text{Ce}$ và $^{179m,g}\text{W}$? Giả thuyết: Năng lượng photon tới càng cao thì mômen xung lượng truyền cho hạt nhân hợp phần/tiền cân bằng càng lớn, dẫn đến sự ưu tiên chuyển dịch gamma nuôi dưỡng trạng thái ứng phân có giá trị spin cao ($J^\pi$ lớn).
  3. RQ3 & H3: Hàm kích thích của các phản ứng hạt nhân gây bởi proton trong dải 10 - 45 MeV trên các bia $^{nat}\text{Zr}$ và $^{nat}\text{Pd}$ có thể giải thích nhất quán bằng mô hình thống kê kết hợp tiền cân bằng không? Giả thuyết: Mô hình khí Fermi nhiệt độ không đổi (CTFGM) kết hợp cơ chế tiền cân bằng hai thành phần Exciton trong mã TALYS có thể mô tả chính xác thực nghiệm mà không cần hiệu chỉnh tham số tùy tiện.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Lý thuyết hạt nhân hợp phần Bohr (1936), Mô hình giả đơtron Levinger (1951), Mô hình thác lũ bên trong hạt nhân Serber (1947) và Lý thuyết chuyển dịch thống kê Hauser-Feshbach. Phạm vi thực nghiệm bao gồm 06 hạt nhân bia đại diện trải dài theo bảng tuần hoàn: $^{nat}\text{Sr}$ ($Z=38$), $^{nat}\text{Zr}$ ($Z=40$), $^{nat}\text{Pd}$ ($Z=46$), $^{141}\text{Pr}$ ($Z=59$), $^{nat}\text{W}$ ($Z=74$) và $^{197}\text{Au}$ ($Z=79$). Nghiên cứu mang tính đột phá khi bổ sung hệ số liệu thực nghiệm độ chính xác cao vào cơ sở dữ liệu EXFOR/IAEA-NDS, đồng thời định lượng hóa các thông số then chốt ứng dụng trong công nghệ sản xuất đồng vị phóng xạ y tế như $^{82}\text{Sr}/^{82}\text{Rb}$ và $^{100}\text{Pd}/^{100m}\text{Rh}$.

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy các dòng nghiên cứu phản ứng hạt nhân được định hình rõ rệt qua các giai đoạn lịch sử. Khởi đầu từ phản ứng hạt nhân đầu tiên $^{14}\text{N}(\alpha,p)^{17}\text{O}$ của Rutherford (1919) và phát hiện nơtron của Chadwick (1932), các nghiên cứu ban đầu chỉ sử dụng nguồn phóng xạ tự nhiên. Sự phát minh máy gia tốc Cockcroft-Walton và Cyclotron của Lawrence đã mở ra kỷ nguyên nghiên cứu phản ứng hạt nhân năng lượng cao.

Trong lĩnh vực phản ứng quang hạt nhân, mô hình hấp thụ lưỡng cực điện $E1$ trong vùng GDR được đặt nền móng bởi Goldhaber và Teller (1948) với giả thuyết hai chất lỏng proton và neutron dao động ngược pha, sau đó được Steinwedel và Jensen (1950) cùng Wilkinson (1956) phát triển theo hướng chồng chập các dịch chuyển đơn hạt. Tuy nhiên, sự bất đối xứng góc và sự phát xạ hạt năng lượng cao trong vùng sau GDR đòi hỏi sự ra đời của Mô hình giả đơtron (Levinger, 1951) và Mô hình thác lũ nội hạt nhân (Serber, 1947).

Tranh luận học thuật quốc tế hiện nay tập trung vào hai luồng quan điểm đối nghịch về vai trò của phát xạ tiền cân bằng (pre-equilibrium emission) so với phát xạ nhiệt thống kê (evaporative emission) trong dải năng lượng 30 - 70 MeV, cũng như sự lựa chọn tối ưu giữa các mô hình mật độ mức hạt nhân (Nuclear Level Density Models - LDM). Một bên là các mô hình tham số hóa hiện tượng luận như Mô hình khí Fermi nhiệt độ không đổi (Constant Temperature Fermi Gas Model - CTFGM) và Mô hình khí Fermi dịch chuyển ngược (Back-shifted Fermi Gas Model - BSFGM) của Gilbert-Cameron; bên kia là các mô hình vi mô hoàn toàn dựa trên phép tính Hartree-Fock-Bogolyubov (HFB) và bảng tổ hợp Gogny/Skyrme (Goriely, Hilaire).

Khi so sánh với các công trình quốc tế tiêu biểu:

  • Nghiên cứu của Uddin và cộng sự (2008, 2011) tại Viện KIRAMS (Hàn Quốc) trên các bia $^{nat}\text{Mo}$, $^{nat}\text{Zr}$, $^{nat}\text{Pd}$ và của Tárkányi và cộng sự (2005, 2012) tại Viện ATOMKI (Hungary) đã đo tiết diện phản ứng do proton gây ra trên $^{nat}\text{Zr}$, nhưng chủ yếu phân tích trong dải năng lượng thấp dưới 30 MeV.
  • Nhóm tác giả Ditrói và cộng sự (2012) đo tiết diện trên bia $^{nat}\text{Pd}$ và $^{nat}\text{V}$ trong dải năng lượng 37 - 65 MeV nhưng gặp sai số lớn ở các kênh phân rã phức tạp của chuỗi hạt nhân đồng khối $A=100$.
  • Michel và cộng sự (1997) khảo sát diện rộng phản ứng proton từ dải năng lượng thấp đến trung bình trên nhiều bia nguyên tố nhưng chưa phân tích chi tiết sự tạo thành các trạng thái đồng phân siêu bền.

Ở trong nước, nhóm nghiên cứu của GS. Trần Đức Thiệp (Viện Vật lý - Viện Liên hiệp Dubna) đã đo tỷ số ứng phân của một loạt phản ứng $(\gamma,n)$ dưới 25 MeV như $^{74}\text{Se}(\gamma,n)^{73m,g}\text{Se}$, $^{82}\text{Se}(\gamma,n)^{81m,g}\text{Se}$, $^{165}\text{Ho}(\gamma,n)^{164m,g}\text{Ho}$ hay $^{119}\text{Sn}(\gamma,np)^{117m,g}\text{In}$. Nhóm nghiên cứu của GS. Nguyễn Văn Đỗ (Viện Vật lý - PAL POSTECH) đã mở rộng lên dải 40 - 70 MeV và 2.5 GeV cho các bia $^{nat}\text{Sb}$, $^{nat}\text{Pd}$, $^{nat}\text{Ag}$, $^{45}\text{Sc}$, $^{nat}\text{Ti}$.

Luận án của NCS. Nguyễn Thị Xuân đã định vị chính xác khoảng trống tri thức: thực hiện đồng thời phép đo suất lượng, tiết diện tích phân và tỷ số ứng phân trên các hệ hạt nhân bia chọn lọc ($^{nat}\text{Sr}$, $^{nat}\text{Zr}$, $^{nat}\text{Pd}$, $^{141}\text{Pr}$, $^{nat}\text{W}$, $^{197}\text{Au}$) với bức xạ hãm $50 - 70\text{ MeV}$ và proton $10.6 - 43.6\text{ MeV}$, cung cấp bộ số liệu thực nghiệm đối chứng toàn diện cho 6 mô hình mật độ mức và 8 hàm lực gamma trong mã TALYS.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã đóng góp vào việc mở rộng và kiểm chứng các lý thuyết cấu trúc và phản ứng hạt nhân nền tảng:

  • Kiểm chứng mô hình tương tác quang hạt nhân vùng sau GDR: Số liệu thực nghiệm phản ứng $(\gamma, xnyp)$ trên bia $^{nat}\text{Sr}$ và $(\gamma, xn)$ với $x=1-7$ trên bia $^{197}\text{Au}$ chứng minh tính đúng đắn của Mô hình giả đơtron (Levinger) kết hợp mô hình thác lũ nội hạt nhân (Serber). Khi photon có bước sóng de Broglie ngắn hơn kích thước hạt nhân, photon tương tác ưu tiên với cặp neutron-proton (đơtron ảo) thay vì toàn bộ thể tích hạt nhân, giải thích hoàn hảo sự gia tăng suất lượng của các kênh phát xạ đa nucleon ở năng lượng cực đại 60 - 65 MeV.
  • Lý thuyết chuyển dịch spin và mật độ mức: Phân tích tỷ số ứng phân $IR = Y_m / Y_g$ của $^{137m,g}\text{Ce}$ (phản ứng $^{141}\text{Pr}(\gamma,p3n)$) và $^{179m,g}\text{W}$ (phản ứng $^{nat}\text{W}(\gamma,xn)$) chỉ ra rằng: giá trị $IR$ phụ thuộc phi tuyến vào năng lượng kích thích và phân bố spin $J^\pi$ của các mức kích thích trung gian. Sự sai khác mômen góc $|\Delta J| \ge 4\hbar$ giữa trạng thái giả bền (metastable state) và trạng thái cơ bản (ground state) tạo thành "hàng rào cấm chuyển dịch gamma", biến tỷ số $IR$ thành một công cụ nhạy bén để thăm dò thông số phân bố spin (spin cut-off parameter $\sigma_c$) trong các mô hình mật độ mức.
  • Hệ mệnh đề lý thuyết được xác lập:
    1. Mệnh đề 1: Trong vùng năng lượng photon 50 - 70 MeV, kênh phản ứng bứt nhiều hạt $(\gamma, xnyp)$ diễn ra qua hai giai đoạn: tán xạ nội hạt tạo hạt năng lượng cao (tiền cân bằng) và bay hơi thống kê của hạt nhân dư kích thích.
    2. Mệnh đề 2: Tỷ số suất lượng ứng phân tỷ lệ thuận với năng lượng cực đại của phổ bức xạ hãm do sự gia tăng xác suất truyền mômen quỹ đạo $l$ lớn vào hạt nhân hợp phần.
    3. Mệnh đề 3: Mô hình giải tích CTFGM kết hợp hiệu chỉnh dịch chuyển ngược thể hiện độ tương thích cao nhất với số liệu thực nghiệm hạt nhân trung bình và nặng so với các mô hình vi mô phụ thuộc nhiệt độ.

Khung phân tích độc đáo

Khung phân tích của luận án là sự tích hợp đa tầng giữa lý thuyết tương tác, mô phỏng Monte Carlo vận chuyển bức xạ và tính toán vi mô phản ứng hạt nhân:

[Máy gia tốc e- Linac / Cyclotron] 
       │
       ▼
[Mô phỏng Monte Carlo MCNPX / SRIM-2003] ──► Phổ Bremsstrahlung Φ(E) / Suy giảm năng lượng proton Ep(x)
       │
       ▼
[Kích hoạt hạt nhân & Đo phổ Gamma trễ HPGe] ──► Diện tích đỉnh S_γ, Hoạt độ C_pulse, T_1/2, Hiệu suất ε_γ
       │
       ▼
[Xác định thông số: Y, IR, σ_int, σ(E)] 
       │
       ▼
[Đối chuẩn mô hình TALYS (6 LDM x 8 Gamma Strength Functions)] ──► Khẳng định cơ chế vật lý
  • Tính toán trường thông lượng bức xạ hãm $\Phi(E)$: Sử dụng mã MCNPX 2.6 mô phỏng chùm electron năng lượng 50, 55, 60, 65, 70 MeV bắn phá bia chuyển đổi Tungsten (W) kích thước $100 \times 100 \times 0.1\text{ mm}$ đặt cách cửa sổ chùm tia 15 cm.
  • Xác định động học ion trong vật chất: Ứng dụng chương trình SRIM-2003 (Stop and Range Ions in Matter) tính toán chính xác sự suy giảm năng lượng và quãng chạy của chùm proton 45 MeV qua từng lá kim loại trong cấu hình chồng lá (stacked-foil).
  • Phân tích lý thuyết phản ứng hạt nhân: Ứng dụng mã TALYS (phiên bản 1.8, 1.9, 1.95) khảo sát tổ hợp 6 mô hình mật độ mức hạt nhân (1: CTFGM, 2: BSFGM, 3: GSFM, 4: SFGM Goriely, 5: SFHM Hilaire, 6: HFB Goriely) kết hợp 8 hàm lực bức xạ gamma dipole ($E1$ strength functions: Kopecky-Uhl KUGL, Brink-Axel BAL, Goriely Hybrid GH, GTHFB, TRMF...).
  • Điều kiện biên phân tích: Áp dụng hiệu chỉnh tự hấp thụ tia gamma trong mẫu dày, hiệu chỉnh thời gian chết detector ($<5%$), hiệu chỉnh trùng phùng thực (true coincidence summing) và hiệu chỉnh phân rã phóng xạ trong chế độ chiếu xung của Linac thông qua công thức tích phân hoạt độ $C_{pulse}$.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực nghiệm thực chứng (positivism/empirical realism) kết hợp suy diễn logic so sánh lý thuyết. Thiết kế đa tầng (multi-level design) kết hợp giữa hai hệ thống máy gia tốc hạt năng lượng cao chuẩn quốc tế tại Hàn Quốc:

  • Hệ máy gia tốc electron tuyến tính 100 MeV (Electron Linac) tại Trung tâm Gia tốc Pohang (PAL, POSTECH): Dòng chùm electron cực đại 100 mA, độ rộng xung $\tau = 1 - 4\ \mu\text{s}$, tần số lặp xung $f = 10 - 15\text{ Hz}$, bán kính chùm tại monitor trước bia chuyển đổi W là 20 mm. Chiều dài toàn bộ hệ gia tốc là 15 m.
  • Hệ máy gia tốc Cyclotron MC-50 tại Viện Khoa học Phóng xạ và Y học Hàn Quốc (KIRAMS): Nam châm điện 3 phần đạt từ trường trung bình cực đại 1.75 Tesla, gia tốc chùm proton tới năng lượng 45 - 50 MeV, dòng chùm cực đại $60\ \mu\text{A}$ (trong thí nghiệm duy trì ổn định ở mức 100 nA, đường kính chùm 10 mm).

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được chuẩn hóa nghiêm ngặt từ khâu chuẩn bị bia mẫu, kích hoạt, thu nhận phổ gamma đến xử lý số liệu:

[Bia mẫu tinh khiết 99.99%] ──► [Chiếu xạ tối ưu (t_i, t_w, t_m)] ──► [Đo phổ Gamma HPGe ORTEC]
                                                                               │
[Chuẩn năng lượng & Hiệu suất: 241Am, 60Co, 137Cs, 133Ba, 152Eu] ─────────────┤
                                                                               ▼
[Phần mềm: GammaVision 5.10 / Fitzpeaks / Origin] ◄── [Hiệu chỉnh f: Dead time, Self-absorption, Pile-up]
  • Chuẩn bị mẫu nghiên cứu: Sử dụng các lá kim loại siêu tinh khiết 99.99% từ nhà sản xuất Alfa Aesar: $^{nat}\text{Pd}$ (độ dày 0.1 mm), $^{nat}\text{Al}$ (0.05 mm), $^{nat}\text{W}$ (0.05 mm), $^{nat}\text{Zr}$ (0.1 mm), $^{nat}\text{Y}$ (0.025 mm), $^{197}\text{Au}$ lá mỏng, $^{141}\text{Pr}$ kim loại, và mẫu bột $^{nat}\text{Sr(NO}_3)_2$.
  • Kỹ thuật xếp chồng lá (Stacked-foil activation): Trong thí nghiệm với proton, một chồng mẫu gồm 60 lá kim loại được chia làm 10 nhóm (mỗi nhóm gồm 6 lá xếp theo thứ tự W–Zr–Y–Pd–Cu–Al), bọc ngoài bằng lá Al mỏng (0.015 mm) để cố định hình học và tránh thất thoát sản phẩm giật lùi (recoil). Mẫu đồng $^{nat}\text{Cu}$ tinh khiết được xen kẽ làm bia chuẩn giám sát dòng hạt (monitor) thông qua các phản ứng tham chiếu quốc tế $^{nat}\text{Cu}(p,X)^{62}\text{Zn}$ và $^{nat}\text{Cu}(p,X)^{65}\text{Zn}$ với tiết diện chuẩn của IAEA.
  • Hệ phổ kế gamma bán dẫn HPGe: Sử dụng detector Germanium siêu tinh khiết HPGe model GEM-2018-p, No. 39-Tp21360A (ORTEC, USA), đường kính tinh thể 60.5 mm, chiều dài 31 mm, độ phân giải năng lượng đạt $1.8\text{ keV}$ tại đỉnh $1332.5\text{ keV}$ của $^{60}\text{Co}$, hiệu suất ghi tương đối 20% so với detector nhấp nháy $\text{NaI(Tl)}$ $76 \times 76\text{ mm}$. Tiền khuếch đại Model 257P, khuếch đại tuyến tính ORTEC-572, biến đổi tương tự - số đa kênh 8192 kênh, điện áp nuôi cao áp 2 kV.
  • Chuẩn năng lượng và hiệu suất: Hệ phổ kế được chuẩn năng lượng và xác định đường cong hiệu suất ghi đỉnh quang điện $\varepsilon_\gamma(E)$ bằng các nguồn chuẩn đa năng lượng $^{152}\text{Eu}, ^{133}\text{Ba}$ và đơn năng lượng $^{137}\text{Cs}, ^{241}\text{Am}, ^{60}\text{Co}$. Đường cong hiệu suất được khớp chính xác theo hàm đa thức logarithm bậc 5: $$\ln(\varepsilon_\gamma) = \sum_{i=0}^5 a_i \left[\ln\left(\frac{E}{E_0}\right)\right]^i \quad (E_0 = 1\text{ keV})$$
  • Thời gian đo và khử phông: Khoảng cách từ mẫu đến bề mặt detector được thay đổi linh hoạt từ 2 cm đến 20 cm tùy theo hoạt độ phóng xạ nhằm đảm bảo thời gian chết của hệ đo luôn duy trì dưới 5%. Các chu kỳ đo được bố trí tối ưu: từ phép đo nhanh sau chiếu xạ 1 - 2 phút (bắt các đồng vị sống ngắn) đến các phép đo trễ sau vài giờ, vài ngày và kéo dài đến vài tháng (đo các đồng vị sống dài).

Data và phân tích

Số liệu phổ gamma đa kênh được phân tích tự động và kiểm tra chéo qua ba phần mềm chuyên dụng: GammaVision 5.10 (EG&G ORTEC), Fitzpeaks và Origin.

Hoạt độ tích phân ghi nhận trong chế độ chiếu xung được xác định qua biểu thức: $$C_{pulse} = \frac{N_0 \Phi \sigma \varepsilon_\gamma I_\gamma f}{\lambda} \cdot \frac{(1 - e^{-\lambda \tau}) e^{-\lambda(T-\tau)} (1 - e^{-\lambda t_i}) e^{-\lambda t_w} (1 - e^{-\lambda t_m})}{(1 - e^{-\lambda T})}$$ trong đó $\tau$ là độ rộng xung, $T$ là chu kỳ xung, $t_i$ là thời gian chiếu, $t_w$ là thời gian chờ (cooling time), $t_m$ là thời gian đo (counting time), $I_\gamma$ là cường độ vạch gamma đặc trưng (trích xuất từ cơ sở dữ liệu NuDat 2.8), và $f$ là tích các hệ số hiệu chỉnh mất số đếm.

Suất lượng phản ứng $Y$ và tiết diện tích phân $\sigma_{int}$ được tích phân số từ phổ bức xạ hãm $\Phi(E)$ và tiết diện vi phân $\sigma(E)$: $$Y(E_{max}) = N \int_{E_{th}}^{E_{max}} \sigma(E)\Phi(E)dE, \quad \sigma_{int} = \int_{E_{th}}^{E_{max}} \sigma(E)dE$$

Các kiểm định độ tin cậy (robustness checks) được thực hiện bằng cách đối chiếu số liệu diện tích đỉnh qua nhiều vạch gamma khác nhau của cùng một đồng vị phân rã và theo dõi quy luật suy giảm hoạt độ theo chu kỳ bán rã thực nghiệm $T_{1/2}$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Phát hiện quy luật suất lượng phản ứng quang hạt nhân sinh nhiều hạt trên $^{nat}\text{Sr}$: Luận án lần đầu tiên xác định đầy đủ suất lượng tạo thành 10 đồng vị phóng xạ: $^{82}\text{Sr}, ^{83(m+g)}\text{Sr}, ^{85m}\text{Sr}, ^{85g}\text{Sr}, ^{87m}\text{Sr}, ^{81(g+0.976m)}\text{Rb}, ^{82m}\text{Rb}, ^{83}\text{Rb}, ^{84(m+g)}\text{Rb}, ^{86(m+g)}\text{Rb}$ ở các mức năng lượng $E_{\gamma max} = 55, 60, 65\text{ MeV}$. Kết quả chứng minh rằng ở vùng năng lượng này, các kênh phát xạ $(\gamma, xn)$ và $(\gamma, pxn)$ cạnh tranh mạnh mẽ, trong đó tính toán TALYS kết hợp mô hình mật độ mức CTFGM cho kết quả tương thích vượt trội so với các mô hình vi mô HFB.
  2. Định lượng tính phụ thuộc năng lượng của tỷ số ứng phân $^{137m,g}\text{Ce}$ và $^{179m,g}\text{W}$:
    • Cặp đồng phân $^{137m,g}\text{Ce}$ ($J^\pi = 11/2^-$ cho mức $137m\text{Ce}$ và $J^\pi = 3/2^+$ cho mức $137g\text{Ce}$) tạo thành từ phản ứng $^{141}\text{Pr}(\gamma, p3n)$ ở năng lượng 50, 60, 70 MeV cho thấy tỷ số $IR = Y_m/Y_g$ tăng đều theo năng lượng bức xạ hãm.
    • Cặp đồng phân $^{179m,g}\text{W}$ ($J^\pi = 1/2^-$ cho mức $179m\text{W}$ và $J^\pi = 7/2^-$ cho mức $179g\text{W}$) tạo thành từ phản ứng $^{nat}\text{W}(\gamma, xn)$ ở năng lượng 50, 55, 60, 65 MeV cũng thể hiện quy luật tương tự, khẳng định sự ưu tiên truyền mômen xung lượng cao khi năng lượng kích thích tăng.
  3. Xác định tiết diện tích phân đa nơtron trên $^{197}\text{Au}$ và $^{110}\text{Pd}$:
    • Đo đạc thành công tiết diện tích phân của chuỗi phản ứng phát xạ tới 7 nơtron $^{197}\text{Au}(\gamma, xn)^{197-x}\text{Au}$ ($x = 1 - 7$) ở $E_{\gamma max} = 60\text{ MeV}$. Kết quả cho thấy tỷ số giữa thực nghiệm và lý thuyết TALYS 1.9 nằm trong khoảng $0.85 - 1.18$, chứng minh mô hình phát xạ thống kê vẫn duy trì tính hợp lý ở các kênh phát xạ đa hạt bậc cao.
    • Xác định chính xác tiết diện tích phân của $^{110}\text{Pd}(\gamma,n)^{109m}\text{Pd}$, $^{110}\text{Pd}(\gamma,n)^{109g}\text{Pd}$ và $^{110}\text{Pd}(\gamma,X)^{108m}\text{Rh}$ tại $E_{\gamma max} = 70\text{ MeV}$, giải quyết sự bất đồng số liệu lịch sử giữa các nhóm tác giả Mayer và cộng sự.
  4. Hàm kích thích và suất lượng bia dày của chùm proton trên $^{nat}\text{Zr}$ và $^{nat}\text{Pd}$:
    • Đo đạc hàm kích thích phản ứng $^{nat}\text{Zr}(p,X)$ tạo thành $^{95}\text{Zr}, ^{95m}\text{Nb}, ^{95g}\text{Nb}$ trong dải năng lượng rộng $10.6 - 43.6\text{ MeV}$, cung cấp số liệu suất lượng bia dày (thick-target yield) chính xác cho ứng dụng công nghệ vật liệu hạt nhân.
    • Xác định hàm kích thích tích lũy của $^{nat}\text{Pd}(p,X)^{100m,g}\text{Rh}$ và $^{nat}\text{Pd}(p,X)^{100}\text{Pd}$ trong dải $21.09 - 42.61\text{ MeV}$, làm rõ sơ đồ phân rã nhánh phức tạp từ $^{100}\text{Ag} \to ^{100}\text{Pd} \to ^{100}\text{Rh}$.

Implications đa chiều

  • Về mặt lý thuyết học thuật: Bổ sung tập dữ liệu thực nghiệm chuẩn xác vào thư viện EXFOR của IAEA; khẳng định mô hình CTFGM và cơ chế tiền cân bằng hai thành phần (two-component exciton model) là công cụ dự đoán tin cậy nhất cho vùng năng lượng trung bình sau GDR.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thực nghiệm kích hoạt chùm tia Bremsstrahlung năng lượng liên tục kết hợp mô phỏng MCNPX và kỹ thuật kích hoạt đa lá (stacked-foil) với proton trên cyclotron MC-50, có khả năng áp dụng rộng rãi cho các nghiên cứu kích hoạt hạt nhân khác.
  • Về mặt ứng dụng thực tiễn và y học hạt nhân: Cung cấp thông số vật lý chính xác (năng lượng ngưỡng, tiết diện, suất lượng cực đại) phục vụ sản xuất hai đồng vị y tế mang tính sống còn:
    1. Đồng vị $^{82}\text{Sr}$ ($T_{1/2} = 25.34\text{ ngày}$) dùng để nạp generator $^{82}\text{Sr}/^{82}\text{Rb}$ cho chụp cắt lớp phát xạ positron (PET) chẩn đoán tưới máu cơ tim.
    2. Đồng vị $^{100}\text{Pd}$ ($T_{1/2} = 3.63\text{ ngày}$) phân rã sang $^{100g}\text{Rh}$ dùng trong xạ trị áp sát (brachytherapy) điều trị ung thư tuyến tiền liệt.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận các giới hạn vật lý và kỹ thuật nội tại:

  1. Phổ bức xạ hãm liên tục: Việc sử dụng bức xạ hãm Bremsstrahlung từ bia chuyển đổi W tạo ra phổ photon liên tục từ 0 đến $E_{max}$, do đó chỉ có thể đo trực tiếp suất lượng tích phân $Y(E_{max})$ hoặc tiết diện tích phân $\sigma_{int}$, không thể thu được trực tiếp tiết diện vi phân đơn năng $\sigma(E)$ như các nguồn photon đơn sắc (Laser Compton Scattering).
  2. Giới hạn của phương pháp phổ kế gamma trễ (off-line): Phương pháp đo trễ không thể phát hiện các trạng thái kích thích có chu kỳ bán rã cực ngắn ($T_{1/2} < \text{vài giây}$) hoặc các sản phẩm hạt nhân bền không phát bức xạ gamma đặc trưng.
  3. Phân rã chuỗi hạt nhân phức tạp: Đối với các phản ứng trên bia đa đồng vị tự nhiên ($^{nat}\text{Zr}, ^{nat}\text{Pd}, ^{nat}\text{Sr}$), sự chồng chập đóng góp từ nhiều đồng vị bia khác nhau vào cùng một hạt nhân sản phẩm gây khó khăn trong việc tách biệt tuyệt đối tiết diện từng kênh riêng lẻ.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra các hướng đột phá:

  • Ứng dụng nguồn bức xạ photon đơn năng chuẩn trực cao dựa trên tán xạ ngược Laser Compton (LCS) để đo trực tiếp hàm kích thích vi phân $\sigma(E_\gamma)$ trong vùng sau GDR.
  • Phát triển hệ đo phổ gamma tức thời trên đường chùm tia (in-beam prompt gamma-ray spectroscopy) kết hợp detector mảng HPGe và detector nhấp nháy nhanh $\text{LaBr}_3\text{(Ce)}$ để bắt các đồng vị sống siêu ngắn cỡ nano-giây.
  • Mở rộng dải năng lượng proton lên 70 - 100 MeV và nghiên cứu thêm chùm hạt mang điện khác như deuteron, hạt alpha và chùm ion nặng trên các bia siêu tinh khiết làm giàu đồng vị (enriched targets).

Tác động và ảnh hưởng

  • Tác động học thuật quốc tế: Các công trình công bố từ luận án trên các tạp chí chuyên ngành quốc tế uy tín (thuộc danh mục ISI/Web of Science như Physical Review C, Nuclear Instruments and Methods in Physics Research, Applied Radiation and Isotopes, European Physical Journal A) đã đóng góp hàng chục bộ số liệu hạt nhân thực nghiệm mới vào cơ sở dữ liệu EXFOR của Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA), phục vụ các nhà vật lý hạt nhân toàn cầu kiểm chuẩn mô hình.
  • Đột phá công nghệ y tế và công nghiệp: Cung cấp cơ sở khoa học để thiết kế bia bắn phá tối ưu và lựa chọn dải năng lượng chùm proton/chùm photon thích hợp cho các máy gia tốc y tế thương mại, giúp nâng cao độ tinh khiết hạt nhân phóng xạ và giảm thiểu tạp chất phóng xạ nguy hại.
  • An toàn bức xạ và thiết kế che chắn: Dữ liệu tiết diện phát xạ đa neutron và proton là tham số đầu vào tối quan trọng cho việc tính toán độ kích hoạt vật liệu che chắn, đánh giá liều bức xạ và an toàn hạt nhân trong các cơ sở máy gia tốc năng lượng cao.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và nhà khoa học vật lý hạt nhân: Tiếp cận phương pháp luận chuẩn mực về kỹ thuật kích hoạt hạt nhân, phân tích phổ gamma phức tạp và phương pháp so sánh đối chuẩn với mã TALYS, MCNPX, SRIM.
  • Chuyên gia vật lý y học và kỹ sư dược chất phóng xạ: Nắm bắt chính xác các cửa sổ năng lượng phản ứng hạt nhân tối ưu để sản xuất đồng vị y học ($^{82}\text{Sr}, ^{100}\text{Pd}, ^{95}\text{Nb}$) trên máy gia tốc Cyclotron hoặc Linac.
  • Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) và Trung tâm Dữ liệu Hạt nhân: Bổ sung nguồn dữ liệu hạt nhân thực nghiệm có độ tin cậy cao phục vụ đánh giá và cập nhật các thư viện quốc tế (ENDF/B, JENDL, JEFF, TENDL).
  • Kỹ sư thiết kế máy gia tốc và an toàn bức xạ: Ứng dụng số liệu tiết diện tương tác để tính toán liều kích hoạt môi trường và thiết kế hệ thống che chắn an toàn cho các cơ sở gia tốc bức xạ.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Trả lời: Luận án đã mở rộng phạm vi kiểm chứng thực nghiệm của Mô hình giả đơtron Levinger (QDM) và cơ chế tiền cân bằng Exciton lên vùng năng lượng sau GDR (50 - 70 MeV), chứng minh mô hình mật độ mức CTFGM có khả năng mô tả chính xác nhất sự hình thành các trạng thái đồng phân có spin cao thông qua tỷ số ứng phân $IR$.
  2. Cải tiến phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế? Trả lời: So với các nghiên cứu của Uddin et al. (chỉ đo dải thấp trên bia Mo, Zr) và Ditrói et al. (gặp sai số lớn ở chuỗi $A=100$), luận án đã thiết lập quy trình chuẩn hóa hình học chiếu 60 lá kim loại xen kẽ monitor $^{nat}\text{Cu}$, kết hợp đồng thời mô phỏng suy giảm năng lượng SRIM-2003, mô phỏng phổ Bremsstrahlung MCNPX 2.6 và hiệu chỉnh hoạt độ xung $C_{pulse}$, triệt tiêu phần lớn các nguồn sai số hệ thống.
  3. Phát hiện gây bất ngờ nhất từ dữ liệu thực nghiệm là gì? Trả lời: Việc đo đạc thành công tiết diện tích phân của phản ứng bứt tới 7 neutron $^{197}\text{Au}(\gamma, 7n)^{190}\text{Au}$ ở $E_{\gamma max} = 60\text{ MeV}$ cho thấy xác suất bứt nhiều neutron vẫn diễn ra mạnh mẽ và tương thích bất ngờ với lý thuyết thống kê, dù năng lượng kích thích trung bình trên mỗi nucleon tham gia là tương đối hạn chế.
  4. Quy trình tái lập thực nghiệm (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Hoàn toàn đầy đủ và chi tiết. Luận án công bố chi tiết thông số kích thước bia, dòng hạt, phổ chùm tia, thông số detector HPGe ORTEC GEM-2018-p, các vạch năng lượng chuẩn, hàm đa thức hiệu suất $\varepsilon_\gamma(E)$, thời gian chiếu ($t_i$), thời gian chờ ($t_w$), thời gian đo ($t_m$) và các từ khóa đầu vào (keywords input) cho mã tính toán TALYS 1.95.
  5. Định hướng nghiên cứu 10 năm tới được vạch ra như thế nào? Trả lời: Tập trung vào: (1) Chuyển dịch từ phổ bức xạ hãm liên tục sang chùm tia photon đơn năng chuẩn trực phân cực dựa trên công nghệ Laser Compton Backscattering; (2) Phép đo phổ gamma tức thời (prompt gamma) trực tiếp trên chùm tia; (3) Mở rộng nghiên cứu cơ chế phân rã sang các chùm ion nặng và vùng hạt nhân giàu neutron/kém neutron xa dải bền.

Kết luận

  1. Luận án đã đo đạc và xác định thành công một hệ số liệu thực nghiệm đồ sộ, có độ chính xác cao về suất lượng, tiết diện tích phân, hàm kích thích và tỷ số ứng phân của các phản ứng quang hạt nhân và phản ứng proton trên 06 nguyên tố tự nhiên: $^{nat}\text{Sr}$, $^{nat}\text{Zr}$, $^{nat}\text{Pd}$, $^{141}\text{Pr}$, $^{nat}\text{W}$, $^{197}\text{Au}$.
  2. Khẳng định vững chắc vai trò chi phối của cơ chế hạt nhân hợp phần kết hợp cơ chế giả đơtron và phát xạ tiền cân bằng trong vùng năng lượng sau GDR ($50 - 70\text{ MeV}$).
  3. Chứng minh mô hình mật độ mức khí Fermi nhiệt độ không đổi (CTFGM) kết hợp hàm lực gamma Lorentzian tổng quát Kopecky-Uhl (KUGL) trong mã TALYS thể hiện sự tương thích tối ưu nhất với số liệu thực nghiệm.
  4. Bổ sung các bộ dữ liệu hạt nhân thực nghiệm quý giá vào thư viện hạt nhân quốc tế EXFOR/IAEA, giải tỏa tình trạng khan hiếm dữ liệu phản ứng sinh nhiều hạt ở năng lượng trung bình.
  5. Cung cấp các thông số vật lý chuẩn mực phục vụ công nghệ sản xuất các đồng vị phóng xạ y học quan trọng trong chẩn đoán PET ($^{82}\text{Sr}/^{82}\text{Rb}$) và điều trị ung thư ($^{100}\text{Pd}/^{100m}\text{Rh}$).
  6. Thiết lập quy trình thực nghiệm và tính toán lý thuyết mẫu mực, mở ra các hướng nghiên cứu hiện đại tiếp theo cho ngành Vật lý hạt nhân thực nghiệm tại Việt Nam và trên trường quốc tế.