Tổng quan về luận án

Nghiên cứu phân tích định lượng đồng thời các nguyên tố phóng xạ họ Actini, đặc biệt là Uran ($\text{U}$) và Thori ($\text{Th}$), giữ vị trí then chốt trong chu trình công nghệ nhiên liệu hạt nhân, chế tạo viên gốm $\text{UO}_2$ cho lò phản ứng, điều tra địa chất khoáng sản phóng xạ và quan trắc an toàn phóng xạ môi trường. Trong tự nhiên, Uran chiếm khoảng $4\cdot 10^{-4}%$ khối lượng vỏ Trái Đất, tồn tại chủ yếu trong các khoáng vật như uraninit ($\text{UO}_2$, hàm lượng $\text{U}$ từ $45-85%$), pitchblend ($\text{UO}_2\cdot n\text{H}_2\text{O}$, $55%$), autunit ($\text{Ca}(\text{UO}_2)_2(\text{PO}_4)_2\cdot n\text{H}_2\text{O}$, $45-55%$), uranophan ($50-53%$), cùng nồng độ vết trong nước sông ($5\cdot 10^{-8}%$) và nước đại dương. Trong khi đó, Thori tập trung chủ yếu trong sa khoáng monazit ($(\text{Th, Ce, Y, La})\text{PO}_4$, $10-30%\ \text{ThO}_2$), thorit ($\text{ThSiO}_4$, $81,5%\ \text{ThO}_2$), thorogummit ($24-58%$) và esgynit ($17%$). Sự phát triển mạnh mẽ của công nghiệp chế biến khoáng sản sa khoáng ven biển cùng nguy cơ phát tán bụi đạn uran nghèo (Depleted Uranium - $\text{DU}$) trong các hoạt động quân sự đặt ra yêu cầu cấp bách về việc xây dựng các quy trình phân tích đồng thời $\text{U}$ và $\text{Th}$ với độ nhạy siêu vết, độ chọn lọc tuyệt đối và khả năng giải quyết các hiệu ứng xen phủ phổ phức tạp.

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ sự hạn chế của các phương pháp phân tích kinh điển:

  1. Trong phương pháp trắc quang với thuốc thử hữu cơ kinh điển Arsenazo III, phức của $\text{U}(\text{VI})$ và $\text{Th}(\text{IV})$ có cực đại hấp thụ quang rất gần nhau, dẫn đến sự xen phủ phổ hấp thụ nghiêm trọng. Các quy trình truyền thống buộc phải chiết tách hóa học bằng dung môi hữu cơ độc hại như tributyl photphat (TBP), axit di-(2-ethylhexyl) photphoric (HDEHP), trioctylphosphin oxit (TOPO) hoặc nhựa trao đổi ion, làm kéo dài thời gian, tốn kém hóa chất và tích lũy sai số ngẫu nhiên.
  2. Phương pháp trắc quang đạo hàm (Derivative Spectrophotometry) tuy loại bỏ được bước tách pha nhưng tín hiệu mật độ quang đo tại điểm cắt zero (zero-crossing) không trùng với cực đại hấp thụ, làm suy giảm độ nhạy nghiêm trọng ở đạo hàm bậc cao.
  3. Phương pháp hồi quy bình phương tối thiểu riêng phần (Partial Least Squares - PLS) và phương pháp Vierordt chỉ xử lý tốt các hệ cộng tính tuyến tính lý tưởng, thất bại hoàn toàn khi xử lý các tương tác liên cấu tử phi tuyến và biến dạng phổ do hiệu ứng nền.
  4. Kỹ thuật khối phổ plasma cảm ứng (ICP-MS) tuy đạt giới hạn phát hiện siêu vết (ppt - ppq) nhưng đối mặt với hiệu ứng nền (matrix effect) và hiệu ứng nhớ (memory effect) khi nồng độ muối khoáng cao.

Nhằm giải quyết toàn diện các thách thức trên, luận án thiết lập hệ thống 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để giải mã và phân tách định lượng phổ hấp thụ quang học xen phủ hoàn toàn của phức $\text{U}$-Arsenazo III, $\text{Th}$-Arsenazo III và $\text{Zr}$-Arsenazo III mà không cần chiết tách hóa học?
  • RQ2: Mô hình mạng nơ-ron nhân tạo Perceptron đa lớp (MLP-ANN) với cấu trúc và thuật toán học nào có khả năng mô hình hóa chính xác các hàm số phi tuyến biểu diễn mối quan hệ giữa bước sóng, mật độ quang và nồng độ các actinide trong hệ mờ?
  • RQ3: Các thông số vận hành của hệ thống ICP-MS (thời gian lấy tích phân tín hiệu IT, nồng độ axit, số khối $m/z$) ảnh hưởng như thế nào đến độ chính xác và độ lặp lại khi phân tích $\text{U}$ và $\text{Th}$ ở nồng độ siêu lệch nhau?
  • RQ4: Kỹ thuật ICP-MS có thể phân biệt chính xác tỷ lệ đồng vị $^{235}\text{U}/^{238}\text{U}$ trong các mẫu uran tự nhiên và uran nghèo $\text{DU}$ phục vụ đánh giá nguồn gốc ô nhiễm hay không?
  • H1: Thuật toán giải chập phổ quang học và mạng MLP-ANN có thể nhận dạng các đặc trưng phi tuyến của phổ hấp thụ, cho phép định lượng trực tiếp $\text{U}$ và $\text{Th}$ với sai số căn quân phương (RMS) dưới $1,5%$ ngay cả khi có mặt nguyên tố cản trở mạnh là Zirconi ($\text{Zr}$).
  • H2: Việc tối ưu hóa ma trận trọng số trong mô hình mạng nơ-ron M6 sẽ loại bỏ hoàn toàn sự phụ thuộc vào điểm cắt zero của trắc quang đạo hàm.
  • H3: Ứng dụng kỹ thuật ICP-MS với phương pháp nội chuẩn và thời gian lấy tín hiệu tối ưu sẽ thiết lập đường chuẩn tuyến tính rộng từ $0,5\text{ ppb}$ đến $50\text{ ppm}$ với độ lệch chuẩn tương đối (RSD) $< 3%$.
  • H4: Phép đo phổ khối ion hóa nhiệt độ cao cho phép xác định trực tiếp tỷ lệ đồng vị $^{235}\text{U}/^{238}\text{U}$ với độ tin cậy tương đương các phương pháp hạt nhân chuyên sâu.

Khung lý thuyết của nghiên cứu tích hợp đa ngành giữa Hóa trắc lượng (Chemometrics), Trí tuệ nhân tạo (Artificial Intelligence), Hóa học phức chất Actini và Vật lý Plasma nhiệt độ cao. Đóng góp đột phá của luận án được lượng hóa qua việc xây dựng thành công phần mềm tách phổ chuyên dụng viết bằng ngôn ngữ Pascal và mô hình mạng nơ-ron lan truyền ngược tối ưu (Model M6) đánh giá trên tập dữ liệu $224$ mẫu thực nghiệm, đồng thời chuẩn hóa quy trình phân tích đồng thời $\text{U}$, $\text{Th}$ và tỷ lệ đồng vị trên hệ ICP-MS Agilent 7500a. Phạm vi nghiên cứu bao quát từ dung dịch chuẩn đơn, hỗn hợp phức đa cấu tử, mẫu giả lập nền khoáng sản phức tạp đến các mẫu thực tế gồm quặng sa khoáng, mẫu chuẩn môi trường và bột/viên gốm $\text{UO}_2$ sạch hạt nhân do Viện Công nghệ Xạ hiếm chế tạo.


Literature Review và Positioning

Nghiên cứu phân tích $\text{U}$ và $\text{Th}$ trong lịch sử hóa phân tích chia thành ba dòng chảy học thuật chính:

Dòng chảy thứ nhất tập trung vào các phương pháp phân tích hạt nhân phóng xạ. Kỹ thuật phân tích kích hoạt nơtron hóa phóng xạ (RNAA) và dụng cụ (INAA) được phát triển bởi Chai, Ebihara và Repinc (2002) cho phép xác định lượng siêu vết $\text{U}$ và $\text{Th}$ trong các mẫu địa chất và hóa vũ trụ với độ nhạy cao. Tuy nhiên, RNAA đòi hỏi nguồn nơtron từ lò phản ứng hạt nhân và quy trình xử lý phóng xạ sau chiếu xạ cực kỳ nghiêm ngặt. Phương pháp phổ kế alpha kết hợp chiết lỏng-lỏng và điện phân làm giàu do Zarki và cộng sự (2001) hay Ayranov và Krähenbühl (2003) nghiên cứu trên mẫu trầm tích biển sử dụng tác nhân chiết TOPO, HDEHP và đo phổ PERALS (Photon-Electron Rejecting Liquid Alpha Spectrometry) đòi hỏi thời gian đếm xung lên tới $500.000\text{ giây}$ đối với mẫu $0,5\text{ g}$ để đạt độ phát hiện $0,34 - 1,15\text{ Bq/kg}$ cho $\text{U}$ và $0,31 - 1,66\text{ Bq/kg}$ cho $\text{Th}$. Tương tự, phương pháp phổ kế gamma phân giải cao của Shriwastwa, Kumar và Raghunath (2004) dùng để phân tích không phá hủy viên oxit $\text{U}$-$\text{Th}$ bị giới hạn bởi độ nhạy khi đo các mẫu môi trường có hoạt độ riêng thấp.

Dòng chảy thứ hai khai thác kỹ thuật quang phổ phát xạ và khối phổ plasma cảm ứng (ICP-AES, ICP-MS). Điển hình là nghiên cứu của Fujino và cộng sự khi xác định $\text{U}$ và $\text{Th}$ trong quặng apatit Florida (Mỹ) bằng ICP-AES; các tác giả đã phải dùng dung môi Diisobutyl ketone và thuốc thử 1-Phenyl-3-methyl-4-trifluoroacetyl-5-pyrazolone để tách nền Canxi photphat, thu được kết quả $\text{U}$ là $103 \pm 2,2\text{ ppm}$ và $\text{Th}$ là $8,84 \pm 0,19\text{ ppm}$, hoàn toàn phù hợp với phép đo ICP-MS. Grinberg, Willie và Sturgeon (2005) xác định vết ppt $\text{U}$ và $\text{Th}$ trong chì siêu tinh khiết sau khi tách nền chì bằng sắc ký chiết dùng nhựa UTEVA. Để đánh giá ô nhiễm vũ khí, Ejnik, Squibb, McDiarmid và Centeno (2001) đã ứng dụng kỹ thuật ICP-DRC-MS (buồng phản ứng động học) để đo tỷ lệ $^{235}\text{U}/^{238}\text{U}$ trong nước tiểu bệnh nhân bị phơi nhiễm mảnh đạn $\text{DU}$.

Dòng chảy thứ ba là sự trỗi dậy của Hóa trắc lượng và Mạng nơ-ron nhân tạo (ANN) trong trắc quang đa cấu tử. Các công trình của Satyanarayana, Kamarajan và Manavalan (2004), Kazemi (2002) hay các nghiên cứu định lượng đồng thời $\text{Cu-Co-Ni}$ bằng thuốc thử PAR cho độ lệch chuẩn tương đối $\text{RSD}$ từ $0,1 - 0,2%$ đã khẳng định tiềm năng vượt trội của ANN so với phương pháp bình phương tối thiểu riêng phần (PLS) trong các hệ có độ tương tác quang phổ phi tuyến phức tạp.

                    TIẾN TRÌNH TIẾP CẬN HỌC THUẬT
┌────────────────────────────────────────────────────────────────────────┐
│ PHƯƠNG PHÁP HẠT NHÂN (INAA, RNAA, PERALS, Phổ Gamma)                   │
│ Nhược điểm: Cần lò phản ứng, thời gian đo dài (500.000s), xử lý mẫu khó│
└──────────────────────────────────┬─────────────────────────────────────┘
                                   ▼
┌────────────────────────────────────────────────────────────────────────┐
│ QUANG PHỔ TRẮC QUANG CỔ ĐIỂN & TRẮC QUANG ĐẠO HÀM (Arsenazo III)       │
│ Nhược điểm: Trùng lấn pic, zero-crossing suy giảm độ nhạy, PLS hạn chế │
└──────────────────────────────────┬─────────────────────────────────────┘
                                   ▼
┌────────────────────────────────────────────────────────────────────────┐
│ ĐỊNH VỊ CỦA LUẬN ÁN: TÍCH HỢP ĐA KỸ THUẬT ĐỘT PHÁ                      │
│ 1. Giải chập phổ quang học bằng phần mềm Pascal chuyên dụng            │
│ 2. Mạng nơ-ron nhân tạo MLP-ANN (Model M6) xử lý phi tuyến            │
│ 3. Khối phổ plasma ICP-MS Agilent 7500a: Tuyến tính rộng, đồng vị DU   │
└────────────────────────────────────────────────────────────────────────┘

Trong bức tranh tổng quan đó, tồn tại hai cuộc tranh luận học thuật lớn:

  1. Chiết tách hóa học bảo thủ đối lập với Xử lý toán học hiện đại: Một trường phái cho rằng cần tinh chế tuyệt đối mẫu qua hệ cột chiết để tránh sai số quang phổ; trường phái ngược lại khẳng định giải thuật toán học thông minh (ANN, deconvolution) có thể xử lý trực tiếp phổ hỗn hợp, giúp tinh giản quy trình xanh hóa phân tích.
  2. Phương pháp nền phù hợp (Matrix Matching) đối lập với Tách pha triệt để trong ICP-MS: Nền phù hợp tiết kiệm thời gian nhưng dễ gây hiệu ứng nhớ (memory effect) trên nón nạp mẫu (sampler/skimmer cones), đòi hỏi kiểm soát tối ưu thời gian tích phân tín hiệu.

Luận án định vị tiên phong tại giao điểm của hai hướng tiếp cận: vừa làm chủ kỹ thuật phần mềm mạng nơ-ron để giải quyết bài toán trắc quang chi phí thấp tại các phòng thí nghiệm thông thường, vừa tối ưu hóa quy trình ICP-MS hiện đại bậc nhất cho các phép đo đòi hỏi độ nhạy siêu vết và phân tích tỷ lệ đồng vị hạt nhân.


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

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

Luận án mang lại những đóng góp nền tảng cho lý thuyết hóa học phân tích và hóa học tính toán:

  1. Mở rộng Định luật Beer-Lambert trong không gian phi tuyến: Định luật cơ bản $A = \sum \varepsilon_i l C_i$ vốn chỉ áp dụng cho các hệ dung dịch lý tưởng, không có tương tác hóa học hay biến dạng quang sai. Luận án đã mở rộng mô hình này sang không gian hàm phi tuyến thông qua Định lý xấp xỉ phổ quát (Universal Approximation Theorem) của mạng Perceptron nhiều lớp, chứng minh rằng mạng ANN có thể tự học các sai lệch quang phổ gây ra bởi sự cạnh tranh tạo phức giữa $\text{UO}_2^{2+}$, $\text{Th}^{4+}$, $\text{Zr}^{4+}$ với Arsenazo III trong môi trường axit mạnh.
  2. Hiện thực hóa Lý thuyết ion hóa Plasma nhiệt động của Saha: Luận án vận dụng phương trình Saha để lượng hóa mức độ ion hóa trong ngọn lửa plasma Argon ở nhiệt độ $6800 - 8000\text{ K}$: $$\frac{n_i^2}{n_n} = 2,4 \cdot 10^{21} \cdot T^{3/2} \cdot \frac{U_i}{k T} \cdot \exp\left(-\frac{U_i}{k T}\right)$$ Trong đó $n_i, n_n$ lần lượt là mật độ nguyên tử bị ion hóa và trung hòa; $k$ là hằng số Boltzmann; $U_i$ là thế ion hóa thứ nhất của nguyên tố. Nghiên cứu đã chứng minh thực nghiệm rằng ở nhiệt độ plasma $6800\text{ K}$, hiệu suất ion hóa tạo ion đơn điện tích $+1$ của $\text{U}$ ($U_i = 6,19\text{ eV}$) và $\text{Th}$ ($U_i = 6,31\text{ eV}$) đạt gần như $100%$, thiết lập cơ sở lý thuyết cho việc đo đồng thời mà không bị suy giảm tín hiệu do phân bố ion hóa phân tán.
                  KHUNG KHÁI NIỆM VÀ MÔ HÌNH LÝ THUYẾT
┌────────────────────────────────────────────────────────────────────────┐
│                        TÍN HIỆU ĐẦU VÀO TOÀN PHẦN                      │
│   (Quang phổ hấp thụ UV-Vis đa cấu tử / Phổ dòng ion khối m/z ICP-MS)  │
└──────────────────────────────────┬─────────────────────────────────────┘
                                   │
         ┌─────────────────────────┴─────────────────────────┐
         ▼                                                   ▼
┌─────────────────────────────────┐ ┌─────────────────────────────────┐
│     KHUNG PHÂN TÍCH QUANG HỌC   │ │      KHUNG PHÂN TÍCH KHỐI PHỔ   │
│   • Phức Actini - Arsenazo III  │ │   • Khí hóa, Nguyên tử hóa,     │
│   • Phổ đạo hàm bậc 1, bậc 2    │ │     Ion hóa plasma Ar (8000 K)  │
│   • Thuật toán giải chập Pascal │ │   • Lọc khối tứ cực (Quadrupole)│
│   • Mạng MLP-ANN (Sigmoid, M6)  │ │   • Detector nhân electron kép  │
└────────────────┬────────────────┘ └────────────────┬────────────────┘
                 │                                   │
                 └─────────────────────────┬─────────┘
                                           ▼
┌────────────────────────────────────────────────────────────────────────┐
│                       KẾT QUẢ ĐẦU RA ĐA MỤC TIÊU                        │
│ 1. Hàm lượng U, Th siêu vết (ng/g - ppm) không cần chiết tách          │
│ 2. Tỷ lệ đồng vị chuẩn xác 235U/238U (Nhận diện Uranium nghèo DU)      │
│ 3. Độc lập với hiệu ứng nền matrix effect và nhiễu quang sai           │
└────────────────────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp chặt chẽ 3 trụ cột lý thuyết:

  • Lý thuyết quang phổ phức chất Actini: Bản chất tương tác phối trí giữa ion kim loại nặng và phối tử Arsenazo III sinh ra phổ hấp thụ có cực đại xen phủ tại vùng khả kiến ($600 - 700\text{ nm}$).
  • Lý thuyết Mạng nơ-ron nhân tạo Perceptron đa lớp (MLP): Sử dụng thuật giải lan truyền ngược sai số (Backpropagation) với hàm kích hoạt phi tuyến Sigmoid $f(x) = \frac{1}{1 + e^{-x}}$, kết hợp hằng số tốc độ học $\eta$ và hằng số động năng (momentum constant $\alpha$) nhằm tránh kẹt tại các điểm cực tiểu cục bộ (local minima).
  • Lý thuyết phân giải khối phổ vô cơ: Cơ chế lọc và hội tụ chùm ion đơn hóa trị $M^+$ qua hệ thấu kính ion và bộ phân giải khối tứ cực (Quadrupole Mass Filter) kết hợp đầu dò nhân electron chế độ kép (Dual Mode Electron Multiplier).

Các điều kiện biên (Boundary Conditions) được xác lập nghiêm ngặt: môi trường axit phân tích $\text{HCl}$ và $\text{HNO}_3$ tinh khiết cao, nồng độ thuốc thử Arsenazo III luôn đảm bảo dư tạo phức hoàn toàn, dải số khối quét chọn lọc tại $m/z = 232\ (\text{Th}), 235\ (\text{U}), 238\ (\text{U})$, loại bỏ hoàn toàn các cấu tử cản trở phổ khối đẳng khối (isobaric interferences).


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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học Thực chứng luận (Positivism) kết hợp Thực tại luận phản biện (Critical Realism), đảm bảo tính khách quan tuyệt đối thông qua dữ liệu thực nghiệm đo đạc định lượng. Thiết kế nghiên cứu đa mức (Multi-level Design) bao gồm:

  • Cấp độ 1 (Mô hình hóa toán học): Xây dựng thuật toán tách phổ viết bằng ngôn ngữ lập trình Pascal và thiết kế cấu trúc mạng nơ-ron MLP trên không gian ma trận trọng số.
  • Cấp độ 2 (Hệ mẫu chuẩn phòng thí nghiệm): Khảo sát hệ dung dịch chuẩn đơn nguyên tố và hệ hỗn hợp hai cấu tử ($\text{U-Th}$), ba cấu tử ($\text{U-Th-Zr}$) ở các tỷ lệ nồng độ biến thiên từ cân bằng đến siêu lệch.
  • Cấp độ 3 (Hệ mẫu giả lập đa nguyên tố): Kiểm chứng tính kháng nhiễu trên các dung dịch mẫu giả chứa đồng thời các ion kim loại chuyển tiếp và tạp chất nền phổ biến ($\text{Fe}^{3+}, \text{Al}^{3+}, \text{Mn}^{2+}, \text{Cr}^{3+}, \text{Cu}^{2+}, \text{Ni}^{2+}, \text{Zn}^{2+}, \text{Pb}^{2+}$).
  • Cấp độ 4 (Hệ mẫu thực địa và chuẩn quốc tế): Phân tích kiểm chứng trên quặng sa khoáng phóng xạ thực tế, mẫu địa chất chuẩn và viên gốm nhiên liệu $\text{UO}_2$.

Quy trình nghiên cứu rigorous

Thiết bị và hóa chất chuẩn độ tinh khiết cao:

  • Hệ thống khối phổ plasma cảm ứng ICP-MS Agilent 7500a điều khiển tự động qua phần mềm ChemStation, sử dụng khí Argon siêu tinh khiết $99,999%$.
  • Máy quang phổ kế hấp thụ phân tử UV-VIS JASCO V-530 kết nối máy tính qua phần mềm JASCO V-500/FP-750 for Windows.
  • Hóa chất tinh khiết phân tích (p.A., Merck): Axit $\text{HCl}\ (d=1,19)$, $\text{HNO}_3\ (d=1,41)$, $\text{HF}\ 47%$, thuốc thử tinh thể Arsenazo III, các dung dịch chuẩn gốc $1000\ \mu\text{g/ml}$ của $\text{U}, \text{Th}, \text{Zr}, \text{Fe}, \text{Mn}, \text{Cr}, \text{Al}, \text{Pb}, \text{Cu}, \text{Ni}, \text{Zn}$.
                 QUY TRÌNH THỰC NGHIỆM ĐA TUYẾN NGHIÊM NGẶT
┌────────────────────────────────────────────────────────────────────────┐
│ MẪU PHÂN TÍCH (Chuẩn đơn / Hỗn hợp U-Th-Zr / Sa khoáng / Mẫu gốm UO2)  │
└──────────────────────────────────┬─────────────────────────────────────┘
                                   │
         ┌─────────────────────────┴─────────────────────────┐
         ▼                                                   ▼
┌─────────────────────────────────┐ ┌─────────────────────────────────┐
│     TUYẾN 1: QUANG PHỔ HẤP THỤ  │ │     TUYẾN 2: KHỐI PHỔ ICP-MS    │
│  Tạo phức Arsenazo III (pH ~1)  │ │   Chuyển mẫu thành sol khí      │
│  Quét phổ UV-Vis (600 - 700 nm) │ │   Ion hóa ngọn lửa Plasma 8000 K│
│  Xử lý 2 nhánh thuật toán:      │ │   Tối ưu hóa thời gian đo IT    │
│   1. Giải chập phổ (Pascal)     │ │   Quét khối m/z = 232, 235, 238 │
│   2. Mạng nơ-ron Model M6       │ │   Nội chuẩn hóa loại trừ matrix │
└────────────────┬────────────────┘ └────────────────┬────────────────┘
                 │                                   │
                 └─────────────────────────┬─────────┘
                                           ▼
┌────────────────────────────────────────────────────────────────────────┐
│ ĐỐI SOÁNH ĐA PHƯƠNG PHÁP & ĐÁNH GIÁ CHUẨN ĐỘ TIN CẬY (RSD, RMS, p-val) │
└────────────────────────────────────────────────────────────────────────┘

Quy trình thu thập dữ liệu và chuẩn hóa độ tin cậy:

  • Thiết lập quy trình chuẩn bị mẫu khép kín, kiểm soát chặt chẽ thể tích sol khí dẫn vào buồng plasma, áp suất chân không buồng phân giải khối và độ ổn định nguồn phát cao tần RF.
  • Thực hiện kiểm chứng chéo tam giác (Triangulation) giữa ba kỹ thuật độc lập: Trắc quang đạo hàm, Tách phổ/Mạng ANN và Khối phổ ICP-MS.
  • Đánh giá độ tin cậy (Reliability) qua độ lệch chuẩn tương đối ($\text{RSD} < 3%$) và tính giá trị nội hàm (Internal Validity) qua hệ số tương quan tuyến tính $R^2 > 0,999$.

Data và phân tích

Phân tích dữ liệu thực nghiệm quy mô lớn:

  • Tập mẫu huấn luyện ANN: Khảo sát hệ dữ liệu trắc quang gồm $224$ mẫu hỗn hợp có tỷ lệ nồng độ $\text{U}$ và $\text{Th}$ biến thiên đa dạng (Bảng 3.12 và 3.13).
  • Phân tích tương quan: Đánh giá ma trận hệ số tương quan giữa hàm lượng $\text{U}, \text{Th}$ với mật độ quang và bước sóng quét (Bảng 3.6, 3.7, 3.8, 3.9) nhằm rút gọn số nút vào của lớp nhập (Input Layer) mà không làm mất thông tin phổ.
  • Tối ưu hóa cấu trúc mạng MLP: So sánh thực nghiệm các mô hình cấu trúc mạng nơ-ron (Bảng 3.10, Hình 3.24) để chọn ra mô hình tối ưu Model M6 (ma trận trọng số chi tiết tại Bảng 3.11 và 3.23) đạt sai số trung bình căn quân phương RMS cực tiểu và tốc độ hội tụ nhanh nhất.
  • Kiểm tra độ vững (Robustness Checks): Kiểm định mô hình mạng trên mẫu giả lập chứa nồng độ cao của các ion cản trở $\text{Zr}^{4+}$ (Bảng 3.15, 3.16, 3.26, 3.27) và các nguyên tố kim loại tạp chất (Bảng 3.21, 3.22, 3.31).

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

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

  1. Phát hiện 1: Giải mã thành công phổ xen phủ phức ba cấu tử $\text{U-Th-Zr}$ bằng phần mềm giải chập Pascal: Chương trình máy tính tách phổ tự lập trình đã phân tách trọn vẹn phổ hấp thụ riêng biệt của phức $\text{U}$ ($5 - 70\ \mu\text{g/25ml}$) và $\text{Th}$ ($20 - 60\ \mu\text{g/25ml}$) từ phổ hấp thụ tổng cộng (Hình 3.10 - 3.15). Ngay cả khi có mặt ion can thiệp mạnh $\text{Zr}$ tạo phức quang học cạnh tranh, thuật toán vẫn tách thành công phổ riêng rẽ của phức $\text{U}, \text{Th}, \text{Zr}$ với độ tái lập hình dạng phổ đạt trên $98%$ so với phổ đo đơn lẻ (Hình 3.18, 3.19).

  2. Phát hiện 2: Mô hình mạng nơ-ron M6 đạt độ chính xác đột phá trên 224 mẫu thực nghiệm: Mạng MLP-ANN cấu trúc tối ưu (Model M6) áp dụng hàm kích hoạt Sigmoid đã chứng minh khả năng định lượng trực tiếp $\text{U}$ và $\text{Th}$ trong 224 mẫu hỗn hợp phức tạp với sai số tương đối cực nhỏ, khắc phục hoàn toàn hiện tượng trôi điểm không và suy giảm độ nhạy của trắc quang đạo hàm (Bảng 3.12, 3.13). Kết quả xác định $\text{U}$ và $\text{Th}$ trong mẫu giả lập đa nguyên tố đạt độ đúng cao, không bị sai lệch bởi nền matrix.

  3. Phát hiện 3: Tối ưu hóa điều kiện đo và thời gian tích phân tín hiệu (IT) trên ICP-MS: Khảo sát động học chùm ion trên hệ Agilent 7500a cho thấy việc lựa chọn thời gian lấy tín hiệu (Integration Time - IT) có tính quyết định đối với độ lệch chuẩn tương đối khi đo nồng độ siêu vết $50\text{ ppt}$ (Hình 3.32). Tăng IT lên mức tối ưu giúp triệt tiêu xung nhiễu thống kê, đảm bảo độ ổn định của detector chế độ xung (Pulse Mode).

  4. Phát hiện 4: Thiết lập đường chuẩn tuyến tính siêu rộng từ $0,5\text{ ppb}$ đến $50\text{ ppm}$: Khảo sát thực nghiệm xác lập khoảng tuyến tính rộng chưa từng có (gấp $10^5$ lần) cho cả $\text{U}$ và $\text{Th}$ trong môi trường axit $\text{HNO}_3$ tối ưu (Hình 3.34, 3.35). Đánh giá ảnh hưởng tương hỗ chỉ ra rằng $\text{Th}$ không gây ảnh hưởng đến phép đo $\text{U}$ và ngược lại ngay cả khi tỷ lệ nồng độ cấu tử này vượt cấu tử kia hàng nghìn lần (Hình 3.36, 3.37, Bảng 3.29, 3.30).

  5. Phát hiện 5: Nhận diện chính xác tỷ lệ đồng vị $^{235}\text{U}/^{238}\text{U}$ và nguồn gốc Uranium nghèo ($\text{DU}$): Phương pháp ICP-MS cho kết quả phân tích tỷ lệ đồng vị trong các mẫu uran tự nhiên (bột và viên gốm $\text{UO}_2$) hoàn toàn trùng khớp với giá trị chuẩn lý thuyết ($^{235}\text{U}/^{238}\text{U} \approx 0,00725$, tương đương $0,72%$) (Bảng 3.32). Đặc biệt, nghiên cứu đã định lượng chính xác các mẫu uran nghèo $\text{DU}$ có tỷ lệ $^{235}\text{U}/^{238}\text{U} < 0,003$ (Bảng 3.33), cung cấp công cụ nhận diện ô nhiễm phóng xạ quân sự.

                    BẢNG TỔNG HỢP SO SÁNH HIỆU NĂNG ĐA KỸ THUẬT
┌────────────────────────┬───────────────────────┬──────────────────────┬──────────────────────┐
│ Tiêu chí đánh giá      │ Trắc quang đạo hàm    │ Mạng nơ-ron MLP (M6) │ Khối phổ ICP-MS      │
├────────────────────────┼───────────────────────┼──────────────────────┼──────────────────────┤
│ Dải nồng độ phân tích  │ 1 - 50 ppm            │ 0,2 - 80 ppm         │ 0,5 ppb - 50 ppm     │
│ Giới hạn phát hiện LOD │ ~ 0,5 ppm             │ ~ 0,05 ppm           │ Cỡ ppt - ppq         │
│ Thời gian phân tích    │ Nhanh (vài phút)      │ Rất nhanh (giây)     │ Rất nhanh, đa mẫu    │
│ Tách pha hóa học       │ Không yêu cầu         │ Không yêu cầu        │ Không yêu cầu        │
│ Phân tích tỷ lệ đồng vị│ Không thể             │ Không thể            │ Rất chính xác (DU)   │
│ Khả năng kháng nhiễu Zr│ Kém (lệch điểm zero)  │ Rất cao (RMS < 1,5%) │ Tuyệt đối (tách m/z) │
└────────────────────────┴───────────────────────┴──────────────────────┴──────────────────────┘

Implications đa chiều

  • Về lý thuyết: Mở ra hướng tiếp cận Hóa trắc lượng phi tuyến sâu, chứng minh mạng nơ-ron là công cụ mô hình hóa tối thượng cho các phản ứng phối trí actinide phức tạp.
  • Về phương pháp luận: Cung cấp giải pháp kỹ thuật kép linh hoạt: phương pháp Trắc quang - ANN chi phí thấp cho các phòng kiểm nghiệm tiêu chuẩn và phương pháp ICP-MS độ nhạy cao cho các viện nghiên cứu hạt nhân chuyên sâu.
  • Về ứng dụng thực tiễn: Quy trình kiểm soát chất lượng trực tiếp viên gốm $\text{UO}_2$ phục vụ chế tạo nhiên liệu lò phản ứng hạt nhân nghiên cứu Đà Lạt và các dự án năng lượng nguyên tử tương lai.
  • Về chính sách và môi trường: Cung cấp dữ liệu quan trắc nền phóng xạ tự nhiên tại các vùng mỏ sa khoáng ven biển (Bình Thuận, Thừa Thiên Huế) và xây dựng quy chuẩn nhận diện ô nhiễm bụi phóng xạ uran nghèo $\text{DU}$.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn biên:

  1. Hiệu ứng nhớ ICP-MS ở nồng độ cao: Khi nồng độ nền kim loại tổng vượt quá ngưỡng cho phép, hiện tượng bám dính sol khí trên buồng phun sương và nón nạp mẫu gây hiệu ứng nhớ, đòi hỏi thời gian rửa giải kéo dài giữa các lượt đo.
  2. Phụ thuộc miền dữ liệu huấn luyện của mạng ANN: Ma trận trọng số của Model M6 được huấn luyện trên hệ thuốc thử Arsenazo III; khi thay đổi thuốc thử tạo phức hoặc mở rộng sang các hệ dung môi hữu cơ mới, mạng cần được tái cấu trúc và huấn luyện lại.
  3. Giới hạn pha mẫu rắn: Quy trình ICP-MS hiện tại vẫn yêu cầu công đoạn phá mẫu vô cơ hóa bằng hỗn hợp axit $\text{HF-HNO}_3$, chưa tích hợp kỹ thuật bắn phá mẫu trực tiếp.

Chương trình nghiên cứu tương lai định hình 4 hướng phát triển:

  • Hướng 1: Mở rộng ứng dụng kiến trúc mạng nơ-ron tích chập (CNN) và học sâu (Deep Learning) để tự động nhận dạng và khử nhiễu toàn phổ khối ICP-MS đa nguyên tố.
  • Hướng 2: Tích hợp hệ thống phân tích dòng chảy nano kết hợp khối phổ cảm ứng cung từ hội tụ đúp (nFI-ICP-SFMS) nhằm hạ giới hạn phát hiện xuống dải femtomol/L đối với các actinide nhân tạo siêu vết ($^{236}\text{U}, ^{237}\text{Np}, ^{239}\text{Pu}$).
  • Hướng 3: Phát triển kỹ thuật quang phổ plasma laser phá hủy kết hợp khối phổ (LA-ICP-MS) phân tích vi cấu trúc trực tiếp trên bề mặt viên gốm nhiên liệu $\text{UO}_2$ mà không cần hòa tan mẫu.
  • Hướng 4: Xây dựng cơ sở dữ liệu tỷ lệ đồng vị phóng xạ môi trường quốc gia phục vụ thanh sát hạt nhân quốc tế.

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

Nghiên cứu tạo lập tầm ảnh hưởng sâu rộng trên nhiều bình diện:

  • Ảnh hưởng học thuật: Công trình đóng góp hệ thống dữ liệu thực nghiệm mẫu mực cho chuyên ngành Hóa phân tích và Hóa học Hạt nhân, tạo tiềm năng trích dẫn cao trên các tạp chí quốc tế chuyên ngành (như Talanta, Analytica Chimica Acta, Journal of Radioanalytical and Nuclear Chemistry).
  • Chuyển dịch công nghiệp chế biến khoáng sản: Giúp các doanh nghiệp khai thác sa khoáng ven biển (titan, ilmenit, zircon) giám sát tức thời hàm lượng $\text{U}$ và $\text{Th}$ đi kèm trong quặng monazit, tối ưu hóa công nghệ thu hồi và giảm thiểu phát thải chất thải nguy hại phóng xạ.
  • Nâng cao năng lực thanh sát và an toàn quốc gia: Cung cấp công cụ định lượng tỷ lệ $^{235}\text{U}/^{238}\text{U}$ độc lập, phục vụ công tác kiểm chứng an toàn hạt nhân và phát hiện sự xâm nhập của các nguồn vật liệu phóng xạ bất hợp pháp hoặc đạn uran nghèo.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một thiết kế nghiên cứu mẫu mực về việc kết hợp giữa trí tuệ nhân tạo (Hóa trắc lượng ANN) và phân tích công cụ hiện đại (ICP-MS, UV-Vis).
  • Chuyên gia R&D Nhiên liệu Hạt nhân: Sở hữu quy trình chuẩn xác để đánh giá độ sạch hạt nhân và tỷ lệ làm giàu đồng vị của các sản phẩm $\text{UO}_2$, $\text{U}_3\text{O}_8$.
  • Cơ quan Quản lý Môi trường và An toàn Bức xạ: Có cơ sở khoa học vững chắc để ban hành quy chuẩn kiểm định hàm lượng $\text{U}$-$\text{Th}$ trong sa khoáng và cảnh báo ô nhiễm phóng xạ nguồn nước.
  • Các phòng thí nghiệm phân tích địa chất - khoáng sản: Ứng dụng giải pháp phần mềm tách phổ và ANN để nâng cao năng lực phân tích trên các máy quang phổ sẵn có mà không cần đầu tư thiết bị mới đắt tiền.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Định luật hấp thụ quang Beer-Lambert sang không gian phi tuyến thông qua việc tích hợp mạng Perceptron đa lớp (MLP-ANN). Luận án chứng minh rằng mạng ANN với hàm kích hoạt Sigmoid có khả năng học và bù trừ tự động các tương tác phi tuyến, sự cạnh tranh tạo phức và biến dạng phổ quang học của hệ đa Actini ($\text{U-Th-Zr}$) với Arsenazo III, giải quyết triệt để vấn đề mà các mô hình tuyến tính cổ điển (Vierordt, PLS) bất lực.

2. Tính đổi mới trong phương pháp luận thể hiện như thế nào khi so sánh với các nghiên cứu quốc tế?

So với nghiên cứu của Fujino và cộng sự (phải dùng dung môi Diisobutyl ketone chiết phức tạp để tách nền Canxi khi đo apatit) hay nghiên cứu của Grinberg, Willie và Sturgeon (dùng nhựa UTEVA tách chì siêu tinh khiết), luận án đã đổi mới bằng cách:

  • Phát triển phần mềm giải chập Pascal và mạng nơ-ron Model M6 đo trực tiếp phổ hỗn hợp không cần chiết tách pha hữu cơ, rút ngắn $80%$ thời gian xử lý mẫu.
  • Tối ưu hóa phương pháp nền phù hợp kết hợp nội chuẩn trên ICP-MS Agilent 7500a, xác định đồng thời $\text{U}$ và $\text{Th}$ trong dải tuyến tính siêu rộng ($0,5\text{ ppb} - 50\text{ ppm}$) với độ lệch chuẩn $\text{RSD} < 3%$.

3. Phát hiện thực nghiệm nào 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à tính bất biến quang phổ của mạng ANN trước sự có mặt của Zirconi ($\text{Zr}$) ở nồng độ cao. Thông thường, $\text{Zr}$ tạo phức cực bền với Arsenazo III gây biến dạng hoàn toàn pic hấp thụ của $\text{U}$ và $\text{Th}$. Tuy nhiên, mạng nơ-ron Model M6 (ma trận trọng số Bảng 3.23) vẫn tách và định lượng chính xác $\text{U}$ và $\text{Th}$ trong hỗn hợp $\text{U-Th-Zr}$ (Bảng 3.26, 3.27) với độ thu hồi đạt $98,5 - 101,2%$, tương đương phép đo trên dung dịch đơn nguyên tố.

4. Luận án có cung cấp đầy đủ giao thức để lặp lại nghiên cứu (Replication Protocol) không?

Luận án cung cấp giao thức thực nghiệm chi tiết và minh bạch tuyệt đối:

  • Công bố đầy đủ cấu trúc mạng, thông số học $\eta$, hằng số động năng $\alpha$ và toàn bộ ma trận trọng số liên kết của mạng Model M6 (Bảng 3.11, 3.23).
  • Chi tiết hóa sơ đồ khối thuật toán giải chập phổ lập trình Pascal (Hình 3.6, 3.21).
  • Mô tả chính xác điều kiện vận hành máy ICP-MS: công suất RF, lưu lượng khí nạp sol khí, các giá trị thời gian tích phân tín hiệu IT và nồng độ môi trường axit $\text{HNO}_3$ tối ưu (Hình 3.32, 3.33).

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghị sự 10 năm tập trung vào 3 mũi nhọn:

  1. Chuẩn hóa bộ công cụ phần mềm Chemometrics-AI tích hợp trực tiếp vào vi điều khiển của các máy quang phổ thế hệ mới.
  2. Ứng dụng kỹ thuật khối phổ plasma laser (LA-ICP-MS) và dòng chảy thể tích nano (nFI-ICP-SFMS) để lập bản đồ phân bố actinide siêu vết ở độ phân giải không gian micro-met.
  3. Hoàn thiện mạng lưới giám sát tỷ lệ đồng vị $^{235}\text{U}/^{238}\text{U}$ phục vụ thanh sát hạt nhân và an ninh môi trường quốc gia.

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc toàn bộ mục tiêu khoa học đặt ra với 5 đóng góp cốt lõi:

  1. Thiết lập thành công cơ sở lý thuyết và phần mềm chuyên dụng giải chập phổ quang học trên ngôn ngữ Pascal, phân tách trọn vẹn phổ hấp thụ riêng phần của $\text{U}$ và $\text{Th}$ từ phổ hỗn hợp đa cấu tử.
  2. Xây dựng và tối ưu hóa mô hình mạng nơ-ron nhân tạo MLP-ANN (Model M6) xử lý triệt để bài toán định lượng đồng thời $\text{U}$ và $\text{Th}$ trong các hệ mẫu phức tạp có mặt nguyên tố cản trở $\text{Zr}$ mà không cần chiết tách hóa học.
  3. Chuẩn hóa quy trình phân tích đồng thời $\text{U}$ và $\text{Th}$ bằng ICP-MS Agilent 7500a trên dải tuyến tính siêu rộng từ $0,5\text{ ppb}$ đến $50\text{ ppm}$, kiểm soát chặt chẽ các yếu tố ảnh hưởng nền và thời gian tích phân tín hiệu IT.
  4. Xác định chuẩn xác tỷ lệ đồng vị $^{235}\text{U}/^{238}\text{U}$ trong các mẫu bột và viên gốm $\text{UO}_2$ sạch hạt nhân, đồng thời nhận diện chính xác các mẫu uran nghèo $\text{DU}$ phục vụ cảnh báo ô nhiễm phóng xạ.
  5. Kiểm chứng đồng bộ và đối so sánh thành công các phương pháp nghiên cứu trên hệ mẫu chuẩn quặng phóng xạ, mẫu môi trường và mẫu giả lập đa nguyên tố.

Công trình đánh dấu bước chuyển đổi mô hình (Paradigm Shift) quan trọng từ phương pháp phân tích hóa tách chiết cổ điển sang phương pháp tích hợp Trí tuệ nhân tạo (Hóa trắc lượng) và Khối phổ Plasma hiện đại, mở ra 3 hướng nghiên cứu mũi nhọn mới cho ngành Hóa phân tích Hạt nhân tại Việt Nam và xác lập vị thế học thuật vững chắc trên trường quốc tế.