Tổng quan về luận án
Nghiên cứu cấu trúc và động lực học hạt nhân thông qua tương tác giữa các nucleon là một trong những bài toán trung tâm của Vật lý Hạt nhân hiện đại. Trong các phản ứng hạt nhân, hạt nhân tàn dư (residual nucleus) có thể được hình thành ở trạng thái cơ bản (ground state - $g$) hoặc trạng thái đồng phân siêu bền (metastable/isomeric state - $m$). Trạng thái đồng phân hạt nhân (nuclear isomer) là trạng thái kích thích có thời gian sống dài (từ nano giây đến hàng triệu năm) do sự kìm hãm mạnh mẽ của các quy tắc chọn lọc phân rã điện từ ($\gamma$-decay) và chuyển dời nội chuyển đổi (internal conversion - IC). Đại lượng tỷ số đồng phân (Isomeric Ratio - IR), được định nghĩa là tỷ số giữa xác suất tạo thành trạng thái đồng phân và trạng thái cơ bản ($IR = \sigma_m / \sigma_g$ đối với chùm hạt đơn năng hoặc $IR = Y_m / Y_g$ đối với phổ hạt liên tục), đóng vai trò như một "đầu dò" (probe) vi mô cực kỳ nhạy cảm để giải mã cấu trúc mức năng lượng, mật độ mức (nuclear level density), sự phân bố spin của hạt nhân kích thích, cũng như cơ chế phản ứng hạt nhân phức hợp.
Mặc dù có vai trò tối quan trọng trong vật lý hạt nhân cơ bản và ứng dụng, dữ liệu thực nghiệm về tỷ số đồng phân hiện nay vẫn tồn tại nhiều khoảng trống nghiên cứu (research gaps) mang tính then chốt. Cụ thể:
- Thiếu hụt dữ liệu thực nghiệm và sự phân tán số liệu: Đối với các đồng vị nặng như $^{196,198}\text{Hg}$ và $^{153}\text{Eu}$ trong phản ứng quang hạt nhân $(\gamma, n)$ ở vùng cộng hưởng khổng lồ lưỡng cực điện (Giant Dipole Resonance - GDR), các dữ liệu thực nghiệm trước đây rất khan hiếm hoặc có sự bất đồng lớn giữa các nhóm nghiên cứu quốc tế.
- Khoảng trống trong phản ứng bắt neutron: Hoàn toàn chưa có dữ liệu đo đạc thực nghiệm về IR cho phản ứng bắt neutron hỗn hợp nhiệt - trên nhiệt (mixed thermal-epithermal) trên các đồng vị $^{108,110}\text{Pd}$ và $^{114,116}\text{Cd}$, cũng như phản ứng bắt neutron cộng hưởng (resonance neutron capture) trên hạt nhân $^{110}\text{Pd}$.
- Thách thức trong việc kiểm chứng mô hình lý thuyết: Các chương trình mô phỏng phản ứng hạt nhân hiện đại như TALYS chỉ tính toán tiết diện với chùm hạt đơn năng, trong khi các thí nghiệm quang hạt nhân chủ yếu sử dụng bức xạ hãm (bremsstrahlung) liên tục; việc tích hợp mô phỏng bức xạ hãm bằng công cụ Monte Carlo GEANT4 với mã TALYS để đánh giá đồng thời 6 mô hình mật độ mức và 8 hàm lực bức xạ gamma chưa được khảo sát một cách có hệ thống.
Luận án tiến sĩ của tác giả Bùi Minh Huệ với tiêu đề "Study of Isomeric Ratio and Related Effects in Photonuclear and Neutron Capture Reactions" (Chuyên ngành: Vật lý nguyên tử và hạt nhân, Mã số: 9440106), dưới sự hướng dẫn khoa học của GS. Trần Đức Thiệp và TS. Sergey Mikhailovich Lukyanov, đã giải quyết trọn vẹn các khoảng trống nói trên thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:
- Câu hỏi nghiên cứu 1 (RQ1): Quy luật biến thiên thực nghiệm của IR theo năng lượng kích thích bức xạ hãm trong toàn bộ dải GDR và theo các dải năng lượng neutron khác nhau (nhiệt, cộng hưởng, hỗn hợp) diễn ra như thế nào?
- Câu hỏi nghiên cứu 2 (RQ2): Các hiệu ứng cấu trúc hạt nhân (độ chênh lệch spin $\Delta J$, cấu hình nucleon, năng lượng kích thích $E^*$, và kênh phản ứng hạt nhân) tác động định lượng ra sao đến giá trị IR?
- Câu hỏi nghiên cứu 3 (RQ3): Tổ hợp mô hình mật độ mức (Level Density - LD) và hàm lực bức xạ gamma ($\gamma$-ray Strength Function - $\gamma\text{SF}$) nào trong mã tính toán vi mô TALYS kết hợp với phổ mô phỏng GEANT4 mô tả chính xác nhất các dữ liệu thực nghiệm thu được?
- Giả thuyết khoa học 1 (H1): Trong phản ứng $(\gamma, n)$ vùng GDR, sự gia tăng năng lượng bức xạ hãm ban đầu sẽ làm tăng xác suất quần cư trạng thái có spin cao, sau đó giá trị IR sẽ đạt trạng thái bão hòa ở cuối dải GDR do sự xuất hiện của các quá trình tiền cân bằng (pre-equilibrium) và tương tác trực tiếp (direct reaction).
- Giả thuyết khoa học 2 (H2): Xung lượng truyền (intake impulse) từ hạt đạn vào hạt nhân bia càng lớn thì xác suất tạo thành trạng thái đồng phân spin cao càng ưu thế, dẫn đến sự phân tách rõ rệt về IR giữa các kênh phản ứng khác nhau ($(n, 2n) > (n, \gamma) \approx (\gamma, n)$).
- Giả thuyết khoa học 3 (H3): Hiệu ứng cấu hình nucleon thể hiện tính quy luật rõ rệt: trong cùng một chuỗi đồng vị (isotopes) hoặc đồng neutron (isotones $N=81$), giá trị IR có xu hướng giảm dần khi số khối $A$ hoặc số proton $Z$ tăng lên.
Khung lý thuyết của luận án được xây dựng trên nền tảng cơ học lượng tử hạt nhân, thuyết phản ứng thống kê Hauser-Feshbach, mô hình Huizenga-Vandebosch (HVM), thuyết khí Fermi dịch chuyển lùi (Back-Shifted Fermi Gas - BSFG), mô hình siêu lỏng suy rộng (Generalized Superfluid Model - GSM), và các mô hình vi mô Hartree-Fock-Bogolyubov (HFB). Về mặt phạm vi, luận án khảo sát toàn diện 6 cặp đồng phân hạt nhân ($^{152m1,m2}\text{Eu}$, $^{195m,g;197m,g}\text{Hg}$, $^{115m,g}\text{Cd}$, $^{109m,g}\text{Pd}$, $^{137m,g}\text{Ce}$, $^{81m,g}\text{Se}$) trên máy gia tốc vi mô MT-25 Microtron tại Phòng thí nghiệm Phản ứng Hạt nhân Flerov (FLNR), Viện Liên hiệp Nghiên cứu Hạt nhân Dubna (JINR, Nga).
Literature Review và Positioning
Lịch sử nghiên cứu đồng phân hạt nhân bắt đầu từ năm 1917 khi Soddy đưa ra dự đoán: "We can have isotopes with the identity of atomic weight, as well as of chemical character, which is different in their stability and mode of breaking up" [1]. Năm 1921, Otto Hahn phát hiện ra trường hợp đồng phân thực nghiệm đầu tiên ở muối Uranium ($^{234}\text{Pa}$ và $^{234m}\text{Pa}$) [2]. Đến năm 1936, von Weizsäcker đã cung cấp lời giải thích lý thuyết đầu tiên, khẳng định rằng sự kìm hãm bức xạ gamma xuất phát từ độ chênh lệch momen động lượng lớn giữa trạng thái kích thích và trạng thái thấp hơn kết hợp với năng lượng chuyển dời nhỏ [3]. Đầu thập niên 1950, Bohr và Mottelson mở rộng lý thuyết sang các hạt nhân biến dạng đối xứng trục ($^{180}\text{Hf}$), định hình khái niệm $K$-isomer dựa trên quy tắc bảo toàn hình chiếu spin trên trục đối xứng [4].
Trong lĩnh vực phản ứng quang hạt nhân, Chadwick và Goldhaber (1934) lần đầu tiên quan sát sự phân rã quang hạt nhân trên bia deuteri [77], tiếp nối bởi phát hiện bức xạ hãm trên Betatron của Baldwin và Klaiber (1947) [79]. Hiện tượng cộng hưởng khổng lồ lưỡng cực điện (GDR) được mô tả thành công qua mô hình thủy động lực học hai chất lưu proton-neutron của Goldhaber - Teller (1948) [87] và Steinwedel - Jensen (1950) [88].
┌────────────────────────────────────────────────────────────────────────────────┐
│ TIẾN TRÌNH LÝ THUYẾT & THỰC NGHIỆM │
└────────────────────────────────────────────────────────────────────────────────┘
1917: Dự đoán của Soddy về các trạng thái phóng xạ đồng khối khác biệt
│
1921: Otto Hahn phát hiện 234m,g Pa (Đồng phân hạt nhân thực nghiệm đầu tiên)
│
1936: von Weizsäcker giải thích cơ chế kìm hãm gamma (Chênh lệch Spin cao)
│
1952: Bohr & Mottelson xây dựng lý thuyết K-isomer trên hạt nhân biến dạng
│
1960: Huizenga & Vandebosch phát triển mô hình thống kê tính IR (HVM)
│
1990s: Mô hình Hauser-Feshbach nâng cao tích hợp Cascade gamma (STAPRE/COMPLET)
│
Hiện đại: TALYS 1.95 + GEANT4 Monte Carlo + Nguồn gia tốc Microtron MT-25
Trong dòng nghiên cứu tính toán lý thuyết, Huizenga và Vandebosch (1960) đã đặt nền móng cho mô hình HVM dựa trên lý thuyết thống kê Hauser-Feshbach [55, 56], tính toán sự phân bố spin qua từng giai đoạn phát xạ hạt và tầng cascade gamma. Tuy nhiên, mô hình HVM cổ điển bộc lộ nhiều hạn chế do giả định thông số cắt spin ($\sigma$) là một hằng số không phụ thuộc năng lượng, đồng thời chỉ xét chuyển dời điện cực $E1$ mà bỏ qua cấu trúc mức rời rạc thực tế.
Trong các tranh luận học thuật quốc tế gần đây, hai luồng quan điểm nổi bật xoay quanh sự phụ thuộc của tỷ số đồng phân:
- Tranh luận về cơ chế phản ứng ở đuôi GDR: Một số tác giả (như Danagulyan et al. [66], Mazur et al. [65]) cho rằng ở vùng năng lượng trên đỉnh GDR, các quá trình phát xạ nhiều hạt $( \gamma, 2n )$, $( \gamma, 3n )$ và tương tác trực tiếp làm thay đổi mạnh mẽ IR; trong khi các nghiên cứu của Tonchev et al. [34] và Thiep et al. [48] chỉ ra rằng đối với kênh đơn hạt $(\gamma, n)$, tỷ số IR đạt giá trị bình nguyên (plateau) ổn định do sự triệt tiêu tương đối giữa các quá trình tiền cân bằng.
- Tranh luận về độ nhạy của các mô hình mật độ mức: Khi áp dụng mã mô phỏng TALYS, các nghiên cứu quốc tế ghi nhận sự sai khác đáng kể giữa mô hình giải tích hiện tượng học (Phenomenological models: CTM, BSFG, GSM) và mô hình vi mô bán cổ điển (Microscopic models: HFB Goriely, HFB Hilaire).
So sánh với các công trình quốc tế tiêu biểu:
- So sánh với Danagulyan et al. (2014, 2017) [60, 66]: Nhóm Danagulyan sử dụng TALYS 1.4 khảo sát IR trên các bia Sn, Te, In, Hf nhưng nhận thấy sự sai lệch lớn ở các trạng thái spin rất cao do TALYS 1.4 đánh giá thấp mật độ mức ở vùng spin ngoại biên. Luận án của Bùi Minh Huệ đã khắc phục điều này bằng cách sử dụng phiên bản TALYS 1.95 tích hợp cơ sở dữ liệu RIPL mới nhất, kết hợp hiệu chỉnh thông số tán xạ neutron.
- So sánh với Makwana et al. (2020) [67]: Makwana thiết lập công thức bán thực nghiệm cho tiết diện $(\gamma, n)$ gần đỉnh GDR cho $Z > 60$. Luận án mở rộng việc kiểm chứng này sang cấu trúc đồng phân phức tạp với 3 trạng thái phân nhánh trên $^{152}\text{Eu}$ ($3^-$, $8^-$, $0^-$), cung cấp bằng chứng thực nghiệm trực tiếp mà công thức của Makwana chưa bao quát tới.
Vị trí của luận án được xác định rõ ràng: Là công trình tiên phong cung cấp số liệu thực nghiệm nguyên gốc (original experimental data) với độ chính xác cao cho các hệ hạt nhân chưa từng được công bố, đồng thời thiết lập khung tính toán kết hợp mô phỏng bức xạ hãm GEANT4 và mô hình vi mô TALYS 1.95 một cách hoàn chỉnh.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đã mở rộng và kiểm chứng các lý thuyết nền tảng trong vật lý hạt nhân:
- Mở rộng mô hình Huizenga-Vandebosch (HVM) và lý thuyết Hauser-Feshbach: Thay vì sử dụng thông số cắt spin cố định $\sigma$, luận án đã chứng minh sự cần thiết phải tích hợp sự phụ thuộc năng lượng của mômen quán tính hạt nhân $\theta(E^)$ và thông số cắt spin $\sigma(E^) = \sqrt{\theta_{\text{rigid}} T / \hbar^2}$. Đồng thời, cascade gamma được mô tả với đầy đủ các chuyển dời đa cực $E1, M1, E2$ dựa trên cấu trúc mức rời rạc thực tế trích xuất từ ENSDF.
- Khung khái niệm về cơ chế hình thành đồng phân qua 3 giai đoạn:
- Giai đoạn 1 - Hấp thụ photon: Hạt nhân bia spin $J_0$ hấp thụ lượng tử $E1$, tạo thành nhân hợp phần ở trạng thái kích thích với spin $J_c = J_0, J_0 \pm 1$, xác suất tỷ lệ với $P(J_c) \approx 2J_c + 1$.
- Giai đoạn 2 - Bốc hơi neutron: Nhân hợp phần phát xạ neutron nhiệt mang mômen động lượng quỹ đạo $l$, hình thành nhân tàn dư ở năng lượng $E^* = U - S_n - \bar{E}n$, với phân bố spin tuân theo:
$$P(J) \approx \sum{S=|J-1/2|}^{|J+1/2|} \sum_{l=|J_c-S|}^{|J_c+S|} \rho(J, E^*) T_l(E_n)$$
- Giai đoạn 3 - Cascade bức xạ gamma và phân nhánh spin: Quá trình de-excitation diễn ra liên tục cho đến khi năng lượng rơi vào "vùng cắt gamma" ($E_u, E_d$). Tại đây, "lượng tử gamma quyết định" (deciding $\gamma$-ray) sẽ phân nhánh về trạng thái đồng phân hoặc trạng thái cơ bản dựa trên tâm spin $COS = (S_h + S_l)/2$. Trạng thái có spin $J > COS$ sẽ ưu tiên nuôi dưỡng mức đồng phân spin cao.
┌─────────────────────────────────────────────────────────────────────────────────┐
│ SƠ ĐỒ KHUNG PHÂN TÍCH QUÁ TRÌNH HÌNH THÀNH ĐỒNG PHÂN │
└─────────────────────────────────────────────────────────────────────────────────┘
[Hạt nhân bia (J0)] + Photon E1 / Neutron s-wave
│
▼
[Hạt nhân hợp phần kích thích (Jc)] ──► E* ≈ Sn (Bắt neutron) hoặc EGDR
│
▼ (Bốc hơi neutron mang momen l - trong phản ứng quang hạt nhân)
[Hạt nhân tàn dư kích thích (E*, J)] ──► Mật độ mức ρ(J, E*)
│
▼ (Cascade gamma E1, M1, E2)
[Vùng năng lượng cắt (Eu, Ed)]
│
┌──────┴─────────────────────────┐
▼ (Spin J > COS) ▼ (Spin J < COS)
[Trạng thái Đồng phân (Isomer - m)] [Trạng thái Cơ bản (Ground - g)]
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đồng thời 3 trụ cột lý thuyết:
- Thuyết phản ứng nhân hợp phần & Thống kê vi mô: Tích hợp các mô hình mật độ mức hiện tượng học (CTM, BSFG, GSM) và giải tích vi mô (HFB Goriely, HFB Hilaire, HFB Koning).
- Lý thuyết tương tác bức xạ điện từ: Sử dụng 8 hàm lực bức xạ gamma bao gồm Mô hình Lorentzian Chuẩn (Standard Lorentzian - SLO), Mô hình Lorentzian Mở rộng (Generalized Lorentzian - GLO), Mô hình SMLO, và các hàm lực vi mô QRPA.
- Mô hình động học hạt nhân: Phân tích sự đóng góp của quá trình tiền cân bằng qua mô hình Exciton hai thành phần (two-component exciton model).
Điều kiện biên: Khung phân tích áp dụng cho vùng kích thích năng lượng bức xạ hãm $E_{\gamma,\text{endpoint}} = 10 - 25\text{ MeV}$ (vùng GDR) và vùng neutron năng lượng thấp ($E_n < 300\text{ eV}$); hạt nhân mục tiêu có khối lượng trung bình và nặng ($A = 74 - 198$).
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ lập trường triết học thực chứng nghiêm ngặt (strict positivism), kết hợp chặt chẽ giữa kỹ thuật thực nghiệm đo phổ gamma kích hoạt độ chính xác cao và mô phỏng số ngẫu nhiên Monte Carlo quy mô lớn.
Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm:
- Tầng 1 - Nguồn chùm hạt: Sử dụng máy gia tốc vi mô MT-25 Microtron (FLNR, JINR Dubna) để tạo chùm electron năng lượng cao ($14 - 25\text{ MeV}$). Chùm electron đập vào bia vonfram (W) dày làm mát bằng nước để sinh bức xạ hãm liên tục, hoặc chiếu xạ vào bia làm chậm paraffin/nước để sinh dòng neutron nhiệt và trên nhiệt.
- Tầng 2 - Chiếu xạ mẫu: Sử dụng kỹ thuật kích hoạt hạt nhân (Activation Technique). Mẫu kim loại có độ sạch hóa học cao ($>99.9%$) của Hg, Eu, Cd, Pd, Ce, Se được chiếu xạ với thời gian kiểm soát chính xác. Đối với phản ứng bắt neutron, phương pháp "bọc lá Cadmi" (Cd-foil cutoff method, chiều dày $0.5 - 1.0\text{ mm}$) được áp dụng để tách biệt tuyệt đối giữa tương tác của neutron nhiệt ($E_n < 0.5\text{ eV}$) và neutron cộng hưởng ($10 - 300\text{ eV}$).
- Tầng 3 - Đo phổ bức xạ gamma ngoại tuyến (Off-line $\gamma$-ray Spectrometry): Sử dụng phổ kế gamma bán dẫn siêu tinh khiết HPGe (Canberra) có độ phân giải năng lượng cao ($FWHM \approx 1.8\text{ keV}$ tại vạch $1332.5\text{ keV}$ của $^{60}\text{Co}$), đặt trong buồng chì che chắn bức xạ phông dày $10\text{ cm}$.
┌─────────────────────────────────────────────────────────────────────────────────┐
│ QUY TRÌNH THỰC NGHIỆM VÀ XỬ LÝ SỐ LIỆU │
└─────────────────────────────────────────────────────────────────────────────────┘
[Máy gia tốc Microtron MT-25] ──► Chùm e- (14 - 25 MeV)
│
├────────────────────────┬────────────────────────┐
▼ ▼ ▼
[Bia biến đổi W] [Bia chậm Paraffin] [Bọc lá Cd (0.5mm)]
│ │ │
Phổ Bremsstrahlung Neutron nhiệt Neutron cộng hưởng
│ │ │
└────────────────────────┴────────────────────────┘
│
▼
[Chiếu xạ Mẫu (Hg, Eu, Cd, Pd...)]
│ (Làm nguội t_cool)
▼
[Đo phổ gamma bằng đầu dò HPGe]
│
▼
[Hiệu chỉnh: Tự hấp thụ + Trùng phùng ngẫu nhiên]
│
▼
[Xác định diện tích đỉnh phổ S_gamma & Tỷ số IR]
Quy trình nghiên cứu rigorous
Nhằm triệt tiêu sai số hệ thống và đảm bảo độ tin cậy tuyệt đối, quy trình đo đạc thực nghiệm áp dụng các hiệu chỉnh vật lý bắt buộc:
-
Hiệu chỉnh tự hấp thụ bức xạ gamma trong mẫu (Self-absorption correction factor - $F_{abs}$):
Do mẫu có chiều dày hữu hạn và mật độ lớn, các bức xạ gamma năng lượng thấp bị suy giảm khi thoát ra khỏi mẫu. Hệ số hiệu chỉnh $F_{abs}$ được tính toán giải tích theo hình học mẫu:
$$F_{abs} = \frac{1 - e^{-\mu(E_\gamma) \rho d}}{\mu(E_\gamma) \rho d}$$
Trong đó $\mu(E_\gamma)$ là hệ số suy giảm khối ($cm^2/g$), $\rho$ là khối lượng riêng ($g/cm^3$), và $d$ là chiều dày mẫu ($cm$).
-
Hiệu chỉnh hiệu ứng trùng phùng nối tiếp thực (True Coincidence Summing - $F_{sum}$):
Khi hạt nhân phân rã phát ra nhiều photon gamma cascade đồng thời, sự ghi nhận cùng lúc hai photon trong thời gian phân giải của đầu dò làm sai lệch diện tích đỉnh hấp thụ toàn phần. Hệ số $F_{sum}$ được hiệu chỉnh chi tiết dựa trên ma trận xác suất chuyển dời và đường cong hiệu suất đỉnh/tổng của đầu dò HPGe.
-
Xác định tỷ số IR thực nghiệm từ diện tích phổ gamma:
Tỷ số diện tích đỉnh phổ thực nghiệm $S_\gamma$ ở các thời điểm đo khác nhau được liên hệ với tốc độ tạo thành đồng phân $Y_m$ và trạng thái cơ bản $Y_g$ thông qua phương trình phân rã phóng xạ tích phân:
$$S_\gamma = N_0 \cdot \epsilon(E_\gamma) \cdot I_\gamma \cdot Y \cdot \frac{(1 - e^{-\lambda t_{irr}}) e^{-\lambda t_{cool}} (1 - e^{-\lambda t_{meas}})}{\lambda}$$
Trong đó $\epsilon(E_\gamma)$ là hiệu suất đỉnh đầu dò, $I_\gamma$ là cường độ bức xạ tuyệt đối, $\lambda$ là hằng số phân rã, $t_{irr}, t_{cool}, t_{meas}$ lần lượt là thời gian chiếu xạ, thời gian làm nguội và thời gian đo phổ.
Data và phân tích
Phần mô phỏng lý thuyết được thực hiện thông qua liên kết hai nền tảng tính toán mạnh mẽ:
- Mô phỏng GEANT4 (CERN): Xây dựng ứng dụng người dùng mô phỏng chính xác cấu trúc hình học của máy gia tốc MT-25, bia chuyển đổi W, bộ chuẩn trực, và vị trí đặt mẫu. Chạy mô phỏng Monte Carlo với $500 \times 10^6$ (500 triệu) hạt electron sơ cấp để thu được phổ thông lượng năng lượng photon hãm $\phi(E_\gamma)$ và thông lượng neutron $E_n$.
- Mã phản ứng hạt nhân TALYS 1.95: Tính toán tiết diện vi sai tạo thành trạng thái đồng phân $\sigma_m(E)$ và trạng thái cơ bản $\sigma_g(E)$ bằng cách quét độc lập 6 mô hình mật độ mức (Level Density: 1-CTM, 2-BSFG, 3-GSM, 4-HFB Goriely, 5-HFB Hilaire, 6-HFB Koning) kết hợp cùng 8 hàm lực bức xạ gamma ($\gamma\text{SF}$: 1-Kopecky-Uhl GLO, 2-Brink-Axel SLO, 3-Modified Lorentzian, 4-Goriely HFB, 5-QRPA, v.v.).
- Mã phân tích CERN ROOT: Lập trình mã phân tích chuyên biệt để tích phân số học tích chập giữa phổ thông lượng GEANT4 và tiết diện TALYS:
$$IR_{theor} = \frac{\int_{E_{th}}^{E_{max}} \sigma_m(E_\gamma) \phi(E_\gamma) dE_\gamma}{\int_{E_{th}}^{E_{max}} \sigma_g(E_\gamma) \phi(E_\gamma) dE_\gamma}$$
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Công bố tập dữ liệu thực nghiệm đầu tiên trên thế giới:
Luận án đã đo đạc và báo cáo lần đầu tiên các giá trị IR cho các cặp đồng phân sau:
- Cặp $^{195m,g}\text{Hg}$ trong dải năng lượng bức xạ hãm $14 - 24\text{ MeV}$ và $^{197m,g}\text{Hg}$ trong dải $18 - 24\text{ MeV}$ từ phản ứng $(\gamma, n)$.
- Bộ ba đồng phân phức tạp $^{152m1,m2,g}\text{Eu}$ ($J^{\pi} = 8^-, 0^-, 3^-$) tại các năng lượng $19, 21, 23\text{ MeV}$ trong phản ứng $(\gamma, n)$.
- Cặp $^{109m,g}\text{Pd}$, $^{111m,g}\text{Pd}$ và $^{115m,g}\text{Cd}$, $^{117m,g}\text{Cd}$ trong phản ứng bắt neutron hỗn hợp nhiệt - cộng hưởng $(n, \gamma)$, cùng phép đo đầu tiên cho $^{111m,g}\text{Pd}$ trong phản ứng bắt neutron cộng hưởng thuần túy.
-
Xác lập quy luật ảnh hưởng của 4 hiệu ứng hạt nhân:
- Hiệu ứng năng lượng kích thích: Trong phản ứng quang hạt nhân $(\gamma, n)$, giá trị IR tăng đơn điệu theo năng lượng bức xạ hãm từ ngưỡng phản ứng đến đỉnh vùng GDR, sau đó đạt trạng thái bão hòa ổn định ở vùng $22 - 25\text{ MeV}$. Điều này chứng minh sự cân bằng giữa phát xạ bốc hơi và quá trình tiền cân bằng.
- Hiệu ứng độ chênh lệch spin ($\Delta J$): Trạng thái kích thích có độ chênh lệch spin $|\Delta J| = |J_m - J_g|$ càng lớn thì xác suất hình thành đồng phân càng giảm mạnh. Ví dụ, đối với các đồng vị Bari ($^{129m,g}\text{Ba}, ^{131m,g}\text{Ba}, ^{133m,g}\text{Ba}$) có cùng spin trạng thái cơ bản $1/2^+$ nhưng spin trạng thái đồng phân tăng dần ($7/2^+, 9/2^-, 11/2^-$), giá trị IR suy giảm theo cấp số nhân.
- Hiệu ứng kênh phản ứng và xung lượng truyền: So sánh cùng một cặp đồng phân tạo ra từ các phản ứng khác nhau cho thấy:
$$IR_{(n, 2n)} > IR_{(n, \gamma)} \approx IR_{(\gamma, n)}$$
Nguyên nhân là do neutron có khối lượng nghỉ khác không, khi tương tác sẽ truyền xung lượng góc lớn hơn nhiều so với photon điện từ ($m_0 = 0$), tạo điều kiện thuận lợi áp đảo cho các mức spin cao.
- Hiệu ứng cấu hình nucleon: Trong chuỗi hạt nhân isotone $N=81$ ($^{137}\text{Ba}, ^{139}\text{Ce}, ^{141}\text{Nd}, ^{143}\text{Sm}$) có cùng spin và cấu trúc vỏ ngoài, giá trị IR của hai hạt nhân nhẹ hơn lớn gấp đôi hai hạt nhân nặng hơn, khẳng định sự phụ thuộc chặt chẽ của IR vào số proton và số khối $A$.
┌────────────────────────────────────────────────────────────────────────────────┐
│ SO SÁNH HIỆU ỨNG KÊNH PHẢN ỨNG TRÊN GIÁ TRỊ IR │
└────────────────────────────────────────────────────────────────────────────────┘
Kênh (n, 2n) : ───► [Xung lượng truyền lớn (Neutron m > 0)] ───► IR rất cao
Kênh (n, γ) : ───► [Hấp thụ s-wave, truyền spin 1/2] ───► IR trung bình
Kênh (γ, n) : ───► [Photon m=0, hấp thụ lưỡng cực E1] ───► IR thấp
- Đánh giá định lượng 48 tổ hợp mô hình trong TALYS 1.95:
Kết quả so sánh giữa lý thuyết và thực nghiệm chỉ ra rằng mô hình khí Fermi dịch chuyển lùi (BSFG - Back-Shifted Fermi Gas) kết hợp với hàm lực Lorentzian suy rộng (GLO - Generalized Lorentzian Model) tái hiện số liệu thực nghiệm tốt nhất với độ lệch tương đối nhỏ hơn $8.5%$ trên toàn dải năng lượng khảo sát.
Implications đa chiều
- Đóng góp lý thuyết: Cung cấp bằng chứng thực nghiệm vững chắc để chuẩn hóa các tham số cắt spin, tham số mật độ mức $a$, và độ biến dạng hạt nhân trong thư viện dữ liệu hạt nhân RIPL-3 của IAEA.
- Đổi mới phương pháp luận: Thiết lập quy trình tính toán chuẩn kết hợp GEANT4 và TALYS 1.95, mở ra phương pháp tiếp cận chính xác để phân tích các phản ứng quang hạt nhân sử dụng nguồn bức xạ hãm liên tục mà không cần đến các nguồn photon đơn sắc đắt tiền.
- Ứng dụng thực tiễn trong Y học hạt nhân: Tối ưu hóa các thông số chiếu xạ để sản xuất các đồng phân phóng xạ y học có độ tinh khiết hạt nhân cao, ví dụ như $^{197m}\text{Hg}$ (chẩn đoán SPECT), $^{115m}\text{Cd}$ và $^{152m}\text{Eu}$ (nghiên cứu đánh dấu phóng xạ và đo liều bức xạ).
- Ứng dụng trong Công nghệ Lò phản ứng & Độc học hạt nhân: Dữ liệu tiết diện và IR của $^{109,111}\text{Pd}$ và $^{115,117}\text{Cd}$ là tham số đầu vào thiết yếu cho việc tính toán cháy nhiên liệu, thiết kế thanh điều khiển hấp thụ neutron, và công nghệ chuyển hóa rác thải hạt nhân (nuclear transmutation).
Limitations và Future Research
Luận án thừa nhận một số hạn chế khách quan mang tính kỹ thuật:
- Đặc tính nguồn bức xạ hãm liên tục: Do sử dụng bức xạ hãm từ bia chuyển đổi W, kết quả thu được là giá trị tích phân theo thông lượng năng lượng photon hãm, chưa thể tách biệt thành tiết diện đơn năng tuyệt đối tại từng điểm năng lượng hẹp nếu không dùng phương pháp nghịch đảo giải tích phức tạp.
- Hạn chế của mô hình TALYS đối với neutron nhiệt: Mã TALYS 1.95 hiện tại chưa giải quyết trọn vẹn cơ chế bắt neutron trong vùng năng lượng nhiệt và trên nhiệt ($E_n < 0.5\text{ eV}$), đòi hỏi phải dựa vào các mô hình giải tích phụ trợ.
- Độ bất định trong số liệu phân rã phổ học: Một số nhánh phân rã gamma yếu của các đồng vị đồng phân có độ bất định tương đối trong bảng dữ liệu ENSDF lên tới $5 - 10%$, ảnh hưởng trực tiếp đến sai số tổng hợp của IR.
Chương trình nghiên cứu tiếp nối trong tương lai:
- Hướng 1: Sử dụng nguồn chùm photon đơn năng phân cực cao tạo từ cơ chế tán xạ ngược Laser-Compton (Laser Compton Backscattering - LCS) tại các cơ sở hiện đại như ELI-NP (Romania) hoặc NewSUBARU (Nhật Bản) để đo trực tiếp tiết diện quang hạt nhân đơn năng vi sai.
- Hướng 2: Mở rộng nghiên cứu sang các đồng phân cực thấp năng lượng (Extremely Low Energy - ELE isomers, ví dụ mức $8.3\text{ eV}$ của $^{229m}\text{Th}$), phục vụ chế tạo đồng hồ hạt nhân siêu chính xác (nuclear clock).
- Hướng 3: Khảo sát sâu cấu trúc vi mô của cộng hưởng Pygmy (Pygmy Dipole Resonance - PDR, Pygmy Quadrupole Resonance - PQR) và cộng hưởng Scissors $M1$ đối với cascade gamma phân nhánh.
Tác động và ảnh hưởng
┌────────────────────────────────────────────────────────────────────────────────┐
│ MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG │
└────────────────────────────────────────────────────────────────────────────────┘
HỌC THUẬT QUỐC TẾ ──► Bổ sung dữ liệu EXFOR (IAEA), chuẩn hóa thư viện RIPL-3
Y HỌC PHÓNG XẠ ──► Tối ưu hóa chu trình sản xuất đồng vị 197m-Hg, 115m-Cd
NĂNG LƯỢNG NGUYÊN TỬ──► Tham số chính xác cho thiết kế lò phản ứng và thanh hấp thụ
ĐÀO TẠO & HỢP TÁC ──► Cầu nối nghiên cứu đỉnh cao giữa VAST (Việt Nam) & JINR Dubna
- Tác động học thuật: Các kết quả đo đạc của luận án đã được công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Physical Review C, European Physical Journal A, Nuclear Physics A), cung cấp số liệu thực nghiệm quý giá cho Ngân hàng Dữ liệu Hạt nhân Quốc tế EXFOR do Cơ quan Năng lượng Nguyên tử Quốc tế (IAEA) quản lý.
- Tác động công nghiệp và Y tế: Cung cấp bộ thông số hạt nhân chuẩn xác để các viện nghiên cứu hạt nhân tối ưu hóa quy trình điều chế đồng vị phóng xạ trên máy gia tốc và lò phản ứng nghiên cứu, giảm thiểu tạp chất phóng xạ không mong muốn.
- Tác động chính sách và xã hội: Khẳng định năng lực làm chủ công nghệ gia tốc và đo phổ hạt nhân đỉnh cao của các nhà khoa học Việt Nam trong mạng lưới hợp tác quốc tế với Viện JINR Dubna, đóng góp vào chiến lược phát triển khoa học công nghệ hạt nhân hòa bình tại Việt Nam.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Nhà nghiên cứu Vật lý hạt nhân: Tiếp cận phương pháp luận mẫu mực trong việc kết hợp giữa kỹ thuật đo phổ gamma thực nghiệm và công cụ mô phỏng GEANT4/TALYS/ROOT.
- Các nhà phát triển mã mô phỏng hạt nhân (TALYS/EMPIRE Developers): Sử dụng dữ liệu thực nghiệm để hiệu chỉnh các tham số vi mô về mật độ mức và hàm lực bức xạ gamma.
- Kỹ sư R&D Y học Hạt nhân & Dược chất phóng xạ: Ứng dụng các tỷ số phản ứng để tính toán liều lượng, hiệu suất thu hồi và độ làm giàu của các đồng vị đồng phân dùng trong chẩn đoán và điều trị ung thư.
- Chuyên gia Vật lý Lò phản ứng & An toàn Bức xạ: Sử dụng dữ liệu tiết diện bắt neutron của Pd và Cd để mô phỏng chính xác sự chuyển hóa vật liệu trong vùng hoạt lò phản ứng hạt nhân.
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ì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc làm sáng tỏ cơ chế bão hòa của tỷ số đồng phân ở đuôi vùng cộng hưởng khổng lồ lưỡng cực điện (GDR) và lượng hóa quy luật suy giảm IR theo spin $\Delta J$ và số khối $A$. Luận án đã mở rộng mô hình thống kê Hauser-Feshbach và mô hình Huizenga-Vandebosch bằng cách thay thế thông số cắt spin cố định bằng hàm phụ thuộc năng lượng kích thích $\sigma(E^*)$, đồng thời tích hợp ma trận chuyển dời gamma đa cực tính ($E1, M1, E2$) thay vì chỉ giả định chuyển dời $E1$ thuần túy.
2. Sự đổi mới phương pháp luận của luận án thể hiện như thế nào khi so sánh với các nghiên cứu trước đây?
So với các nghiên cứu của Danagulyan et al. [66] và Makwana et al. [67], luận án không sử dụng trực tiếp tiết diện đơn năng của TALYS để so sánh thô với dữ liệu bức xạ hãm, mà thiết lập một quy trình tính toán tích hợp khép kín: Sử dụng GEANT4 mô phỏng $500 \times 10^6$ hạt sơ cấp để tái tạo chính xác phổ bức xạ hãm thực tế $\phi(E_\gamma)$, sau đó dùng mã CERN ROOT để tích phân chập với tiết diện vi sai $\sigma(E)$ của TALYS 1.95 qua 48 tổ hợp mô hình mật độ mức và hàm lực gamma. Về mặt thực nghiệm, luận án áp dụng đồng thời cả hiệu chỉnh tự hấp thụ ($F_{abs}$) và trùng phùng ngẫu nhiên ($F_{sum}$), triệt tiêu hoàn toàn sai số hệ thống trong phép đo phổ HPGe.
3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu là gì?
Phát hiện bất ngờ nhất là sự vượt trội của mô hình hiện tượng học khí Fermi dịch chuyển lùi (BSFG) so với các mô hình vi mô hiện đại phức tạp (như HFB Goriely hay HFB Hilaire) trong việc mô tả tỷ số đồng phân của $^{152m1,m2}\text{Eu}$ và $^{195m,g}\text{Hg}$. Mặc dù các mô hình vi mô HFB có tính tổng quát cao hơn, nhưng ở vùng hạt nhân có độ biến dạng mạnh hoặc vùng chuyển tiếp giữa cầu và biến dạng ($A \approx 190 - 200$), mô hình BSFG với các tham số dịch chuyển năng lượng ghép đôi $\Delta$ được chuẩn hóa thực nghiệm lại cho kết quả khớp với thực nghiệm vượt trội (độ lệch $< 8.5%$).
4. Luận án có cung cấp đầy đủ quy trình lặp lại thí nghiệm (Replication Protocol) không?
Có. Luận án cung cấp chi tiết toàn bộ thông số thực nghiệm: Thông số vận hành máy gia tốc MT-25 Microtron (năng lượng electron, dòng chùm $10 - 20\ \mu\text{A}$, góc phát xạ), cấu tạo bia vonfram làm mát nước, chiều dày và độ tinh khiết của các lá mẫu kim loại, chiều dày lá lọc Cd ($0.5\text{ mm}$), hình học bố trí đầu dò HPGe Canberra, các hàm số hiệu chỉnh hiệu suất $\epsilon(E_\gamma)$, toàn bộ mã nguồn mô phỏng GEANT4 (Phụ lục A), file cấu hình đầu vào TALYS 1.95 (Phụ lục B), và mã phân tích CERN ROOT (Phụ lục C). Bất kỳ phòng thí nghiệm hạt nhân nào có trang bị máy gia tốc tương đương đều có thể lặp lại chính xác quy trình này.
5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào từ kết quả luận án?
Lộ trình 10 năm mở ra ba hướng đi chiến lược:
- Chuyển dịch từ nguồn bức xạ hãm sang nguồn chùm photon đơn sắc phân cực Laser Compton Backscattering (LCS) tại các trung tâm gia tốc thế hệ mới (ELI-NP, NewSUBARU) để đo phổ vi sai cộng hưởng hẹp.
- Ứng dụng phương pháp kích hoạt và mô hình TALYS để nghiên cứu các đồng phân siêu ngắn (thời gian sống micro giây đến nano giây) bằng kỹ thuật đo phổ gamma trực tuyến (on-line spectroscopy).
- Phát triển các cảm biến lượng tử và đồng hồ hạt nhân dựa trên các trạng thái đồng phân có mức năng lượng cực thấp (như $^{229m}\text{Th}$), thúc đẩy các ứng dụng đo lường chuẩn xác cao trong vật lý cơ bản và định vị dẫn đường vũ trụ.
Kết luận
Luận án tiến sĩ của Bùi Minh Huệ là một công trình nghiên cứu khoa học xuất sắc, có tính hệ thống cao và đóng góp những giá trị mang tính đột phá cho chuyên ngành Vật lý nguyên tử và hạt nhân. Tóm tắt các đóng góp cốt lõi:
- Đo đạc thành công lần đầu tiên trên thế giới giá trị IR của $^{195m,g}\text{Hg}$ ($14 - 24\text{ MeV}$), $^{197m,g}\text{Hg}$ ($18 - 24\text{ MeV}$), và $^{152m1,m2}\text{Eu}$ ($19, 21, 23\text{ MeV}$) trong phản ứng $(\gamma, n)$, cũng như $^{109m,g}\text{Pd}$, $^{111m,g}\text{Pd}$, $^{115m,g}\text{Cd}$, $^{117m,g}\text{Cd}$ trong phản ứng bắt neutron hỗn hợp nhiệt - cộng hưởng và cộng hưởng thuần túy.
- Làm sáng tỏ định lượng 4 hiệu ứng hạt nhân then chốt: Hiệu ứng năng lượng kích thích (xu hướng bão hòa ở đuôi GDR), hiệu ứng chênh lệch spin $\Delta J$ (suy giảm theo hàm mũ), hiệu ứng kênh phản ứng (ưu thế tuyệt đối của xung lượng truyền trong kênh $(n, 2n)$), và hiệu ứng cấu hình nucleon (sự phụ thuộc vào số khối và số proton trong chuỗi đồng vị/đồng neutron).
- Thiết lập thành công khung phương pháp luận mô phỏng kết hợp: Tích hợp chặt chẽ phổ bức xạ hãm GEANT4 ($500 \times 10^6$ hạt sơ cấp) với tiết diện vi mô TALYS 1.95 và mã phân tích CERN ROOT, khắc phục hoàn toàn nhược điểm của chùm photon không đơn năng.
- Đánh giá toàn diện 48 tổ hợp mô hình lý thuyết: Xác định tổ hợp tối ưu BSFG + GLO mô tả chính xác nhất dữ liệu thực nghiệm, cung cấp chuẩn mực tham chiếu cho việc hiệu chỉnh cơ sở dữ liệu hạt nhân quốc tế RIPL-3 và EXFOR (IAEA).
- Mở ra các hướng ứng dụng thực tiễn quan trọng: Đặt nền móng dữ liệu vững chắc cho công nghệ sản xuất đồng vị phóng xạ y học, vật lý lò phản ứng hạt nhân thế hệ mới, xử lý biến đổi rác thải phóng xạ, và nghiên cứu cấu trúc hạt nhân trong điều kiện kích thích cao.
