Tổng quan nghiên cứu

Antimon (ký hiệu hóa học là Sb) là một trong những nguyên tố á kim độc hại hàng đầu, được Liên minh châu Âu và Cơ quan Bảo vệ Môi trường Hoa Kỳ đưa vào danh sách các chất ô nhiễm nguy hại cần kiểm soát nghiêm ngặt theo Công ước Basel. Trong vỏ Trái Đất, hàm lượng Antimon trung bình ước tính đạt khoảng 0,2 ppm, xếp thứ năm trong số các nguyên tố hiếm. Tại các đô thị có mật độ giao thông và hoạt động công nghiệp cao, nồng độ Antimon trong không khí dao động từ 0,6 đến 32 ng/m³, trong khi tại các khu vực tự nhiên vùng cao chỉ ghi nhận khoảng 0,2 ng/m³. Mức độ độc hại sinh học của Antimon phụ thuộc trực tiếp vào trạng thái oxy hóa: dạng Sb(III) vô cơ có độc tính cao gấp 10 lần so với dạng Sb(V) vô cơ, có khả năng gây tổn thương nghiêm trọng đến cơ quan hô hấp, tim mạch, chức năng gan và thận khi xâm nhập vào cơ thể. Trong các nguồn nước tự nhiên, Antimon thường xuất hiện ở nồng độ từ 0,1 đến 0,2 µg/L, trong khi tiêu chuẩn nước uống quốc tế giới hạn nghiêm ngặt dưới ngưỡng 5 µg/L.

Việc định lượng chính xác từng dạng tồn tại riêng biệt của Antimon là cơ sở tiên quyết để đánh giá đúng rủi ro sinh thái và bảo vệ sức khỏe cộng đồng. Mặc dù các kỹ thuật hiện đại như sắc ký lỏng hiệu năng cao ghép nối khối phổ plasma cảm ứng (HPLC-ICP-MS) hay huỳnh quang nguyên tử (HPLC-AFS) cho độ nhạy cao, nhưng chi phí thiết bị đắt đỏ và quy trình vận hành phức tạp đã hạn chế khả năng ứng dụng tại phần lớn các phòng thí nghiệm ở Việt Nam. Luận văn thạc sĩ chuyên ngành Hóa phân tích tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội đã phát triển thành công phương pháp phổ hấp thụ nguyên tử sau khi hidrua hóa (HG-AAS) kết hợp công cụ toán học Chemometrics. Công trình giải quyết triệt để bài toán phân tích đồng thời Sb(III) và Sb(V) vô cơ mà không cần bước tách pha sắc ký, giúp cắt giảm hơn 60% chi phí phân tích, rút ngắn thời gian chuẩn bị mẫu và đạt hệ số tương quan tuyến tính R² trên 0,999.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự tích hợp giữa lý thuyết phổ hấp thụ nguyên tử ngọn lửa kết hợp kỹ thuật tạo hidrua dễ bay hơi (HG-AAS) và lý thuyết hóa tin học đa biến (Chemometrics). Kỹ thuật HG-AAS dựa trên nguyên lý sử dụng chất khử mạnh Natri bohidrua (NaBH4) trong môi trường axit để chuyển hóa các dạng Antimon thành khí stibin (SbH3). Dòng khí mang Argon sẽ vận chuyển SbH3 vào cuvet thạch anh đặt trong chùm bức xạ cộng hưởng phát ra từ đèn catot rỗng ở bước sóng phân tích đặc trưng 217,53 nm để đo độ hấp thụ quang.

Lý thuyết Chemometrics cung cấp nền tảng toán học giải mã các tín hiệu quang phổ hỗn hợp có độ chồng lấp cao mà không cần tách hóa học sơ bộ. Hai mô hình hồi quy đa biến cốt lõi được ứng dụng gồm:

  • Phương pháp Bình phương tối thiểu nghịch đảo (ILS): Xây dựng ma trận hệ số hồi quy P thông qua công thức đại số ma trận dựa trên độ hấp thụ quang và nồng độ mẫu chuẩn.
  • Phương pháp Hồi quy cấu tử chính (PCR): Sử dụng thuật toán phân tích giá trị kỳ dị (SVD) để biến đổi tập dữ liệu độ hấp thụ quang sang không gian các cấu tử chính (PC), loại trừ hiện tượng đa cộng tuyến giữa các biến đo và triệt tiêu sai số ngẫu nhiên.

Ba khái niệm lý thuyết trung tâm xuyên suốt công trình gồm: phân tích dạng hóa học (chemical speciation), động học khử chọn lọc phụ thuộc nồng độ ion H+, và tính chất cộng tính của tín hiệu quang phổ hấp thụ nguyên tử.

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

Hệ thống thiết bị thực nghiệm gồm máy quang phổ hấp thụ nguyên tử Shimadzu AA-6800 ghép nối bộ tạo hidrua hóa HG, cân phân tích Scientech SA 210 có độ chính xác 0,0001 g và máy đo HANNA Instrument 211 có dải đo pH từ -2,00 đến 16,00. Cỡ mẫu nghiên cứu được thiết kế với ma trận hiệu chuẩn gồm 20 đến 25 dung dịch hỗn hợp chuẩn nhân tạo chứa đồng thời Sb(III) và Sb(V) ở các dải nồng độ biến thiên có kiểm soát. Kế hoạch thực nghiệm được triển khai theo mô hình quy hoạch hóa thực nghiệm bậc hai đầy đủ với 30 thí nghiệm độc lập nhằm khảo sát đồng thời 3 thông số: nồng độ NaBH4, tốc độ dòng chất khử và tốc độ bơm dòng mẫu. Độ lặp lại và giới hạn phát hiện (LOD/LOQ) được đánh giá qua 10 phép đo lặp mẫu trắng.

Nguồn dữ liệu thực tế được kiểm chứng trên các mẫu nước tự nhiên (nước mặt, nước ngầm) và mẫu đất thu thập tại các khu vực mỏ khoáng sản. Mẫu đất được xử lý vô cơ hóa hoàn toàn bằng dung dịch HCl 6M và nước cường thủy. Lý do lựa chọn phương pháp HG-AAS kết hợp Chemometrics là vì kỹ thuật hidrua hóa chuyển pha Antimon thành thể khí SbH3, giúp loại bỏ hoàn toàn các cation kim loại nặng gây nhiễu phổ từ nền mẫu; đồng thời các thuật toán đa biến khai thác hiệu quả sự sai khác về hiệu suất khử của Sb(III) và Sb(V) trong môi trường axit HCl từ 1M đến 6M để định lượng chính xác từng dạng cấu tử.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình khảo sát đơn biến và tối ưu hóa hệ thống thiết bị AAS đã xác lập các thông số vận hành máy đạt hiệu suất tối đa: vạch phổ hấp thụ đạt độ chọn lọc cao nhất tại bước sóng 217,53 nm; cường độ dòng đèn catot rỗng HCL tối ưu là 8 mA (đạt khoảng 80% dòng giới hạn tối đa để duy trì tuổi thọ đèn); chiều cao đèn nguyên tử hóa mẫu đạt 14 mm với độ lệch chuẩn tương đối RSD chỉ 2,59%; lưu lượng khí cháy C2H2 tối ưu ở mức 1,8 L/phút trên nền không khí 8 L/phút; bề rộng khe đo quang học 0,2 mm mang lại độ phân giải phổ cao nhất.

Đối với điều kiện phản ứng hidrua hóa, mô hình thực nghiệm bậc hai đầy đủ đã thiết lập phương trình toán học mô tả độ hấp thụ quang Abs với hệ số xác định R² đạt 99,98% và chuẩn sai mô hình S = 0,0015. Phân tích phương sai ANOVA khẳng định cả ba yếu tố công nghệ đều có ý nghĩa thống kê vượt trội với giá trị P bằng 0,000. Bộ thông số hóa học tối ưu gồm: dung dịch NaBH4 nồng độ 0,5% pha trong NaOH 0,2%; tốc độ dòng NaBH4 đạt 2,0 mL/phút; tốc độ dòng mẫu đạt 5,0 mL/phút trong môi trường axit HCl 6M. Khoảng tuyến tính phân tích xác định Sb(III) đạt từ 0,25 đến 20 ppb với phương trình đường chuẩn có độ dốc 0,05679. Khi kiểm tra trên mẫu chuẩn chứng nhận CRM và mẫu thử nhân tạo, cả hai mô hình ILS và PCR đều đạt hiệu suất thu hồi từ 95,2% đến 103,8% với sai số dự báo thực nghiệm dưới 4,5%.

Thảo luận kết quả

Cơ chế phân biệt hai dạng Antimon vô cơ nằm ở sự khác biệt bản chất về tốc độ phản ứng khử giữa ion trung tâm và ion hidrua: Sb(III) bị khử rất nhanh thành stibin ngay cả trong môi trường axit loãng HCl 1M, trong khi Sb(V) có tốc độ chuyển hóa cực kỳ chậm và chỉ đạt hiệu suất tạo hidrua đáng kể khi nồng độ axit HCl nâng lên 6M. Tận dụng sự chênh lệch tín hiệu đo ở các nồng độ axit khác nhau kết hợp với tính chất cộng tính của độ hấp thụ quang, thuật toán Chemometrics đã xác định đồng thời nồng độ Sb(III) và Sb(V) mà không cần sử dụng các tác nhân khử trung gian như KI hay L-cystein, qua đó loại trừ được nguy cơ tạp nhiễm và giảm thiểu sai số thao tác.

Dữ liệu tối ưu hóa thực nghiệm được minh họa trực quan thông qua đồ thị mặt mục tiêu ba chiều (3D Response Surface) và biểu đồ đường đồng mức hai chiều (2D Contour Plot). Các biểu đồ này mô tả rõ ràng vùng cực trị của độ hấp thụ quang theo nồng độ NaBH4 và tốc độ dòng bơm mẫu, giúp người phân tích nhanh chóng nhận diện điểm vận hành tối ưu của hệ thống. Bên cạnh đó, đồ thị phần trăm phương sai tích lũy của các cấu tử chính (Scree Plot) trong thuật toán PCR chỉ ra rằng chỉ cần 2 cấu tử chính đầu tiên (PC1 và PC2) đã nắm giữ hơn 98% tổng phương sai của toàn bộ ma trận dữ liệu phổ, chứng minh khả năng nén thông tin và loại bỏ nhiễu xuất sắc của mô hình. So với các công bố quốc tế sử dụng HPLC-ICP-MS, giải pháp này đạt độ đúng tương đương trên các mẫu chuẩn môi trường nhưng rút ngắn thời gian phân tích từ vài giờ xuống còn dưới 15 phút cho mỗi chu kỳ đo.

Đề xuất và khuyến nghị

  1. Ban hành quy trình chuẩn quốc gia về phân tích dạng Sb(III) và Sb(V) bằng hệ HG-AAS kết hợp Chemometrics cho các trung tâm quan trắc môi trường đạt chuẩn ISO/IEC 17025. Mục tiêu là cắt giảm tối thiểu 50% chi phí hóa chất thử nghiệm và rút ngắn 70% thời gian phân tích, thực hiện trong vòng 6 tháng bởi Tổng cục Môi trường phối hợp cùng các viện nghiên cứu phân tích.
  2. Xây dựng và tích hợp module phần mềm giao diện đồ họa (GUI) tự động hóa các thuật toán tính ma trận ILS và PCR trên nền tảng máy tính phòng thí nghiệm. Mục tiêu là tự động tính toán kết quả nồng độ từ ma trận phổ trong thời gian dưới 30 giây mỗi mẫu, hoàn thiện trong thời hạn 3 tháng do các chuyên gia hóa tin học chủ trì.
  3. Triển khai chương trình điều tra diện rộng hàm lượng các dạng Antimon vô cơ trong đất trồng trọt và nước ngầm xung quanh các khu công nghiệp luyện kim và mỏ quặng. Mục tiêu là kiểm soát 100% các nguồn phát thải có nguy cơ vượt ngưỡng an toàn 5 µg/L, thiết lập bản đồ cảnh báo ô nhiễm với lộ trình 12 tháng dưới sự chỉ đạo của Sở Tài nguyên và Môi trường các địa phương.
  4. Mở rộng hướng nghiên cứu ứng dụng Chemometrics để xác định đồng thời các dạng Antimon hữu cơ độc hại như metylstibin và trimetylstibin trong các đối tượng mẫu sinh học phức tạp. Mục tiêu là đạt giới hạn phát hiện dưới 0,05 ppb cho các dẫn xuất hữu cơ trong thời gian 18 đến 24 tháng, thực hiện bởi các nhóm nghiên cứu sau đại học tại các trường đại học trọng điểm.

Đối tượng nên tham khảo luận văn

  1. Học viên cao học, nghiên cứu sinh và giảng viên thuộc các chuyên ngành Hóa phân tích, Hóa môi trường và Hóa tin học cần nguồn tài liệu học thuật bài bản về quy hoạch thực nghiệm bậc hai, phổ hấp thụ nguyên tử và kỹ thuật hồi quy ma trận đa biến.
  2. Kỹ thuật viên và chuyên viên phân tích tại các trung tâm kiểm nghiệm chất lượng môi trường, thực phẩm và y tế công cộng mong muốn ứng dụng một giải pháp phân tích dạng kim loại vi lượng hiệu quả cao, chi phí thấp trên nền thiết bị AAS hiện có.
  3. Cán bộ quản lý môi trường và hoạch định chính sách sinh thái tại các cơ quan quản lý nhà nước cần cơ sở dữ liệu khoa học về độc tính, hành vi phát tán và ngưỡng giới hạn an toàn của Antimon để hoàn thiện các quy chuẩn kỹ thuật môi trường.
  4. Các nhà nghiên cứu và lập trình viên trong lĩnh vực thiết bị phân tích đang tìm kiếm các mô hình toán học thực chứng (SVD, ILS, PCR) để tích hợp vào các phần mềm điều khiển thiết bị đo lường hóa học tự động.

Câu hỏi thường gặp

  1. Tại sao cần phân tích riêng biệt các dạng Sb(III) và Sb(V) thay vì chỉ xác định Antimon tổng số? Việc phân tích riêng từng dạng là bắt buộc vì độc tính sinh học của Sb(III) cao gấp 10 lần so với Sb(V). Nếu chỉ đo tổng lượng Antimon, cơ quan quản lý sẽ không thể đánh giá chính xác mức độ nguy hại thực tế đối với hệ sinh thái và sức khỏe con người, đặc biệt khi nồng độ trong nước tiệm cận ngưỡng giới hạn 5 µg/L.

  2. Phương pháp HG-AAS kết hợp Chemometrics có ưu thế kinh tế vượt trội nào so với HPLC-ICP-MS? Phương pháp này tận dụng được dòng máy AAS phổ biến trong hầu hết các phòng thí nghiệm tiêu chuẩn, giúp giảm hơn 60% chi phí đầu tư thiết bị và vật tư tiêu hao so với hệ thống HPLC-ICP-MS đắt tiền. Đồng thời, kỹ thuật này không sử dụng cột sắc ký tách đắt đỏ và loại bỏ hoàn toàn các dung môi hữu cơ độc hại.

  3. Những thông số kỹ thuật nào giữ vai trò quyết định đến độ nhạy của hệ đo HG-AAS? Độ nhạy của phép đo phụ thuộc chặt chẽ vào bước sóng cộng hưởng 217,53 nm, cường độ dòng đèn 8 mA, chiều cao ngọn lửa 14 mm, lưu lượng khí cháy C2H2 1,8 L/phút và bề rộng khe đo 0,2 mm. Kết hợp với nồng độ NaBH4 0,5% và axit HCl 6M sẽ tạo ra lưu lượng khí stibin ổn định với độ hấp thụ quang cao nhất.

  4. Giữa hai mô hình toán học ILS và PCR, thuật toán nào cho kết quả xử lý dữ liệu tin cậy hơn? Thuật toán PCR thể hiện tính ưu việt vượt trội nhờ kỹ thuật phân rã ma trận SVD, cho phép tập trung hơn 98% phương sai dữ liệu vào 2 cấu tử chính đầu tiên. Điều này giúp PCR triệt tiêu hoàn toàn hiện tượng đa cộng tuyến giữa các biến phổ, hạ sai số dự báo thực nghiệm xuống dưới mức 4,5%.

  5. Giới hạn làm việc tuyến tính và độ tin cậy thực nghiệm của phương pháp đạt mức bao nhiêu? Phương pháp xác lập khoảng làm việc tuyến tính rộng từ 0,25 đến 20 ppb đối với Sb(III) với hệ số tương quan R² đạt 99,98%. Khi kiểm chứng trên mẫu chuẩn quốc tế CRM và mẫu phân tích thực tế, hiệu suất thu hồi mẫu đạt từ 95,2% đến 103,8%, khẳng định độ tin cậy khoa học tuyệt đối.

Kết luận

  • Xây dựng thành công quy trình phân tích đồng thời hai dạng Antimon vô cơ Sb(III) và Sb(V) bằng phương pháp HG-AAS kết hợp Chemometrics với sai số dự báo thực nghiệm dưới 4,5%.
  • Tối ưu hóa toàn diện các thông số thiết bị quang phổ và động học phản ứng hidrua hóa thông qua mô hình quy hoạch thực nghiệm bậc hai đầy đủ đạt độ tương quan R² = 99,98%.
  • Ứng dụng xuất sắc các thuật toán hồi quy đa biến ILS và PCR để định lượng trực tiếp các dạng chất mà không cần qua giai đoạn tách sắc ký hay khử hóa trung gian phức tạp.
  • Xác lập dải đo tuyến tính từ 0,25 đến 20 ppb và kiểm chứng thành công trên mẫu chuẩn CRM cùng mẫu đất, nước thực tế với hiệu suất thu hồi đạt 95,2% đến 103,8%.
  • Đem lại giải pháp phân tích kinh tế cao giúp cắt giảm hơn 60% chi phí vận hành, mở ra tiềm năng ứng dụng rộng rãi tại các phòng kiểm nghiệm môi trường trong lộ trình 6 đến 12 tháng tới.

Các đơn vị kiểm nghiệm, trung tâm quan trắc và viện nghiên cứu môi trường hãy chủ động liên hệ chuyển giao và áp dụng quy trình phân tích tiên tiến này nhằm nâng cao hiệu quả giám sát độc chất và bảo vệ sức khỏe cộng đồng.