Nghiên cứu xử lý DDT & HCH trên kim loại/oxide mang trên g-C3N4 bằng DFT-D

Tổng quan về luận án

Ô nhiễm môi trường bởi các hợp chất hữu cơ khó phân hủy (Persistent Organic Pollutants – POPs) là một trong những thách thức sinh thái và y tế cộng đồng nghiêm trọng nhất trên phạm vi toàn cầu. Theo ghi nhận trong tài liệu nghiên cứu, "Công ước quốc tế Stockholm với mục đích loại bỏ hoặc hạn chế sản xuất và sử dụng các chất ô nhiễm hữu cơ khó phân hủy tính đến tháng 5/2013 đã có sự tham gia của 178 quốc gia và liên minh châu âu". Tại Việt Nam – một quốc gia nông nghiệp có lịch sử sử dụng hóa chất bảo vệ thực vật (BVTV) lâu đời và hiện tồn tại nhiều kho lưu trữ thuốc BVTV quá hạn xuống cấp – dư lượng của dichlorodiphenyltrichloroethane (DDT, công thức $C_{14}H_9Cl_5$) và hexachlorocyclohexane ($\gamma$-HCH hay lindane, công thức $C_6H_6Cl_6$) trong đất và nguồn nước ngầm vẫn duy trì ở mức báo động. Luận án nhấn mạnh: "DDT là một thuốc BVTV rất bền vững do nó bền, không phản ứng với oxygen trong không khí. DDT có thời gian bán phân hủy từ 10-15 năm", đồng thời được Cơ quan Nghiên cứu Ung thư Quốc tế (IARC) xếp vào nhóm 2B và phân loại độc học thủy sinh nhóm 1 ($TL_m \le 1\text{ mg/l}$).

graph TD
    A["POPs Ô nhiễm: DDT ($C_{14}H_9Cl_5$) & $\gamma$-HCH ($C_6H_6Cl_6$)"] --> B["Thách thức: Tồn lưu 10-15 năm, Bền nhiệt, Kháng phân hủy sinh học"]
    B --> C["Hạn chế của g-C3N4 nguyên khối ($E_g \approx 2.7\text{ eV}$)"]
    C --> D["Tái tổ hợp e-/h+ nhanh, Diện tích bề mặt bé, Hấp phụ vật lý yếu"]
    D --> E["Giải pháp Biến tính: Doping Kim loại (Fe, Ni) & Oxide ((TiO2)7)"]
    E --> F["Khảo sát Lượng tử: GFN2-xTB, CREST, RP (14 images), TD-DFT, MD"]
    F --> G["Cơ chế Đột phá: Hấp phụ hóa học + Chuyển dịch điện tích cắt liên kết C-Cl"]

Trước thực trạng các công nghệ xử lý truyền thống bộc lộ nhiều hạn chế về hiệu suất ở nồng độ vết và chi phí vận hành, công nghệ hấp phụ kết hợp quang xúc tác trên nền graphitic carbon nitride ($g\text{-C}_3\text{N}_4$) đã nổi lên như một giải pháp xanh giàu tiềm năng. Luận án xác định: "g-C3N4 – một chất bán dẫn polymer nổi lên như một ứng viên tiềm năng trong công nghệ quang xúc tác do sở hữu nhiều đặc tính hứa hẹn và hấp dẫn như: ổn định hóa học và bền nhiệt, năng lượng vùng cấm vừa phải (khoảng 2,7 eV) thể hiện khả năng quang xúc tác trong vùng ánh sáng nhìn thấy". Tuy nhiên, $g\text{-C}_3\text{N}_4$ nguyên khối tồn tại nhược điểm chí mạng: tốc độ tái tổ hợp cặp electron - lỗ trống ($e^-/h^+$) quang sinh cao, diện tích bề mặt riêng nhỏ và hiệu suất lượng tử thấp.

Khoảng trống nghiên cứu then chốt (Research Gap) được luận án chỉ ra là: Hầu hết các công trình trước đây chỉ dừng lại ở khảo sát thực nghiệm vĩ mô, thiếu hoàn toàn các nghiên cứu lý thuyết chuyên sâu ở cấp độ phân tử và lượng tử nhằm giải thích bản chất tương tác hấp phụ, cơ chế truyền điện tích và con đường phân hủy trực tiếp của DDT và $\gamma$-HCH trên bề mặt $g\text{-C}_3\text{N}_4$ được biến tính. Để giải quyết khoảng trống này, luận án thiết lập 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học tương ứng:

• Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc hình học, phân bố electron và tính chất quang của $g\text{-C}_3\text{N}_4$ ở dạng gợn sóng (cGN) và dạng phẳng (pGN) biến đổi như thế nào khi được pha tạp bởi 6 nguyên tử kim loại ($K, Ca, Ga, Fe, Ni, Cu$) và 2 cụm oxide kim loại ($(ZnO)_3, (TiO_2)_7$)?
• Câu hỏi nghiên cứu 2 (RQ2): Bản chất tương tác hấp phụ của DDT và $\gamma$-HCH trên $g\text{-C}_3\text{N}_4$ nguyên bản và biến tính là hấp phụ vật lý hay hấp phụ hóa học?
• Câu hỏi nghiên cứu 3 (RQ3): Năng lượng vùng cấm ($E_g$), vị trí dải năng lượng ($E_{VB}, E_{CB}$) và phổ hấp thụ quang thay đổi ra sao để thúc đẩy phân tách điện tích quang sinh?
• Câu hỏi nghiên cứu 4 (RQ4): Có tồn tại một cơ chế phân hủy quang xúc tác trực tiếp thông qua chuyển điện tích làm suy yếu liên kết $C\text{-}Cl$ thay vì cơ chế gián tiếp qua gốc tự do thông thường?

Tương ứng với đó là các giả thuyết khoa học:

• Giả thuyết 1 (H1): Sự kết hợp của kim loại chuyển tiếp ($Fe, Ni$) và oxide bán dẫn ($(TiO_2)_7$) lên $g\text{-C}_3\text{N}_4$ sẽ tạo ra các trạng thái xen kẽ trong vùng cấm, giảm $E_g$ và phân tách không gian giữa HOMO và LUMO.
• Giả thuyết 2 (H2): DDT và $\gamma$-HCH chuyển dịch từ tương tác phân tán van der Waals yếu trên $g\text{-C}_3\text{N}_4$ nguyên bản sang liên kết phối trí/hóa học mạnh trên các tâm biến tính $Fe, Ni, TiO_2$.
• Giả thuyết 3 (H3): Dung môi có độ phân cực khác nhau (nước, ethanol, acetonitrile, benzene) và nhiệt độ dao động nhiệt sẽ điều biến năng lượng hấp phụ nhưng không phá vỡ cấu hình bền vững của phức hệ xúc tác.
• Giả thuyết 4 (H4): Sự chuyển dịch electron trực tiếp từ vùng dẫn/LUMO của xúc tác sang orbital phản liên kết $\sigma^*(C\text{-}Cl)$ của phân tử POPs là nguyên nhân cốt lõi gây kéo dài và cắt đứt liên kết $C\text{-}Cl$.

Khung lý thuyết của công trình được định hình vững chắc trên nền tảng cơ học lượng tử, bao gồm Lý thuyết phiếm hàm mật độ liên kết chặt mở rộng thế hệ 2 (GFN2-xTB của Grimme), Lý thuyết phiếm hàm mật độ phụ thuộc thời gian giản lược (sTD-DFT), phương pháp lấy mẫu cấu hình phân tử tự động (CREST), phương pháp tìm đường phản ứng (Reaction Path - RP) và Động lực học phân tử lượng tử (QMD). Phạm vi nghiên cứu bao gồm mô hình hóa $g\text{-C}_3\text{N}_4$ đơn lớp triazine, đánh giá 8 hệ vật liệu biến tính, 4 hệ dung môi liên tục và mô phỏng động học nhiệt động ở 298 K cùng các mức nhiệt độ nâng cao trên quy mô dữ liệu gồm 41 bảng thông số và 83 đồ thị phổ/hình học phân tử chi tiết.

Literature Review và Positioning

Lịch sử nghiên cứu xúc tác quang bán dẫn ghi nhận bước ngoặt lớn từ công trình tiên phong của Wang và cộng sự (2009) khi phát hiện $g\text{-C}_3\text{N}_4$ có khả năng phân tách nước quang hóa tạo $H_2$ và $O_2$ dưới ánh sáng khả kiến. Kể từ đó, nhiều hướng nghiên cứu quốc tế đã tập trung tối ưu hóa vật liệu này. Xu và cộng sự (2013) chứng minh việc gia tăng diện tích bề mặt riêng giúp cải thiện đáng kể hoạt tính phân hủy phẩm màu methylene blue. Song song đó, Li và cộng sự (2011) chỉ ra sự hình thành dị thể $g\text{-C}_3\text{N}_4/ZnO$ nâng cao hoạt tính quang xúc tác nhờ tăng cường khả năng lưu giữ chất ô nhiễm trên bề mặt. Hướng biến tính pha tạp cation kim loại chuyển tiếp ($Fe, Co, Ni$) được Durairaj và cộng sự (2017) xác nhận có khả năng ức chế tái tổ hợp $e^-/h^+$ và thúc đẩy phân hủy các chất hữu cơ bền vững như reactive red 120 và 2-amino-4-chlorophenol.

graph LR
    subgraph "Trường phái Cổ điển (Empirical AOPs)"
        CP1["Zeleska et al. (2010)<br/>TiO2 bột + UV: Hấp phụ yếu, tách 50% HCH / 85% DDT"]
        CP2["Zhao et al. (2012)<br/>TiO2/Montmorillonite: Xúc tác gián tiếp qua gốc tự do OH•/O2•-"]
    end
    subgraph "Trường phái Lượng tử Đột phá (Luận án Phạm Thị Bé)"
        AP1["Tổ hợp Đa tầng: GFN2-xTB + CREST + RP + MD"]
        AP2["Làm sáng tỏ Bản chất: Physisorption -> Chemisorption trên Fe, Ni, TiO2"]
        AP3["Cơ chế Mới: Direct Electron Transfer kích hoạt giãn dài & đứt liên kết C-Cl"]
    end
    CP1 -.->|Chuyển giao Paradigm| AP1
    CP2 -.->|Khắc phục khoảng trống lý thuyết| AP3

Trong bức tranh tổng quan đó, tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

• Quan điểm truyền thống (Empirical Advanced Oxidation Processes): Cho rằng sự phân hủy quang xúc tác của các chất hữu cơ khó phân hủy POPs diễn ra chủ yếu thông qua con đường gián tiếp, phụ thuộc hoàn toàn vào việc sinh ra các gốc oxy hóa hoạt tính cao như hydroxyl ($\cdot\text{OH}$) và superoxide ($\cdot\text{O}_2^-$) trong dung dịch, trong khi bản thân chất xúc tác chỉ đóng vai trò trơ về mặt phối trí trực tiếp với cơ chất.
• Quan điểm lượng tử tương tác bề mặt (Luận án đề xuất): Khẳng định rằng đối với các phân tử hữu cơ chứa liên kết halogen cực kỳ bền như DDT và HCH, tương tác phối trí hóa học trực tiếp giữa cơ chất và tâm kim loại/oxide kim loại trên bề mặt xúc tác đóng vai trò quyết định; trong đó, sự chuyển dịch electron quang sinh trực tiếp vào orbital phản liên kết của cơ chất sẽ kích hoạt sự phân cắt liên kết $C\text{-}Cl$ với hàng rào thế năng thấp hơn nhiều so với việc tấn công ngẫu nhiên của các gốc tự do.

So sánh với hai nghiên cứu quốc tế điển hình trong xử lý POPs:

1. Nghiên cứu của Zeleska và cộng sự (2010): Khảo sát phân hủy HCH, DDT và methoxychlor trong nước bằng bột $TiO_2$ dưới tia UV, đạt hiệu suất loại bỏ 50% HCH và 85% DDT sau 100 phút. Tuy nhiên, công trình này bị giới hạn ở bức xạ tử ngoại năng lượng cao ($<387\text{ nm}$) và không giải thích được cấu trúc trạng thái chuyển tiếp tại bề mặt phân chia pha.
2. Nghiên cứu của Ye và cộng sự (2018): Ứng dụng nanocomposite $g\text{-C}_3\text{N}_4/TiO_2$ để phân hủy polybromodiphenyl ethers (PBDEs) dưới ánh sáng khả kiến đạt hiệu suất 90% trong 2 giờ, xác nhận sản phẩm tạo thành có ít nguyên tử halogen hơn qua GC-MS nhưng để ngỏ cơ chế phân bố orbital phân tử tương tác.

Luận án của Phạm Thị Bé đã định vị chính xác vị trí khoa học của mình: bắc cây cầu nối giữa hóa học tính toán lượng tử chính xác cao và hóa học môi trường thực nghiệm. Bằng việc áp dụng phương pháp GFN2-xTB kết hợp mô phỏng động lực học phân tử MD, luận án không chỉ giải mã tường tận bản chất cấu trúc electron mà còn tiên phong chứng minh cơ chế phân hủy quang xúc tác trực tiếp – một đóng góp vượt trội so với các công bố thực nghiệm thuần túy trước đây.

Đó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à làm sâu sắc thêm các lý thuyết nền tảng của hóa học lượng tử và quang xúc tác:

• Phát triển ứng dụng Lý thuyết Phiếm hàm Mật độ (Hohenberg & Kohn, 1964; Kohn & Sham, 1965): Mở rộng giới hạn áp dụng của phương pháp phiếm hàm mật độ từ các hệ phân tử nhỏ sang các siêu phân tử polymer 2D phức tạp chứa hàng trăm nguyên tử thông qua phương pháp liên kết chặt bán kinh nghiệm GFN2-xTB (Grimme et al., 2019).
• Làm rõ vai trò của tương tác phân tán tầm xa (Grimme, 2017): Ứng dụng thành công mô hình hiệu chỉnh phân tán thế hệ mới DFT-D4 với hàm giảm chấn Becke-Johnson ($f_{damp}$), chứng minh rằng tương tác van der Waals chi phối hoàn toàn quá trình hấp phụ vật lý trên $g\text{-C}_3\text{N}_4$ nguyên bản, nhưng bị lấn át bởi tương tác trao đổi điện tích cục bộ khi có mặt các tâm kim loại $Fe, Ni$.
• Định lượng trạng thái kích thích bằng Lý thuyết Runge-Gross (1984): Thông qua sTD-DFT, luận án đã lượng hóa chính xác các bước chuyển dời electron đơn phân tử, năng lượng kích thích đầu tiên và sự hình thành các trạng thái trung gian quang sinh.
• Chuyển dịch Paradigm (Paradigm Shift): Chuyển dịch mô hình nhận thức từ việc xem $g\text{-C}_3\text{N}_4$ như một chất mang thụ động sang một hệ vật liệu xúc tác quang chủ động có khả năng chuyển dịch điện tích trực tiếp vào cơ chất POPs, làm biến dạng hình học và suy yếu liên kết $C\text{-}Cl$.

classDiagram
    class TheoreticalFramework {
        +Hohenberg_Kohn_Theorems (1964)
        +Kohn_Sham_Equations (1965)
        +Runge_Gross_TDDFT (1984)
        +Grimme_GFN2_xTB (2019)
        +Verlet_Berendsen_MD (1967/1984)
    }
    class AnalyticalWorkflow {
        +Tầng 1: Tối ưu hóa cấu trúc & Electron (IP, EA, GEI, FOD)
        +Tầng 2: Quét bề mặt thế năng & Đường phản ứng (CREST, RP 14-images)
        +Tầng 3: Khảo sát Ổn định Động học & Môi trường (MD, Solvation Models)
    }
    TheoreticalFramework --> AnalyticalWorkflow : Hiện thực hóa thành

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

Khung phân tích của luận án tích hợp đa tầng lý thuyết lượng tử thành một quy trình tính toán liền mạch:

$$\text{Khung Phân Tích} = \text{GFN2-xTB} \xrightarrow{\text{Cấu trúc & FOD}} \text{CREST/RP} \xrightarrow{\text{Đường phản ứng}} \text{sTD-DFT} \xrightarrow{\text{Tính chất quang}} \text{QMD} \xrightarrow{\text{Độ bền nhiệt/dung môi}}$$

1. Khái niệm và Định nghĩa Cốt lõi:
   • Chỉ số Electrophile toàn phần (Global Electrophilicity Index - GEI): $\text{GEI} = \frac{(\text{IP}+\text{EA})^2}{8(\text{IP}-\text{EA})}$, đặc trưng cho tính acid Lewis và khả năng thu nhận cặp electron của hệ xúc tác.
   • Mật độ bị chiếm từng phần (Fractional Occupation Density - FOD): Định lượng mức độ bất ổn định electron và các electron không ghép đôi bề mặt: $\text{FOD}(r) = \sum_i^{N} |1 - f_i| |\psi_i(r)|^2$.
   • Năng lượng tương tác và hấp phụ ($E_{int}, E_{ads}$): Xác định biến thiên năng lượng tự do khi hình thành liên kết phân tử.
2. Điều kiện biên (Boundary Conditions): Mô hình hóa vật liệu trên cụm đơn lớp gợn sóng (cGN) và phẳng (pGN) gồm các đơn vị triazine tuần hoàn; chuẩn hóa mức năng lượng electron tự do $E_e = 4.500\text{ eV}$; nhiệt độ cân bằng electron $T_{el} = 300\text{ K}$; tiêu chuẩn hội tụ năng lượng nghiêm ngặt $0.1 \times 10^{-5}\text{ Hartree}$ và hàm sóng $0.1 \times 10^{-3}\text{ e}$.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi triết lý thực chứng lượng tử (Quantum Positivism / Computational Realism), xây dựng trên quy trình mô phỏng đa cấp độ (Multi-level Computational Design):

• Cấp độ 1 (Điện tử cục bộ): Khảo sát cấu hình orbital biên (HOMO/LUMO), mật độ trạng thái riêng (PDOS) và bản chất liên kết hóa học.
• Cấp độ 2 (Tọa độ phản ứng): Khảo sát bề mặt thế năng tương tác giữa phối tử POPs và chất mang biến tính.
• Cấp độ 3 (Động lực học hệ thống): Đánh giá tính ổn định cấu trúc dưới tác động của nhiệt độ và tương tác tĩnh điện từ các dung môi khác nhau.

flowchart TD
    subgraph S1["1. Chuẩn bị Mô hình"]
        M1["g-C3N4: cGN & pGN"]
        M2["Kim loại: K, Ca, Ga, Fe, Ni, Cu"]
        M3["Oxide: (ZnO)3, (TiO2)7"]
        M4["Cơ chất: DDT & γ-HCH"]
    end
    subgraph S2["2. Tính toán Lượng tử (GFN2-xTB / sTD-DFT)"]
        C1["Tối ưu hóa: Hội tụ 0.1x10^-5 Hartree, 0.1x10^-3 e"]
        C2["Tính toán: IP, EA, GEI, FOD, HOMO-LUMO"]
        C3["Quang phổ: TD-DFT, Eg, EVB, ECB (Ee = 4.500 eV)"]
    end
    subgraph S3["3. Khảo sát Hấp phụ & Động học"]
        A1["CREST: Tìm kiếm cấu hình tối ưu toàn cục"]
        A2["Reaction Path (RP): 14 ảnh nội suy MEP"]
        A3["Mô phỏng MD: Thuật toán Verlet, Berendsen (298 K+)"]
        A4["Dung môi: Implicit Solvation (Nước, EtOH, ACN, Benzene)"]
    end
    S1 --> S2 --> S3

Quy trình nghiên cứu rigorous

Luận án tuân thủ nghiêm ngặt các giao thức tính toán chuẩn mực quốc tế:

• Tiêu chuẩn tính toán lượng tử: Năng lượng toàn phần $E$ trong GFN2-xTB được phân rã chi tiết: $$E = E_{rep} + E_{disp} + E_{EHT} + E_{IES+IXC} + E_{AES} + E_{AXC} + G_{Fermi}$$ đảm bảo tính toán đồng thời lực đẩy lõi, hiệu chỉnh phân tán D4, tương tác tĩnh điện và trao đổi đẳng hướng/dị hướng.
• Quy trình tối ưu hóa cấu hình toàn cục bằng CREST: Sử dụng thuật toán Conformer-Rotamer Ensemble Sampling Tool kết hợp metadynamics để quét toàn bộ không gian pha cấu hình của DDT và HCH trên bề mặt $g\text{-C}_3\text{N}_4$, loại bỏ triệt để nguy cơ rơi vào các cực tiểu cục bộ giả định.
• Nghiên cứu đường phản ứng (Reaction Path - RP): Sử dụng 14 điểm ảnh (images) rời rạc dọc theo tọa độ phản ứng để theo dõi liên tục sự biến thiên năng lượng tương đối ($E_{rel}$) và độ dài liên kết khi phân tử POPs tiến từ vô cùng vào bề mặt xúc tác.
• Kiểm định độ tin cậy (Triangulation & Reliability): Đối chiếu chéo kết quả giữa các chỉ số năng lượng orbital (Koopmans), hiệu ứng phổ tử ngoại - khả kiến (sTD-DFT) và hình ảnh mật độ electron bị chiếm từng phần (FOD) tại đẳng giá trị $0.03\text{ e \AA}^{-3}$.

Data và phân tích

• Mô phỏng Động lực học Phân tử (MD): Tích phân phương trình chuyển động Newton thông qua thuật toán Verlet với bước thời gian siêu mịn $\Delta t = 1.0\text{ fs}$: $$R(t+\Delta t) = 2R(t) - R(t-\Delta t) + \frac{f(R(t))}{m}\Delta t^2 + O(\Delta t^4)$$ Hệ thống được cân bằng nhiệt động thông qua bể nhiệt Berendsen và thuật toán velocity scaling ($\mathbf{v}{mới} = \mathbf{v}{cũ}\sqrt{T_{nghiên cứu}/T(t)}$) ở 298 K và dải nhiệt độ mở rộng để kiểm tra độ bền liên kết.
• Khảo sát hiệu ứng dung môi: Đánh giá độ bền phức hệ trong 4 môi trường dung môi liên tục có hằng số điện môi trải rộng: Nước ($\varepsilon = 78.4$), Ethanol ($\varepsilon = 24.5$), Acetonitrile ($\varepsilon = 35.9$), và Benzene ($\varepsilon = 2.27$).

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

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

graph TB
    subgraph "So sánh Tương tác Cơ bản vs Biến tính"
        N1["g-C3N4 Nguyên bản (cGN / pGN)"] -->|"Tương tác van der Waals"| R1["Hấp phụ Vật lý Yếu<br/>d(C-Cl) không đổi (~1.78 Å)<br/>Không kích hoạt phản ứng"]
        N2["Me/g-C3N4 (Fe, Ni) & TiO2/g-C3N4"] -->|"Chuyển dịch điện tích mạnh"| R2["Hấp phụ Hóa học Bền vững<br/>d(C-Cl) giãn dài > 1.95 Å<br/>Phân cắt trực tiếp liên kết C-Cl"]
    end

1. Bản chất biến tính mạng $g\text{-C}_3\text{N}_4$ bởi kim loại và oxide: Khảo sát 6 kim loại ($K, Ca, Ga, Fe, Ni, Cu$) và cụm $(ZnO)_3, (TiO_2)_7$ chỉ ra rằng các kim loại chuyển tiếp $Fe, Ni$ và cluster $(TiO_2)_7$ tạo liên kết phối trí bền vững nhất trong các hốc tam-triazine của cGN và pGN. Sự kết hợp này làm biến đổi sâu sắc mật độ trạng thái (PDOS), hạ thấp đáng kể năng lượng vùng cấm ($E_g$), đồng thời chuyển dịch phổ hấp thụ UV-Vis mạnh mẽ về phía vùng ánh sáng nhìn thấy ($\lambda > 400\text{ nm}$).
2. Sự chuyển pha cơ chế hấp phụ từ Vật lý sang Hóa học:
   • Trên $g\text{-C}_3\text{N}_4$ nguyên khối (cả cGN và pGN), năng lượng hấp phụ của DDT và HCH rất nhỏ, khoảng cách tương tác lớn ($>3.2\text{ \AA}$), độ dài liên kết $C\text{-}Cl$ không thay đổi so với phân tử tự do $\rightarrow$ Bản chất thuần túy là hấp phụ vật lý.
   • Trên các hệ $Fe/g\text{-C}_3\text{N}_4$, $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$, năng lượng hấp phụ tăng vọt, xuất hiện sự dịch chuyển điện tích rõ rệt từ phối tử sang tâm kim loại $\rightarrow$ Bản chất chuyển hóa hoàn toàn thành hấp phụ hóa học.
3. Hiện tượng phân tách không gian electron và lỗ trống ($e^-/h^+$): Phân tích hình ảnh orbital biên tại đẳng giá trị $0.03\text{ e \AA}^{-3}$ chứng minh: HOMO tập trung chủ yếu trên khung $g\text{-C}_3\text{N}_4$ trong khi LUMO định cư tập trung tại các tâm kim loại $Fe, Ni$ hoặc cluster $TiO_2$. Sự cô lập không gian này kéo dài thời gian sống của hạt tải điện quang sinh, giải quyết triệt để nhược điểm tái tổ hợp nhanh của vật liệu nguyên bản.
4. Cơ chế phân hủy trực tiếp thông qua kéo dài và phân cắt liên kết $C\text{-}Cl$: Khảo sát 14 điểm ảnh trên đường hấp phụ (Reaction Path) phát hiện sự giãn dài bất thường của liên kết $C\text{-}Cl$ trong phân tử DDT (nhóm $-CCl_3$ và vòng chlorobenzene) và vòng $\gamma$-HCH từ giá trị chuẩn $1.78\text{ \AA}$ lên đến $1.95 - 2.15\text{ \AA}$ khi tiếp xúc với tâm xúc tác $Fe/Ni$. Mật độ điện tích âm tích tụ lớn trên nguyên tử $Cl$, chứng minh electron quang sinh từ vùng dẫn ($E_{CB}$) đã truyền thẳng vào orbital $\sigma^*(C\text{-}Cl)$, dẫn đến phản ứng cắt đứt liên kết $C\text{-}Cl$ mà không cần chất oxy hóa trung gian.
5. Độ bền nhiệt động và tính thích ứng dung môi: Mô phỏng MD xác nhận các cấu hình hấp phụ hóa học của DDT và HCH trên $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$ duy trì cấu trúc liên kết ổn định ở 298 K và không bị giải hấp ở nhiệt độ cao. Hiệu ứng dung môi phân cực (nước, ethanol) tăng cường khả năng ổn định phức hợp nhờ hiệu ứng solvat hóa tĩnh điện.

Bảng tổng hợp so sánh tính chất các hệ vật liệu then chốt

Implications đa chiều

• Ý nghĩa Lý thuyết: Khẳng định tính đúng đắn của việc tích hợp lý thuyết liên kết chặt GFN2-xTB trong giải quyết các bài toán xúc tác dị thể phức tạp; thiết lập mô hình truyền electron cảm ứng cho các phản ứng khử dehalogen hóa quang hóa.
• Đột phá Phương pháp luận: Thiết lập quy trình tính toán chuẩn từ sàng lọc đơn phân tử đến mô phỏng động học chất lỏng, có thể chuyển giao hoàn toàn để nghiên cứu các nhóm độc chất POPs khác (Dioxin, Furan, PCBs).
• Ứng dụng Thực tiễn & Định hướng Tổng hợp: Cung cấp thông số cấu trúc chính xác để các nhà thực nghiệm định hướng tổng hợp vật liệu composite $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$ bằng các phương pháp nhiệt phân hoặc sol-gel, tránh việc thử nghiệm mò mẫm tốn kém hóa chất độc hại.
• Khuyến nghị Chính sách: Cung cấp cơ sở khoa học để Bộ Tài nguyên và Môi trường cùng các cơ quan quản lý ban hành quy chuẩn kỹ thuật mới về xử lý triệt để các điểm nóng tồn lưu thuốc BVTV POPs theo cam kết Công ước Stockholm.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án duy trì thái độ học thuật khách quan khi nhận diện 4 giới hạn cụ thể:

1. Mô hình hóa cụm hữu hạn (Cluster Model Approximation): Luận án sử dụng các mô hình phiến phẳng và gợn sóng kích thước hữu hạn thay vì mô hình tuần hoàn 3 chiều điều kiện biên (PBC), điều này có thể dẫn đến hiệu ứng biên nhất định ở các rìa phân tử.
2. Gần đúng bán kinh nghiệm trong GFN2-xTB: Dù đã hiệu chỉnh tham số cho các nguyên tố đến $Z=86$, phép gần đúng "f-in-core" và việc bán kinh nghiệm hóa các tích phân trao đổi electron có thể tạo ra dung sai nhỏ đối với các trạng thái kích thích Rydberg bậc cao so với các phép tính ab initio cấp độ CASSCF/MRCI.
3. Mô hình dung môi liên tục (Implicit Solvation): Chưa phản ánh được tương tác vi mô của các cầu liên kết hydro động học từ từng phân tử nước cụ thể bao quanh tâm phản ứng.
4. Bỏ qua ảnh hưởng của các ion vô cơ cạnh tranh: Chưa tính toán đến sự có mặt đồng thời của các anion phổ biến trong nước tự nhiên ($Cl^-, NO_3^-, SO_4^{2-}, HCO_3^-$) vốn có thể cạnh tranh tâm hấp phụ.

Chương trình nghiên cứu tương lai (Future Agenda) được vạch ra với 5 định hướng:

• Hướng 1: Mở rộng tính toán DFT tuần hoàn (Periodic DFT) với phiếm hàm lai chất lượng cao HSE06 để kiểm chứng cấu trúc vùng năng lượng tinh thể khối.
• Hướng 2: Thực hiện mô phỏng Động lực học phân tử lượng tử với mô hình dung môi tường minh (Explicit Solvation QM/MM) để quan sát sự chuyển dịch proton đồng thời trong pha lỏng.
• Hướng 3: Triển khai thực nghiệm tổng hợp vật liệu $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$ trong phòng thí nghiệm, kiểm chứng phổ FTIR, XRD, HR-TEM và đo hoạt tính quang xúc tác thực tế với DDT/HCH.
• Hướng 4: Mở rộng nghiên cứu cơ chế sang các nhóm hợp chất ô nhiễm mới nổi (PFAS, thuốc bảo vệ thực vật gốc lân hữu cơ và dược phẩm tồn dư).
• Hướng 5: Thiết kế mô hình phản ứng dòng chảy liên tục (Continuous-flow photoreactor) ứng dụng tấm phủ quang xúc tác $g\text{-C}_3\text{N}_4$ biến tính năng lượng mặt trời.

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

graph TD
    Impact["Tác động của Luận án Tiến sĩ"]
    Impact --> AC["Học thuật (Academic Impact)"]
    Impact --> IND["Công nghiệp & Công nghệ (Industry/R&D)"]
    Impact --> POL["Môi trường & Xã hội (Societal/Policy)"]
    
    AC --> AC1["Mô hình GFN2-xTB + CREST tiên phong"]
    AC --> AC2["07 bài báo công bố trong nước & quốc tế"]
    
    IND --> IND1["Công nghệ Hấp phụ - Quang xúc tác phối hợp"]
    IND --> IND2["Thiết kế Module xử lý nước thải chứa POPs"]
    
    POL --> POL1["Hiện thực hóa cam kết Công ước Stockholm"]
    POL --> POL2["Xử lý điểm nóng ô nhiễm kho thuốc BVTV quá hạn"]

• Tác động Học thuật: Công trình là tài liệu tham khảo chuẩn mực trong lĩnh vực Hóa lý thuyết và Hóa lý tại Việt Nam. Phương pháp tiếp cận kết hợp GFN2-xTB, CREST và RP mở ra tiêu chuẩn mới cho các nghiên cứu mô phỏng xúc tác môi trường, dự kiến thu hút trích dẫn cao trong các tạp chí chuyên ngành xúc tác quang và hóa học tính toán.
• Chuyển dịch Công nghệ & R&D: Đặt nền móng kỹ thuật cho việc phát triển các modul xử lý nước thải di động sử dụng màng quang xúc tác kích hoạt bằng ánh sáng mặt trời, giúp các doanh nghiệp công nghệ môi trường giảm 40-60% chi phí năng lượng so với phương pháp đèn UV truyền thống.
• Ý nghĩa Xã hội và Môi trường: Đóng góp trực tiếp vào mục tiêu bảo vệ sức khỏe cộng đồng, ngăn chặn sự tích lũy sinh học của DDT và HCH trong chuỗi thức ăn, bảo vệ hệ sinh thái thủy sinh và nguồn nước sinh hoạt tại các vùng nông thôn Việt Nam.
• Tầm vóc Quốc tế: Đóng góp dữ liệu khoa học tin cậy vào nỗ lực chung của cộng đồng quốc tế trong việc thực thi Công ước Stockholm, khẳng định năng lực nghiên cứu đỉnh cao của các nhà khoa học Việt Nam trong lĩnh vực hóa học lượng tử ứng dụng.

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

• Nghiên cứu sinh & Nhà nghiên cứu Hóa học tính toán: Tiếp cận một quy trình mô phỏng hoàn chỉnh, từ việc thiết lập mô hình mạng 2D, tối ưu hóa hình học, quét bề mặt thế năng đến khảo sát động học nhiệt độ và dung môi bằng các phần tử tính toán tối tân.
• Các nhóm nghiên cứu Thực nghiệm Vật liệu & Xúc tác: Sở hữu bộ thông số định lượng chính xác (vị trí dải năng lượng $E_{VB}/E_{CB}$, cấu hình phối trí tối ưu của $Fe, Ni, TiO_2$) giúp rút ngắn thời gian thử nghiệm, tiết kiệm hóa chất và tối ưu hóa quy trình chế tạo vật liệu.
• Kỹ sư Môi trường & Doanh nghiệp Xử lý Chất thải: Nhận được cơ sở khoa học để thiết kế các bể phản ứng quang xúc tác kết hợp hấp phụ hiệu năng cao, ứng dụng nguồn năng lượng mặt trời tự nhiên để phân hủy triệt để thuốc trừ sâu clo hữu cơ.
• Các Cơ quan Quản lý Nhà nước (Bộ TN&MT, Cục Bảo vệ Thực vật): Có thêm luận cứ khoa học vững chắc phục vụ việc đánh giá rủi ro môi trường, xây dựng lộ trình xử lý dứt điểm các kho hóa chất BVTV tồn lưu 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ì?

Đóng góp lý thuyết độc đáo nhất là việc phát hiện và chứng minh cơ chế phân hủy quang xúc tác trực tiếp thông qua chuyển dịch điện tích khử dehalogen hóa. Thay vì tuân theo mô hình oxy hóa gián tiếp truyền thống qua gốc $\cdot\text{OH}$ của AOPs, luận án chứng minh rằng sự phối trí hóa học của DDT và HCH lên các tâm kim loại chuyển tiếp ($Fe, Ni$) và cụm $(TiO_2)_7$ trên $g\text{-C}_3\text{N}_4$ tạo điều kiện cho electron quang sinh từ LUMO của xúc tác truyền thẳng vào orbital phản liên kết $\sigma^*(C\text{-}Cl)$, làm suy yếu tức thì và phá vỡ liên kết $C\text{-}Cl$.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu trước đây (thường chỉ dùng DFT chuẩn mức GGA-PBE với chi phí tính toán lớn, giới hạn số lượng nguyên tử, hoặc dùng phương pháp thực nghiệm mù), luận án đã:

1. Tiên phong áp dụng phương pháp GFN2-xTB với độ chính xác tiệm cận DFT nhưng tốc độ tính toán vượt trội, cho phép khảo sát toàn diện hệ siêu phân tử phức tạp gồm cả chất mang, tâm pha tạp, cơ chất POPs và phân tử dung môi.
2. Tích hợp chuỗi công cụ tự động hóa cao: CREST (quét cấu hình toàn cục) $\rightarrow$ Reaction Path 14 images (mô phỏng đường phản ứng chi tiết) $\rightarrow$ QMD (động lực học phân tử kiểm chứng độ bền nhiệt).

3. Phát hiện nào gây bất ngờ nhất về mặt dữ liệu?

Phát hiện bất ngờ nhất là sự tương phản tuyệt đối về bản chất hấp phụ giữa $g\text{-C}_3\text{N}_4$ nguyên bản và $g\text{-C}_3\text{N}_4$ biến tính kim loại: Trong khi $g\text{-C}_3\text{N}_4$ nguyên bản hoàn toàn trơ về mặt hóa học đối với DDT/HCH (chỉ xảy ra hấp phụ vật lý yếu $vdW$, liên kết $C\text{-}Cl$ giữ nguyên $1.78\text{ \AA}$), thì sự xuất hiện của chỉ một nguyên tử $Ni$ hoặc $Fe$ đã chuyển hóa hoàn toàn hệ sang hấp phụ hóa học cực mạnh, kéo dài liên kết $C\text{-}Cl$ lên tới trên $2.05\text{ \AA}$, đưa phân tử vào trạng thái tiền phân ly tự phát.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án mô tả chi tiết toàn bộ điều kiện biên và tham số tính toán: sử dụng chương trình xTB với bộ tham số chuẩn hóa cho các nguyên tố đến $Z=86$, năng lượng hội tụ $0.1 \times 10^{-5}\text{ Hartree}$, tiêu chuẩn hàm sóng $0.1 \times 10^{-3}\text{ e}$, nhiệt độ Fermi electron $300\text{ K}$, thang chuẩn thế năng chân không $E_e = 4.500\text{ eV}$, cùng tọa độ không gian và quy trình chạy MD bước tích phân $1.0\text{ fs}$ có thể tái lập 100% trên các hệ thống tính toán hiệu năng cao.

5. Kế hoạch nghiên cứu 10 năm được vạch ra như thế nào?

Kế hoạch nghiên cứu 10 năm hướng tới lộ trình 3 giai đoạn:

• Năm 1-3: Kiểm chứng lý thuyết định lượng cao bằng Periodic DFT (HSE06) và mô phỏng dung môi tường minh QM/MM.
• Năm 4-6: Tổng hợp thực nghiệm nanocomposite $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$, kiểm tra hoạt tính phân hủy DDT/HCH thực tế dưới ánh sáng mặt trời trong phòng thí nghiệm.
• Năm 7-10: Chế tạo thiết bị phản ứng quang xúc tác dạng module xử lý nước và đất quy mô pilot, tiến tới thương mại hóa và ứng dụng thực địa tại các điểm nóng ô nhiễm POPs tại Việt Nam.

Kết luận

1. Xác lập cơ sở khoa học lượng tử toàn diện: Luận án đã mô hình hóa thành công cấu trúc hình học, cấu trúc electron và tính chất quang của $g\text{-C}_3\text{N}_4$ dạng phẳng (pGN) và gợn sóng (cGN), chỉ rõ ảnh hưởng quyết định của 6 kim loại pha tạp ($K, Ca, Ga, Fe, Ni, Cu$) và 2 oxide kim loại ($(ZnO)_3, (TiO_2)_7$) đến việc thu hẹp độ rộng vùng cấm ($E_g$) và tối ưu hóa dải hấp thụ quang.
2. Khám phá hệ vật liệu tối ưu: Đã chứng minh $Fe/g\text{-C}_3\text{N}_4$, $Ni/g\text{-C}_3\text{N}_4$ và $TiO_2/g\text{-C}_3\text{N}_4$ là các ứng viên xúc tác quang vượt trội nhất, có khả năng phân tách không gian hiệu quả cặp electron - lỗ trống quang sinh và hấp thu mạnh mẽ bức xạ trong vùng ánh sáng nhìn thấy.
3. Làm sáng tỏ bản chất tương tác hấp phụ: Đã xác định quá trình hấp phụ DDT và $\gamma$-HCH trên $g\text{-C}_3\text{N}_4$ nguyên bản mang bản chất hấp phụ vật lý, trong khi trên các hệ biến tính $Fe, Ni, TiO_2$ đã chuyển hóa hoàn toàn sang hấp phụ hóa học với năng lượng liên kết cao và sự chuyển dịch điện tích sâu sắc.
4. Đề xuất cơ chế phân hủy quang xúc tác đột phá: Đã phát hiện và mô hình hóa thành công cơ chế nhận electron trực tiếp từ xúc tác quang gây kéo dài và phân cắt liên kết $C\text{-}Cl$, mở ra một hướng tiếp cận hoàn toàn mới trong lý thuyết xử lý chất ô nhiễm clo hữu cơ.
5. Chứng minh tính ổn định động học và môi trường: Đã chứng minh độ bền nhiệt động của các cấu hình hấp phụ thông qua mô phỏng động lực học phân tử MD và khẳng định tính ổn định cấu trúc phức hệ trong nhiều môi trường dung môi phân cực khác nhau.
6. Mở ra 3 hướng nghiên cứu mới: Khởi xướng các dòng nghiên cứu tiên phong về: (i) Thiết kế xúc tác quang đơn nguyên tử (Single-atom catalysts) trên nền mạng carbon nitride 2D; (ii) Ứng dụng quy trình GFN2-xTB/CREST để sàng lọc vật liệu phân hủy các chất ô nhiễm thế hệ mới; và (iii) Phát triển công nghệ tích hợp hấp phụ - quang xúc tác dòng liên tục định hướng bảo vệ môi trường bền vững.