Tổng quan về luận án
Nghiên cứu hạt nhân và kỹ thuật đo bức xạ hiện đại đóng vai trò then chốt trong việc kiểm soát an toàn lò phản ứng, nghiên cứu cấu trúc vật liệu, y học hạt nhân và an ninh phóng xạ. Trong môi trường trường bức xạ hỗn hợp, việc ghi đo nơtron luôn đối mặt với thách thức kỹ thuật phức tạp do sự hiện diện đồng thời không thể tránh khỏi của bức xạ gamma nền. Luận án tiến sĩ của tác giả Phan Văn Chuân, thực hiện dưới sự hướng dẫn của PGS. TS Nguyễn Đức Hòa và TS Nguyễn Xuân Hải tại Trường Đại học Đà Lạt phối hợp với Viện Nghiên cứu Hạt nhân Đà Lạt, đại diện cho một bước tiến tiên phong trong việc làm chủ công nghệ thiết kế, chế tạo phần cứng đetectơ nhấp nháy lỏng kết hợp kỹ thuật xử lý tín hiệu số (DSP) và công nghệ mảng cổng lập trình được (FPGA) tại Việt Nam.
Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án giải quyết xuất phát từ giới hạn vật lý của các phương pháp phân biệt dạng xung nơtron/gamma (Pulse Shape Discrimination - PSD) truyền thống. Trong các nghiên cứu kinh điển về PSD tương tự (Brooks, 1959; Knoll, 2010), khả năng phân loại nơtron và gamma bị suy giảm nghiêm trọng ở vùng năng lượng thấp dưới 200 keVee do tỷ số tín hiệu trên nhiễu (SNR) thấp và hiện tượng chồng chập xung khi tốc độ đếm vượt quá 200 kcps. Ngược lại, các thuật toán số hóa hiện đại hoạt động trong miền tần số như biến đổi Fourier rời rạc DFTM (Safari et al., 2016) hay biến đổi Wavelet DWT (Singh et al., 2019) tuy cải thiện được độ phân biệt nhưng lại đòi hỏi năng lực tính toán cực lớn, tài nguyên phần cứng phức tạp và khó triển khai thời gian thực (online) trên các hệ đo nhỏ gọn.
Để giải quyết triệt để vấn đề này, luận án tập trung vào các câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để tối ưu hóa thiết kế hình học và điện tử đầu dò nhấp nháy lỏng EJ-301 nhằm thu được xung có thời gian tăng phù hợp cho việc truyền dẫn và số hóa mà không làm suy giảm các đặc trưng phân rã chậm nhạy cảm với nơtron?
- Câu hỏi nghiên cứu 2 (RQ2): Liệu một phương pháp phân biệt dạng xung mới dựa trên xung tham khảo (Reference-Pulse - RP) trong miền thời gian có thể vượt trội hơn các phương pháp tiêu chuẩn (DCI, CPR, PGA, TCT) về hệ số phẩm chất (Figure-of-Merit - FoM) ở vùng năng lượng dưới 200 keVee mà vẫn duy trì thuật toán đơn giản, tối ưu tài nguyên phần cứng?
- Câu hỏi nghiên cứu 3 (RQ3): Hệ đo số hóa tích hợp chip lấy mẫu nhanh DRS4 và phần mềm MCA_DRS4 có khả năng nhận dạng xung chồng chập và tái tạo phổ nơtron/gamma đồng thời trong điều kiện thực nghiệm khắc nghiệt tại kênh dẫn nơtron lò phản ứng hạt nhân hay không?
Các giả thuyết nghiên cứu tương ứng:
- Giả thuyết 1 (H1): Tích hợp mạch tiền khuếch đại (TKĐ) nhạy điện áp sử dụng vi mạch băng thông siêu cao THS3202 với thời gian tăng $t_{rise} > 10\text{ ns}$ sẽ loại bỏ hoàn toàn hiện tượng phản xạ xung trên đường truyền cáp đồng trục mà vẫn bảo toàn tỷ lệ tích phân đuôi phục vụ thuật toán PSD.
- Giả thuyết 2 (H2): Phương pháp xung tham khảo (RP) dựa trên tối ưu hóa thống kê độ tương đồng dạng sóng trung bình chuẩn hóa sẽ duy trì giá trị $\text{FoM} \ge 1.0$ tại ngưỡng năng lượng 50 keVee, vượt trội so với phương pháp tích phân điện tích số (DCI).
- Giả thuyết 3 (H3): Thuật toán phân tích số hóa trực tiếp cho phép loại trừ chính xác các xung chồng chập thời gian thực, nâng cao độ tin cậy của phổ nơtron 148 keV tại Kênh số 4 Lò phản ứng hạt nhân Đà Lạt.
Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng động học phân rã phát quang hữu cơ của Birks (1964), lý thuyết tương tác bức xạ với vật chất, và kỹ thuật chuyển đổi tương tự - số tốc độ siêu cao dựa trên mảng điện dung đóng ngắt DRS4 (Domino Ring Sampler) do Viện Paul Scherrer (PSI - Thụy Sĩ) phát triển. Đóng góp đột phá của luận án là chế tạo thành công đầu dò EJ-301 thể tích hoạt động 51,73 ml với cấu trúc quang - điện tử tối ưu, phát triển thuật toán RP độc quyền có hiệu năng vượt trội ở vùng năng lượng thấp (50 - 200 keVee), và hiện thực hóa thành công hệ đo MCA_DRS4 trực tuyến tại Lò phản ứng hạt nhân Đà Lạt (LPƯHN Đà Lạt) cũng như nguồn chuẩn $^{252}\text{Cf}$, $^{60}\text{Co}$, $^{137}\text{Cs}$, $^{22}\text{Na}$.
Literature Review và Positioning
Nghiên cứu về phân biệt bức xạ nơtron và gamma bằng đetectơ nhấp nháy hữu cơ bắt đầu từ công trình nền tảng của Brooks (1959), người đầu tiên phát hiện ra sự khác biệt về thành phần phát xạ chậm khi các hạt mang điện khác nhau ion hóa môi trường nhấp nháy. Trong nhiều thập kỷ, các hệ thống NIM tương tự sử dụng kỹ thuật giao điểm không (Zero-Crossing Method) hoặc kỹ thuật tích phân tỷ số dòng điện (Charge Comparison Method) là tiêu chuẩn công nghiệp (Knoll, 2010; Perkins et al., 1999). Tuy nhiên, các hệ thống tương tự bộc lộ nhược điểm cố hữu: thời gian chết (dead-time) lớn, độ trôi tham số theo nhiệt độ của các mạch khuếch đại tương tự, và sự bất lực trong việc xử lý các xung méo do chồng chập ở tốc độ đếm cao.
Sự bùng nổ của kỹ thuật số hóa tốc độ cao (>100 MSPS đến hàng GSPS) đã mở ra một trường phái nghiên cứu mới về DSP trong đo bức xạ hạt nhân. Các hướng tiếp cận hiện đại trong y văn quốc tế phân chia thành hai dòng chính:
┌─────────────────────────────────────────────────────────────┐
│ CÁC PHƯƠNG PHÁP PHÂN BIỆT DẠNG XUNG NƠTRON/GAMMA (PSD) │
└──────────────────────────────┬──────────────────────────────┘
│
┌───────────────────────────────┴───────────────────────────────┐
▼ ▼
┌───────────────────────────────────────────┐ ┌───────────────────────────────────────────┐
│ PHÂN TÍCH TRONG MIỀN THỜI GIAN │ │ PHÂN TÍCH TRONG MIỀN TẦN SỐ │
└─────────────────────┬─────────────────────┘ └─────────────────────┬─────────────────────┘
│ │
┌─────────────────┬───────────┴───────────┬──────────────────┐ ┌───────────────┴───────────────┐
▼ ▼ ▼ ▼ ▼ ▼
┌─────────┐ ┌───────────┐ ┌───────────┐ ┌───────────┐ ┌───────────┐ ┌───────────┐
│Phương │ │Phương │ │Phương │ │Phương │ │Biến đổi │ │Biến đổi │
│pháp TCT │ │pháp PGA │ │pháp DCI │ │pháp CPR │ │Fourier │ │Wavelet │
│(Vượt │ │(Độ dốc │ │(Tích phân │ │(Nhận dạng │ │DFTM │ │DWT / WPTM │
│ngưỡng) │ │sườn sau) │ │điện tích) │ │mẫu góc) │ │(Safari │ │(Singh │
│ │ │ │ │ │ │(Takaku, │ │et al., │ │et al., │
│ │ │ │ │ │ │Kamada) │ │2016) │ │2019) │
└─────────┘ └───────────┘ └─────┬─────┘ └───────────┘ └───────────┘ └───────────┘
│
▼
┌───────────────────────────┐
│ PHƯƠNG PHÁP MỚI CỦA │
│ LUẬN ÁN NÀY: │
│ Xung tham khảo tối ưu RP │
│ (Reference-Pulse Method) │
│ (FoM cao ở <200 keVee, │
│ tiết kiệm tài nguyên) │
└───────────────────────────┘
- Phân tích trong miền tần số và không gian phi tuyến: Safari et al. (2016) đã đề xuất phương pháp biến đổi Fourier rời rạc (DFTM) để trích xuất phổ biên độ pha, chứng minh tính bất biến cao đối với nhiễu trắng so với DCI. Singh et al. (2019) ứng dụng biến đổi Wavelet rời rạc (DWT) và Wavelet Packet Transform (WPTM), cho phép lọc đa độ phân giải để tối đa hóa khoảng cách giữa cụm nơtron và gamma. Gần đây, Lotfi et al. (2016) áp dụng thuật toán phân cụm lượng tử (Discrimination based on Quantum Clustering - DBQC) đạt giá trị $\text{FoM} = 1.6$ tại ngưỡng 100 keVee. Tuy nhiên, các kỹ thuật này đều gặp rào cản nghiêm trọng khi đưa vào chip xử lý nhúng do khối lượng phép tính ma trận phức tạp và độ trễ thuật toán lớn.
- Phân tích hình học và thống kê trong miền thời gian: Phương pháp tích phân điện tích số (DCI) vẫn là chuẩn quy chiếu toàn cầu (Flaska & Pozzi, 2007). Takaku et al. (2001) và Kamada et al. (2017) đã phát triển phương pháp nhận dạng mẫu tương quan (CPR) sử dụng tích vô hướng vector để tính góc cosin giữa xung đo và xung mẫu. Dù CPR cải thiện độ phân giải ở dải năng lượng cao, nó suy giảm nhanh chóng ở vùng năng lượng thấp dưới 200 keVee do sự biến dạng phi đối xứng của đỉnh xung bị chi phối bởi dao động thống kê photon.
Luận án của NCS Phan Văn Chuân định vị chính xác vào điểm giao thoa giữa hai trường phái: phát triển phương pháp xung tham khảo (RP) vận hành hoàn toàn trong miền thời gian nhưng khai thác thông tin hình học toàn diện của dạng sóng chuẩn hóa. Luận án tiến hành so sánh đối chuẩn trực tiếp với hai nghiên cứu quốc tế tiêu biểu:
- So với nghiên cứu của Safari et al. (2016) về DFTM: Phương pháp RP của luận án đạt hệ số phẩm chất tương đương trong dải $100 - 1000\text{ keVee}$ nhưng giảm hơn 70% tài nguyên logic DSP, loại bỏ hoàn toàn yêu cầu tính toán biến đổi Fourier phức.
- So với nghiên cứu của Kamada et al. (2017) về CPR với đetectơ NE213: Thuật toán RP giải quyết triệt để hiện tượng nhạy cảm nhiễu đường cơ bản, nâng cao FoM tại ngưỡng 50 keVee từ mức không thể phân tách lên mức tách biệt rõ nét ($\text{FoM} \approx 0.95 - 1.05$).
Đó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à cụ thể hóa mô hình động học phát quang của chất nhấp nháy hữu cơ EJ-301 dưới tác động của các hạt ion hóa có mật độ truyền năng lượng tuyến tính (LET) khác nhau. Quá trình phát quang trong chất nhấp nháy hữu cơ tuân theo phương trình tổng quát mô tả cường độ ánh sáng $I(t)$ do sự kích hoạt đồng thời các trạng thái đơn tuyến (singlet) và tam tuyến (triplet):
$$I(t) = A e^{-t/\tau_f} + B e^{-t/\tau_s}$$
Trong đó, $\tau_f$ là thời hằng phân rã nhanh ($3.2\text{ ns}$ đối với EJ-301), $\tau_s$ là thời hằng phân rã chậm ($32.2\text{ ns}$). Đối với tương tác của tia gamma, các điện tử tán xạ Compton có LET thấp chủ yếu kích thích trạng thái singlet, làm cho tỷ số biên độ thành phần chậm trên nhanh $B/A$ rất nhỏ. Ngược lại, nơtron nhanh tương tác tán xạ đàn hồi $(n, p)$ tạo ra các proton giật lùi có LET cao, làm gia tăng xác suất tương tác giữa các trạng thái triplet kích thích, dẫn đến sự gia tăng vượt trội của thành phần phân rã chậm $B$.
Mối quan hệ phi tuyến giữa năng lượng proton giật lùi $E_p$ và lượng ánh sáng phát ra $L_p$ (tính bằng đơn vị MeVee hoặc keVee) được luận án mô hình hóa chặt chẽ theo công thức Birks:
$$L_p = \int_0^E \frac{S , dE}{1 + kB \left(\frac{dE}{dx}\right)}$$
Đóng góp lý thuyết cốt lõi của luận án thể hiện qua các mệnh đề được kiểm chứng bằng thực nghiệm:
- Mệnh đề 1: Sự phân bố thống kê của tham số PSD không đơn thuần tuân theo phân bố chuẩn đơn (single Gaussian) mà là một hàm kết hợp của năng lượng photon thu gom tại photocathode và nhiễu nhiệt dynode.
- Mệnh đề 2: Tồn tại một không gian vector chuẩn hóa trong đó sai số bình phương trung bình giữa xung thực nghiệm $MP(t)$ và xung tham khảo nơtron $RP_n(t)$ hoặc gamma $RP_g(t)$ là một đại lượng bất biến tương đối đối với sự thăng giáng biên độ ở vùng năng lượng thấp.
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 trúc lý thuyết:
- Lý thuyết dập tắt phát quang Birks (Scintillation Quenching Theory): Giải thích sự sai lệch biên độ quang năng giữa điện tử và ion nặng.
- Định lý lấy mẫu Nyquist-Shannon mở rộng trên mảng điện dung chuyển mạch DRS4: Cho phép tái tạo dạng sóng pico-giây mà không làm méo dạng sườn sau.
- Lý thuyết kiểm định giả thuyết thống kê trong không gian vector đa chiều (Statistical Decision Theory): Tối ưu hóa việc phân loại nhị phân (neutron vs gamma) và nhận diện xung chồng chập (pile-up rejection).
┌────────────────────────────────────────────────────────────┐
│ KHUNG PHÂN TÍCH ĐA TẦNG CỦA HỆ ĐO │
└─────────────────────────────┬──────────────────────────────┘
│
┌─────────────────────────────────────────────┼─────────────────────────────────────────────┐
▼ ▼ ▼
┌──────────────────────────┐ ┌──────────────────────────┐ ┌──────────────────────────┐
│ TẦNG 1: QUANG - ĐIỆN │ │ TẦNG 2: SỐ HÓA DRS4 │ │ TẦNG 3: THUẬT TOÁN DSP │
│ - Nhấp nháy EJ-301 │ ──────► │ - Mảng tụ điện ngắt │ ───────────────► │ - Trừ đường cơ bản │
│ - PMT Hamamatsu R9420 │ Tín hiệu│ - Tần số 1.0 - 2.5 GSPS │ Dữ liệu số hóa │ - Chuẩn hóa biên độ │
│ - TKĐ THS3202 (4 tầng) │ tương tự│ - Hiệu chuẩn điện áp │ 1024 điểm/xung │ - Đối sánh vector RP │
│ - Phân cực DY6-8 tối ưu │ │ và thời gian pico-giây│ │ - Loại trừ chồng chập │
└──────────────────────────┘ └──────────────────────────┘ └────────────┬─────────────┘
│
▼
┌──────────────────────────┐
│ TẦNG 4: PHÂN TÍCH VÀ PHỔ │
│ - Tách phổ Nơtron │
│ - Tách phổ Gamma │
│ - Đánh giá FoM đa dải │
└──────────────────────────┘
Điều kiện biên (boundary conditions) của khung phân tích được xác lập chặt chẽ: Dải năng lượng khảo sát từ 50 keVee đến 5000 keVee; dải tốc độ đếm lên tới $10^5\text{ cps}$; hệ số khuếch đại điện áp của TKĐ được chuẩn định tối ưu ở hai mức $K = 8.9$ và $K = 13.2$.
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ủ nghiêm ngặt thế giới quan thực chứng (positivism) kết hợp phương pháp luận thực nghiệm định lượng cao cấp. Thiết kế nghiên cứu bao gồm sự đồng thiết kế giữa kỹ thuật phần cứng (Hardware Co-design), thiết kế vi mạch lai tương tự/số, lập trình logic phần cứng FPGA và phát triển phần mềm xử lý dữ liệu chuyên dụng MCA_DRS4.
┌─────────────────────────────────────────────────────────────────────────────────────────────┐
│ SƠ ĐỒ PHẦN CỨNG CHI TIẾT CỦA HỆ ĐO │
└────────────────────────────────────────────────┬────────────────────────────────────────────┘
│
┌──────────────────────┐ ┌──────────┴──────────┐ ┌──────────────────┐
│ ĐẦU NHẤP NHÁY │ Quang │ PMT HAMAMATSU │ Tín hiệu │ TIỀN KHUẾCH ĐẠI │
│ EJ-301 LỎNG │ ──────────► │ R9420 FAST PMT │ ──────────► │ THS3202 OP-AMP │
│ - Thể tích 51,73 ml │ Kính UV │ - Cao thế -1500V │ Anode │ - 4 tầng khuếch │
│ - D = 34mm, L = 60mm│ Silicone │ - Mạch chia thế │ ghép trực │ đại & tạo xung│
│ - Vỏ nhôm đánh bóng │ grease │ RY=50Ω, C=10nF │ tiếp │ - trise > 10 ns │
└──────────────────────┘ └─────────────────────┘ └────────┬─────────┘
│
┌────────────────────────────────────────────┘ Cáp đồng trục 50Ω (>2m)
▼
┌─────────────────────────────────┐
│ BẢN MẠCH SỐ HÓA DRS4 V5 │
│ - Chip lấy mẫu DRS4 (1-2.5GSPS)│
│ - FPGA Altera điều khiển logic │
│ - Vi điều khiển truyền thông │
│ - Giao tiếp USB tốc độ cao │
└───────────────┬─────────────────┘
│
▼
┌─────────────────────────────────┐
│ PHẦN MỀM MÁY TÍNH MCA_DRS4 │
│ - Hiển thị trực tuyến phổ 2D/1D│
│ - Thực thi thuật toán RP/DCI │
│ - Tách phổ đồng thời n & gamma │
└─────────────────────────────────┘
Quy trình nghiên cứu rigorous
Quy trình chế tạo phần cứng và thiết lập thực nghiệm bao gồm các thông số kỹ thuật chính xác:
- Chế tạo cell chứa chất nhấp nháy EJ-301: Vỏ đầu nhấp nháy được gia công bằng nhôm tinh khiết dạng hình trụ đường kính trong 34 mm, chiều dài 60 mm. Mặt trong được đánh bóng quang học để đạt hệ số phản xạ cực đại. Cửa sổ quang sử dụng kính truyền qua tử ngoại (UV-glass) ghép nối quang học bằng mỡ silicone chuyên dụng. Thể tích chất lỏng được định lượng chính xác:
"thể tích nhấp nháy EJ-301 được bơm vào đầu nhấp nháy EJ-301 là $V = 3,14 \times 0,017^2 \times 0,06 \times 0,95 = 51,73\text{ ml}$"
(để lại 5% thể tích giãn nở nhiệt an toàn).
- Phân cực và phối hợp trở kháng PMT R9420: Ống nhân quang Hamamatsu R9420 có thời gian tăng siêu nhanh 1,6 ns, độ phân tán thời gian truyền FWHM đạt 550 ps, hệ số khuếch đại $5 \times 10^5$ tại điện áp $-1500\text{ V}$. Để triệt tiêu hoàn toàn dao động ký sinh do sườn tăng cực ngắn kích thích, mạch chia áp được cải tiến đặc biệt:
"Khảo sát các xung ra từ anode của PMT đã xác định được giá trị điện trở và tụ phân thế là: $R_Y = 50\ \Omega$, $C_1=C_2=C_3=10\text{ nF}$. Với các giá trị này, xung ra ít bị dao động do ảnh hưởng của sườn tăng của xung quá ngắn."
- Thiết kế mạch tiền khuếch đại (TKĐ) băng thông rộng: Sử dụng IC khuếch đại thuật toán siêu cao tần THS3202. Để loại bỏ hiện tượng méo tín hiệu trên đường truyền cáp đồng trục dài:
"hình thành lại xung có thời gian tăng của sườn trước lớn hơn 10 ns để tránh hiệu ứng phản xạ xung khi truyền trên cáp đồng trục có chiều dài trên 2 m."
- Số hóa và hiệu chuẩn hệ đo DRS4: Tín hiệu từ TKĐ được đưa vào bản mạch DRS4 V5 do PSI thiết kế. Bản mạch thực hiện lấy mẫu tương tự trên mảng 1024 tụ điện đóng ngắt với tốc độ cấu hình từ 1.0 GSPS đến 2.5 GSPS. Thuật toán hiệu chuẩn điện áp từng cell (voltage calibration) và hiệu chuẩn phi tuyến thời gian (time calibration) được thực hiện tự động trước mỗi phiên đo, đưa mức nhiễu nền RMS của hệ xuống dưới 1.2 mV.
┌──────────────────────────────────────────────────────────────┐
│ LƯU ĐỒ THUẬT TOÁN XUNG THAM KHẢO (RP) VÀ LỌC │
└──────────────────────────────┬───────────────────────────────┘
│
▼
┌──────────────────────────┐
│ Thu nhận xung số hóa │
│ MP(i), i=1..N │
└────────────┬─────────────┘
│
▼
┌──────────────────────────┐
│ Trừ đường cơ bản (BL) │
│ tính từ 30 mẫu đầu tiên │
└────────────┬─────────────┘
│
▼
┌──────────────────────────┐
│ Chuẩn hóa biên độ cực │
│ đại về 1.0 đơn vị │
└────────────┬─────────────┘
│
▼
┌──────────────────────────┐
│ Tính tổng sai số vector │
│ với RPn và RPg: │
│ Sn = Σ |MP(i) - RPn(i)| │
│ Sg = Σ |MP(i) - RPg(i)| │
└────────────┬─────────────┘
│
┌─────────────────────┴─────────────────────┐
▼ ▼
┌─────────────────────┐ ┌─────────────────────┐
│ Sn/Sg > Thresh_PU │ │ Phân loại PSD │
│ (Xung chồng chập) │ │ Tham số = Sn - Sg │
└──────────┬──────────┘ └──────────┬──────────┘
│ │
▼ ▼
┌─────────────────────┐ ┌─────────────────────┐
│ LOẠI BỎ SỰ KIỆN │ │ > 0: Phổ Gamma │
│ (Pile-up Rejection)│ │ < 0: Phổ Nơtron │
└─────────────────────┘ └─────────────────────┘
Data và phân tích
Số liệu thực nghiệm được thu thập từ các chiến dịch đo đa dạng:
- Nguồn chuẩn gamma đơn năng: $^{137}\text{Cs}$ ($662\text{ keV}$), $^{60}\text{Co}$ ($1.17\text{ MeV}$ và $1.33\text{ MeV}$), $^{22}\text{Na}$ ($511\text{ keV}$ và $1274\text{ keV}$) để chuẩn hóa thang năng lượng electron tương đương (keVee).
- Nguồn nơtron phân hạch tự phát $^{252}\text{Cf}$ đặt trong điều kiện có che chắn và đo phông chì/paraffin.
- Chùm nơtron chuẩn 148 keV trích xuất từ Kênh ngang số 4 của Lò phản ứng hạt nhân Đà Lạt.
Các thuật toán PSD tích hợp trong phần mềm MCA_DRS4 bao gồm:
- Thời gian vượt ngưỡng (TCT): Tính khoảng thời gian tín hiệu duy trì trên 10% đỉnh xung.
- Phân tích độ dốc (PGA): Tính tỷ số vi phân biên độ tại hai thời điểm $t_1 = t_{peak} + 15\text{ ns}$ và $t_2 = t_{peak} + 45\text{ ns}$.
- Tích phân điện tích số (DCI): Tính tỷ số $Q_{tail} / Q_{total}$ với khoảng tích phân đuôi từ $30\text{ ns}$ sau đỉnh đến hết cửa sổ xung $200\text{ ns}$.
- Nhận dạng mẫu (CPR): Tính góc tương quan $\theta = \arccos\left(\frac{\mathbf{X} \cdot \mathbf{Y}}{|\mathbf{X}| |\mathbf{Y}|}\right)$.
- Phương pháp xung tham khảo (RP): Tối ưu hóa ma trận chênh lệch dạng sóng chuẩn hóa đối với hai xung mẫu nơtron $RP_n$ và gamma $RP_g$.
Hệ số phẩm chất phân biệt (FoM) được tính toán định lượng cho từng dải năng lượng:
$$\text{FoM} = \frac{|ch_n - ch_\gamma|}{\text{FWHM}n + \text{FWHM}\gamma}$$
Trong đó, $ch_n, ch_\gamma$ là tọa độ đỉnh phân bố Gaussian của nơtron và gamma; $\text{FWHM}n, \text{FWHM}\gamma$ là độ rộng nửa chiều cao tương ứng của từng phân bố.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Vượt ngưỡng giới hạn phân biệt năng lượng thấp (50 - 200 keVee): Trong khi các phương pháp tiêu chuẩn DCI, PGA và TCT đều suy giảm nghiêm trọng khi năng lượng giảm dưới 200 keVee (với $\text{FoM}_{\text{DCI}} < 0.7$ tại 50 keVee, hai đỉnh nơtron và gamma chồng lấn hoàn toàn), phương pháp xung tham khảo (RP) duy trì $\text{FoM} \ge 1.02$ tại 50 keVee và đạt $\text{FoM} = 1.48$ tại 100 keVee. Điều này cho phép mở rộng dải phổ nơtron hữu ích xuống vùng năng lượng nhiệt hóa và trung gian thấp chưa từng đạt được trước đây trên các hệ đo nội địa.
- Đặc tính tuyến tính quang năng tuyệt vời của đetectơ EJ-301 tự chế tạo: Kết quả chuẩn năng lượng với các đỉnh tán xạ Compton của $^{137}\text{Cs}$ (477 keVee), $^{22}\text{Na}$ (341 keVee và 1062 keVee), và $^{60}\text{Co}$ (963 keVee và 1118 keVee) chứng minh hệ số tương quan tuyến tính $R^2 > 0.999$, khẳng định kết cấu quang học và mạch TKĐ THS3202 không gây méo phi tuyến trong dải từ 0 đến 1.5 V.
- Cơ chế loại trừ xung chồng chập thời gian thực chính xác: Dựa trên phân bố tham số nhận dạng bất thường của vector độ lệch chuẩn $S$, thuật toán RP nhận diện và loại bỏ tự động trên 98.5% các xung chồng chập ngẫu nhiên ở tốc độ đếm cao, ngăn chặn hoàn toàn hiện tượng tạo đỉnh giả trong phổ nơtron.
- Tách thành công vạch phổ nơtron 148 keV tại Lò phản ứng hạt nhân Đà Lạt: Thực nghiệm trên Kênh số 4 LPƯHN Đà Lạt trong điều kiện mở kênh và đóng kênh (shutter open/close) đã chứng minh khả năng bóc tách hoàn toàn phổ gamma nền cực mạnh của lò phản ứng, thu được phổ phân bố nơtron 148 keV rõ nét, phù hợp hoàn hảo với các tính toán mô phỏng Monte Carlo MCNP chuẩn tắc.
| Phương pháp PSD |
FoM tại 50 keVee |
FoM tại 100 keVee |
FoM tại 500 keVee |
FoM tại 1000 keVee |
Nhu cầu tài nguyên phần cứng |
| TCT (Thời gian vượt ngưỡng) |
0.32 |
0.54 |
0.91 |
1.15 |
Rất thấp |
| PGA (Độ dốc xung) |
0.28 |
0.48 |
0.88 |
1.20 |
Rất thấp |
| DCI (Tích phân điện tích) |
0.68 |
0.92 |
1.35 |
1.62 |
Thấp |
| CPR (Nhận dạng mẫu) |
0.72 |
1.05 |
1.42 |
1.70 |
Trung bình - Cao |
| RP (Luận án phát triển) |
1.02 |
1.48 |
1.76 |
2.05 |
Thấp - Trung bình |
Implications đa chiều
- Về mặt học thuật và lý thuyết: Chứng minh rằng việc khai thác toàn bộ hình học vector xung số hóa trong miền thời gian có thể đạt hiệu quả tương đương hoặc vượt trội so với các phép biến đổi miền tần số phức tạp (DFTM, Wavelet), thiết lập chuẩn so sánh mới cho các nghiên cứu phát triển thuật toán PSD.
- Về mặt kỹ thuật và phương pháp luận: Khẳng định tính khả thi của việc sử dụng các mảng lấy mẫu điện dung chuyển mạch tốc độ cao (DRS4) kết hợp vi mạch op-amp thương mại giá thành thấp để thay thế các khối ADC Flash đắt tiền và cồng kềnh trong thiết kế thiết bị hạt nhân.
- Về mặt ứng dụng thực tiễn: Tạo ra một thiết bị hoàn chỉnh có thể ứng dụng ngay lập tức trong giám sát an toàn bức xạ nơtron tại Viện Nghiên cứu Hạt nhân Đà Lạt, đo liều bức xạ hỗn hợp trong y tế xạ trị nơtron (BNCT), và kiểm tra vật liệu không phá hủy trong công nghiệp.
- Về mặt chính sách và tự chủ công nghệ: Đóng góp trực tiếp vào chiến lược nội địa hóa trang thiết bị bức xạ hạt nhân kỹ thuật số của Việt Nam, giảm sự phụ thuộc vào các khối NIM/VME nhập khẩu đắt đỏ từ các tập đoàn phương Tây.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án cũng thẳng thắn thừa nhận các giới hạn kỹ thuật:
- Tính chất dễ cháy và độc tính của chất nhấp nháy lỏng: EJ-301 là dung môi hữu cơ gốc xylene dễ bay hơi và dễ cháy, đòi hỏi quy trình gia công vỏ kim loại cực kỳ nghiêm ngặt và hạn chế khả năng ứng dụng trong các môi trường di động dã ngoại khắc nghiệt so với các chất nhấp nháy chất dẻo mới xuất hiện như EJ-299-33A.
- Độ nhạy nhiệt của ống nhân quang PMT: Hệ số nhân của PMT R9420 và điện áp phân cực dynode phụ thuộc vào nhiệt độ môi trường, dẫn đến hiện tượng trôi nhẹ đường cơ bản khi hoạt động liên tục nhiều giờ mà không có mạch bù nhiệt phần cứng chủ động.
- Giới hạn băng thông truyền dữ liệu USB: Bản mạch DRS4 V5 truyền dữ liệu số hóa về máy tính qua chuẩn USB 2.0, tạo ra nút thắt cổ chai về tốc độ đếm tối đa của toàn hệ thống khi hoạt động ở chế độ ghi nhận từng xung đầy đủ (list-mode).
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Hướng 1: Nghiên cứu tích hợp toàn bộ thuật toán RP vào lõi IP-Core trên FPGA (sử dụng ngôn ngữ VHDL/Verilog) để thực hiện tính toán PSD trực tiếp trên phần cứng mà không cần chuyển dữ liệu thô về PC, nâng tốc độ đếm lên $> 500\text{ kcps}$.
- Hướng 2: Ứng dụng đầu dò chất nhấp nháy dẻo PSD (như EJ-276, EJ-299) thay thế nhấp nháy lỏng nhằm tăng tính cơ động, độ bền cơ học và an toàn môi trường.
- Hướng 3: Mở rộng hệ đo thành mảng nhiều đầu dò (detector array) kết hợp kỹ thuật trùng phùng thời gian bay (ToF) để xây dựng hệ phổ kế tán xạ nơtron góc phục vụ nghiên cứu cấu trúc vật liệu tiên tiến trên lò phản ứng nghiên cứu mới.
Tác động và ảnh hưởng
- Ảnh hưởng học thuật: Luận án cung cấp hệ phương pháp luận thực nghiệm và bộ mã thuật toán chuẩn hóa, dự kiến đóng góp lớn vào các trích dẫn trong lĩnh vực thiết bị hạt nhân (Nuclear Instruments and Methods in Physics Research Section A; IEEE Transactions on Nuclear Science).
- Chuyển đổi công nghiệp & An ninh: Cung cấp giải pháp đo nơtron nhanh/nhiệt chính xác cao cho ngành thăm dò địa vật lý dầu khí (giếng khoan nơtron), kiểm tra an ninh biên giới phát hiện vật liệu phân hạch hạt nhân bất hợp pháp ($^{235}\text{U}, ^{239}\text{Pu}$).
- Lợi ích xã hội và y tế: Ứng dụng trực tiếp trong đo liều sinh học nơtron và bảo vệ bức xạ cho nhân viên y tế tại các cơ sở xạ trị gia tốc và lò phản ứng nghiên cứu.
Đối tượng hưởng lợi
- Nghiên cứu sinh và học viên cao học: Cung cấp tài liệu tham khảo mẫu mực về quy trình thiết kế đetectơ bức xạ từ cấp độ vật lý linh kiện đến lập trình DSP và phân tích phổ hạt nhân.
- Các nhà khoa học vật lý hạt nhân thực nghiệm: Sở hữu một công cụ đo phổ nơtron nhanh hiệu năng cao, chi phí thấp phục vụ các thí nghiệm tán xạ và phản ứng hạt nhân.
- Kỹ sư R&D thiết bị y tế và bức xạ: Nắm bắt được phương pháp ứng dụng chip DRS4 và thuật toán RP để chế tạo các máy phân tích đa kênh MCA chuyên dụng thế hệ mới.
- Cơ quan quản lý an toàn bức xạ và hạt nhân: Có thêm công cụ đo kiểm độc lập, tin cậy để thẩm định trường bức xạ hỗn hợp tại các cơ sở hạt nhân 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ì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc xây dựng Mô hình phân biệt dạng xung không gian vector tham khảo (Reference-Pulse Vector Discrimination Model), mở rộng lý thuyết động học phát quang Birks và lý thuyết dập tắt ion hóa trong môi trường nhấp nháy hữu cơ. Luận án đã chứng minh toán học và thực nghiệm rằng: sự khác biệt hình dạng giữa xung proton giật lùi và xung điện tử Compton được bảo toàn trọn vẹn dưới dạng độ lệch chuẩn hình học khi xung đo được chuẩn hóa biên độ cực đại và đối sánh vi phân trực tiếp với hai dạng sóng tham khảo trung bình $RP_n$ và $RP_g$. Điều này giải quyết bài toán suy giảm độ phân giải ở dải năng lượng thấp mà không cần viện dẫn đến các phép gần đúng tích phân đuôi cắt cụt của lý thuyết Brooks truyền thống.
2. Đột phá về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
So sánh với:
- Nghiên cứu của Safari et al. (2016) sử dụng phương pháp DFTM (miền tần số): Phương pháp RP của luận án vận hành trực tiếp trong miền thời gian, loại bỏ hoàn toàn các bước biến đổi ma trận phức tạp nhưng đạt hệ số phẩm chất $\text{FoM} = 1.48$ tại 100 keVee (vượt mức $\text{FoM} \approx 1.35$ của DFTM), đồng thời giảm độ phức tạp tính toán xuống hàng chục lần.
- Nghiên cứu của Kamada et al. (2017) sử dụng phương pháp CPR (nhận dạng mẫu góc cosin): Luận án đã khắc phục điểm yếu chí mạng của CPR ở vùng năng lượng thấp ($< 200\text{ keVee}$) bằng cách kết hợp thuật toán khử trôi đường cơ bản 30 mẫu đầu và tối ưu hóa ma trận trọng số đuôi xung, nâng FoM tại 50 keVee từ $0.72$ lên $1.02$.
3. Phát hiện bất ngờ nhất trong quá trình thực nghiệm kèm minh chứng dữ liệu?
Phát hiện bất ngờ nhất là sự suy giảm phi tuyến của độ phân biệt khi tăng hệ số khuếch đại của tiền khuếch đại THS3202 lên mức quá cao ($K = 13.2$) ở vùng năng lượng thấp. Về mặt trực giác, hệ số khuếch đại lớn sẽ làm tăng biên độ tín hiệu nhỏ của vùng 50 keVee, kỳ vọng làm tăng SNR. Tuy nhiên, dữ liệu thực nghiệm chỉ ra rằng ở $K = 13.2$, tỷ số tín hiệu trên nhiễu SNR thực tế bị suy giảm do băng thông của THS3202 bị co hẹp theo định luật Gain-Bandwidth Product, làm biến dạng sườn tăng và tăng mức nhiễu nền RMS từ 0.8 mV lên 1.9 mV. Mức khuếch đại $K = 8.9$ mới là điểm tối ưu tuyệt đối mang lại giá trị FoM cao nhất cho phương pháp RP.
4. Luận án có cung cấp đầy đủ quy trình lặp lại thực nghiệm (Replication Protocol) không?
Có. Luận án cung cấp quy trình công nghệ chi tiết từ bản vẽ cơ khí vỏ đetectơ (đường kính 34 mm, dài 60 mm, dung tích 51,73 ml), sơ đồ nguyên lý mạch chia thế PMT R9420 ($R_Y = 50\ \Omega, C_1=C_2=C_3=10\text{ nF}$), sơ đồ nguyên lý 4 tầng TKĐ THS3202, quy trình nạp mẫu EJ-301 trong môi trường trơ, lưu đồ thuật toán chi tiết từng bước của chương trình MCA_DRS4 và các tham số khớp đường cong Gaussian trên phần mềm phân tích.
5. Kế hoạch nghiên cứu 10 năm được vạch ra như thế nào?
Kế hoạch 10 năm phát triển công nghệ đo bức xạ nơtron được luận án định hình qua 3 giai đoạn:
- Năm 1 - 3: Hoàn thiện đóng gói lõi IP-Core thuật toán RP trên chip SoC FPGA (Zynq UltraScale+), chế tạo máy đo nơtron/gamma cầm tay thương mại tích hợp GPS và truyền thông không dây.
- Năm 4 - 6: Nghiên cứu mảng đetectơ nhấp nháy đa kênh (16 - 64 kênh) kết hợp chất nhấp nháy dẻo PSD thế hệ mới phục vụ hệ đo tán xạ nơtron góc tại Trung tâm Nghiên cứu Khoa học Công nghệ Hạt nhân mới (CNST).
- Năm 7 - 10: Phát triển hệ đo phổ liều nơtron thời gian thực trong y học xạ trị proton/nơtron và hệ thống chụp ảnh cắt lớp nơtron nhanh (Fast Neutron Tomography) cho công nghiệp hàng không vũ trụ.
Kết luận
- Chế tạo thành công đetectơ nhấp nháy lỏng EJ-301 hoàn chỉnh: Tối ưu hóa kích thước hình học 51,73 ml, ghép nối quang học hoàn hảo với PMT R9420 và mạch tiền khuếch đại siêu nhanh THS3202, tạo ra các xung điện áp có $t_{rise} > 10\text{ ns}$ loại trừ phản xạ cáp và bảo toàn đặc trưng phân rã hạt nhân.
- Phát triển thành công phương pháp phân biệt xung tham khảo (RP) đột phá: Xác lập kỷ lục về hiệu năng phân biệt ở vùng năng lượng thấp với $\text{FoM} = 1.02$ tại 50 keVee và $\text{FoM} = 1.48$ tại 100 keVee, vượt trội so với các phương pháp chuẩn quốc tế DCI, CPR, PGA và TCT.
- Xây dựng hệ đo số hóa đồng thời nơtron - gamma MCA_DRS4: Làm chủ công nghệ tích hợp chip lấy mẫu mảng tụ điện DRS4 tốc độ 1.0 - 2.5 GSPS với máy tính, tích hợp thuật toán khử chồng chập xung trực tuyến đạt hiệu suất loại bỏ $> 98.5%$.
- Kiểm chứng thực nghiệm thành công tại Lò phản ứng hạt nhân Đà Lạt: Tách và chuẩn hóa chính xác phổ nơtron lọc 148 keV tại Kênh số 4 trong trường bức xạ gamma nền công suất cao, khẳng định độ tin cậy tuyệt đối của hệ đo.
- Mở ra 3 nhánh nghiên cứu công nghệ bức xạ mũi nhọn: Thiết bị đo liều nơtron số hóa di động; Phổ kế nơtron thời gian bay đa kênh; và Hệ thống kiểm soát an ninh phóng xạ thông minh.
- Đóng góp di sản khoa học định lượng: Đặt nền móng kỹ thuật vững chắc cho việc tự chủ thiết kế, chế tạo các thiết bị ghi đo bức xạ hạt nhân kỹ thuật số cao cấp tại Việt Nam, phục vụ trực tiếp cho chương trình phát triển ứng dụng năng lượng nguyên tử vì mục đích hòa bình.
