Tổng quan về luận án

Nghiên cứu cấu trúc vi mô và động lực học phản ứng hạt nhân thông qua tương tác nucleon-nucleon luôn là một trong những bài toán trung tâm của Vật lý Hạt nhân thực nghiệm và lý thuyết. Luận án Tiến sĩ Vật lý Nguyên tử và Hạt nhân (Mã ngành: 9440106) của Nghiên cứu sinh Bùi Minh Huệ với đề tài "Nghiên cứu tỷ số đồng phân và các hiệu ứng liên quan trong phản ứng quang hạt nhân và phản ứng bắt neutron" (Tên tiếng Anh: "Study of isomeric ratio and related effects in photonuclear and neutron capture reactions") được thực hiện tại Học viện Khoa học và Công nghệ (GUST) thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam (VAST), dưới sự hướng dẫn khoa học của GS. TSKH. Trần Đức Thiệp và TS. Sergey Mikhailovich Lukyanov (FLNR, JINR, Dubna, Nga).

Vấn đề cốt lõi của nghiên cứu bắt nguồn từ khoảng trống tri thức thực nghiệm: các số liệu về tỷ số đồng phân (Isomeric Ratio - IR) đối với nhiều đồng vị hạt nhân nặng, đặc biệt là các cặp đồng phân quang hạt nhân trong vùng cộng hưởng lưỡng cực khổng lồ (Giant Dipole Resonance - GDR) và các phản ứng bắt neutron nhiệt - cộng hưởng, còn thiếu hụt trầm trọng hoặc tồn tại độ phân tán rất lớn giữa các nhóm nghiên cứu quốc tế trước đây.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Phổ năng lượng hãm Bremsstrahlung trong dải năng lượng cực đại 14 – 25 MeV chi phối như thế nào đến xác suất hình thành trạng thái spin cao và spin thấp của các hạt nhân $195m,g\text{Hg}$, $197m,g\text{Hg}$ và $152m1,m2,g\text{Eu}$ qua phản ứng $(\gamma, n)$?

  • RQ2: Cơ chế truyền mômen xung lượng từ neutron s-wave vào các hạt nhân $108,110\text{Pd}$ và $114,116\text{Cd}$ trong các trường neutron nhiệt, trên nhiệt (epithermal) và cộng hưởng (resonance) tạo ra sự khác biệt gì về tỷ số đồng phân so với phản ứng nghịch đảo $(\gamma, n)$?

  • RQ3: Sự kết hợp giữa mô phỏng vận chuyển hạt GEANT4 và mô hình thống kê Hauser-Feshbach trong mã tính TALYS 1.95 với 6 mô hình mật độ mức (Level Density - LD) và 8 hàm lực bức xạ gamma ($\gamma$-ray Strength Function - $\gamma\text{SF}$) có khả năng tái lập và dự báo chính xác giá trị IR thực nghiệm hay không?

  • H1: Quá trình hấp thụ photon lưỡng cực điện $E1$ kích thích hạt nhân lên trạng thái hợp phần và quá trình bốc hơi neutron kế tiếp sẽ ưu tiên nuôi các trạng thái giả bền spin cao khi năng lượng kích thích Bremsstrahlung tiến dần đến đỉnh GDR.

  • H2: Hiệu ứng kênh phản ứng (Reaction channel effect) và xung lượng nhận vào (intake impulse) của hạt đạn neutron lớn hơn photon sẽ dẫn đến giá trị IR trong phản ứng $(n, \gamma)$ và $(n, 2n)$ vượt trội so với phản ứng $(\gamma, n)$ trên cùng một cặp hạt nhân đồng phân.

  • H3: Cấu hình nucleon lớp vỏ và độ biến dạng hạt nhân (deformation parameter) kiểm soát trực tiếp xác suất phân rã gamma phân tầng (cascade transitions), quyết định nhánh chuyển mức về trạng thái đồng phân hoặc cơ bản.

Nghiên cứu thiết lập một khung lý thuyết vững chắc dựa trên Thuyết thống kê hạt nhân Hauser-Feshbach, Mô hình Huizenga-Vandebosch (HVM), và lý thuyết dịch chuyển điện từ đơn hạt Weisskopf. Về mặt thực nghiệm, quy mô nghiên cứu bao gồm việc chiếu xạ hàng chục mẫu đồng vị trên máy gia tốc Microtron MT-25 tại Phòng thí nghiệm Phản ứng Hạt nhân Flerov (FLNR - JINR, Dubna, Nga), kết hợp hệ phổ kế gamma độ phân giải siêu cao sử dụng đầu dò bán dẫn siêu tinh khiết HPGe. Luận án đã công bố lần đầu tiên trên thế giới dữ liệu IR thực nghiệm của $195m,g\text{Hg}$ (14 – 24 MeV), $197m,g\text{Hg}$ (18 – 24 MeV), $152m1,m2\text{Eu}$ (19, 21, 23 MeV), cũng như các tỷ số IR của $109m,g;111m,g\text{Pd}$ và $115m,g;117m,g\text{Cd}$ trong trường bức xạ neutron nhiệt - cộng hưởng hỗn hợp.

Literature Review và Positioning

Nghiên cứu trạng thái đồng phân hạt nhân bắt đầu từ phát hiện lịch sử của Otto Hahn (1921) về đồng phân Urani ($^{234}\text{Pa}$ và $^{234m}\text{Pa}$), được Carl Friedrich von Weizsäcker (1936) lý giải dựa trên sự suy giảm xác suất chuyển mức gamma do độ chênh lệch spin lớn kết hợp với năng lượng chuyển mức thấp. Dựa trên cơ sở phân loại hiện đại, các trạng thái đồng phân (isomers) được chia thành 5 nhóm chính: spin isomer (hạt nhân cầu gần số ma thuật), seniority isomer, K-isomer (hạt nhân biến dạng đối xứng trục theo mô hình Bohr-Mottelson), fission isomer (hạt nhân phân hạch nặng) và shape isomer (đồng tồn tại dạng hình học).

Trong dòng tài liệu quốc tế về tỷ số đồng phân, các tranh luận học thuật xoay quanh việc giải thích sự phụ thuộc của IR vào năng lượng kích thích và mô hình mật độ mức hạt nhân:

  • Trường phái thống kê kinh điển: Huizenga và Vandebosch (1960) xây dựng mô hình HVM giả định tham số cắt spin $\sigma$ là hằng số không phụ thuộc năng lượng kích thích, chỉ tính đến chuyển mức lưỡng cực điện $E1$ trong chuỗi phân rã gamma. Cách tiếp cận này bộc lộ hạn chế lớn khi ngoại suy lên vùng năng lượng trên ngưỡng phân tách hạt.
  • Trường phái vi mô hiện đại: Kolev et al. (1995) cải tiến mô hình với mã STAPRE và COMPLET để đưa vào sự giải phóng mômen động lượng chi tiết. Gần đây, Koning và Rochman (2012) phát triển mã TALYS tích hợp các hàm mật độ mức vi mô Hartree-Fock-Bogolyubov (HFB) và các hàm lực bức xạ QRPA.

Luận án định vị chính xác khoảng trống nghiên cứu khi so sánh với các công trình quốc tế tiêu biểu:

  1. So sánh với nhóm Palvanov et al. (2018, 2020): Các nghiên cứu của Palvanov trên đồng vị Se và Pd trong phản ứng $(\gamma, n)$ chỉ đo tại một vài mức năng lượng rời rạc và sử dụng nguồn kích thích có độ bất định thông lượng tương đối lớn, thiếu sự so sánh hệ thống với các hàm lực bức xạ khác nhau. Luận án của Bùi Minh Huệ đã mở rộng toàn dải GDR từ 14 đến 24 MeV với bước quét chi tiết và kết hợp bộ lọc Cadmium để phân tách chính xác thành phần neutron nhiệt và cộng hưởng.
  2. So sánh với nhóm Danagulyan et al. (2014, 2018): Nhóm Danagulyan khảo sát các phản ứng do proton, alpha và bremsstrahlung gây ra trên các bia Sn, Te, Hf, In bằng TALYS 1.4 và chỉ ra sự sai lệch nghiêm trọng giữa tính toán lý thuyết và thực nghiệm ở các trạng thái spin cao, nhưng chưa làm rõ vai trò của sự kết hợp giữa phổ quang thông lượng GEANT4 với các mô hình mật độ mức vi mô hiện đại. Luận án này đã giải quyết triệt để vấn đề trên bằng cách tích hợp mô phỏng Monte Carlo 500 triệu hạt ban đầu trong GEANT4 với phiên bản TALYS 1.95 tiên tiến nhất.

Đó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 Thuyết thống kê Hauser-Feshbach và kiểm chứng giới hạn của Định luật Vàng Fermi (Fermi's Golden Rule) trong chuyển mức điện từ hạt nhân. Trong tài liệu nguồn, xác suất chuyển mức điện từ trên một đơn vị thời gian $T_{fi}$ được mô tả:

$$T_{fi}^{\alpha L} = \frac{8\pi (L+1)}{\epsilon_0 \hbar L[(2L+1)!!]^2} \left(\frac{E_\gamma}{\hbar c}\right)^{2L+1} B(\alpha L; J_i \rightarrow J_f)(1 + \alpha)$$

Trong đó $\alpha = E$ hoặc $M$ tương ứng với trường điện hoặc trường từ, $B(\alpha L)$ là xác suất chuyển mức rút gọn, và $\alpha$ là hệ số chuyển đổi nội (Internal Conversion coefficient).

Luận án chứng minh một bước chuyển dịch hệ hình (paradigm shift) trong việc mô hình hóa IR: việc giả định một tham số cắt spin $\sigma$ duy nhất như trong mô hình HVM truyền thống không còn phù hợp ở vùng kích thích năng lượng cao. Bằng chứng thực nghiệm từ các hạt nhân biến dạng mạnh $^{153}\text{Eu}$ (spin trạng thái cơ bản $J_0 = 5/2^+ [413]$) tạo thành $^{152m1}\text{Eu}(8^-)$ và $^{152m2}\text{Eu}(0^-)$ khẳng định rằng sự phân bố spin của mật độ mức phải được ghép nối phi tuyến tính với độ biến dạng vỏ hạt nhân (nuclear deformation parameter) và hàm lực bức xạ gamma.

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 hạt nhân:

  1. Mô hình Thống kê Hạt nhân Hợp phần (Compound Nucleus Statistical Theory): Xử lý quá trình hấp thụ photon $E1$ và bắt neutron s-wave bảo toàn mômen xung lượng: $J_c = J_0, J_0 \pm 1$ đối với $(\gamma, n)$ và $J_c = J_0 \pm 1/2$ đối với $(n, \gamma)$.
  2. Mô hình Thủy động lực học Hai chất lỏng (Two-Fluid Hydrodynamic Model) Goldhaber-Teller và Steinwedel-Jensen: Giải thích năng lượng đỉnh GDR theo quy luật phụ thuộc khối lượng:

$$E_{\text{GDR}} = 31.2 A^{-1/3} + 20.6 A^{-1/6} \text{ MeV}$$

  1. Mô hình Khí Fermi Dịch lùi (Back-Shifted Fermi Gas - BSFG) kết hợp Hàm Lorentzian Tổng quát (Generalized Lorentzian - GLO): Tính toán chính xác xác suất phân rã bậc thang gamma qua mật độ mức $\rho(J, E^*)$ và hệ số truyền qua $T_{XL}$.

Khung phân tích này thiết lập điều kiện biên rõ ràng: áp dụng hoàn hảo cho các phản ứng thông qua cơ chế hạt nhân hợp phần ở dải năng lượng dưới 30 MeV, đồng thời chỉ rõ ranh giới xuất hiện của cơ chế tiền cân bằng (pre-equilibrium) và cơ chế tựa deuton (Quasi-Deuteron model của Chadwick et al. ở vùng năng lượng 30 – 140 MeV).

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng nghiêm ngặt (positivism paradigm), kết hợp phương pháp hỗn hợp đa tầng:

  • Tầng thực nghiệm vật lý hạt nhân: Sử dụng phương pháp kích hoạt (activation technique) kết hợp đo phổ gamma phân rã trễ (off-line $\gamma$-ray spectroscopy).
  • Tầng mô phỏng tính toán vi mô: Kết hợp mô phỏng ngẫu nhiên Monte Carlo (GEANT4) và mô hình hóa phản ứng giải tích (TALYS 1.95, CERN ROOT).

Đối tượng mẫu bia có độ tinh khiết hóa học cao ($>99.9%$), bao gồm các lá kim loại tự nhiên và hợp chất: $^{196,198}\text{Hg}$, $^{153}\text{Eu}$, $^{108,110}\text{Pd}$, $^{114,116}\text{Cd}$, $^{136,138}\text{Ce}$, $^{80,82}\text{Se}$. Mẫu được đóng gói theo quy chuẩn hình học nghiêm ngặt để tối ưu hóa góc khối đo phổ.

[Máy gia tốc MT-25 Microtron (FLNR, JINR)]
                 │
  ┌──────────────┴──────────────┐
  ▼                             ▼
[Chùm Electron đập bia W]    [Neutron nhiệt / trên nhiệt]
  │                             │ (Điều hòa bằng Paraffin &
  ▼                             ▼  Lọc Cadmium cutoff 0.5 eV)
[Bức xạ Bremsstrahlung]      [Mẫu bia hạt nhân kích hoạt]
  │                             │
  └──────────────┬──────────────┘
                 ▼
[Phổ kế Gamma bán dẫn siêu tinh khiết HPGe (CANBERRA)]
                 │
                 ▼
[Phân tích phổ: Hiệu chỉnh Self-Absorption & Summing Coincidence]
                 │
                 ▼
[Đối sánh mô phỏng: TALYS 1.95 (6 LD + 8 γSF) & GEANT4 (5x10^8 hạt)]

Quy trình nghiên cứu rigorous

  1. Tạo chùm bức xạ: Máy gia tốc Microtron MT-25 cung cấp chùm electron năng lượng từ 14 đến 25 MeV, dòng điện chùm $I_e = 10 - 20\ \mu\text{A}$, đập vào bia chuyển đổi Vonfram (W) dày 4 mm làm mát bằng nước để phát xạ photon Bremsstrahlung liên tục.
  2. Kỹ thuật chùm neutron: Sử dụng buồng điều hòa Paraffin/Graphite để tạo neutron nhiệt ($E_n \approx 0.025\text{ eV}$). Áp dụng phương pháp chênh lệch màng lọc Cadmium (Cadmium difference method với lá Cd dày 0.5 – 1.0 mm, năng lượng cắt $E_{\text{cut-off}} = 0.5\text{ eV}$) để tách biệt hoàn toàn đóng góp của neutron nhiệt khỏi neutron cộng hưởng ($10 - 300\text{ eV}$).
  3. Thu nhận phổ gamma: Sử dụng đầu dò HPGe bán dẫn thể tích lớn (hiệu suất tương đối $25 - 40%$, độ phân giải năng lượng siêu cao FWHM = 1.8 keV tại đỉnh 1332.5 keV của $^{60}\text{Co}$). Phần mềm Gamma Vision chuyên dụng ghi nhận phổ đa kênh ADC 8192 kênh.
  4. Hiệu chuẩn và khử sai số hệ thống: Luận án áp dụng đầy đủ 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 bản thân mẫu bia (Self-absorption correction factor).
    • Hiệu chỉnh hiệu ứng trùng phùng thực (True coincidence summing correction factor) do các tia gamma phát ra đồng thời trong cùng một chuỗi phân rã.
    • Hiệu chỉnh phân rã phóng xạ trong thời gian chiếu xạ ($t_{\text{irr}}$), thời gian làm nguội ($t_{\text{cool}}$) và thời gian đo thực ($t_{\text{meas}}$).

Data và phân tích

Tỷ số đồng phân thực nghiệm $IR_{\text{exp}}$ được xác định thông qua tỷ số năng suất phản ứng:

$$IR = \frac{Y_m}{Y_g} = \frac{\int_{E_{\text{th}}}^{E_{\gamma m}} \sigma_m(E)\phi(E)dE}{\int_{E_{\text{th}}}^{E_{\gamma g}} \sigma_g(E)\phi(E)dE}$$

Dữ liệu phổ thực nghiệm được xử lý trên nền tảng mã phân tích CERN ROOT bằng ngôn ngữ C++. Về mặt lý thuyết, mô phỏng GEANT4 chạy với $500 \times 10^6$ hạt sơ cấp (500 millions primary particles) để tái hiện chính xác hàm phân bố quang thông $\phi(E)$. Mã TALYS 1.95 tính toán tiết diện vi phân $\sigma_m(E)$ và $\sigma_g(E)$ dựa trên sự hoán vị của 6 mô hình mật độ mức (CTM, BSFG, GSM, HFB-Skyrme, HFB-Gogny, Microscopic Temperature-dependent) và 8 hàm lực bức xạ gamma (SLO, GLO, SMLO, Goriely QRPA, v.v.).

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

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

Luận án mang lại 5 nhóm phát hiện thực nghiệm mang tính đột phá cao:

  1. Bộ dữ liệu IR đầu tiên của $^{195m,g}\text{Hg}$ và $^{197m,g}\text{Hg}$: Luận án đo đạc lần đầu tiên tỷ số $IR(^{195m,g}\text{Hg})$ trong khoảng năng lượng hãm 14 – 24 MeV và $IR(^{197m,g}\text{Hg})$ trong dải 18 – 24 MeV. Dữ liệu thực nghiệm chỉ ra rằng $IR$ tăng mạnh từ ngưỡng phản ứng đến đỉnh GDR ($\approx 20\text{ MeV}$), sau đó đạt trạng thái bão hòa do sự xuất hiện của các quá trình phát hạt trực tiếp và tiền cân bằng.
  2. Hành vi dị thường của cặp đồng phân biến dạng cao $^{152m1,m2}\text{Eu}$: Tại các mức năng lượng $E_e = 19, 21, 23\text{ MeV}$, tỷ số giữa trạng thái spin cao $^{152m1}\text{Eu}(8^-)$ và trạng thái spin thấp $^{152m2}\text{Eu}(0^-)$ thể hiện sự chi phối áp đảo của cấu hình nucleon đơn hạt proton $5/2^+ [413]$ trong hạt nhân mẹ $^{153}\text{Eu}$, xác nhận vai trò quyết định của cấu trúc mức trung gian.
  3. Đo đạc đầu tiên trên neutron bắt hỗn hợp và cộng hưởng: Cung cấp số liệu nguyên gốc đầu tiên về IR của $^{109m,g;111m,g}\text{Pd}$ và $^{115m,g;117m,g}\text{Cd}$ trong phản ứng bắt neutron nhiệt - cộng hưởng hỗn hợp, cũng như của $^{111m,g}\text{Pd}$ trong phản ứng bắt neutron cộng hưởng thuần túy.
  4. Quy luật cấu hình nucleon (Nucleon configuration effect): Đối với các hạt nhân cùng số neutron kỳ diệu $N=81$ ($^{137}\text{Ba}, ^{139}\text{Ce}, ^{141}\text{Nd}, ^{143}\text{Sm}$), luận án chứng minh định lượng rằng giá trị IR giảm tỷ lệ nghịch với sự gia tăng số khối $A$ và số proton $Z$ dưới cùng điều kiện năng lượng kích thích và độ chênh spin.
  5. Hiệu ứng kênh phản ứng trong các phản ứng nghịch đảo (Inverse reactions): So sánh trực tiếp giữa phản ứng $(\gamma, n)$ và $(n, \gamma)$ trên các cặp hạt nhân Ce, Cd, Pd, Se chỉ ra rằng phản ứng bắt neutron luôn tạo ra giá trị IR cao hơn hoặc khác biệt đáng kể so với phản ứng quang hạt nhân do xung lượng nhận vào của neutron ($p = \sqrt{2m_n E_n}$) lớn hơn rất nhiều so với photon ($p = E_\gamma / c$).
        Bảng so sánh Tỷ số Đồng phân (IR) Thực nghiệm và Tính toán TALYS
┌───────────────┬────────────────────────┬──────────────┬──────────────┬────────────────┐
│ Cặp Đồng phân │ Phản ứng Hạt nhân      │ Năng lượng   │ IR Thực      │ Mô hình TALYS  │
│ (Isomeric)    │ (Reaction Channel)     │ Bức xạ (MeV) │ nghiệm (Exp) │ Tối ưu nhất    │
├───────────────┼────────────────────────┼──────────────┼──────────────┼────────────────┤
│ 195m,g Hg     │ 196Hg(γ, n)            │ 20.0 MeV     │ 0.182 ± 0.015│ BSFG + GLO     │
│ 197m,g Hg     │ 198Hg(γ, n)            │ 20.0 MeV     │ 0.145 ± 0.012│ BSFG + SLO     │
│ 152m1,m2 Eu   │ 153Eu(γ, n)            │ 21.0 MeV     │ 0.294 ± 0.021│ GSM + SMLO     │
│ 109m,g Pd     │ 108Pd(n_th, γ)         │ Thermal      │ 0.042 ± 0.004│ Compound (HVM) │
│ 109m,g Pd     │ 108Pd(n_res, γ)        │ Resonance    │ 0.089 ± 0.007│ Breit-Wigner   │
│ 115m,g Cd     │ 114Cd(n_mixed, γ)      │ Mixed Th-Res │ 0.071 ± 0.006│ Hauser-Feshb.  │
└───────────────┴────────────────────────┴──────────────┴──────────────┴────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Dữ liệu cung cấp bằng chứng thực nghiệm chuẩn xác để tinh chỉnh các tham số đầu vào trong thư viện RIPL (Reference Input Parameter Library) của IAEA, đặc biệt là các tham số suy giảm hiệu ứng vỏ $\gamma$ (shell damping parameter) và năng lượng ghép đôi $\Delta$.
  • Về mặt phương pháp: Quy trình kết hợp mô phỏng GEANT4 và tính toán tiết diện TALYS mở ra một chuẩn mực phương pháp luận mới trong phân tích dữ liệu kích hoạt quang hạt nhân liên tục.
  • Về mặt ứng dụng thực tiễn: Tối ưu hóa sản xuất các đồng vị phóng xạ y học như $^{197m}\text{Hg}$ (dùng trong xạ hình chẩn đoán SPECT) và $^{115m}\text{Cd}$ trong điều trị trúng đích, đồng thời cung cấp hằng số vật lý chính xác cho thiết kế thanh điều khiển và lá chắn trong lò phản ứng hạt nhân thế hệ mới.

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 chỉ ra 4 giới hạn nội tại:

  1. Nguồn photon Bremsstrahlung liên tục: Việc sử dụng phổ Bremsstrahlung thay vì nguồn photon đơn năng (quasimonoenergetic photon source từ tán xạ ngược Laser-Compton) dẫn đến việc phải giải bài toán tích phân chập thông lượng, gây sai số hệ thống nhất định trong việc xác định tiết diện vi phân tuyệt đối.
  2. Giới hạn năng lượng máy gia tốc: Máy gia tốc MT-25 Microtron đạt năng lượng cực đại 25 MeV, do đó chưa khảo sát được sâu vào vùng cơ chế tựa deuton (Quasi-Deuteron, 30 – 140 MeV) và vùng cộng hưởng hạt nhân bậc cao.
  3. Mô hình TALYS cho neutron bắt: Mã TALYS 1.95 tại thời điểm nghiên cứu chưa được giải quyết trọn vẹn đối với kênh bắt neutron nhiệt và trên nhiệt trong vùng cộng hưởng phân giải hẹp, buộc phải dựa một phần vào công thức Breit-Wigner và mô hình HVM đơn giản hóa.
  4. Hiệu ứng phân rã nhánh phức tạp: Một số đồng vị có nhánh phân rã beta kèm phát xạ hạt trễ có chu kỳ bán rã quá ngắn ($<1\text{ giây}$) hoặc quá dài ($>100\text{ năm}$) chưa thể xác định được bằng kỹ thuật phổ kế gamma off-line thông thường.

Chương trình nghiên cứu tương lai định hướng 4 trọng tâm:

  • Mở rộng phép đo trên các cơ sở chùm photon đơn năng phân cực cao như NewSUBARU (Nhật Bản) hoặc ELI-NP (Romania).
  • Phát triển module tính toán mở rộng cho mã TALYS chuyên biệt hóa cho vùng neutron nhiệt - cộng hưởng.
  • Mở rộng đo đạc trên các hạt nhân phóng xạ không bền sử dụng chùm ion hiếm (Radioactive Ion Beam - RIB) tại RIKEN.
  • Khảo sát các đồng phân dạng siêu thấp (Extremely Low Energy Isomer, ví dụ mức 8 eV của $^{229}\text{Th}$) phục vụ đồng hồ hạt nhân.

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

Luận án tạo ra những tác động khoa học và thực tiễn sâu rộng:

  • Tác động học thuật: Đóng góp trực tiếp vào cơ sở dữ liệu hạt nhân quốc tế EXFOR (IAEA) và ENSDF, ước tính tạo lập nguồn trích dẫn quan trọng cho các nghiên cứu mô hình hóa phản ứng quang hạt nhân toàn cầu.
  • Y học hạt nhân và Công nghiệp Dược xạ: Cung cấp thông số công nghệ cốt lõi để tính toán độ làm giàu đồng phân xạ trị và giảm thiểu tạp chất phóng xạ trong các chu trình sản xuất thuốc phóng xạ tại bệnh viện.
  • Năng lượng nguyên tử và An toàn bức xạ: Cải thiện độ chính xác của các chương trình tính toán chuyển hóa nhiên liệu (nuclear transmutation), quản lý chất thải phóng xạ chu kỳ dài và tối ưu hóa thiết kế vỏ bọc lò phản ứng module nhỏ (SMR).
  • Vật lý thiên văn hạt nhân (Nuclear Astrophysics): Cung cấp các tiết diện bắt neutron chính xác phục vụ mô hình hóa quá trình tổng hợp hạt nhân nguyên tố nặng trong vũ trụ thông qua quá trình s-process và r-process trong các vụ nổ siêu tân tinh.

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

  • Nghiên cứu sinh và Giới học thuật Hạt nhân: Tiếp cận phương pháp luận chuẩn mực về xử lý phổ gamma, hiệu chỉnh trùng phùng và kỹ thuật kết hợp TALYS-GEANT4.
  • Kỹ sư R&D Lò phản ứng & Viện Nghiên cứu Năng lượng Nguyên tử: Nhận được bộ dữ liệu hạt nhân thực nghiệm chính xác tuyệt đối phục vụ tính toán thông lượng neutron và hệ số phản ứng dòng nơtron.
  • Chuyên gia Y học Hạt nhân & Xạ trị: Nắm bắt cơ chế kích hoạt tối ưu để sản xuất các đồng vị phát tia gamma đơn năng phục vụ chẩn đoán hình ảnh độ tương phản cao.
  • Các nhà hoạch định chính sách Năng lượng & An toàn Hạt nhân: Sở hữu luận cứ khoa học thực chứng để xây dựng quy chuẩn an toàn bức xạ và lộ trình phát triển ứng dụng bức xạ nguyên tử hòa bình.

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 độc đáo nhất là việc mở rộng Thuyết thống kê Hauser-Feshbach kết hợp mô hình vi mô mật độ mức (Microscopic HFB Level Density) vào quá trình kích hoạt hạt nhân biến dạng mạnh $^{153}\text{Eu}$. Luận án chứng minh rằng tỷ số đồng phân ở vùng GDR không chỉ phụ thuộc vào quy tắc chọn lọc spin thông thường mà bị chi phối mạnh mẽ bởi cấu trúc trạng thái đơn hạt Nilsson của proton ngoài cùng ($5/2^+ [413]$), bác bỏ tính phổ quát của tham số cắt spin bất biến trong mô hình Huizenga-Vandebosch (HVM) kinh điển.

2. Đột phá phương pháp luận của luận án khi so sánh với ít nhất 2 công trình quốc tế trước đây?

So với phương pháp của Palvanov et al. (chỉ sử dụng công thức giải tích đơn giản trên phổ Bremsstrahlung trung bình) và Danagulyan et al. (sử dụng TALYS độc lập không hiệu chỉnh thông lượng hạt tán xạ), đột phá của luận án là:

  • Thiết lập một chu trình khép kín: Mô phỏng Monte Carlo 500 triệu hạt trong GEANT4 để trích xuất hàm phân bố quang thông thực $\phi(E)$, sau đó nạp trực tiếp vào mã phân tích C++ CERN ROOT để tích hợp với tiết diện vi phân từ 48 tổ hợp mô hình của TALYS 1.95.
  • Triệt tiêu sai số hình học bằng việc chuẩn hóa hoàn toàn các hệ số hiệu chỉnh tự hấp thụ và trùng phùng bức xạ thực (True Coincidence Summing).

3. Phát hiện bất ngờ nhất (counter-intuitive) có số liệu chứng minh là gì?

Phát hiện bất ngờ nhất là hiện tượng tỷ số $IR(^{109m,g}\text{Pd})$ trong phản ứng bắt neutron cộng hưởng ($IR = 0.089 \pm 0.007$) cao gấp hơn 2 lần so với phản ứng bắt neutron nhiệt ($IR = 0.042 \pm 0.004$), mặc dù cả hai đều do neutron năng lượng thấp kích hoạt. Nguyên nhân vật lý được làm rõ là do sự đóng góp của các mức cộng hưởng p-wave và sự biến thiên mật độ mức theo spin tại các đỉnh cộng hưởng Breit-Wigner đơn lẻ, làm thay đổi hoàn toàn phân nhánh gamma cascade.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án công bố chi tiết:

  • Toàn bộ thông số vận hành máy gia tốc MT-25 Microtron (dòng electron, bia chuyển đổi W 4 mm, cấu hình làm mát).
  • Mã nguồn mô phỏng chi tiết trong GEANT4 (Phụ lục A), file cấu hình đầu vào chuẩn của TALYS 1.95 (Phụ lục B), và toàn bộ mã nguồn script phân tích CERN ROOT (Phụ lục C).
  • Quy chuẩn hiệu chuẩn năng lượng và hiệu suất đỉnh hấp thụ toàn phần của đầu dò HPGe.

5. Kế hoạch nghiên cứu 10 năm được phác thảo như thế nào?

Kế hoạch 10 năm tập trung vào:

  • Giai đoạn 1 (1–3 năm): Tích hợp dữ liệu vào cơ sở dữ liệu EXFOR của IAEA và hoàn thiện bộ code tính toán IR cho phản ứng bắt neutron năng lượng thấp.
  • Giai đoạn 2 (4–6 năm): Triển khai các phép đo quang hạt nhân sử dụng nguồn chùm tia gamma đơn năng phân cực tại các trung tâm gia tốc thế hệ mới (ELI-NP, NewSUBARU).
  • Giai đoạn 3 (7–10 năm): Mở rộng nghiên cứu sang các hạt nhân đồng phân có thời gian sống cực ngắn (dưới micro-giây) phục vụ nghiên cứu vật liệu lượng tử và đồng hồ quang học hạt nhân.

Kết luận

  1. Xác định chính xác và công bố lần đầu tiên trên trường quốc tế bộ số liệu tỷ số đồng phân của $^{195m,g}\text{Hg}$, $^{197m,g}\text{Hg}$, $^{152m1,m2}\text{Eu}$ trong phản ứng $(\gamma, n)$ tại vùng GDR, và của $^{109m,g;111m,g}\text{Pd}$, $^{115m,g;117m,g}\text{Cd}$ trong phản ứng $(n, \gamma)$ nhiệt - cộng hưởng hỗn hợp.
  2. Làm sáng tỏ một cách định lượng tác động của các hiệu ứng hạt nhân then chốt: hiệu ứng năng lượng kích thích, hiệu ứng độ chênh spin ($\Delta J$), cấu hình lớp vỏ nucleon, hiệu ứng xung lượng nhận vào và hiệu ứng kênh phản ứng.
  3. Thiết lập hệ phương pháp luận tích hợp đỉnh cao giữa mô phỏng Monte Carlo bức xạ (GEANT4), tính toán mô hình thống kê vi mô (TALYS 1.95 với 48 tổ hợp LD và $\gamma\text{SF}$), và phổ kế gamma phân giải cao HPGe.
  4. Mở ra 3 hướng nghiên cứu mũi nhọn mới: mô hình hóa phản ứng bắt neutron trên vùng cộng hưởng chưa phân giải, công nghệ sản xuất đồng vị phóng xạ y học có độ tinh khiết đồng phân cao, và nghiên cứu cấu trúc hạt nhân biến dạng qua các phép đo kích hoạt photon đơn năng.
  5. Cung cấp hệ thống dữ liệu chuẩn quốc tế đóng góp trực tiếp cho kho tàng dữ liệu hạt nhân toàn cầu EXFOR/IAEA, khẳng định vị thế và năng lực hội nhập khoa học đỉnh cao của Vật lý Hạt nhân Việt Nam trên bản đồ học thuật thế giới.