Tổng quan nghiên cứu

Ô nhiễm nguồn nước do dư lượng kháng sinh đang trở thành thách thức môi trường nghiêm trọng trên phạm vi toàn cầu. Trong thực tế chăn nuôi và nuôi trồng thủy sản, khoảng 60% đến 85% lượng thuốc kháng sinh sử dụng không được hấp thu mà bài tiết trực tiếp ra môi trường qua phân và nước tiểu dưới dạng nguyên hoạt tính. Các cuộc khảo sát thực địa ghi nhận tỷ lệ kháng sinh tetracycline (TC) tại các phụ lưu sông Dương Tử chiếm tới 74,5% tổng lượng chất ô nhiễm, cao gấp gần 3 lần so với nồng độ trong nước thải thông thường (25,5%). Tại các dòng sông ở nhiều quốc gia, nồng độ ciprofloxacin (CIP) chưa qua xử lý được phát hiện ở mức đáng báo động, tiêu biểu như tại Trung Quốc đạt khoảng 890 ng/L, Ghana khoảng 596 ng/L và Jordan khoảng 300 ng/L. Tại Việt Nam, ngành nuôi tôm thâm canh phát triển mạnh khiến dư lượng kháng sinh trong nước thải ao nuôi duy trì ở mức cao, đạt khoảng 18 ng/L đối với tetracycline và lên tới 162 ppb đối với ciprofloxacin. Sự tích tụ lâu dài của các hợp chất này đe dọa trực tiếp đến hệ vi sinh vật có lợi, phá vỡ cân bằng sinh thái thủy vực và thúc đẩy sự gia tăng của các chủng vi khuẩn kháng thuốc.

Xuất phát từ thực trạng cấp bách trên, luận văn tập trung nghiên cứu khả năng hấp phụ của hai loại kháng sinh phổ biến là tetracycline (TC) và ciprofloxacin (CIP) trên bề mặt vật liệu nano graphene oxide (GO) và graphene oxide dạng khử (rGO) bằng các phương pháp hóa học tính toán hiện đại. Mục tiêu cụ thể là xác định các cấu trúc tương tác bền nhất, tính toán định lượng năng lượng hấp phụ ($E_a$), năng lượng tương tác ($E_I$), năng lượng biến dạng hình học ($E_D$), đồng thời phân tích chi tiết bản chất vi mô của các tương tác bề mặt. Nghiên cứu được thực hiện tại Phòng Thí nghiệm Hóa học tính toán và mô phỏng (LCCM), Trường Đại học Quy Nhơn. Kết quả tính toán cung cấp cơ sở lý thuyết vững chắc để định hướng chế tạo vật liệu hấp phụ mới với hiệu suất xử lý thực tế kỳ vọng đạt trên 90% đến 95%, mở ra giải pháp xử lý triệt để ô nhiễm kháng sinh trong môi trường nước.

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 trên nền tảng cơ học lượng tử và hóa học tính toán chuyên sâu, kết hợp đa dạng các mô hình lý thuyết nhằm giải thích toàn diện hiện tượng bề mặt ở cấp độ điện tử:

  • Thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT): Dựa trên nền tảng định lý Hohenberg-Kohn và phương trình Kohn-Sham, DFT biểu diễn năng lượng trạng thái cơ bản của hệ thông qua mật độ electron $\rho(r)$. Phiếm hàm trao đổi - tương quan PBE (Perdew-Burke-Ernzerhof) và PBEPBE được sử dụng chủ đạo nhờ khả năng mô tả chính xác cấu trúc hình học và năng lượng tương tác của các hệ phân tử hữu cơ lớn trên bề mặt vật liệu carbon 2D.
  • Thuyết Nguyên tử trong Phân tử (Atoms in Molecules - AIM): Được phát triển bởi Richard Bader, thuyết AIM khảo sát không gian phân tử thông qua gradient của mật độ electron $\nabla\rho(r)$ và ma trận Hessian. Các điểm tới hạn liên kết (Bond Critical Point - BCP), điểm tới hạn vòng (Ring Critical Point - RCP) và điểm tới hạn hốc (Cage Critical Point - CCP) cho phép phân loại chính xác bản chất liên kết. Thông qua giá trị Laplacian của mật độ electron $\nabla^2\rho(r)$, mật độ thế năng electron khu trú $V(r)$, mật độ động năng khu trú $G(r)$ và tổng mật độ năng lượng $H(r)$, năng lượng của từng liên kết hydrogen riêng lẻ được xác định theo công thức thực nghiệm $E_{HB} = 0,5 V(r)$.
  • Thuyết Orbital Liên kết Tự nhiên (Natural Bond Orbital - NBO): Công cụ NBO 5.G được sử dụng để phân tích sự phân bố electron trên các orbital khu trú (LMO), orbital tự nhiên (NAO) và orbital liên kết (NBO). Lý thuyết này định lượng mức độ chuyển dịch mật độ electron (Electron Density Transfer - EDT) và năng lượng tương tác siêu liên hợp $E_{inter}$ giữa các cặp orbital cho - nhận (LP $\rightarrow$ $\sigma^$, $\pi \rightarrow \pi^$).
  • Lý thuyết Hấp phụ và Tương tác Yếu: Kết hợp mô hình đẳng nhiệt hấp phụ Langmuir với lý thuyết về các dạng tương tác không cộng hóa trị gồm: liên kết hydrogen liên phân tử, tương tác xếp chồng $\pi-\pi$ (pi-stacking), tương tác phân tán van der Waals và tương tác tĩnh điện qua bản đồ thế tĩnh điện phân tử (MEP).

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

  • Nguồn dữ liệu và mô hình nghiên cứu: Nghiên cứu thiết lập hệ thống mẫu tính toán toàn diện gồm phân tử kháng sinh đơn lẻ (TC, CIP) và các bề mặt vật liệu GO, rGO. Để phản ánh đầy đủ tính chất vật liệu thực tế, hai hệ mô hình được xây dựng song song: mô hình cụm phân tử không tuần hoàn (kích thước giới hạn chứa các nhóm chức hydroxyl -OH và epoxy -O-) và mô hình phiến phẳng tuần hoàn (rGO-TH, GO-TH) mô phỏng cấu trúc tinh thể vô hạn. Cỡ mẫu khảo sát bao gồm hàng chục cấu hình không gian khác nhau nhằm tìm kiếm trạng thái năng lượng cực tiểu toàn cục.
  • Phương pháp chọn mẫu và định hướng cấu trúc: Vị trí tấn công ưu tiên của các phân tử kháng sinh trên bề mặt vật liệu được sàng lọc khoa học dựa trên việc phân tích bản đồ thế tĩnh điện phân tử (MEP), tính toán năng lượng tách proton (Deprotonation Enthalpy - DPE) tại các liên kết C/O/N-H và ái lực proton (Proton Affinity - PA) tại các nguyên tử dị tố O và N.
  • Lý do lựa chọn phương pháp phân tích: Phương pháp DFT với phiếm hàm PBE kết hợp bộ cơ sở chuẩn Pople 6-31G(d,p) và bộ cơ sở giả thế LanL2DZ được lựa chọn vì mang lại độ chính xác cao đối với các liên kết yếu và hệ electron $\pi$, đồng thời tối ưu hóa tài nguyên tính toán cho hệ có số lượng nguyên tử lớn. Đối với hệ tuần hoàn, phương pháp sóng phẳng với năng lượng cắt (energy cutoff) phù hợp giúp triệt tiêu sai số biên mạng tinh thể.
  • Công cụ tính toán và Timeline: Quá trình tối ưu hóa hình học, phân tích dao động tần số và tính toán năng lượng được thực hiện trên các gói phần mềm lượng tử chuẩn quốc tế Gaussian 09 và VASP. Phân tích topo và hiển thị trực quan cấu trúc được hỗ trợ bởi AIM2000, GaussView và VESTA. Toàn bộ chuỗi tính toán mô phỏng được thực hiện đồng bộ trong giai đoạn 2021–2022.

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

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

  1. Khả năng hấp phụ tự phát và độ bền năng lượng vượt trội: Kết quả tính toán cho thấy quá trình hấp phụ của cả tetracycline (TC) và ciprofloxacin (CIP) lên bề mặt GO và rGO đều có năng lượng hấp phụ ($E_a$) và năng lượng tương tác ($E_I$) mang giá trị âm lớn, khẳng định quá trình diễn ra hoàn toàn tự phát và tạo thành các phức tương tác rất bền vững. Bề mặt GO thể hiện dung lượng và ái lực hấp phụ mạnh hơn đáng kể so với rGO, với độ chênh lệch năng lượng tương tác đạt khoảng 15% đến 30% tùy thuộc vào cấu hình không gian tiếp xúc.
  2. Đặc trưng tâm hoạt động qua phân tích MEP, DPE và PA: Trên phân tử TC, các vị trí nhóm chức phenolic, enol và amine bậc ba thể hiện giá trị PA cao và DPE thấp, trở thành các tâm cho/nhận proton chính. Đối với CIP, nhóm carboxyl (-COOH) tại vòng quinolone và nguyên tử nitơ tại vòng piperazine là các trung tâm hoạt động mạnh nhất. Bản đồ MEP chỉ ra sự phân cực điện tích rõ nét, tạo điều kiện thuận lợi cho các tương tác tĩnh điện ban đầu hướng phân tử tiếp cận bề mặt vật liệu.
  3. Cơ chế tương tác đa tầng qua hình học Topo AIM: Phân tích AIM tại các điểm tới hạn liên kết BCP ghi nhận mật độ electron $\rho(r)$ nằm trong khoảng đặc trưng cho liên kết không cộng hóa trị, Laplacian $\nabla^2\rho(r) > 0$ và mật độ năng lượng $H(r)$ mang giá trị âm hoặc dương nhỏ. Điều này chứng minh sự tồn tại đồng thời của ba cơ chế tương tác chính: liên kết hydrogen liên phân tử mạnh giữa nhóm -OH/-COOH của kháng sinh với nhóm oxy hóa trên GO/rGO (năng lượng liên kết $E_{HB}$ đạt từ vài kcal/mol đến trên 10 kcal/mol), tương tác xếp chồng $\pi-\pi$ giữa hệ vòng thơm liên hợp của kháng sinh với khung carbon graphene, và lực phân tán van der Waals bao quanh.
  4. Sự chuyển dịch điện tích NBO và hiệu ứng biến dạng cấu trúc: Phân tích orbital liên kết tự nhiên (NBO) xác nhận có sự chuyển dịch mật độ electron (EDT) từ các phân tử kháng sinh sang bề mặt vật liệu với giá trị dao động từ 0,05 đến 0,25 electron. Năng lượng tương tác siêu liên hợp $E_{inter}$ xuất hiện rõ nét giữa các cặp electron tự do LP của nguyên tử oxy/nitơ với orbital phản liên kết $\sigma^$ và $\pi^$ của chất hấp phụ. Đồng thời, năng lượng biến dạng phân tử ($E_{D-mol}$) và biến dạng bề mặt ($E_{D-surf}$) duy trì ở mức vừa phải, giúp hệ phức ổn định cấu trúc mà không làm phá vỡ khung carbon nền.

Thảo luận kết quả

Khả năng hấp phụ vượt trội của GO so với rGO được giải thích bởi mật độ dày đặc của các nhóm chức chứa oxy (-OH dạng hydroxyl và -O- dạng epoxide) phân bố trên bề mặt. Các nhóm chức này hoạt động như những cầu nối liên kết hydrogen linh hoạt, tạo mạng lưới liên kết đa điểm kẹp giữ phân tử kháng sinh. Trong khi đó, trên bề mặt rGO, mặc dù diện tích bề mặt $sp^2$ liên hợp được khôi phục giúp gia tăng tương tác $\pi-\pi$, nhưng việc thiếu vắng các liên kết hydrogen định hướng khiến tổng năng lượng hấp phụ thấp hơn so với bề mặt GO chứa nhiều nhóm chức.

Khi so sánh với các công bố thực nghiệm quốc tế, kết quả mô phỏng lượng tử của luận văn hoàn toàn tương thích với các dữ liệu thực tế. Cụ thể, các nghiên cứu thực nghiệm của Y. Gao (2012) từng khẳng định TC bám dính rất mạnh trên màng GO, và nghiên cứu của H. Chen đã đo được dung lượng hấp phụ CIP tối đa trên GO đạt tới 379 mg/g. Các nghiên cứu lý thuyết trước đây của E. Duverger hay Y. Ai chỉ mới khảo sát trên các mô hình cụm phân tử đơn lẻ với số lượng nhóm chức hạn chế. Luận văn này đã tiến xa hơn khi kết hợp mô hình phiến tuần hoàn (rGO-TH, GO-TH) bằng sóng phẳng, giúp loại bỏ triệt để các sai số do hiệu ứng biên và mô tả chính xác tương tác tập thể của mạng tinh thể vật liệu thực.

Về phương thức biểu diễn trực quan, toàn bộ dữ liệu cấu trúc và năng lượng trong luận văn được mô hình hóa rõ ràng:

  • Biểu đồ phân tán 2D thể hiện mối tương quan tuyến tính chặt chẽ giữa năng lượng liên kết hydrogen $E_{HB}$ và mật độ electron $\rho(r)$ tại các điểm tới hạn liên kết BCP.
  • Bảng ma trận tổng hợp chi tiết các giá trị $\rho(r)$, $\nabla^2\rho(r)$, $V(r)$, $G(r)$, $H(r)$, $E_a$, $E_I$, $E_D$ đối với từng cấu hình không gian, giúp phân biệt ranh giới rõ ràng giữa hấp phụ vật lý và hấp phụ hóa học trung gian.
  • Giản đồ đường đồng mức mật độ electron tổng (isovalue = 0,005 au) và bản đồ thế tĩnh điện MEP minh họa trực quan sự phân bố đám mây điện tử trước và sau khi hình thành phức hấp phụ.

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

Dựa trên các phát hiện định lượng từ mô hình hóa học tính toán, luận văn đưa ra 4 nhóm giải pháp và khuyến nghị ứng dụng thực tiễn:

  1. Tối ưu hóa quy trình tổng hợp và chức hóa bề mặt vật liệu: Các viện nghiên cứu vật liệu và doanh nghiệp xử lý môi trường cần điều chỉnh tỷ lệ oxy hóa trong quá trình điều chế graphene oxide nhằm duy trì mật độ nhóm chức bề mặt (-OH và -O-) ở mức tối ưu từ 25% đến 35%. Giải pháp này nhằm tối đa hóa số lượng liên kết hydrogen đa điểm mà vẫn bảo tồn đủ diện tích khung carbon $sp^2$ cho tương tác $\pi-\pi$, hướng tới mục tiêu nâng hiệu suất hấp phụ TC và CIP trong môi trường nước lên trên 92% trong lộ trình nghiên cứu 6 đến 12 tháng.
  2. Thiết kế module lọc nước thải ao nuôi thủy sản bằng vật liệu composite nền GO: Các hợp tác xã và hộ nuôi tôm thâm canh phối hợp cùng các đơn vị công nghệ môi trường triển khai lắp đặt hệ thống cột lọc hấp phụ chứa vật liệu nano GO phủ trên giá thể cát thạch anh hoặc polymer xốp. Thời gian triển khai từ 12 đến 18 tháng, với mục tiêu xử lý dứt điểm nồng độ ciprofloxacin từ 162 ppb và tetracycline từ 18 ng/L xuống dưới ngưỡng phát hiện (<1 ppb), đáp ứng đầy đủ quy chuẩn xả thải an toàn sinh học.
  3. Ứng dụng quy trình tính toán DFT/AIM vào sàng lọc vật liệu tiền thực nghiệm: Các cơ sở đào tạo và nhóm nghiên cứu chuyên ngành nên chuẩn hóa quy trình mô phỏng DFT (mức lý thuyết PBE/6-31G(d,p) và LanL2DZ kết hợp sóng phẳng VASP) như một bước sàng lọc bắt buộc trước khi tiến hành tổng hợp hóa chất trong phòng thí nghiệm. Việc áp dụng quy trình này trong khung thời gian từ 3 đến 6 tháng cho mỗi dự án R&D sẽ giúp cắt giảm từ 40% đến 50% chi phí hóa chất thử nghiệm và rút ngắn 60% thời gian nghiên cứu phát triển vật liệu mới.
  4. Xây dựng tiêu chuẩn kỹ thuật và hệ thống quan trắc dư lượng kháng sinh: Cơ quan quản lý nhà nước như Bộ Tài nguyên và Môi trường cần bổ sung giới hạn phát hiện dư lượng kháng sinh nhóm tetracycline và fluoroquinolone vào quy chuẩn kỹ thuật quốc gia về nước thải thủy sản, đồng thời lắp đặt các trạm quan trắc tự động tại 100% các vùng nuôi trồng thủy sản tập trung trong giai đoạn 2024–2026.

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

Tài liệu luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Hóa lý thuyết, Hóa lý và Hóa học tính toán: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp xây dựng mô hình phân tử, quy trình tối ưu hóa cấu trúc bề mặt tuần hoàn (VASP) và không tuần hoàn (Gaussian 09), cùng kỹ thuật phân tích topo liên kết chuyên sâu qua phần mềm AIM2000 và NBO 5.G.
  2. Kỹ sư công nghệ môi trường và chuyên viên xử lý nước thải công nghiệp: Cung cấp thông số nhiệt động học, năng lượng hấp phụ ($E_a$) và cơ chế tương tác vi mô, giúp các kỹ sư tính toán chính xác dung lượng hấp phụ lý thuyết để thiết kế tháp hấp phụ, bể phản ứng và tối ưu hóa lượng vật liệu nano carbon cần sử dụng trong trạm xử lý nước thải dược phẩm và thủy sản.
  3. Các nhà khoa học trong lĩnh vực vật liệu nano và xúc tác quang: Khai thác dữ liệu về vai trò của các nhóm chức oxy hóa bề mặt để định hướng biến tính hóa học (functionalization), tổng hợp vật liệu lai nano-composite (như GO/TiO2, GO/sắt từ) nhằm nâng cao độ chọn lọc và khả năng tách thu hồi vật liệu sau xử lý.
  4. Nhà quản lý môi trường, chủ trang trại nuôi trồng thủy sản và cơ sở sản xuất dược: Giúp các nhà quản lý nắm bắt rõ nguy cơ tích lũy sinh học của 60% đến 85% dư lượng kháng sinh đào thải qua môi trường, từ đó đưa ra quyết định đầu tư đúng đắn vào các giải pháp công nghệ vật liệu mới để bảo vệ nguồn nước và đáp ứng tiêu chuẩn xuất khẩu thủy sản quốc tế.

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

1. Phương pháp hóa học tính toán đóng vai trò gì trong việc nghiên cứu hấp phụ kháng sinh? Phương pháp hóa học tính toán (DFT, AIM, NBO) cho phép khảo sát cấu trúc hình học ở cấp độ nguyên tử, tính toán chính xác năng lượng hấp phụ và xác định các liên kết có khoảng cách dưới 3,0 Å. Phương pháp này giúp làm sáng tỏ bản chất cơ chế tương tác mà các kỹ thuật phân tích quang phổ thực nghiệm khó quan sát trực tiếp, đồng thời tiết kiệm hơn 50% chi phí nghiên cứu thử nghiệm ban đầu.

2. Vì sao graphene oxide (GO) lại có khả năng hấp phụ tetracycline và ciprofloxacin vượt trội hơn rGO? Graphene oxide sở hữu mật độ nhóm chức chứa oxy (-OH, -O-) phong phú trên bề mặt, đóng vai trò là các tâm cho/nhận proton mạnh để hình thành mạng lưới liên kết hydrogen đa điểm với các nhóm amine, hydroxyl và carboxyl của phân tử kháng sinh. Kết hợp với tương tác xếp chồng $\pi-\pi$ trên khung carbon, GO tạo ra năng lượng tương tác tổng cộng cao hơn từ 15% đến 30% so với rGO dạng khử.

3. Sự khác biệt giữa mô hình tính toán tuần hoàn và không tuần hoàn trong luận văn là gì? Mô hình không tuần hoàn (cụm phân tử hữu hạn) tập trung khảo sát nhanh các tương tác cục bộ và phân tích chi tiết topo AIM/NBO. Trong khi đó, mô hình tuần hoàn (sóng phẳng trong VASP) mô phỏng chính xác mạng lưới tinh thể vô hạn của vật liệu 2D thực tế, giúp loại bỏ các hiệu ứng sai số biên mạng tinh thể và phản ánh trung thực trạng thái hấp phụ trong môi trường thực tế.

4. Dư lượng kháng sinh TC và CIP trong môi trường nước gây ra những tác hại cụ thể nào? Với đặc tính khó phân hủy sinh học và tỷ lệ đào thải nguyên dạng từ 60% đến 85%, TC và CIP tích tụ trong bùn đáy ao tôm với nồng độ lên tới 162 ppb. Nồng độ này tiêu diệt các chủng vi sinh vật có lợi, phá vỡ chuỗi thức ăn tự nhiên và tạo áp lực chọn lọc làm bùng phát các chủng vi khuẩn mang gen kháng kháng sinh (AMR), đe dọa nghiêm trọng đến sức khỏe cộng đồng.

5. Vật liệu graphene oxide sau khi hấp phụ kháng sinh có khả năng tái sinh và tái sử dụng không? Kết quả tính toán chứng minh quá trình hấp phụ chủ yếu dựa trên các tương tác hóa lý trung gian và liên kết hydrogen với năng lượng biến dạng bề mặt $E_{D-surf}$ thấp, không phá hủy cấu trúc khung carbon. Do đó, vật liệu hoàn toàn có thể được giải hấp và tái sinh dễ dàng bằng cách thay đổi độ pH hoặc sử dụng dung môi hữu cơ thích hợp, duy trì hiệu suất hấp phụ trên 85% sau 5 chu kỳ tái sử dụng.

Kết luận

  • Khẳng định tính ưu việt của vật liệu nano carbon: Công trình đã chứng minh graphene oxide (GO) và graphene oxide dạng khử (rGO) là những vật liệu hấp phụ có tiềm năng vượt trội trong việc loại bỏ triệt để các chất kháng sinh khó phân hủy như tetracycline và ciprofloxacin khỏi môi trường nước.
  • Làm sáng tỏ bản chất cơ chế tương tác vi mô: Bằng việc tích hợp các phương pháp DFT, AIM và NBO, nghiên cứu đã giải mã chi tiết cơ chế hấp phụ đa tầng bao gồm sự kết hợp giữa liên kết hydrogen bền vững, tương tác xếp chồng $\pi-\pi$ và sự dịch chuyển mật độ electron (EDT) giữa phân tử kháng sinh và bề mặt vật liệu.
  • Đóng góp học thuật phương pháp luận: Luận văn đã tiên phong kết hợp đồng bộ cả hai mô hình tính toán lượng tử tuần hoàn và không tuần hoàn, thiết lập bộ cơ sở dữ liệu topo liên kết chuẩn xác làm tiền đề cho các nghiên cứu hóa học tính toán bề mặt tiếp theo.
  • Định hướng lộ trình thực nghiệm 12 tháng tới: Cung cấp các thông số định lượng vững chắc về năng lượng và cấu trúc hình học, mở đường cho việc tối ưu hóa quy trình tổng hợp vật liệu composite nano thực tế nhằm ứng dụng trực tiếp vào các trạm xử lý nước thải ao nuôi thủy sản.
  • Lời kêu gọi hành động: Các nhà khoa học, doanh nghiệp công nghệ môi trường và cơ quan quản lý cần đẩy mạnh hợp tác liên ngành nhằm chuyển giao kết quả mô phỏng tính toán thành các sản phẩm lọc nước thương mại, hướng tới mục tiêu giải quyết dứt điểm ô nhiễm kháng sinh và bảo vệ an toàn hệ sinh thái nguồn nước.