Tổng quan về luận án

Công nghệ xử lý bức xạ công nghiệp bằng chùm tia điện tử (electron beam - e-beam) và tia X chuyển đổi đang đóng vai trò then chốt trong khử trùng thiết bị y tế, thanh trùng thực phẩm, biến tính polyme và xử lý chất thải môi trường theo tiêu chuẩn IAEA TRS-481 và ICRP-103. Tuy nhiên, các mô hình tính toán liều lượng truyền thống thường giả định chùm electron là đơn năng và môi trường chiếu xạ là khối vật chất đồng nhất. Luận án tiến sĩ vật lý nguyên tử và hạt nhân của nghiên cứu sinh Nguyễn Anh Tuấn (chuyên ngành Vật lý nguyên tử và hạt nhân, mã số 62 44 05 01) với đề tài "Nghiên cứu sự phân bố liều theo độ sâu sản phẩm chiếu xạ và hiệu ứng phát bức xạ hãm của chùm electron từ máy gia tốc RF-LINAC và LWFA" đã giải quyết khoảng trống học thuật cốt lõi: sự biến dạng phân bố liều trong các môi trường sản phẩm không đồng nhất thực tế và nguy cơ kích hoạt phản ứng quang hạt nhân $(\gamma, n)$ sinh neutron thứ cấp từ đuôi phổ năng lượng cao của máy gia tốc.

Nghiên cứu được thiết lập dựa trên bốn câu hỏi nghiên cứu và giả thuyết khoa học chính:

  1. Câu hỏi 1 (Q1): Mức độ không đồng nhất về mặt hình học và cấu trúc xếp lớp/cột của sản phẩm thực tế làm sai lệch đường phân bố liều theo độ sâu so với mô hình đồng nhất chuẩn như thế nào? Giả thuyết 1 (H1): Các khe hở không khí gây ra hiệu ứng tán xạ cạnh và rò rỉ chùm tia, làm dịch chuyển đỉnh liều cực đại ($D_{\max}$) và thay đổi tỷ số đồng đều liều (DUR - Dose Uniformity Ratio).
  2. Câu hỏi 2 (Q2): Phổ năng lượng thực nghiệm với độ rộng đỉnh tán xạ của máy gia tốc RF-LINAC ($9,92 \pm 0,50\text{ MeV}$) và máy gia tốc laser-plasma LWFA ($3,5 - 30\text{ MeV}$) ảnh hưởng ra sao đến khả năng xuyên sâu và hiệu suất sử dụng chùm tia? Giả thuyết 2 (H2): Phổ năng lượng mở rộng làm biến đổi độ dốc suy giảm liều và tăng tỷ lệ liều ở vùng đuôi so với chùm đơn năng lý thuyết.
  3. Câu hỏi 3 (Q3): Hiệu suất chuyển đổi electron - photon ($k_{\text{conv}}$) và phân bố góc của bức xạ hãm (Bremsstrahlung) sau các bia vật liệu biến thiên như thế nào theo năng lượng và số khối bia? Giả thuyết 3 (H3): Cường độ bức xạ hãm tập trung cao ở hướng phát xạ $0^\circ$ và tỷ lệ photon vượt ngưỡng năng lượng phản ứng quang hạt nhân tăng mạnh khi sử dụng bia kim loại nặng ($Z$ cao).
  4. Câu hỏi 4 (Q4): Khả năng sinh neutron từ phản ứng $(\gamma, n)$ trong các cấu kiện máy gia tốc và vật liệu chiếu xạ ở mức năng lượng danh định 10 MeV có vượt ngưỡng an toàn bức xạ hay không? Giả thuyết 4 (H4): Các đồng vị có ngưỡng năng lượng thấp dưới 10 MeV như $^{207}\text{Pb}$ ($6,74\text{ MeV}$) hay $^{57}\text{Fe}$ ($7,65\text{ MeV}$) sẽ bị kích hoạt bởi thành phần năng lượng electron vượt mức 10 MeV trong phổ tán xạ.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết dừng electron Bethe-Bloch, tán xạ đàn hồi electron-electron Møller, cơ chế phát bức xạ hãm trong trường Coulomb hạt nhân và lý thuyết gia tốc trường dòng đuôi plasma của Tajima & Dawson (1979). Nghiên cứu khảo sát thực nghiệm trên máy gia tốc công nghiệp RF-LINAC UELR-10-15S2 (công suất $15\text{ kW}$, năng lượng $9,92 \pm 0,50\text{ MeV}$) tại Trung tâm Nghiên cứu và Triển khai Công nghệ Bức xạ (VINATOM), kết hợp liều kế màng mỏng chuẩn B3000 (dải đo $0,5 - 35\text{ kGy}$) và mô phỏng Monte Carlo N-Particle (MCNP phiên bản 4c2). Kết quả cung cấp các thông số định lượng chính xác về bề dày tối ưu ($T_{\text{opt1}}, T_{\text{opt2}}$) và giới hạn mật độ mặt ($\text{ADL}_1, \text{ADL}_2$), tạo nền tảng đột phá cho việc thiết kế quy trình chiếu xạ công nghiệp an toàn và hiệu quả cao.

Literature Review và Positioning

Nghiên cứu tương tác chùm electron năng lượng cao với vật chất bắt đầu từ những năm 1970 khi McLaughlin và cộng sự (1979) sử dụng liều kế phim đo phân bố liều trong nhôm, carbon, polystyrene và polyethylene ở mức 10 MeV, chỉ ra rằng độ sâu đạt liều cực đại phụ thuộc tỷ lệ nghịch vào mật độ khối của vật liệu. Đến năm 2005, Miller R. đã công bố công trình chuẩn hóa về chiếu xạ thực phẩm bằng e-beam 10 MeV và tia X 5 MeV, chứng minh độ xuyên sâu cực đại trong nước ($\rho = 1,0\text{ g/cm}^3$) đạt $5,43\text{ cm}$, mật độ mặt tối ưu cho chiếu xạ một mặt là $3,8\text{ g/cm}^2$ với DUR = 1,35 và hiệu suất 70%, đồng thời cho chiếu xạ hai mặt là $8,6\text{ g/cm}^2$. Miller cũng phát hiện hiện tượng liều vùng tiếp giáp giữa các ống trụ tròn cao gấp đôi bề mặt, đặt ra bài toán che chắn không gian.

Tuy nhiên, các tranh luận học thuật lớn trong y văn tập trung vào hai quan điểm đối lập:

  • Quan điểm giải tích đồng nhất hóa: Đa số các hướng dẫn tiêu chuẩn của IAEA (TRS-188, TRS-481) và ASTM giả định sản phẩm là môi trường đồng nhất tương đương nước để đơn giản hóa thuật toán lập kế hoạch chiếu xạ (Treatment Planning Systems/Industrial Dosimetry).
  • Quan điểm bất đồng nhất cấu trúc: Ponomarev (2001) và Chmielewski cùng cộng sự (2007) lập luận rằng trong các sản phẩm dạng ống rỗng đa lớp (PVC-Al-PVC-PE) hoặc bao gói công nghiệp có khoảng trống không khí, sự tán xạ đa góc Coulomb tại các mặt phân giới làm thay đổi cấu trúc trường liều, khiến thuật toán đồng nhất hóa gây sai số liều cục bộ nghiêm trọng từ 15% đến 30%.

Về mặt bức xạ hãm và an toàn quang hạt nhân, Kazuaki và cộng sự (2001) đã đo đạc phổ Bremsstrahlung trên bia Cu, W ở dải năng lượng 18, 28, 38 MeV, chỉ ra tính bất đối xứng không gian cực mạnh khi cường độ hướng $0^\circ$ lớn hơn $10^4$ lần so với hướng $90^\circ$ theo khuyến cáo TRS-188. Tại Việt Nam, Bùi Văn Loát và cộng sự (2007) khảo sát chùm bức xạ hãm trên máy gia tốc y tế Primus Siemens (6 - 15 MeV), trong khi Dương Thành Tài và cộng sự (2018) mô phỏng liều xạ trị bằng BEAMnrc/DOSXYZnrc với sai số 0,7%. Dẫu vậy, các nghiên cứu này chưa tích hợp phổ năng lượng thực tế với cấu trúc hình học sản phẩm phức tạp trong chiếu xạ công nghiệp quy mô lớn.

Luận án của Nguyễn Anh Tuấn định vị chính xác vào khoảng trống này bằng cách:

  1. Kết hợp đồng thời mô phỏng MCNP4c2 với thực nghiệm đo liều kế màng mỏng B3000 trên máy gia tốc RF-LINAC công nghiệp UELR-10-15S2 ($15\text{ kW}$, $9,92\text{ MeV}$) trong các cấu trúc hình học bất đồng nhất (dạng cột $\phi = 2,5\text{ cm}; 3,5\text{ cm}$ và dạng lớp cách khoảng $d = 1,0 - 5,0\text{ cm}$).
  2. Mở rộng so sánh với máy gia tốc plasma thế hệ mới LWFA dựa trên dữ liệu phổ thực nghiệm của Nakanii và cộng sự (2008) tại Viện Kỹ thuật Laser, Đại học Osaka ($3,5 - 30\text{ MeV}$), xác định trực tiếp tỷ lệ phát photon trên ngưỡng $(\gamma, n)$ trên các bia polypropylene (PP), Al, Fe, Pb.

Đó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 dừng hạt mang điện trong môi trường không gian giới hạn:

  • Mở rộng lý thuyết dừng Bethe-Bloch: Công thức Bethe-Bloch chuẩn xác định độ mất năng lượng ion hóa riêng $(-\frac{dE}{dx}){\text{ion}}$: $$-\left(\frac{dE}{dx}\right){\text{ion}} = 4\pi N_A r_e^2 m_e c^2 \frac{\rho Z}{A \beta^2} \left[ \ln\left(\frac{2m_e c^2 \beta^2 \gamma^2 T_0}{I^2}\right) + F(\tau) \right]$$ Nghiên cứu đã chỉ ra rằng trong môi trường bất đồng nhất có cấu trúc lớp và cột, mật độ electron hiệu dụng $n_e$ bị gián đoạn, làm cho tích phân mất năng lượng không còn là hàm liên tục theo chiều sâu hình học, dẫn đến hiện tượng tích tụ liều thứ cấp (build-up effect) tại các mặt phân giới vật liệu - không khí.
  • Bổ chính lý thuyết tán xạ Møller và Bremsstrahlung: Đóng góp làm rõ sự thay đổi tiết diện tán xạ vi phân $\frac{d\sigma_{\text{Møller}}}{dE}$ khi chùm electron tán xạ liên tiếp qua các khe rãnh không khí, giải thích nguyên nhân gây tăng liều cục bộ tại biên tiếp xúc giữa các cột sản phẩm hình trụ.
  • Chuyển dịch mô hình đánh giá kích hoạt hạt nhân: Chứng minh sự cần thiết phải chuyển từ mô hình kích hoạt đơn năng sang mô hình tích phân chập giữa phổ Bremsstrahlung liên tục với hàm tiết diện phản ứng $\sigma_{(\gamma, n)}(E_\gamma)$, đặc biệt khi đuôi phổ electron của RF-LINAC và LWFA vượt qua năng lượng ngưỡng $Q_{\gamma-n}$ của các cấu tử che chắn.
+-----------------------------------------------------------------------------------+
|                            KHUNG PHÂN TÍCH ĐỘC ĐÁO                                |
+-----------------------------------------------------------------------------------+
|  [Đặc tính nguồn electron]            [Tương tác môi trường phức hợp]             |
|  - RF-LINAC: Phổ Gauss                - Khối đồng nhất PP (rho = 0.95 g/cm3)       |
|    (9.92 +- 0.50 MeV)                 - Khối bất đồng nhất:                       |
|  - LWFA: Phổ rộng (3.5 - 30 MeV)        + Dạng cột (phi = 2.5 cm; 3.5 cm)         |
|                                         + Dạng lớp (d = 1.0 - 5.0 cm)             |
+----------------------------------+------------------------------------------------+
                                   |
                                   v
+-----------------------------------------------------------------------------------+
|              MÔ HÌNH TOÁN HỌC & MÔ PHỎNG MONTE CARLO (MCNP4c2)                    |
|  - Tán xạ ion hóa (Bethe-Bloch & Moller) -> Phân bố liều theo độ sâu Dn           |
|  - Phát xạ hãm hạt nhân (Bethe-Heitler) -> Phổ photon Bremsstrahlung              |
|  - Tích phân thông lượng Tally F4 -> Chuyển đổi liều hấp thụ (kGy)                |
+----------------------------------+------------------------------------------------+
                                   |
                                   v
+-----------------------------------------------------------------------------------+
|                       KẾT QUẢ ĐẦU RA & THỰC NGHIỆM ĐỐI CHỨNG                      |
|  - Kiểm chứng liều kế màng B3000 (0.5 - 35 kGy) trên RF-LINAC UELR-10-15S2        |
|  - Tối ưu hóa bề dày chiếu xạ: T_opt1, T_opt2 & Giới hạn mật độ mặt ADL           |
|  - Xác định tỷ lệ kích hoạt quang hạt nhân delta_(gamma, n) trên bia Al, Fe, Pb   |
+-----------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng: lý thuyết tương tác bức xạ ion hóa Bethe-Bloch/Møller, lý thuyết điện động lực học phát bức xạ hãm Coulomb và lý thuyết trường gia tốc sóng plasma phi tuyến Tajima-Dawson. Điểm mới nằm ở việc phát triển thuật toán liên kết giữa thông lượng hạt tính toán từ Tally F4 trong MCNP với công thức chuyển đổi liều hấp thụ thực tế: $$D(\text{kGy}) = \frac{D_e(\text{MeV}\cdot\text{cm}^2/\text{g}) \times I(\text{mA}) \times \eta}{10 \times S_w(\text{m}) \times v(\text{m/s})}$$ trong đó $D_e$ là năng lượng electron để lại trên đơn vị mật độ mặt, $I$ là dòng chùm electron ($15\text{ kW}$ hệ thống), $S_w$ là bề rộng quét ($40\text{ cm}$), $v$ là vận tốc băng tải ($1,0\text{ m/phút}$) và $\eta$ là hệ số sử dụng chùm tia. Điều kiện biên của khung phân tích được xác định rõ: dải năng lượng electron khảo sát từ 3,5 đến 30 MeV, vật liệu hấp thụ chuẩn có mật độ khối từ 0,1 đến $1,0\text{ g/cm}^3$, và khoảng cách hình học khe hở không khí từ 1,0 đến $5,0\text{ cm}$.

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ý thực chứng (positivism) nghiêm ngặt với phương pháp luận kết hợp chặt chẽ giữa mô phỏng tính toán Monte Carlo ngẫu nhiên và thực nghiệm đo liều đo quang phổ hạt nhân. Thiết kế nghiên cứu đa tầng được tổ chức theo sơ đồ:

  • Tầng 1 (Vi mô): Tương tác mức nguyên tử/hạt nhân, tính toán tiết diện vi phân của tán xạ Møller, quá trình sinh delta-ray và phổ photon phát xạ hãm bằng MCNP4c2.
  • Tầng 2 (Trung mô): Phân bố liều hấp thụ trong các cell hình cầu bán kính $r = 0,25\text{ cm}$ phân bố dọc theo trục đối xứng và các mặt phân giới của dummy box.
  • Tầng 3 (Vĩ mô): Trường bức xạ không gian buồng máy, thông lượng Bremsstrahlung theo các góc $0^\circ, 30^\circ, 45^\circ, 60^\circ, 90^\circ, 120^\circ$ và liều tích lũy trên băng tải hàng công nghiệp.

Mô hình giả lập được thực hiện chính xác theo thông số máy gia tốc RF-LINAC UELR-10-15S2: chùm electron quét xung tần số $141\text{ Hz}$, độ rộng xung $17\text{ }\mu\text{s}$, năng lượng trung bình $9,92 \pm 0,50\text{ MeV}$, bề rộng quét $40\text{ cm}$, liều bề mặt danh định cài đặt $13,5\text{ kGy}$.

+-----------------------------------------------------------------------------------+
|              THIẾT KẾ MÔ PHỎNG VÀ ĐO ĐẠC LIỀU THEO CHIỀU SÂU                      |
+-----------------------------------------------------------------------------------+
|  Chùm Electron (9.92 MeV) ---> [ Đầu quét 40 cm ]                                 |
|                                       |                                           |
|                                       v                                           |
|       +-------------------------------------------------------------------+       |
|       | THÙNG GIẢ HÀNG PP (24 x 20 x 10 cm3, rho = 0.95 g/cm3)            |       |
|       |                                                                   |       |
|       |  Mặt trên (d = 0 cm):   [ Liều kế B3000 ] -> Dn = 1.00 +- 0.05    |       |
|       |  Độ sâu d = 1.0 cm:     [ Liều kế B3000 ] -> Dn = 1.17 +- 0.06    |       |
|       |  Độ sâu d = 2.0 cm:     [ Liều kế B3000 ] -> Dn = 1.28 +- 0.06    |       |
|       |  Đỉnh liều d = 2.2 cm:  [ ĐỈNH LIỀU CỰC ĐẠI ] -> Dn = 1.32 (MCNP) |       |
|       |  Độ sâu d = 3.0 cm:     [ Liều kế B3000 ] -> Dn = 1.19 +- 0.06    |       |
|       |  Độ sâu d = 4.0 cm:     [ Liều kế B3000 ] -> Dn = 0.58 +- 0.04    |       |
|       |  Độ sâu d = 5.0 cm:     [ Liều kế B3000 ] -> Dn = 0.13 +- 0.01    |       |
|       |  Độ sâu d = 5.5 cm:     [ TẮT CHÙM ELECTRON HOÀN TOÀN ]          |       |
|       +-------------------------------------------------------------------+       |
+-----------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm và mô phỏng được chuẩn hóa theo các bước nghiêm ngặt:

  1. Chế tạo thùng giả hàng (Dummy Phantoms):
    • Thùng đồng nhất: Xếp chồng khít các tấm nhựa polypropylene (PP, mật độ $\rho = 0,95\text{ g/cm}^3$) kích thước $24 \times 20 \times 1\text{ cm}^3$.
    • Thùng dạng cột: Sử dụng các ống trụ tròn PP đặc đường kính $\phi = 2,5\text{ cm}$ và $\phi = 3,5\text{ cm}$ xếp sát nhau trong thùng carton $20 \times 20 \times 10\text{ cm}^3$, cài đặt liều kế dọc theo trục đối xứng cách nhau đúng $1,0\text{ cm}$.
    • Thùng dạng lớp: Bố trí các tấm PP dày $1,0\text{ cm}$ với khoảng cách khe hở không khí thay đổi có kiểm soát $d = 1,0\text{ cm}; 2,0\text{ cm}; 3,5\text{ cm}; 5,0\text{ cm}$.
  2. Quy trình đo liều chuẩn sắc ký:
    • Sử dụng hệ liều kế màng mỏng polymer B3000 của hãng GEX Corporation (Mỹ), độ dày đồng nhất, dải liều tuyến tính $0,5 - 35\text{ kGy}$, sai số đo chuẩn $\pm 5%$.
    • Chiếu xạ trên băng chuyền chuyển động đều tốc độ $1,0\text{ m/phút}$ qua đầu quét máy RF-LINAC.
    • Ổn định màu bằng phương pháp nhiệt luyện (annealing) trong tủ sấy chuyên dụng ở nhiệt độ $60^\circ\text{C}$ trong đúng 15 phút.
    • Đọc mật độ quang học (Optical Density) trên máy quang phổ phân tích GENESYS 20 (Thermo Scientific, Mỹ) tại bước sóng hấp thụ cực đại $\lambda = 512\text{ nm}$.
  3. Kiểm chuẩn độ tin cậy và giá trị:
    • Tam giác đạc phương pháp (Methodological Triangulation): Đối chiếu trực tiếp giữa mô phỏng vi hạt Monte Carlo MCNP4c2 và kết quả đo vật lý liều kế B3000.
    • Độ lặp lại (Reliability): Mỗi điểm đo thực nghiệm được lặp lại 3 lần, tính giá trị trung bình và độ lệch chuẩn ($\sigma_N$). Chuẩn hóa liều theo liều bề mặt $D_n = D(z)/D(0)$ để triệt tiêu dao động dòng tức thời của máy phát.

Data và phân tích

Dữ liệu mô phỏng MCNP được thiết lập với số lượng lịch sử hạt $N = 10^7$ hạt electron để đảm bảo sai số thống kê tương đối của Tally F4 luôn nhỏ hơn $3%$ ở hầu hết các cell tính liều. Kết quả so sánh giữa mô phỏng MCNP và thực nghiệm đo liều kế B3000 trong sản phẩm đồng nhất chiếu xạ một mặt được tổng hợp chi tiết trong bảng số liệu sau:

Độ sâu $z$ (cm) Liều chuẩn hóa mô phỏng MCNP ($D_n$) Liều chuẩn hóa đo phim B3000 ($D_n$) Độ sai lệch tuyệt đối (%)
0 (Bề mặt) $1,00 \pm 0,03$ $1,00 \pm 0,05$ 0,00
1 $1,20 \pm 0,04$ $1,17 \pm 0,06$ 2,56
2 $1,32 \pm 0,04$ $1,28 \pm 0,06$ 3,13
3 $1,21 \pm 0,04$ $1,19 \pm 0,06$ 1,68
4 $0,63 \pm 0,05$ $0,58 \pm 0,04$ 8,70
5 $0,12 \pm 0,01$ $0,13 \pm 0,01$ 7,69

Độ sai lệch giữa mô phỏng và thực nghiệm luôn duy trì dưới $3,2%$ trong vùng liều hiệu dụng ($0 - 3\text{ cm}$), và đạt cực đại $8,7%$ tại độ sâu $4\text{ cm}$. Sai lệch này hoàn toàn nằm trong giới hạn cho phép của đo liều công nghiệp, xuất phát từ nguyên nhân chùm tia thực nghiệm có góc mở quét nhất định trong khi mô phỏng giả định chùm song song lý tưởng chiếu thẳng góc.

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

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

  1. Quy luật biến dạng đường cong phân bố liều trong môi trường bất đồng nhất:

    • Trong sản phẩm đồng nhất ($\rho = 0,95\text{ g/cm}^3$), liều chuẩn hóa đạt đỉnh $D_{N\max} = 1,32$ tại độ sâu $2,2\text{ cm}$, sau đó suy giảm nhanh và tắt hoàn toàn tại $5,5\text{ cm}$. Bề dày tối ưu cho chiếu xạ một mặt là $T_{\text{opt1}} = 3,45\text{ cm}$ ứng với giới hạn mật độ mặt $\text{ADL}1 = 3,28\text{ g/cm}^2$; chiếu xạ hai mặt đạt $T{\text{opt2}} = 7,60\text{ cm}$ ứng với $\text{ADL}_2 = 7,22\text{ g/cm}^2$.
    • Trong sản phẩm dạng cột ($\phi = 2,5\text{ cm}$ và $3,5\text{ cm}$), do mật độ trung bình thể tích giảm ($\rho_{\text{eff}} < \rho_{\text{thực}}$), đường phân bố liều kéo dài sâu hơn vào trong thùng hàng. Tuy nhiên, sự tán xạ electron qua khe hở biên làm giảm liều cực đại tại trục cột xuống $D_{N\max} \approx 1,22 - 1,25$, đồng thời làm tăng độ không đồng đều liều DUR lên $1,48 - 1,55$.
    • Trong sản phẩm dạng lớp với khe hở không khí $d = 1,0 - 5,0\text{ cm}$, đường phân bố liều bị đứt gãy cục bộ. Tại mỗi bề mặt tiếp giáp lớp PP kế tiếp, xuất hiện hiện tượng tích tụ liều thứ cấp làm tăng vọt liều bề mặt lớp từ 8% đến 14% so với tính toán giải tích liên tục.

  2. Đặc trưng góc phát xạ và cấu trúc phổ bức xạ hãm:

    • Kết quả khảo sát phổ Bremsstrahlung trên các bia PP, Al, Fe, Pb cho thấy hiệu suất chuyển đổi $k_{\text{conv}}$ tỷ lệ thuận với nguyên tử số $Z$ và bề dày bia.
    • Chùm photon bức xạ hãm có tính định hướng cực cao: thông lượng photon tại góc phát xạ $0^\circ$ (dọc trục tới của electron) lớn hơn $10^4$ lần so với hướng ngang $90^\circ$ và hướng ngược $120^\circ$. Điều này khẳng định che chắn an toàn bức xạ phòng máy phải tập trung tối đa vào hướng trực diện của chùm tia (beam stop).
    • Thiết kế bia chuyển đổi tia X hỗn hợp ba lớp $\text{Ti}-\text{H}_2\text{O}-\text{Pb}$ đạt hiệu suất chuyển đổi tối ưu mà vẫn đảm bảo khả năng tản nhiệt bằng dòng nước làm mát trung gian mà không làm suy giảm cường độ photon phát ra.

  3. Nguy cơ kích hoạt phản ứng quang hạt nhân $(\gamma, n)$ sinh neutron thứ cấp:

    • Năng lượng ngưỡng tạo phản ứng $Q_{\gamma-n}$ của các vật liệu phổ biến: $^{12}\text{C}$ ($18,72\text{ MeV}$), $^{16}\text{O}$ ($15,66\text{ MeV}$), $^{27}\text{Al}$ ($13,06\text{ MeV}$), $^{56}\text{Fe}$ ($11,20\text{ MeV}$), $^{57}\text{Fe}$ ($7,65\text{ MeV}$), $^{206}\text{Pb}$ ($8,09\text{ MeV}$), $^{207}\text{Pb}$ ($6,74\text{ MeV}$), $^{208}\text{Pb}$ ($7,37\text{ MeV}$).
    • Trên máy gia tốc RF-LINAC ($9,92 \pm 0,50\text{ MeV}$), do đuôi phổ năng lượng mở rộng đạt tới $10,5\text{ MeV}$, tỷ lệ chùm photon sinh ra có năng lượng vượt ngưỡng $Q_{\gamma-n}$ trên bia Pb là $\delta_{\gamma-n} = 0,08%$ đối với $^{207}\text{Pb}$ và $0,02%$ đối với $^{208}\text{Pb}$. Mặc dù tỷ lệ nhỏ, nhưng trong các máy công nghiệp công suất cao ($15\text{ kW}$), lượng neutron sinh ra đủ để gây kích hoạt phóng xạ các cấu kiện che chắn nếu sử dụng chì làm bia hãm trực tiếp.
    • Trên máy gia tốc LWFA với phổ năng lượng trải rộng $3,5 - 30\text{ MeV}$, tỷ lệ $\delta_{\gamma-n}$ tăng vọt: đạt $3,45%$ trên bia Al, $5,12%$ trên bia Fe và $8,67%$ trên bia Pb. Phổ LWFA không qua lọc từ trường sẽ tạo ra thông lượng neutron thứ cấp nguy hiểm, đe dọa nghiêm trọng an toàn bức xạ.
+-----------------------------------------------------------------------------------+
|               PHỔ BỨC XẠ HÃM VÀ NGƯỠNG PHẢN ỨNG QUANG HẠT NHÂN                     |
+-----------------------------------------------------------------------------------+
|  Năng lượng (MeV)                                                                 |
|   30 MeV +-------------------------------------------------------+ [Đỉnh LWFA]   |
|          |                                                       |                |
|   18.72  | - - - - - - - - Ngưỡng (gamma, n) của C-12 - - - - -  |                |
|   15.66  | - - - - - - - - Ngưỡng (gamma, n) của O-16 - - - - -  |                |
|   13.06  | - - - - - - - - Ngưỡng (gamma, n) của Al-27 - - - - - |                |
|   11.20  | - - - - - - - - Ngưỡng (gamma, n) của Fe-56 - - - - - |                |
|   10.00  | ================= Năng lượng danh định RF-LINAC ===== |                |
|    8.09  | * * * * * * * * Ngưỡng (gamma, n) của Pb-206 * * * * *| [VÙNG NGUY CƠ  |
|    7.65  | * * * * * * * * Ngưỡng (gamma, n) của Fe-57  * * * * *|  KÍCH HOẠT     |
|    7.37  | * * * * * * * * Ngưỡng (gamma, n) của Pb-208 * * * * *|  NEUTRON TRÊN  |
|    6.74  | * * * * * * * * Ngưỡng (gamma, n) của Pb-207 * * * * *|  BIA NẶNG]     |
|    0.00  +-------------------------------------------------------+                |
+-----------------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết học thuật: Luận án cung cấp bộ dữ liệu thực nghiệm - mô phỏng chuẩn tắc về sự tương tác của chùm electron phổ thực trong môi trường dị thể, lấp đầy khiếm khuyết của các mô hình giải tích đơn năng cổ điển.
  • Về mặt thiết kế công nghệ và thực tiễn: Cung cấp công thức hiệu chỉnh bề dày đóng gói sản phẩm và bảng thông số $T_{\text{opt}}, \text{ADL}$ cho các nhà máy chiếu xạ thực phẩm và y tế, giúp triệt tiêu hiện tượng thiếu liều cục bộ (under-dosing gây mất vô trùng) hoặc quá liều (over-dosing gây hỏng phẩm chất cơ lý của vật liệu polyme và nông sản).
  • Về mặt chính sách và an toàn hạt nhân: Đưa ra khuyến cáo thực tiễn gửi Cục An toàn Bức xạ và Hạt nhân:
    1. Tuyệt đối không sử dụng các tấm chắn bằng chì (Pb) hoặc thép hợp kim chứa $^{57}\text{Fe}$ trực tiếp dưới đầu chiếu chùm electron 10 MeV công suất lớn mà không có lớp hấp thụ trung gian.
    2. Bắt buộc lắp đặt hệ thống nam châm lọc năng lượng (magnetic chicane / energy filter) khi ứng dụng máy gia tốc laser-plasma LWFA trong công nghệ bức xạ để cắt bỏ phần đuôi phổ năng lượng cao trên 10 MeV.

Limitations và Future Research

Nghiên cứu có các giới hạn khoa học và điều kiện biên cụ thể:

  1. Giới hạn hình học Dummy: Mô hình thực nghiệm mới khảo sát các dạng hình học cơ bản (hình trụ tròn xếp cột và tấm phẳng xếp lớp). Sản phẩm công nghiệp thực tế có thể có hình dạng bất định hình phức tạp hoặc dạng bột rời lỏng có độ xốp biến thiên liên tục.
  2. Giới hạn phổ LWFA thực nghiệm: Dữ liệu phổ electron LWFA được lấy từ bộ số liệu thực nghiệm công bố của hệ Petawatt Laser tại Viện Kỹ thuật Laser Osaka (Nakanii et al., 2008), chưa được đo đạc trực tiếp trên hệ máy đặt tại Việt Nam do hạn chế hạ tầng laser công suất siêu cao trong nước.
  3. Giới hạn đo thông lượng neutron trực tiếp: Xác suất phản ứng $(\gamma, n)$ được đánh giá gián tiếp thông qua tích phân phổ Bremsstrahlung vượt ngưỡng $Q_{\gamma-n}$ kết hợp thư viện số liệu hạt nhân NNDC/IAEA, chưa tiến hành kích hoạt vàng ($^{197}\text{Au}$) hoặc dùng ống đếm tỉ lệ $^{3}\text{He}$ đo phổ neutron thực tế tại phòng máy.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng thuật toán mô phỏng MCNP tự động hóa liên kết với chụp cắt lớp vi tính công nghiệp (Industrial CT Scan) để tái tạo bản đồ phân bố liều 3D thời gian thực cho kiện hàng có cấu trúc ngẫu nhiên.
  • Nghiên cứu thực nghiệm đo thông lượng và phổ năng lượng neutron thứ cấp sinh ra xung quanh đầu quét máy RF-LINAC $15\text{ kW}$ bằng phương pháp kích hoạt lá kim loại đa nguyên tố (Activation Foil Method).
  • Thiết kế hệ quang học từ trường (beamline transport system) tối ưu hóa việc đơn sắc hóa chùm electron từ nguồn LWFA nhằm ứng dụng cho các thiết bị xạ trị di động thế hệ mới.

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

Luận án mang lại các giá trị khoa học và kinh tế - xã hội có thể lượng hóa rõ ràng:

  • Tác động học thuật: Đóng góp trực tiếp vào hệ thống dữ liệu vật lý hạt nhân ứng dụng của Đại học Quốc gia TP.HCM và Viện Năng lượng Nguyên tử Việt Nam (VINATOM). Các bài báo công bố từ luận án trên các tạp chí chuyên ngành uy tín trong và ngoài nước cung cấp tài liệu tham khảo chuẩn cho nghiên cứu tính liều và che chắn bức xạ.
  • Tác động công nghiệp và chuyển giao công nghệ: Ứng dụng trực tiếp tại Trung tâm Nghiên cứu và Triển khai Công nghệ Bức xạ trong việc tối ưu hóa quy trình chiếu xạ thanh trùng thủy hải sản xuất khẩu và găng tay y tế. Tăng hiệu suất sử dụng chùm tia electron từ 60% lên xấp xỉ 70%, giúp tiết kiệm hàng trăm triệu đồng chi phí điện năng vận hành máy gia tốc mỗi năm và giảm thời gian lưu kho của hàng hóa.
  • Tác động xã hội và an toàn bức xạ: Cung cấp cơ sở khoa học để cơ quan quản lý nhà nước rà soát quy chuẩn an toàn bức xạ máy gia tốc công nghiệp, ngăn ngừa nguy cơ phơi nhiễm neutron thứ cấp cho nhân viên vận hành và bảo vệ môi trường sinh thái xung quanh cơ sở chiếu xạ.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Vật lý Hạt nhân/Kỹ thuật Hạt nhân: Tiếp cận phương pháp luận chuẩn tắc kết hợp giữa mô phỏng MCNP chuyên sâu và kỹ thuật đo liều màng phim B3000 thực nghiệm.
  • Kỹ sư vận hành và chuyên gia an toàn bức xạ tại các trung tâm chiếu xạ công nghiệp: Nắm vững phương pháp tính toán bề dày tối ưu ($T_{\text{opt}}$), xác định giới hạn mật độ mặt ($\text{ADL}$) và quy tắc thiết kế che chắn beam dump an toàn.
  • Các doanh nghiệp chế biến xuất khẩu nông sản, thực phẩm và vật tư y tế: Được hưởng lợi từ quy trình chiếu xạ đồng đều, đảm bảo liều hấp thụ đạt chuẩn kiểm dịch quốc tế (FDA, USDA, ISO 11137) mà không làm suy giảm chất lượng cảm quan của sản phẩm.
  • Cơ quan quản lý an toàn bức xạ hạt nhân: Sở hữu luận cứ khoa học tin cậy để thẩm định hồ sơ an toàn, cấp phép vận hành các máy gia tốc chùm electron năng lượng cao trên toàn quốc.

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 lý thuyết mất năng lượng Bethe-Bloch và tán xạ Møller sang môi trường bất đồng nhất có cấu trúc hình học gián đoạn (dạng cột và lớp). Nghiên cứu chỉ ra rằng các khe hở không khí gây ra hiệu ứng tán xạ cạnh làm suy giảm liều cực đại tại tâm trục cột từ $D_{N\max} = 1,32$ xuống $1,22 - 1,25$, đồng thời hiện tượng tích tụ liều thứ cấp tại mặt phân giới làm tăng liều tức thời bề mặt lớp kế tiếp từ $8%$ đến $14%$.

  2. Cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây như thế nào? Trả lời: So với nghiên cứu của Miller R. (2005) (chỉ khảo sát môi trường nước đồng nhất bằng code Tiger) và Kazuaki K. et al. (2001) (chỉ đo phổ Bremsstrahlung trên bia kim loại đơn chất ở góc cố định), luận án đã:

    • Xây dựng mô hình vi mô đa cell trong MCNP4c2 mô phỏng chính xác cấu trúc hình học dummy 3D bất đồng nhất.
    • Đối chứng thực nghiệm với hệ liều kế B3000 đạt độ sai lệch chặt chẽ $< 3,2%$ ở vùng hiệu dụng trên máy gia tốc RF-LINAC công nghiệp UELR-10-15S2 ($15\text{ kW}$).
    • Tích hợp phổ liên tục thực nghiệm của máy gia tốc laser-plasma LWFA Osaka để đánh giá định lượng nguy cơ kích hoạt hạt nhân.

  3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm và vật lý hạt nhân? Trả lời: Phát hiện bất ngờ nhất là mặc dù máy gia tốc RF-LINAC có năng lượng danh định là 10 MeV, nhưng thành phần đuôi phổ năng lượng cao ($9,92 \pm 0,50\text{ MeV}$, vươn tới $10,5\text{ MeV}$) vẫn tạo ra tỷ lệ photon bức xạ hãm $\delta_{\gamma-n} = 0,08%$ vượt ngưỡng kích hoạt quang hạt nhân của $^{207}\text{Pb}$ ($Q_{\gamma-n} = 6,74\text{ MeV}$) và $^{57}\text{Fe}$ ($7,65\text{ MeV}$). Điều này cảnh báo chì không phải là vật liệu che chắn an toàn tuyệt đối cho chùm electron 10 MeV công suất lớn nếu không có giải pháp triệt tiêu neutron.

  4. Luận án có cung cấp quy trình nhân bản thực nghiệm (Replication Protocol) không? Trả lời: Toàn bộ quy trình được chuẩn hóa chi tiết: từ thông số hình học dummy ($24 \times 20 \times 1\text{ cm}^3$, $\rho = 0,95\text{ g/cm}^3$), quy cách cài đặt liều kế B3000, thông số máy RF-LINAC ($141\text{ Hz}, 17\text{ }\mu\text{s}, 40\text{ cm}, 1,0\text{ m/phút}$), quy trình ủ nhiệt ổn định màu ($60^\circ\text{C}$ trong 15 phút) và đo mật độ quang học ở bước sóng $512\text{ nm}$ trên máy GENESYS 20. Bất kỳ phòng thí nghiệm bức xạ nào cũng có thể tái lập chính xác nghiên cứu này.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Trả lời: Định hướng phát triển tập trung vào ba mũi nhọn: (1) Tích hợp trí tuệ nhân tạo và xử lý ảnh CT Scan để tối ưu hóa tự động trường chiếu xạ theo thời gian thực cho kiện hàng phi đồng nhất; (2) Chế tạo hệ lọc chùm tia từ trường nhỏ gọn ứng dụng cho máy gia tốc plasma LWFA hướng tới xạ trị ung thư tại giường bệnh; (3) Đo đạc phổ neutron tán xạ bằng phương pháp phổ kế khối hạt nhân phục vụ an toàn bức xạ thế hệ mới.

Kết luận

Luận án tiến sĩ của Nguyễn Anh Tuấn là một công trình nghiên cứu thực nghiệm và tính toán mô phỏng xuất sắc trong lĩnh vực Vật lý nguyên tử và hạt nhân tại Việt Nam, mang lại các đóng góp cốt lõi:

  1. Xác lập hoàn chỉnh bản đồ phân bố liều theo độ sâu của chùm electron 10 MeV trong các mô hình sản phẩm bất đồng nhất dạng cột và dạng lớp, chỉ ra quy luật biến dạng liều và độ lệch chuẩn so với môi trường đồng nhất.
  2. Xác định chính xác các đại lượng công nghệ then chốt cho máy gia tốc RF-LINAC UELR-10-15S2: bề dày tối ưu chiếu một mặt $T_{\text{opt1}} = 3,45\text{ cm}$ ($\text{ADL}1 = 3,28\text{ g/cm}^2$), chiếu hai mặt $T{\text{opt2}} = 7,60\text{ cm}$ ($\text{ADL}_2 = 7,22\text{ g/cm}^2$).
  3. Khám phá đặc trưng định hướng cực mạnh của phổ bức xạ hãm Bremsstrahlung tại góc $0^\circ$ (lớn hơn $10^4$ lần so với góc $90^\circ$) và phát triển thành công cấu trúc bia chuyển đổi tia X ba lớp $\text{Ti}-\text{H}_2\text{O}-\text{Pb}$ tối ưu tản nhiệt.
  4. Phát hiện và định lượng nguy cơ sinh neutron thứ cấp từ phản ứng quang hạt nhân $(\gamma, n)$ do đuôi phổ năng lượng cao trên máy RF-LINAC và phổ năng lượng rộng ($3,5 - 30\text{ MeV}$) của máy gia tốc laser-plasma LWFA, đưa ra cảnh báo bắt buộc về an toàn vật liệu che chắn.
  5. Xây dựng cầu nối phương pháp luận vững chắc giữa tính toán Monte Carlo hạt nhân (MCNP4c2) với thực nghiệm đo liều màng phim chuẩn (B3000/GENESYS 20), đóng góp thiết thực cho ngành công nghệ bức xạ và bảo vệ môi trường hạt nhân tại Việt Nam và khu vực.