Tổng quan nghiên cứu

Theo thống kê từ Cơ quan Bảo vệ Môi trường Hoa Kỳ, khí carbon dioxide chiếm tới 82% tổng lượng phát thải khí nhà kính và 85% nguồn năng lượng toàn cầu vẫn đang phụ thuộc vào việc đốt nhiên liệu hóa thạch. Cùng với đó, sulfur dioxide là loại khí ô nhiễm nguy hại hàng đầu với 73% lượng phát sinh từ hoạt động năng lượng và 20% từ các quy trình sản xuất công nghiệp, trực tiếp gây ra mưa axit và các bệnh lý hô hấp nghiêm trọng. Vấn đề cấp bách hiện nay là tìm kiếm giải pháp xử lý đồng thời hai loại khí này tại các nguồn xả thải cố định.

Luận văn thạc sĩ chuyên ngành Vật lý Kỹ thuật tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh thực hiện từ tháng 01/2015 đến tháng 11/2015 đã tập trung giải quyết bài toán trên. Mục tiêu trọng tâm của đề tài là khảo sát chi tiết cơ chế hấp phụ đơn lẻ và hấp phụ đồng thời hai khí SO2 và CO2 bên trong vật liệu khung hữu cơ kim loại Ni(bdc)(ted)0,5 bằng phương pháp mô phỏng lượng tử. Nghiên cứu mang ý nghĩa then chốt trong việc giải thích bản chất tương tác vi mô, hỗ trợ định hướng thiết kế vật liệu xốp có khả năng bắt giữ SO2 đạt mức 9,97 mmol/g và lưu trữ CO2 lên tới 60 wt% ở áp suất 25 bar, góp phần tối ưu hóa hiệu quả bảo vệ môi trường và giảm chi phí xử lý khí thải công nghiệp.

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 lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) với hai định lý Hohenberg-Kohn và hệ phương trình Kohn-Sham tự hợp để giải quyết bài toán đa hạt trong cơ học lượng tử. Để mô tả chính xác tương tác tầm xa giữa các phân tử khí và khung vật liệu, nghiên cứu áp dụng xấp xỉ gradient tổng quát (GGA) kết hợp phiếm hàm mật độ có hiệu chỉnh tương tác van der Waals (vdW-DF) dựa trên thế trao đổi tương quan revPBE.

Ba khái niệm lý thuyết cốt lõi được sử dụng bao gồm: mật độ trạng thái điện tử (DOS), sự biến thiên mật độ điện tích (CDD) và phân tích điện tích định lượng Bader. Mô hình cấu trúc mạng tinh thể ba chiều của Ni(bdc)(ted)0,5 bao gồm các cụm oxit kim loại nickel liên kết với cầu nối hữu cơ benzene dicarbonxylate (bdc) và triethylenediamine (ted), tạo thành mạng lưới mao quản có thể tích rỗng tự do chiếm tới 90% tổng thể tích vật liệu.

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

Nghiên cứu ứng dụng gói phần mềm mô phỏng Vienna Ab initio Simulation Package (VASP) sử dụng phương pháp sóng chiếu tăng cường (PAW) kết hợp bộ cơ sở sóng phẳng. Ô cơ sở mô phỏng có kích thước hình học 11,15 × 11,15 × 9,53 ų với cỡ mẫu tính toán chuẩn hóa gồm đúng 54 nguyên tử. Phương pháp chọn mẫu cấu hình được tiến hành bằng cách quét toàn diện bốn vị trí hấp phụ khả dĩ: cụm oxit kim loại, vòng benzen, cầu nối ted và vị trí trung gian benzen-ted.

Không gian xung lượng vùng Brillouin được lấy tích phân trên lưới điểm k Monkhorst-Pack kích thước 5 × 5 × 5. Mức năng lượng cắt sóng phẳng được thiết lập ở 700 eV nhằm đảm bảo tiêu chuẩn hội tụ năng lượng tuyệt đối. Lý do lựa chọn phương pháp DFT trên VASP là nhờ khả năng tính toán chuẩn xác từ nguyên lý ban đầu, cho phép phân tích sâu các orbital phân tử mà thực nghiệm không thể quan sát trực tiếp, đồng thời tiết kiệm hơn 80% thời gian và chi phí so với việc thử nghiệm trong phòng lab.

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

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

Thứ nhất, năng lượng hấp phụ đơn lẻ của SO2 trong khung mạng Ni(bdc)(ted)0,5 dao động từ -0,351 eV đến -1,010 eV, cao hơn từ 64% đến 165% so với mức năng lượng của CO2 (dao động từ -0,213 eV đến -0,380 eV), khẳng định ái lực của SO2 mạnh vượt trội tại mọi vị trí.

Thứ hai, vị trí hấp phụ bền vững nhất của cả hai chất khí đều tập trung tại cụm oxit kim loại nickel (-1,010 eV đối với SO2 và -0,380 eV đối với CO2), trong khi vị trí trên vòng benzen là kém bền vững nhất (-0,351 eV cho SO2 và -0,213 eV cho CO2).

Thứ ba, phân tích điện tích Bader xác nhận nguyên tử carbon của CO2 nhường từ -0,019e đến -0,066e, còn hai nguyên tử oxy nhận từ +0,035e đến +0,087e; khoảng cách tương tác O(CO2)---H(benzen) đo được là 3,46 Å và C(CO2)---C(benzen) là 3,47 Å.

Thứ tư, khi hấp phụ đồng thời, sự xuất hiện của CO2 làm thay đổi không đáng kể cấu trúc hình học và năng lượng liên kết của SO2, chứng minh SO2 luôn chiếm giữ ưu tiên các tâm hoạt tính bền nhất.

Thảo luận kết quả

Các dữ liệu tính toán được biểu diễn trực quan thông qua bảng đối chiếu năng lượng liên kết và đồ thị mật độ trạng thái điện tử (DOS) với mức Fermi quy chuẩn về 0 eV. Bản chất khác biệt về năng lượng hấp phụ giữa hai khí bắt nguồn từ cấu trúc điện tử: phân tử SO2 có mômen lưỡng cực lớn, tạo liên kết tĩnh điện mạnh với các orbital của nickel và oxy trên khung mạng. Ngược lại, phân tử CO2 có phân bố tứ cực nên chỉ tương tác yếu qua liên kết hydro dạng (CO2)O---H(ted) và liên kết tĩnh điện nhẹ.

So sánh với các nghiên cứu khác, dung lượng hấp phụ SO2 của Ni(bdc)(ted)0,5 đạt mức 9,97 mmol/g tại 298 K và 1,13 bar, vượt trội hơn vật liệu NOTT-300 (8,1 mmol/g) và Mg-MOF-74 (năng lượng hấp phụ đạt 73 kJ/mol). Khi quan sát trên biểu đồ phân bố mật độ điện tích (CDD) với mặt đẳng trị 5 × 10⁻³ e/Bohr³, sự dịch chuyển trạng thái 2σu/1πu tại vùng năng lượng -8 eV đã lý giải chính xác cơ chế cạnh tranh hấp phụ mà không gây biến dạng cấu trúc khung tinh thể ban đầu.

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

  1. Tối ưu hóa cấu trúc vật liệu: Các nhà khoa học thuộc nhóm nghiên cứu vật liệu cần tiến hành biến tính bề mặt khung hữu cơ bằng các nhóm chức phân cực nhằm tăng dung lượng lưu trữ CO2 từ mức 10 wt% lên trên 35 wt% ở điều kiện 1 bar trong giai đoạn 2026-2027.

  2. Mở rộng quy mô mô phỏng: Đội ngũ kỹ sư tính toán lượng tử cần thiết lập mô hình mô phỏng dòng khí đa cấu tử gồm SO2, CO2, hơi nước và N2 dưới dải áp suất thay đổi từ 1 đến 50 bar, hoàn thành trước quý 4 năm 2026.

  3. Chế tạo màng lọc composite: Các doanh nghiệp công nghệ môi trường cần hợp tác cùng viện nghiên cứu để gia công thử nghiệm màng lọc Ni(bdc)(ted)0,5 có độ bền đạt trên 50 chu kỳ hấp phụ và giải hấp nhiệt độ thấp vào năm 2027.

  4. Ứng dụng xử lý khí thải nguồn điểm: Ban quản lý các nhà máy nhiệt điện đốt than cần lắp đặt thử nghiệm hệ thống hấp phụ khí sử dụng vật liệu MOF để thu hồi trên 95% lượng SO2 và 40% lượng CO2 từ năm 2028.

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

  1. Giảng viên và học viên cao học chuyên ngành Vật lý Kỹ thuật, Khoa học Vật liệu: Khai thác chi tiết phương pháp thiết lập ô cơ sở 54 nguyên tử và các thông số kỹ thuật VASP với mức năng lượng cắt 700 eV phục vụ nghiên cứu cấu trúc điện tử.

  2. Kỹ sư môi trường và chuyên gia xử lý khí thải: Sử dụng số liệu nhiệt động học của hệ MOF để thiết kế tháp hấp phụ tách lọc hỗn hợp khí chứa hơn 70% nồng độ SO2 phát thải từ các lò đốt công nghiệp.

  3. Nhà khoa học tổng hợp vật liệu nano xốp: Vận dụng cơ chế tương tác orbital giữa cụm kim loại nickel và các linker bdc/ted để tổng hợp các thế hệ khung hữu cơ kim loại mới có dung lượng hấp phụ vượt mốc 10 mmol/g.

  4. Doanh nghiệp sản xuất thiết bị lọc khí công nghệ cao: Tham khảo dữ liệu về tính chọn lọc và độ bền nhiệt của Ni(bdc)(ted)0,5 nhằm xây dựng phương án thương mại hóa vật liệu trong lộ trình 2026-2030.

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

Vật liệu MOF Ni(bdc)(ted)0,5 có ưu điểm gì nổi bật so với các vật liệu xốp truyền thống? Vật liệu sở hữu cấu trúc tinh thể ba chiều với thể tích lỗ xốp tự do chiếm tới 90% và diện tích bề mặt lớn hơn nhiều so với zeolite (500-800 m²/g). Ni(bdc)(ted)0,5 đạt dung lượng hấp phụ SO2 ấn tượng lên tới 9,97 mmol/g ở điều kiện 298 K, đồng thời có khả năng tái sinh cấu trúc thuận nghịch dễ dàng mà không làm tiêu hao nhiều năng lượng.

Tại sao năng lượng hấp phụ của SO2 lại âm hơn đáng kể so với CO2? Năng lượng hấp phụ của SO2 đạt mức -1,010 eV tại tâm oxit kim loại, cao hơn nhiều so với mức -0,380 eV của CO2. Nguyên nhân là do phân tử SO2 có mômen lưỡng cực vĩnh viễn và mật độ điện tích tập trung, tạo ra tương tác tĩnh điện mạnh với các orbital của nickel, trong khi CO2 chủ yếu tương tác tứ cực yếu hơn.

Phương pháp DFT đóng vai trò như thế nào trong việc nghiên cứu tương tác khí và MOF? Phương pháp DFT kết hợp hiệu chỉnh vdW-DF và năng lượng cắt 700 eV cho phép mô tả chính xác tương tác lượng tử ở cấp độ orbital. Tính toán mô phỏng giúp xác định chính xác vị trí hấp phụ tối ưu, sự chuyển dịch điện tích Bader và mật độ trạng thái điện tử DOS, giúp tiết kiệm hơn 80% thời gian nghiên cứu thực nghiệm.

Khí CO2 có làm suy giảm khả năng bắt giữ SO2 khi cả hai cùng xuất hiện không? Kết quả mô phỏng cho thấy sự xuất hiện của CO2 không làm thay đổi trật tự vị trí ưu tiên của SO2 trên cụm oxit kim loại. Do năng lượng hấp phụ của SO2 vượt trội hơn CO2 từ 64% đến 165%, SO2 luôn chiếm giữ các vị trí liên kết bền nhất, đảm bảo hiệu suất chọn lọc cao của màng vật liệu.

Dung lượng lưu trữ khí CO2 của vật liệu biến thiên như thế nào theo điều kiện áp suất? Ở điều kiện tiêu chuẩn 298 K và áp suất 1 bar, Ni(bdc)(ted)0,5 đạt dung lượng hấp phụ CO2 khoảng 10 wt%. Khi tăng áp suất lên mức 25 bar, dung lượng hấp phụ tăng mạnh đạt 60 wt%, tương đương với các dòng vật liệu xốp hàng đầu như MOF-177 (60,8 wt%) và MIL-101 (56,9 wt%).

Kết luận

  • Luận văn đã xác lập hoàn chỉnh mô hình tính toán lượng tử DFT mô tả tương tác hấp phụ của SO2 và CO2 trong mạng tinh thể Ni(bdc)(ted)0,5 với ô cơ sở 54 nguyên tử.
  • Chứng minh độ bền hấp phụ vượt trội của SO2 so với CO2 với mức chênh lệch năng lượng liên kết cực đại từ -0,380 eV lên -1,010 eV tại tâm oxit kim loại.
  • Làm sáng tỏ bản chất chuyển dịch điện tích Bader và sự biến thiên orbital tại vùng trạng thái 2σu/1πu ở mức năng lượng -8 eV.
  • Khẳng định tính chọn lọc bền vững của vật liệu đối với SO2 trong điều kiện hấp phụ đồng thời cùng CO2.
  • Cung cấp bộ thông số vi mô chuẩn xác phục vụ thiết kế hệ thống màng lọc xử lý khí thải đạt hiệu suất bắt giữ trên 95%.

Lộ trình tiếp theo hướng tới việc mở rộng tính toán động học phân tử và thử nghiệm pilot trên dòng khí thải thực tế trong giai đoạn 2026-2028. Các nhà nghiên cứu và doanh nghiệp quan tâm có thể khai thác toàn văn luận văn thạc sĩ để ứng dụng các thông số tính toán vào quy trình phát triển vật liệu xử lý môi trường.