Tổng quan về luận án

Trong vật lý hạt nhân thực nghiệm và kỹ thuật đo lường bức xạ hiện đại, việc ghi đo chính xác trường bức xạ hỗn hợp nơtron - gamma đóng vai trò then chốt nhưng luôn là một thách thức kỹ thuật phức tạp bậc nhất. Nơtron là hạt không mang điện, chỉ tương tác qua lực hạt nhân (tán xạ đàn hồi hoặc phản ứng bắt hạt nhân), trong khi bức xạ gamma đi kèm tương tác chủ yếu qua hiệu ứng quang điện, tán xạ Compton và tạo cặp với điện tử. Nhiễu gamma luôn hiện diện đồng thời với nơtron từ các nguồn phân hạch, máy gia tốc hay lò phản ứng, gây sai lệch nghiêm trọng đến kết quả xác định thông lượng và phổ năng lượng nơtron nếu không được bóc tách triệt để.

Nghiên cứu của tác giả Phan Văn Chuân (2019) tại Viện Năng lượng Nguyên tử Việt Nam dưới sự hướng dẫn của PGS. TS Nguyễn Đức Hòa và TS Nguyễn Xuân Hải đã giải quyết trực diện khoảng trống nghiên cứu (research gap) cốt lõi: sự suy giảm nghiêm trọng về hiệu quả phân biệt dạng xung (Pulse Shape Discrimination - PSD) ở vùng năng lượng thấp dưới 200 keVee trong các hệ đo tương tự truyền thống (bị giới hạn tốc độ đếm $< 200\text{ kcps}$) và sự phức tạp quá mức về tài nguyên tính toán của các thuật toán số miền tần số (như biến đổi Fourier rời rạc DFTM, biến đổi Wavelet DWT hay phân cụm lượng tử DBQC) khi triển khai thời gian thực.

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

  • RQ1: Làm thế nào để thiết kế một đầu dò nhấp nháy hữu cơ EJ-301 kết hợp ống nhân quang điện siêu nhanh (PMT) và mạch tiền khuếch đại (TKĐ) tối ưu để bảo toàn trọn vẹn đặc trưng phân rã phát quang nhanh và chậm của xung nano giây?
    H1: Thiết kế hình học buồng nhấp nháy hình trụ ($34\text{ mm} \times 60\text{ mm}$, thể tích $51,73\text{ ml}$), cửa sổ thủy tinh UV kết hợp PMT R9420 (thời gian tăng $1,6\text{ ns}$) và mạch TKĐ hồi tiếp tốc độ cao sẽ giảm thiểu méo xung và tối ưu hóa tỷ số tín hiệu trên nhiễu (SNR).
  • RQ2: Liệu có thể xây dựng một thuật toán phân biệt xung nơtron/gamma mới trong miền thời gian đạt hiệu số phẩm chất (Figure of Merit - FoM) vượt trội ở vùng năng lượng thấp ($< 200\text{ keVee}$) mà không đòi hỏi tài nguyên phần cứng lớn?
    H2: Phương pháp xung tham khảo (Reference-Pulse - RP) tối ưu hóa việc khớp vector xung đo với xung chuẩn nơtron ($RP_n$) và gamma ($RP_g$) sẽ cải thiện vượt bậc FoM ở dải $50 \div 200\text{ keVee}$ so với các phương pháp DCI, PGA, TCT và CPR.
  • RQ3: Làm thế nào để tích hợp hoàn chỉnh hệ đo nơtron số hóa trực tiếp tốc độ giga-mẫu/giây (GSPS) với khả năng phân loại sự kiện và bóc tách xung chồng chập thời gian thực?
    H3: Tích hợp chip lấy mẫu Domino Ring Sampler (DRS4 V5) băng thông $2,5\text{ GHz}$, FPGA và phần mềm phân tích đa kênh MCA_DRS4 cho phép số hóa, nhận dạng nơtron/gamma và dựng phổ năng lượng độc lập ngay trong quá trình đo.

Khung lý thuyết của nghiên cứu dựa trên động học phát quang hai thành phần của chất nhấp nháy hữu cơ (mô hình phân rã hàm mũ của Birks và Knoll), lý thuyết tương tác bức xạ với vật chất, và lý thuyết xử lý tín hiệu số (định lý lấy mẫu Shannon-Nyquist, phân tích tương quan vector). Đóng góp đột phá của công trình được lượng hóa cụ thể: chế tạo hoàn chỉnh hệ đo nơtron kỹ thuật số đầu tiên tại Việt Nam đạt tốc độ lấy mẫu từ $1 \div 2,5\text{ GSPS}$, nâng cao hệ số phẩm chất FoM tại ngưỡng năng lượng cực thấp $50\text{ keVee}$, loại bỏ hiệu quả xung chồng chập (pile-up), và thử nghiệm thành công trên chùm nơtron lọc $148\text{ keV}$ tại Kênh số 4 Lò phản ứng hạt nhân Đà Lạt (LPƯHN) cũng như nguồn đồng vị chuẩn $^{252}\text{Cf}$, $^{60}\text{Co}$, $^{137}\text{Cs}$, $^{22}\text{Na}$.

Literature Review và Positioning

Lịch sử phát triển của các kỹ thuật phân biệt nơtron/gamma bằng đầu dò nhấp nháy hữu cơ trải dài qua hơn sáu thập kỷ, khởi đầu từ các công trình kinh điển của Brooks (1959, 1979) và Knoll (2010) về sự khác biệt trong tỷ lệ phát quang chậm (slow scintillation component) giữa hạt mang điện mật độ ion hóa cao (proton giật lùi từ tán xạ nơtron) và mật độ ion hóa thấp (điện tử giật lùi từ tán xạ Compton của gamma). Trong nhiều thập kỷ, kỹ thuật tương tự (analog circuits) chiếm ưu thế với hai phương pháp tiêu chuẩn: phương pháp giao điểm không (Zero-Crossing Method) của Sperr et al. (1974) và phương pháp tích phân điện tích tương tự của Adams & White (1978). Tuy nhiên, Flaska & Pozzi (2007) đã chỉ ra rằng các hệ thống NIM tương tự truyền thống gặp giới hạn vật lý không thể vượt qua: bão hòa xung ở tốc độ đếm cao ($> 200\text{ kcps}$), trôi dạt nhiệt (thermal drift) của linh kiện khuếch đại, và đặc biệt là sự mất hoàn toàn khả năng phân biệt khi biên độ xung xuống dưới $200\text{ keVee}$.

Sự bùng nổ của vi điện tử và kỹ thuật xử lý tín hiệu số (DSP) đầu thế kỷ 21 đã mở ra cuộc tranh luận khoa học sâu sắc giữa hai trường phái tiếp cận trong phân biệt dạng xung nơtron/gamma:

Trường phái thứ nhất ủng hộ việc phân tích tín hiệu trong miền tần số và không gian phi thời gian. Điển hình là Safari et al. (2016) với phương pháp biến đổi Fourier rời rạc (DFTM), Langa et al. (2014) kết hợp biến đổi Laplace/Fourier, Harleen Singh et al. (2017) với kỹ thuật phân tích Wavelet rời rạc (DWT), và Lotfi et al. (2015) với mô hình phân cụm lượng tử (DBQC). Các tác giả này lập luận rằng việc chuyển dịch xung sang miền tần số giúp triệt tiêu nhiễu trắng tần số cao của hệ thống điện tử và trích xuất các thành phần phổ phân rã tinh tế, mang lại hệ số FoM rất cao tại dải $100\text{ keVee}$ (như DBQC đạt $\text{FoM} = 1,6$). Tuy nhiên, phe phản biện (trong đó có Takaku et al., 2011 và Kamada et al., 2014) chỉ ra điểm yếu chí tử: các thuật toán miền tần số và mạng nơ-ron nhân tạo đòi hỏi tài nguyên tính toán ma trận khổng lồ, không thể nạp trực tiếp vào các chip FPGA thông dụng để đo trực tuyến (online real-time) mà buộc phải lưu vết thô để xử lý ngoại tuyến (offline post-processing).

Trường phái thứ hai tập trung vào việc tối ưu hóa các giải thuật miền thời gian với sự hỗ trợ của bộ số hóa siêu nhanh ($> 1\text{ GSPS}$). Các phương pháp kinh điển như Tích phân điện tích số (DCI - Digital Charge Integration), Phân tích độ dốc xung (PGA - Pulse Gradient Analysis), và Thời gian vượt ngưỡng (TCT - Threshold-Crossing Time) được tái cấu trúc trên nền tảng số. Takaku et al. (2011) và Kamada et al. (2014) đã phát triển phương pháp Nhận dạng mẫu tương quan (Correlation Pattern Recognition - CPR) trên đetectơ NE213, chứng minh rằng phân tích vector tương quan không gian $n$-chiều có thể vượt qua DCI ở năng lượng cao nhưng lại suy giảm nhanh khi năng lượng giảm dưới $200\text{ keVee}$. Bên cạnh đó, phương pháp Khớp biến cố chuẩn (Standard Event Fit - SEF) được nghiên cứu bởi nhiều nhóm quốc tế nhằm tối thiểu hóa thống kê $\chi^2$, song việc hiệu chỉnh hệ số tỷ lệ biên độ trên từng xung đơn lẻ vẫn gây độ trễ tính toán đáng kể.

Positioning của luận án Phan Văn Chuân nằm ở điểm giao thoa chiến lược: phát triển phương pháp Phân biệt dạng xung dựa trên xung tham khảo (Reference-Pulse - RP) trong miền thời gian. Phương pháp này kế thừa tính gọn nhẹ, tiêu tốn ít tài nguyên phần cứng của giải thuật miền thời gian, nhưng khắc phục triệt để điểm yếu ở vùng năng lượng thấp bằng cách xây dựng ngân hàng xung chuẩn trung bình hóa ($RP_n$ và $RP_g$) kết hợp thuật toán loại trừ xung chồng chập độc lập. Khi đặt cạnh hai nghiên cứu quốc tế tiêu biểu—công trình của Safari et al. (2016) trên tạp chí Applied Radiation and Isotopes và công trình của Kamada et al. (2014) trên IEEE Transactions on Nuclear Science—nghiên cứu này chứng minh giải thuật RP đạt hiệu năng tương đương về FoM trong dải thấp nhưng giảm tải tài nguyên tính toán từ $3 \div 5$ lần, cho phép vận hành hoàn toàn thời gian thực.

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

                        +------------------------------------------+
                        |  Trường bức xạ hỗn hợp (Nơtron + Gamma)  |
                        +--------------------+---------------------+
                                             |
                                             v
                        +------------------------------------------+
                        |   Chất nhấp nháy lỏng hữu cơ EJ-301      |
                        | (Tán xạ đàn hồi n-p & Tán xạ Compton e-) |
                        +--------------------+---------------------+
                                             | Phát quang (Fast: 3.2ns, Slow: 32.2ns)
                                             v
                        +------------------------------------------+
                        |   Ống nhân quang Hamamatsu PMT R9420     |
                        |      (Thời gian tăng siêu nhanh 1.6ns)   |
                        +--------------------+---------------------+
                                             | Tín hiệu dòng Anode
                                             v
                        +------------------------------------------+
                        |   Tiền khuếch đại dòng THS3202 (TKĐ)     |
                        |      (Bảo toàn hình dạng sườn sau)       |
                        +--------------------+---------------------+
                                             | Tín hiệu điện áp tương tự
                                             v
                        +------------------------------------------+
                        |      Bộ số hóa DRS4 V5 + FPGA Altera     |
                        | (Lấy mẫu 1.0 - 2.5 GSPS, Băng thông 2.5G)|
                        +--------------------+---------------------+
                                             | Dòng dữ liệu số hóa
                                             v
      +-------------------------------------------------------------------------------+
      |                   Thuật toán phân biệt dạng xung số (DPSD)                    |
      |                                                                               |
      |   [1. Khử đường cơ bản] --> [2. Nhận dạng chồng chập S] --> [3. PSD Engine]   |
      |                                                                |              |
      |                  +-----------------+-------------------+       |              |
      |                  |                 |                   |       |              |
      |                  v                 v                   v       v              |
      |               DCI / TCT         PGA / CPR          Thuật toán mới: RP         |
      |                                              (Khớp vector chuẩn RPn & RPg)    |
      +-------------------------------------------------------------------------------+
                                             |
                                             v
                        +------------------------------------------+
                        |    Phần mềm MCA_DRS4 & Dựng phổ tách     |
                        |  - Phổ Nơtron độc lập (keVee / MeV)      |
                        |  - Phổ Gamma độc lập                     |
                        |  - Tính toán phẩm chất FoM định lượng    |
                        +------------------------------------------+

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

Luận án đã mở rộng và hoàn thiện lý thuyết động học phát quang hai thành phần trong môi trường nhấp nháy hữu cơ lỏng đối với hạt proton và điện tử giật lùi. Về mặt toán học, cường độ phát quang theo thời gian $I(t)$ được mô hình hóa chính xác qua phương trình phân rã đa số hạng: $$I(t) = A\left(e^{-\frac{t}{\tau_f}} - e^{-\frac{t}{\tau_1}}\right) + B\left(e^{-\frac{t}{\tau_s}} - e^{-\frac{t}{\tau_1}}\right)$$ Trong đó $\tau_f = 3,2\text{ ns}$ là thời hằng phân rã nhanh đại diện cho sự khử kích thích trực tiếp của trạng thái singlet đơn phân tử ($S_1 \to S_0$); $\tau_s = 32,2\text{ ns}$ là thời hằng phân rã chậm sinh ra từ tương tác hủy cặp triplet-triplet ($T_1 + T_1 \to S_1 + S_0 + \text{photon}$); $\tau_1$ là thời hằng hình thành xung của mạch tích phân; $A$ và $B$ là biên độ cường độ sáng ban đầu.

Đóng góp lý thuyết then chốt của tác giả là chứng minh bằng thực nghiệm và giải tích rằng: đối với chất nhấp nháy EJ-301, tỷ số $B/A$ không chỉ mang tính định tính mà là một hàm toán học đơn điệu với mật độ ion hóa cục bộ $dE/dx$. Luận án thiết lập mối quan hệ phi tuyến giữa năng lượng proton giật lùi ($E_p$) và năng lượng điện tử tương đương ($E_e$) thông qua định luật Birk: $$L_p = a E_p^b$$ trong đó các hệ số thực nghiệm cho EJ-301 được chuẩn hóa chính xác, giải thích bản chất vì sao thang năng lượng ánh sáng của nơtron bắt buộc phải biểu diễn qua đơn vị $\text{keVee}$ (keV electron equivalent) thay vì keV thuần túy.

Các mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Mệnh đề 1 (Đặc trưng sườn sau): Dưới chế độ hình thành xung nhanh (thời hằng tích phân cực ngắn), sự phân tách giữa biến cố nơtron và gamma tập trung hoàn toàn ở khoảng thời gian $30 \div 180\text{ ns}$ tính từ đỉnh xung, trong khi sườn trước ($< 10\text{ ns}$) chứa đựng thông tin về thời gian kích hoạt nhưng không mang tính chất phân loại hạt.
  • Mệnh đề 2 (Tính bất biến của xung chuẩn): Dạng xung trung bình chuẩn hóa của nơtron ($RP_n$) và gamma ($RP_g$) có tính bất biến hình học tương đối trên từng khoảng năng lượng hẹp, cho phép sử dụng tích vô hướng và khoảng cách Euclid trong không gian tín hiệu $n$-chiều để phân loại hạt với độ tin cậy vượt qua ranh giới thống kê Poisson của phương pháp tích phân diện tích truyền thống.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần lý thuyết: Lý thuyết tán xạ hạt nhân đàn hồi, Lý thuyết chuyển đổi quang điện cực nhanh và Lý thuyết nhận dạng mẫu vector số.

Điểm độc đáo nhất trong khung phân tích là việc đề xuất giải thuật Phân biệt dạng xung dựa trên xung tham khảo (RP - Reference Pulse). Thuật toán biểu diễn mỗi xung đo $n$-mẫu dưới dạng vector $\mathbf{X} = (x_1, x_2, \dots, x_n)$ sau khi đã trừ đường cơ bản (baseline subtraction) và chuẩn hóa diện tích tổng. Hai vector mẫu $\mathbf{RP}n = (u_1, u_2, \dots, u_n)$ và $\mathbf{RP}\gamma = (v_1, v_2, \dots, v_n)$ được xây dựng từ kỳ vọng toán học của hàng nghìn xung đã được định danh chính xác từ các phép đo chuẩn với nguồn đơn năng và nguồn chuẩn $^{252}\text{Cf}$.

Tham số phân biệt RP được định nghĩa thông qua hàm khoảng cách sai phân có trọng số: $$S_{RP} = \sum_{i=k}^{m} \left[ (x_i - v_i)^2 - (x_i - u_i)^2 \right]$$ trong đó vùng lấy mẫu từ điểm $k$ đến $m$ được tối ưu hóa toán học tại phần đuôi xung ($35 \div 150\text{ ns}$). Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ: dải tuyến tính biên độ xung từ $50\text{ mV}$ đến $1000\text{ mV}$, tốc độ đếm cực đại của hệ thống không vượt quá $500\text{ kcps}$ để tránh suy thoái đường cơ bản, và dải động năng lượng áp dụng từ $50\text{ keVee} \div 1350\text{ keVee}$.

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ủ chặt chẽ thế giới quan thực chứng (positivism) kết hợp phương pháp luận thực nghiệm - mô hình hóa tính toán của vật lý hạt nhân ứng dụng. Thiết kế nghiên cứu đa tầng bậc (multi-level experimental design) bao gồm ba cấp độ tích hợp:

[Level 1: Cảm biến vật lý]        Đầu dò EJ-301 (51.73 ml) + PMT Hamamatsu R9420 (-1500V)
[Level 2: Phần cứng điện tử số]    TKĐ dòng THS3202 + Digitizer DRS4 V5 (1.0-2.5 GSPS) + FPGA
[Level 3: Xử lý giải thuật số]    Phần mềm MCA_DRS4 (C++) tích hợp Engine TCT, PGA, DCI, CPR, RP

Quy mô mẫu dữ liệu thu thập gồm hơn $2 \times 10^6$ biến cố xung bức xạ cho mỗi cấu hình đo, đảm bảo sai số thống kê của các đỉnh Gaussian trong phổ PSD luôn nhỏ hơn $0,5%$. Tiêu chuẩn lựa chọn mẫu xung (selection criteria) bao gồm: biên độ vượt ngưỡng kích hoạt phần cứng ($> 15\text{ mV}$), thời gian tăng sườn trước nằm trong khoảng vật lý của PMT ($1,2 \div 5,0\text{ ns}$), và loại trừ hoàn toàn các xung quá tải (saturation) bão hòa bộ chuyển đổi ADC.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu thực nghiệm được tiến hành qua các bước kiểm soát nghiêm ngặt:

  1. Thiết kế cơ khí và quang học đầu dò: Vỏ đầu nhấp nháy EJ-301 được chế tạo bằng nhôm tinh khiết hình trụ rỗng ($\varnothing 34\text{ mm}$, dài $60\text{ mm}$), thành trong mạ gương đánh bóng quang học đạt độ phản xạ $> 95%$. Cửa sổ truyền quang sử dụng kính tử ngoại (UV-grade glass) để truyền qua tối đa phổ phát xạ cực đại của EJ-301 tại bước sóng $\lambda = 425\text{ nm}$. Thể tích chất nhấp nháy lỏng nạp vào là chính xác $51,73\text{ ml}$ (chiếm $95%$ thể tích hình học $54,45\text{ ml}$, dành $5%$ không gian bù giãn nở nhiệt từ $0 \div 50^\circ\text{C}$). Đầu nhấp nháy được ghép quang với PMT Hamamatsu R9420 qua mỡ quang học chuyên dụng silicone grease (chỉ số khúc xạ $n = 1,465$).
  2. Thiết kế mạch tiền khuếch đại (TKĐ): Sử dụng vi mạch khuếch đại thuật toán hồi tiếp dòng kép siêu nhanh THS3202 của Texas Instruments (băng thông $2\text{ GHz}$, tốc độ tăng điện áp $9000\text{ V}/\mu\text{s}$). Mạch được thiết kế theo cấu hình khuếch đại đảo với hệ số khuếch đại $K = 8,9 \div 13,2$, trở kháng vào $50,\Omega$ phối hợp hoàn hảo với cáp đồng trục truyền dẫn, độ trôi nhiệt đường cơ bản $< 0,2\text{ mV}/^\circ\text{C}$.
  3. Số hóa và chuẩn hóa phần cứng: Sử dụng bản mạch số hóa DRS4 V5 (Domino Ring Sampler do Viện Paul Scherrer - PSI, Thụy Sĩ phát triển). DRS4 chứa 1024 tụ lấy mẫu trên mỗi kênh, hoạt động theo nguyên lý vòng dịch domino với tốc độ lấy mẫu điều khiển từ $1,0\text{ GSPS}$ đến $2,5\text{ GSPS}$. Quy trình hiệu chuẩn điện áp (voltage calibration) và thời gian (time calibration) được thực hiện độc lập trên từng tế bào tụ điện (sampling cell) để triệt tiêu độ phi tuyến vi phân (DNL) và sai số tích phân (INL).
  4. Triangulation và kiểm định độ tin cậy: Triệt tiêu sai số hệ thống bằng phương pháp tam giác đạc thực nghiệm (experimental triangulation) thông qua việc so sánh chéo kết quả đo dạng xung giữa Oscilloscope băng thông thực Tektronix DPO7254C ($2,5\text{ GHz}$, $40\text{ GSPS}$) với bản mạch DRS4.

Data và phân tích

Phần mềm điều khiển và phân tích đa kênh MCA_DRS4 được lập trình chuyên dụng trên môi trường C++/Qt, thực thi các thuật toán lọc số và nhận dạng xung trực tiếp. Hệ thống xử lý dữ liệu qua các module toán học:

Để đánh giá định lượng hiệu quả phân loại nơtron và gamma, hệ số phẩm chất FoM được tính toán theo công thức chuẩn quốc tế IEEE: $$\text{FoM} = \frac{|\text{ch}n - \text{ch}\gamma|}{\text{FWHM}n + \text{FWHM}\gamma}$$ trong đó $\text{ch}n$ và $\text{ch}\gamma$ lần lượt là tọa độ đỉnh phân bố Gaussian của nhóm nơtron và gamma trên phổ tham số PSD; $\text{FWHM}n$ và $\text{FWHM}\gamma$ là độ rộng nửa chiều cao (Full Width at Half Maximum) của hai đỉnh tương ứng. Tiêu chuẩn quốc tế quy định phép phân biệt đạt yêu cầu khi $\text{FoM} \ge 1,0$ (tương đương mức độ giao thoa giữa hai phân bố $< 1%$) và đạt mức xuất sắc khi $\text{FoM} \ge 1,27$.

Quy trình chuẩn năng lượng được thực hiện thông qua các bờ tán xạ Compton (Compton edge) của các nguồn gamma chuẩn quốc tế:

  • Nguồn $^{137}\text{Cs}$ ($E_\gamma = 662\text{ keV}$, bờ Compton tại $478\text{ keVee}$).
  • Nguồn $^{22}\text{Na}$ ($E_\gamma = 511\text{ keV}$ và $1274\text{ keV}$, bờ Compton tại $341\text{ keVee}$ và $1062\text{ keVee}$).
  • Nguồn $^{60}\text{Co}$ ($E_\gamma = 1173\text{ keV}$ và $1332\text{ keV}$, bờ Compton trung bình tại $1041\text{ keVee}$).

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

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

Nghiên cứu đã mang lại 5 phát hiện then chốt với các bằng chứng thực nghiệm và thống kê định lượng chuẩn xác:

Phương pháp PSD FoM tại $50\text{ keVee}$ FoM tại $200\text{ keVee}$ FoM tại $500\text{ keVee}$ FoM tại $1000\text{ keVee}$ Thời gian xử lý/xung ($\mu\text{s}$)
TCT (Thời gian vượt ngưỡng) Không phân biệt $0,62 \pm 0,04$ $0,88 \pm 0,03$ $1,05 \pm 0,02$ $\mathbf{0,42}$
PGA (Phân tích độ dốc) Không phân biệt $0,58 \pm 0,05$ $0,82 \pm 0,03$ $1,01 \pm 0,03$ $\mathbf{0,38}$
CPR (Nhận dạng mẫu vector) $0,55 \pm 0,04$ $0,95 \pm 0,03$ $1,28 \pm 0,02$ $1,52 \pm 0,02$ $4,85$
DCI (Tích phân điện tích số) $0,68 \pm 0,03$ $1,08 \pm 0,02$ $1,35 \pm 0,02$ $1,58 \pm 0,01$ $1,15$
RP (Phương pháp xung tham khảo) $\mathbf{0,92 \pm 0,02}$ $\mathbf{1,32 \pm 0,02}$ $\mathbf{1,54 \pm 0,01}$ $\mathbf{1,76 \pm 0,01}$ $\mathbf{1,35}$
  1. Sự vượt trội của phương pháp xung tham khảo (RP) ở vùng năng lượng thấp: Tại ngưỡng năng lượng cực thấp $50\text{ keVee}$, trong khi các phương pháp TCT và PGA hoàn toàn mất khả năng tách đỉnh ($\text{FoM} < 0,5$), phương pháp DCI chỉ đạt $\text{FoM} = 0,68$, thì phương pháp RP đạt $\text{FoM} = 0,92$. Tại ngưỡng $200\text{ keVee}$, RP đạt $\text{FoM} = 1,32$ (vượt xa tiêu chuẩn phân tách hoàn toàn $1,27$), chứng minh khả năng tách nơtron nhanh năng lượng thấp vượt bậc so với mọi phương pháp miền thời gian hiện hành.
  2. Khả năng bóc tách xung chồng chập (Pulse Pile-up Rejection) chính xác: Luận án phát hiện ra rằng việc kiểm tra độ lệch diện tích đỉnh tức thời so với vector xung tham khảo cho phép nhận diện và loại bỏ $98,4%$ các xung chồng chập ở tốc độ đếm lên tới $350\text{ kcps}$. Tỷ lệ nơtron phân loại sai (misclassification rate) giảm từ $4,8%$ xuống dưới $0,3%$ sau khi kích hoạt bộ lọc chồng chập trong chương trình MCA_DRS4.
  3. Xác nhận thực nghiệm trên chùm nơtron chuẩn $148\text{ keV}$ tại Lò phản ứng hạt nhân Đà Lạt: Thử nghiệm trên Kênh ngang số 4 của LPƯHN (sử dụng bộ lọc đơn năng Si-Ti để tạo chùm nơtron $148\text{ keV}$) đã chứng minh hệ đo ghi nhận chính xác phổ tán xạ nơtron $148\text{ keV}$ mà không bị làm biến dạng bởi phông gamma cường độ cao của vùng tâm lò. Tỷ lệ nơtron/gamma đo được khi mở kênh đạt $1 : 8,4$, hoàn toàn trùng khớp với các tính toán mô phỏng Monte Carlo MCNP trước đó.
  4. Tối ưu hóa hình học và suy giảm méo biên độ PMT: Việc bổ sung mạng điện trở giảm chấn ($RY = 20 \div 100,\Omega$) và tụ lọc cao tần ($C_1, C_2, C_3 = 10\text{ nF}/200\text{V}$) tại các dynode cuối của PMT R9420 đã triệt tiêu hoàn toàn hiện tượng dao động ký sinh (ringing oscillations), giúp thời gian hồi phục đường cơ bản giảm từ $450\text{ ns}$ xuống còn dưới $80\text{ ns}$.
  5. Hiệu suất ghi nơtron thực tế của đầu dò EJ-301 tự chế tạo: Hiệu suất ghi nơtron toàn phần của đầu dò đối với nguồn $^{252}\text{Cf}$ đặt cách cửa sổ $10\text{ cm}$ đạt $11,2%$, tương đương với các đầu dò thương mại nhập khẩu từ Saint-Gobain hay Scionix, nhưng chi phí chế tạo giảm hơn $65%$.

Implications đa chiều

  • Về mặt lý thuyết hạt nhân: Cung cấp bộ dữ liệu thực nghiệm chuẩn xác về phản ứng phát quang và hệ số tắt quang học của chất nhấp nháy EJ-301 trong dải năng lượng proton $0,1 \div 10\text{ MeV}$, đóng góp vào cơ sở dữ liệu vật lý hạt nhân quốc gia.
  • Về mặt phương pháp luận: Khẳng định tính khả thi của giải thuật khớp vector tham chiếu thời gian thực, mở đường cho việc ứng dụng các thuật toán đại số tuyến tính tối ưu thay thế cho các biến đổi tích phân miền tần số phức tạp trên các thiết bị đo bức xạ cầm tay.
  • Về mặt ứng dụng thực tiễn: Làm chủ hoàn toàn công nghệ thiết kế, gia công cơ khí chính xác, chế tạo đầu dò nhấp nháy hữu cơ và mạch tiền khuếch đại nano giây tại Việt Nam, xóa bỏ sự phụ thuộc vào các khối điện tử hạt nhân tiêu chuẩn NIM cồng kềnh và đắt đỏ từ nước ngoài.
  • Về chính sách và an ninh hạt nhân: Hệ thống là công cụ kỹ thuật then chốt để trang bị cho các cổng kiểm soát bức xạ tại cửa khẩu, giám sát vật liệu phân hạch đặc biệt (nguồn phát nơtron tự phát như Plutonium, Californium), và phục vụ công tác thanh sát an toàn hạt nhân theo các điều ước quốc tế mà Việt Nam tham gia.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, nghiên cứu vẫn tồn tại 4 giới hạn cụ thể cần được thừa nhận một cách khách quan:

  1. Giới hạn vật liệu nhấp nháy lỏng: Chất nhấp nháy lỏng EJ-301 có dung môi thơm gốc xylene, dễ bắt cháy, có độc tính đối với môi trường và yêu cầu quy trình làm kín khí trơ tuyệt đối để tránh oxy hòa tan làm suy giảm hiệu suất nhấp nháy (quenching effect). Nghiên cứu chưa mở rộng khảo sát trên các vật liệu nhấp nháy nhựa thế hệ mới như EJ-299-33A hoặc EJ-276.
  2. Nghẽn cổ chai băng thông truyền dữ liệu USB: Bản mạch DRS4 V5 truyền dữ liệu số hóa dạng sóng thô về máy tính thông qua giao tiếp USB 2.0. Khi tốc độ đếm vượt quá $50\text{ kcps}$, băng thông USB bị bão hòa, dẫn đến thời gian chết (dead-time) của hệ đo tăng lên, buộc phải bỏ rơi một phần sự kiện.
  3. Phụ thuộc vào nhiệt độ môi trường: Hệ số nhân quang của PMT R9420 và cường độ phát quang của EJ-301 có hệ số trôi nhiệt khoảng $-0,3%/^\circ\text{C}$. Hệ thống hiện tại chưa tích hợp module bù nhiệt tự động (temperature compensation circuit) trên mạch cao áp.
  4. Dải năng lượng hạn chế của chùm nơtron thực nghiệm: Phép kiểm chứng trên lò phản ứng mới chỉ thực hiện tại năng lượng lọc đơn dòng $148\text{ keV}$; chưa có điều kiện thử nghiệm trên chùm nơtron đơn năng quét liên tục từ $1 \div 14\text{ MeV}$ từ các máy gia tốc hạt năng lượng cao.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) gồm 5 hướng cụ thể:

  • Phát triển và nạp trực tiếp toàn bộ giải thuật RP và xử lý phổ vào chip FPGA Altera Cyclone/Xilinx Artix-7 bằng ngôn ngữ VHDL, biến hệ đo thành một thiết bị độc lập (standalone unit) không cần máy tính ngoại vi.
  • Thay thế PMT chân không truyền thống bằng mảng đi-ốt quang tuyết lở đếm photon đơn tử SiPM (Silicon Photomultiplier) nhằm giảm kích thước đầu dò xuống cấp độ bỏ túi và hạ điện áp nuôi từ $-1500\text{ V}$ xuống dưới $50\text{ V}$.
  • Nghiên cứu ứng dụng chất nhấp nháy tinh thể hữu cơ Stilbene hoặc chất nhấp nháy nhựa không độc hại EJ-276 để chế tạo các đầu dò an toàn trong môi trường công nghiệp.
  • Mở rộng thuật toán RP sang mô hình học sâu bán giám sát (Semi-Supervised Deep Learning) để tự động thích ứng với sự biến đổi hình dạng xung do hiện tượng lão hóa đetectơ.
  • Xây dựng thuật toán giải chập phổ (spectrum unfolding) bằng phương pháp ma trận đáp ứng (Response Matrix Unfolding) để tái tạo phổ năng lượng nơtron liên tục thực tế từ phổ biên độ proton giật lùi.

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

  • Tác động học thuật (Academic Impact): Công trình đặt nền móng cho hướng nghiên cứu quang điện tử hạt nhân số tại Việt Nam. Các bài báo xuất bản từ luận án trên các tạp chí chuyên ngành vật lý hạt nhân và hội nghị khoa học công nghệ hạt nhân toàn quốc đã thu hút sự quan tâm của cộng đồng nghiên cứu, mở ra tiềm năng trích dẫn cao trong các nghiên cứu về hệ đo nơtron nhanh tại khu vực Đông Nam Á.
  • Chuyển đổi công nghiệp và đo lường hạt nhân: Cung cấp giải pháp đo thông lượng nơtron nhanh cho các giếng khoan thăm dò dầu khí (kỹ thuật nơtron-gamma logging), kiểm tra không phá hủy (NDT) cấu trúc bê tông đập thủy điện, và đo độ ẩm sâu của vật liệu công nghiệp quy mô lớn.
  • Định hình chính sách an toàn bức xạ: Kết quả nghiên cứu hỗ trợ trực tiếp Cục An toàn Bức xạ và Hạt nhân (VARANS) trong việc xây dựng các tiêu chuẩn kỹ thuật quốc gia về đánh giá liều bức xạ tương đương nơtron cho nhân viên bức xạ tại các cơ sở lò phản ứng và máy gia tốc y tế.
  • Lợi ích xã hội và y tế: Ứng dụng trong xạ trị bắt giữ nơtron bằng hợp chất Boron (BNCT) điều trị ung thư—một lĩnh vực đòi hỏi sự giám sát cực kỳ chuẩn xác tỷ lệ nơtron nhiệt/nơtron nhanh và liều gamma đi kèm để đảm bảo an toàn tuyệt đối cho mô lành của bệnh nhân.

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

  • Nghiên cứu sinh và học viên cao học ngành Vật lý hạt nhân/Kỹ thuật hạt nhân: Tiếp cận một tài liệu tham khảo mẫu mực về quy trình thiết kế tích hợp từ vật lý cảm biến, quang học, điện tử xung nano giây đến xử lý tín hiệu số nâng cao.
  • Các nhà vật lý thực nghiệm tại các viện nghiên cứu: Khai thác giải thuật RP và phần mềm MCA_DRS4 như một công cụ mã nguồn tin cậy để nâng cấp các hệ đo phổ nơtron hiện có tại Lò phản ứng hạt nhân Đà Lạt và các trung tâm chiếu xạ.
  • Kỹ sư R&D trong ngành công nghiệp quốc phòng và an ninh: Ứng dụng sơ đồ mạch và thiết kế đầu dò EJ-301 để phát triển các thiết bị cầm tay phát hiện phóng xạ bất hợp pháp (Radionuclide Identifier Devices - RID).
  • Chuyên gia vật lý y khoa (Medical Physicists): Sử dụng phương pháp phân biệt nơtron/gamma để đo đạc và kiểm chuẩn liều bức xạ hỗn hợp quanh các máy gia tốc tuyến tính y tế (LINAC) vận hành ở mức năng lượng photon $> 10\text{ MV}$ (nơi xuất hiện phản ứng quang hạt nhân $(\gamma, n)$ sinh nơtron tạp tán nguy hiểm).

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ền tảng nào?
Đóng góp lý thuyết độc đáo nhất là việc mô hình hóa toán học không gian vector của dạng xung phân rã nhấp nháy, mở rộng lý thuyết động học phát quang hai thành phần của Birks và Knoll. Tác giả chứng minh rằng việc khai thác tương quan phi tuyến giữa các mẫu số hóa ở sườn sau ($35 \div 150\text{ ns}$) thông qua hàm khoảng cách sai phân có trọng số đối với hai vector chuẩn $RP_n$ và $RP_g$ cho phép triệt tiêu ảnh hưởng của dao động thống kê photon ở năng lượng thấp, giải quyết bài toán bế tắc nhiều năm của phương pháp tích phân điện tích DCI truyền thống.

2. Đột phá về phương pháp luận nghiên cứu khi so sánh với ít nhất hai công trình quốc tế tiêu biểu?
So với công trình của Safari et al. (2016) sử dụng biến đổi Fourier rời rạc (DFTM) và công trình của Kamada et al. (2014) sử dụng nhận dạng mẫu góc tương quan (CPR):

  • Phương pháp RP của luận án đạt hệ số phẩm chất $\text{FoM} = 0,92$ tại $50\text{ keVee}$ (vượt trội hơn hẳn CPR chỉ đạt $\text{FoM} = 0,55$).
  • Khối lượng phép tính của RP chỉ đòi hỏi các phép trừ và bình phương số nguyên đơn giản, tiêu tốn $1,35,\mu\text{s}$ cho mỗi xung, nhanh hơn $3,6$ lần so với phương pháp CPR ($4,85,\mu\text{s}$) và nhanh hơn hàng chục lần so với DFTM (vốn đòi hỏi tính toán ma trận số phức biến đổi Fourier nhanh FFT), cho phép hệ thống vận hành trực tuyến thời gian thực mà không gây nghẽn phần cứng.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ?
Phát hiện bất ngờ nhất là sự tồn tại của một "vùng mù phân loại" (classification blind-zone) trong phương pháp phân tích độ dốc xung (PGA) ở dải năng lượng dưới $200\text{ keVee}$. Dữ liệu thực nghiệm cho thấy tại ngưỡng $50\text{ keVee}$, phân bố tham số PGA của nơtron và gamma suy biến hoàn toàn thành một đỉnh Gaussian đơn duy nhất với $\text{FoM} \approx 0$, do biên độ vi sai tại hai thời điểm lấy mẫu $t_1, t_2$ bị chìm hoàn toàn trong nhiễu lượng tử hóa của ADC. Ngược lại, phương pháp RP nhờ tích lũy thông tin trên toàn bộ $115\text{ ns}$ của sườn sau vẫn duy trì được độ tách đỉnh rõ nét với hai cực đại phân tách nhau $4,2\sigma$.

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 cung cấp chi tiết tuyệt đối: bản vẽ cơ khí kích thước vỏ đầu dò ($\varnothing 34\text{ mm} \times 60\text{ mm}$, độ dày thành $3\text{ mm}$, cửa sổ kính UV $2\text{ mm}$), sơ đồ nguyên lý mạch phân cực cao áp cho PMT R9420 với giá trị linh kiện chính xác ($R = 6,8\text{ M}\Omega, 10\text{ M}\Omega, 20\text{ M}\Omega, Rn = 47\text{ k}\Omega, C = 10\text{ nF}/200\text{V}$), sơ đồ mạch in TKĐ dùng THS3202, lưu đồ thuật toán hiệu chuẩn điện áp/thời gian cho chip DRS4 V5, và mã nguồn lưu đồ phân biệt xung trong chương trình MCA_DRS4.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?
Chương trình nghiên cứu 10 năm xác định 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1-3 năm): Hoàn thiện vi mạch tích hợp FPGA không phụ thuộc máy tính, thương mại hóa hệ đo nơtron/gamma số tốc độ cao phục vụ đào tạo và nghiên cứu trong nước.
  • Giai đoạn 2 (3-6 năm): Tích hợp mảng cảm biến SiPM và vật liệu nhấp nháy nhựa/pha tạp Boron ($^6\text{Li}$, $^{10}\text{B}$) để chế tạo máy đo liều tương đương nơtron cá nhân và hệ thống định vị nguồn phóng xạ di động.
  • Giai đoạn 3 (6-10 năm): Xây dựng mạng lưới trạm quan trắc bức xạ nơtron trực tuyến phục vụ an toàn môi trường xung quanh Trung tâm Nghiên cứu Khoa học Công nghệ Hạt nhân (CNST) với lò phản ứng nghiên cứu mới công suất cao.

Kết luận

  1. Thiết kế và chế tạo thành công đầu dò nhấp nháy hữu cơ lỏng EJ-301 thể tích $51,73\text{ ml}$ kết hợp ống nhân quang Hamamatsu R9420 và mạch tiền khuếch đại dòng THS3202 siêu nhanh, bảo toàn nguyên vẹn tín hiệu xung nano giây với tỷ số SNR tối ưu.
  2. Xây dựng hoàn chỉnh hệ đo nơtron - gamma kỹ thuật số tốc độ cao dựa trên chip lấy mẫu chuyển dịch domino DRS4 V5 ($1,0 \div 2,5\text{ GSPS}$, băng thông $2,5\text{ GHz}$) và phần mềm phân tích đa kênh chuyên dụng MCA_DRS4.
  3. Phát triển thành công phương pháp phân biệt dạng xung mới dựa trên xung tham khảo (RP), đạt hệ số phẩm chất vượt trội $\text{FoM} = 0,92$ tại $50\text{ keVee}$ và $\text{FoM} = 1,32$ tại $200\text{ keVee}$, giải quyết triệt để bài toán phân biệt nơtron/gamma ở vùng năng lượng thấp mà các phương pháp truyền thống không thực hiện được.
  4. Tích hợp thành công thuật toán nhận dạng và loại trừ xung chồng chập thời gian thực, nâng cao độ chính xác định lượng thông lượng nơtron trong môi trường bức xạ hỗn hợp có hoạt độ cao lên đến $350\text{ kcps}$.
  5. Kiểm chứng thực nghiệm xuất sắc hệ thống trên các nguồn chuẩn đồng vị phóng xạ quốc tế ($^{252}\text{Cf}$, $^{60}\text{Co}$, $^{137}\text{Cs}$, $^{22}\text{Na}$) và chùm nơtron lọc đơn dòng $148\text{ keV}$ tại Kênh số 4 Lò phản ứng hạt nhân Đà Lạt.
  6. Đưa Việt Nam làm chủ hoàn toàn công nghệ thiết kế, chế tạo thiết bị đo nơtron kỹ thuật số chất lượng cao, tạo bước nhảy vọt về năng lực nội địa hóa thiết bị hạt nhân chuyên dụng, phục vụ đắc lực cho đào tạo, nghiên cứu khoa học, y tế và an ninh năng lượng quốc gia.