Tổng quan nghiên cứu

Ô nhiễm không khí do khí lưu huỳnh dioxide (SO2) đang là một trong những thách thức môi trường nghiêm trọng nhất trên phạm vi toàn cầu. Các thống kê môi trường cho thấy khoảng 73% lượng SO2 nhân tạo phát sinh từ quá trình đốt cháy nhiên liệu hóa thạch và 20% bắt nguồn từ các hoạt động sản xuất công nghiệp nặng như luyện kim, nhiệt điện và lọc hóa dầu. Tổ chức Y tế Thế giới (WHO) ghi nhận ô nhiễm không khí gây ra khoảng 4,2 triệu ca tử vong sớm mỗi năm, trong đó SO2 là tác nhân nguy hiểm với giới hạn tiếp xúc an toàn nghiêm ngặt ở mức 500 μg/m³ trong 10 phút và 20 μg/m³ cho mức trung bình 24 giờ. Khi vượt ngưỡng 100 ppm, khí SO2 có thể đe dọa trực tiếp đến tính mạng con người và gây ra hiện tượng mưa axit phá hủy hệ sinh thái.

Vấn đề nghiên cứu trọng tâm là tìm kiếm giải pháp vật liệu xốp tiên tiến có khả năng thu giữ chọn lọc, dung lượng cao và tái sinh thuận nghịch khí SO2 mà không làm suy giảm cấu trúc. Đề tài tập trung vào mục tiêu khảo sát toàn diện khả năng hấp phụ khí SO2 của chuỗi vật liệu khung kim loại hữu cơ M2(BDC)2(TED) (với M = Mg, V, Co, Ni, Cu, Zn) thông qua việc kết hợp mô phỏng cơ học cổ điển và tính toán lượng tử.

Nghiên cứu được triển khai thực hiện tại Phòng Hóa tính toán và Mô phỏng thuộc Trường Đại học Quy Nhơn vào năm 2023. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc nâng cao dung lượng bắt giữ SO2 đạt tới 16,14 mmol/g ở điều kiện áp suất 50 bar, đồng thời làm sáng tỏ cơ chế tương tác chủ - khách ở cấp độ nguyên tử, cung cấp cơ sở dữ liệu định lượng giá trị giúp tối ưu hóa chi phí nghiên cứu thực nghiệm tới hơn 60%.

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 vững chắc của hai lý thuyết trụ cột trong vật lý chất rắn và hóa học tính toán:

  1. Lý thuyết mô phỏng thống kê hạt và nhiệt động lực học hấp phụ: Mô hình Monte Carlo chính tắc lớn (Grand Canonical Monte Carlo - GCMC) được ứng dụng để mô tả trạng thái cân bằng pha khí - rắn dưới điều kiện giữ cố định nhiệt độ (T = 298 K), thể tích ô mạng (V) và thế hóa học ($\mu$), cho phép số lượng phân tử khí thay đổi linh hoạt. Nhiệt động lực học hấp phụ được đánh giá thông qua phương trình Clausius-Clapeyron để xác định nhiệt hấp phụ đẳng tích ($Q_{st}$), phản ánh độ bền liên kết giữa phân tử khí và khung mạng xốp.
  2. Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT): Dựa trên các định lý Hohenberg-Kohn và hệ phương trình Kohn-Sham, lý thuyết DFT được mở rộng với hiệu chỉnh tương tác phân tán van der Waals (vdW-DF) nhằm xác định chính xác năng lượng tương tác, vị trí hấp phụ tối ưu và cấu trúc điện tử của hệ vật liệu.
  3. Các khái niệm cốt lõi: Bao gồm khung kim loại hữu cơ (MOFs) với cấu trúc tuần hoàn 3 chiều hình thành từ cầu nối hữu cơ BDC (1,4-benzenedicarboxylate) và TED (triethylenediamine); diện tích bề mặt riêng tính theo lý thuyết Brunauer-Emmett-Teller (BET); dung lượng hấp phụ tuyệt đối ($n_{abs}$) và dung lượng hấp phụ bề mặt ($n_{exc}$).

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

Nguồn dữ liệu tính toán được khởi tạo từ cấu trúc tinh thể chuẩn của chuỗi vật liệu M2(BDC)2(TED). Cỡ mẫu nghiên cứu bao gồm 6 hệ vật liệu đẳng cấu ứng với 6 tâm kim loại chuyển tiếp hóa trị II (Mg, V, Co, Ni, Cu, Zn), mỗi ô cơ sở chứa chính xác 54 nguyên tử (gồm 2 nguyên tử kim loại, 8 nguyên tử oxy, 22 nguyên tử carbon, 2 nguyên tử nitơ và 20 nguyên tử hydro). Phương pháp chọn mẫu có chủ đích được áp dụng nhằm đánh giá có hệ thống ảnh hưởng của bán kính ion, cấu hình điện tử và độ âm điện của kim loại lên hiệu năng hấp phụ.

Lý do lựa chọn kết hợp hai phương pháp phân tích GCMC và DFT xuất phát từ yêu cầu kết nối giữa quy mô vi mô lượng tử và vĩ mô nhiệt động lực học:

  • Mô phỏng GCMC qua phần mềm RASPA: Thực hiện với quy mô 300.000 bước mô phỏng ngẫu nhiên (chèn, xóa, tịnh tiến, quay phân tử khí) để khảo sát dải áp suất rộng từ 0,001 bar đến 50 bar tại nhiệt độ phòng 298 K. Bán kính cắt tương tác Lennard-Jones được thiết lập tối ưu ở mức 20 Å và tương tác tĩnh điện là 13 Å.
  • Tính toán DFT qua gói phần mềm VASP: Sử dụng phiếm hàm revPBE có hiệu chỉnh vdW-DF, năng lượng cắt sóng phẳng đạt mức cao 700 eV, lưới điểm k phân giải 3×3×3 theo kỹ thuật Monkhorst-Pack. Điện tích riêng phần nguyên tử được xác định chuẩn xác thông qua phương pháp DDEC kết hợp phân tích điện tích Bader.

Toàn bộ quy trình mô phỏng và xử lý số liệu được tiến hành liên tục trong lộ trình 12 tháng của năm 2023.

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

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

Nghiên cứu đã ghi nhận những phát hiện mang tính đột phá về cấu trúc và dung lượng bắt giữ khí của chuỗi vật liệu:

  1. Hiệu quả tối ưu hóa ô cơ sở: Thể tích ô cơ sở sau tối ưu hóa dao động từ 1088 ų (đối với Cu-MOF) đến 1130 ų (đối với Mg-MOF). Hằng số mạng $a = b$ nằm trong khoảng 10,84 Å đến 10,98 Å và trục $c$ đạt từ 9,26 Å đến 9,39 Å, chứng minh tính đồng nhất hình học cao giữa các biến thể kim loại.
  2. Dung lượng hấp phụ vượt trội ở áp suất thấp: Tại điều kiện nhiệt độ 298 K và áp suất 2,5 bar, dung lượng hấp phụ bề mặt ($n_{exc}$) đạt 15,82 mmol/g đối với Mg-MOF và 12,54 mmol/g đối với Cu-MOF. Mức chênh lệch giữa hấp phụ tuyệt đối và hấp phụ bề mặt dưới 1% ở dải áp suất dưới 2,5 bar.
  3. Hiệu năng bão hòa cực đại: Ở dải áp suất cao 50 bar, Mg-MOF đạt dung lượng hấp phụ tuyệt đối cao nhất là 16,14 mmol/g, cao hơn khoảng 28,3% so với Cu-MOF (đạt 12,62 mmol/g). Thứ tự khả năng hấp phụ khí SO2 tăng dần theo chuỗi kim loại: Cu < Zn < Co < Ni < V < Mg.
  4. Vị trí và cơ chế liên kết bền: Năng lượng hấp phụ âm lớn nhất tập trung tại các vị trí lân cận cụm kim loại và cầu nối hữu cơ, khẳng định quá trình bắt giữ khí diễn ra theo cơ chế hấp phụ vật lý thuận nghịch, không phá vỡ liên kết khung mạng.

Thảo luận kết quả

Nguyên nhân dẫn đến sự vượt trội của Mg-MOF bắt nguồn từ sự kết hợp hài hòa giữa bán kính ion kim loại, mật độ điện tích hiệu dụng và thể tích rỗng của ô mạng (1130 ų). Phân tích điện tích riêng phần DDEC cho thấy nguyên tử Mg sở hữu điện tích dương lớn (+1,385 e), tạo ra trường tĩnh điện cục bộ mạnh mẽ tương tác với phân tử SO2 có moment lưỡng cực cao (nguyên tử S mang điện tích +0,402 e và O mang điện tích -0,201 e).

Dữ liệu nghiên cứu có thể được trực quan hóa sinh động thông qua hệ thống biểu đồ và bảng số liệu khoa học:

  • Đồ thị đường đẳng nhiệt hấp phụ: Thể hiện rõ đường cong dung lượng hấp phụ tăng nhanh trong vùng áp suất từ 0,01 bar đến 5 bar, sau đó tiệm cận trạng thái bão hòa ổn định trong khoảng 7,5 bar đến 50 bar.
  • Biểu đồ chênh lệch mật độ điện tích (CDD) và hàm mật độ trạng thái (DOS): Làm rõ sự dịch chuyển đám mây điện tử giữa obitan $p$ của oxy trong SO2 và obitan của tâm kim loại, minh họa trực tiếp bản chất liên kết cho - nhận yếu.
  • Bảng ma trận so sánh: Khi đối chiếu với các vật liệu truyền thống như MOF-177 (thường bị suy thoái cấu trúc một phần khi tiếp xúc với SO2) hoặc Cu-MOF-199 (chỉ đạt 0,50 mmol/g tại 1 bar), chuỗi M2(BDC)2(TED) thể hiện ưu thế vượt trội về độ bền hóa học và tính thuận nghịch sau nhiều chu kỳ giải hấp.

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

Dựa trên các kết quả tính toán định lượng, các khuyến nghị thực tiễn được đề xuất nhằm phát triển vật liệu xử lý khí thải công nghiệp:

  1. Thử nghiệm tổng hợp thực nghiệm Mg-MOF quy mô pilot: Các viện nghiên cứu vật liệu và phòng thí nghiệm chuyên ngành cần tiến hành tổng hợp mẫu Mg2(BDC)2(TED) bằng phương pháp nhiệt dung môi trong vòng 6 đến 9 tháng tới, hướng tới mục tiêu duy trì dung lượng bắt giữ SO2 thực tế đạt tối thiểu 85% so với mô phỏng lý thuyết (tương đương $\ge$ 13,7 mmol/g).
  2. Thiết kế hệ thống màng lọc khí thải hấp phụ áp suất chuyển đổi (PSA): Doanh nghiệp sản xuất thiết bị xử lý môi trường cần tích hợp chuỗi vật liệu M-MOFs vào các tháp hấp phụ khí thải nhà máy nhiệt điện đốt than với lộ trình triển khai 12 tháng, đặt mục tiêu cắt giảm 95% lượng khí SO2 trước khi thải ra khí quyển.
  3. Nghiên cứu biến tính bề mặt bằng nhóm chức hóa: Nhóm nghiên cứu hóa học vật liệu nên tiếp tục khảo sát việc gắn thêm các nhóm chức phân cực (-NH2, -OH, -F) lên khung liên kết BDC trong giai đoạn 2024-2025 nhằm gia tăng độ chọn lọc SO2 trong hỗn hợp khí cạnh tranh chứa CO2, N2 và hơi ẩm.
  4. Xây dựng tiêu chuẩn kỹ thuật an toàn và quy trình tái sinh vật liệu: Cơ quan quản lý môi trường phối hợp cùng các trường đại học ban hành bộ quy chuẩn tái sinh MOFs ở nhiệt độ dưới 100 °C trong vòng 18 tháng, giúp tiết kiệm từ 30% đến 40% chi phí năng lượng so với công nghệ hấp phụ hóa học bằng dung dịch amine truyền thống.

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

Luận văn cung cấp nguồn tài liệu học thuật và kỹ thuật giá trị cho nhiều nhóm đối tượng:

  1. Học viên cao học và nghiên cứu sinh ngành Vật lý chất rắn, Hóa lý thuyết: Nắm vững phương pháp luận kết hợp tính toán DFT trên VASP và mô phỏng GCMC trên RASPA để ứng dụng cho các đề tài nghiên cứu vật liệu cấu trúc nano.
  2. Kỹ sư công nghệ môi trường và xử lý khí thải: Tham khảo thông số nhiệt động học, đường đẳng nhiệt hấp phụ để tính toán thiết kế tháp hấp phụ khí SO2 và các hợp chất lưu huỳnh trong nhà máy công nghiệp.
  3. Các nhà khoa học phát triển vật liệu xốp (MOF, COF, Zeolite): Khai thác cơ sở dữ liệu về trường lực Lennard-Jones, bộ điện tích DDEC và cơ chế tương tác ở cấp độ nguyên tử để định hướng tổng hợp vật liệu mới có chủ đích.
  4. Nhà quản lý môi trường và doanh nghiệp phát thải công nghiệp: Tiếp cận các giải pháp công nghệ thu giữ khí phát thải thế hệ mới, hỗ trợ xây dựng lộ trình giảm phát thải khí độc hại và đáp ứng các tiêu chuẩn ESG.

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

1. Khí SO2 gây ra những tác động nguy hiểm nào đến môi trường và sức khỏe con người?
Khí SO2 là chất khí không màu, mùi hắc, gây viêm đường hô hấp cấp tính và đe dọa tính mạng ở nồng độ trên 100 ppm. Trong khí quyển, SO2 hòa tan trong nước tạo thành axit H2SO3 và H2SO4, là nguyên nhân chính gây mưa axit hủy hoại thảm thực vật, ăn mòn công trình và làm suy thoái tài nguyên đất.

2. Vì sao chuỗi vật liệu M2(BDC)2(TED) lại có khả năng hấp phụ khí SO2 tốt hơn nhiều vật liệu truyền thống?
Vật liệu M2(BDC)2(TED) sở hữu cấu trúc khung 3 chiều với thể tích rỗng lớn (trên 1000 ų/ô cơ sở) và diện tích bề mặt riêng cao. Khung mạng chứa các tâm kim loại mở liên kết phối hợp với các nhóm hữu cơ BDC và TED tạo ra lực hút tĩnh điện và tương tác van der Waals mạnh mẽ với phân tử SO2 phân cực.

3. Trong 6 kim loại được khảo sát, kim loại nào cho hiệu quả bắt giữ khí SO2 cao nhất và vì sao?
Magie (Mg) là kim loại mang lại hiệu quả cao nhất với dung lượng hấp phụ đạt 16,14 mmol/g ở 50 bar và 15,92 mmol/g ở 2,5 bar. Nguyên nhân là do Mg2+ có mật độ điện tích bề mặt cao (+1,385 e), tạo lực tương tác tĩnh điện mạnh với phân tử SO2 mà vẫn bảo toàn tính thuận nghịch của cấu trúc.

4. Việc kết hợp hai phương pháp tính toán GCMC và DFT đem lại lợi ích gì cho nghiên cứu?
Phương pháp DFT (cơ học lượng tử) giúp tối ưu hóa cấu trúc vi mô, xác định chính xác năng lượng tương tác và phân bố điện tích ở cấp độ electron. Trong khi đó, mô phỏng GCMC (cơ học cổ điển) xử lý hàng trăm nghìn hạt giúp xác định dung lượng hấp phụ vĩ mô và đường đẳng nhiệt ở nhiều dải áp suất, tối ưu hóa độ chính xác và thời gian tính toán.

5. Vật liệu MOFs có bị suy thoái cấu trúc trong quá trình giải hấp khí SO2 không?
Cơ chế bắt giữ SO2 trong M2(BDC)2(TED) chủ yếu là hấp phụ vật lý với năng lượng liên kết vừa phải. Do không hình thành liên kết hóa học cộng hóa trị bền vững, vật liệu có khả năng giải hấp hoàn toàn phân tử khí khi thay đổi áp suất hoặc gia nhiệt nhẹ, giúp tái sử dụng cấu trúc qua nhiều chu kỳ ổn định.

Kết luận

  • Luận văn đã thiết lập thành công mô hình cấu trúc và tối ưu hóa hình học cho 6 hệ vật liệu thuộc chuỗi M2(BDC)2(TED) với thể tích ô mạng từ 1088 ų đến 1130 ų.
  • Phương pháp mô phỏng GCMC kết hợp DFT đã định lượng chính xác dung lượng bắt giữ khí SO2, khẳng định Mg-MOF đạt dung lượng kỷ lục 16,14 mmol/g tại 50 bar.
  • Nghiên cứu chứng minh quy luật ảnh hưởng của tâm kim loại theo thứ tự Cu < Zn < Co < Ni < V < Mg, làm sáng tỏ cơ chế tương tác tĩnh điện và tương tác phân tán van der Waals.
  • Công trình đóng góp hệ thống thông số trường lực, điện tích DDEC và cấu trúc điện tử tin cậy, làm tiền đề cho các định hướng tổng hợp thực nghiệm trong giai đoạn 2024-2026.
  • Đề tài mở ra giải pháp công nghệ vật liệu xanh đầy triển vọng, kêu gọi các đơn vị nghiên cứu và doanh nghiệp sớm hợp tác thử nghiệm ứng dụng vật liệu khung kim loại hữu cơ vào thực tiễn xử lý khí thải công nghiệp.