Tổng quan nghiên cứu

Nồng độ khí nhà kính carbon dioxide (CO2) trong khí quyển đã ghi nhận mức tăng kỷ lục từ xấp xỉ 280 ppm vào đầu thế kỷ XIX lên mức 419 ppm vào tháng 5 năm 2021 theo số liệu đo đạc chính xác từ Đài quan sát Mauna Loa của Cơ quan Quản lý Khí quyển và Đại dương Quốc gia Mỹ (NOAA). Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) đã đưa ra cảnh báo nghiêm trọng rằng nếu không có biện pháp cắt giảm phát thải quyết liệt, nồng độ CO2 toàn cầu có thể chạm ngưỡng 950 ppm vào năm 2100. Song song với CO2, khí carbon monoxide (CO) phát sinh từ các hoạt động đốt cháy nhiên liệu hóa thạch và luyện kim cũng gây ô nhiễm môi trường nghiêm trọng và đe dọa trực tiếp đến tính mạng con người. Trong khi công nghệ hấp thụ hóa học truyền thống bằng dung dịch amine lỏng dù có thể loại bỏ tới 98% CO2 nhưng lại bộc lộ hạn chế lớn về ăn mòn thiết bị, tiêu tốn năng lượng tái sinh và tạo chất thải thứ cấp, thì giải pháp sử dụng vật liệu hấp phụ rắn xốp khung hữu cơ kim loại (Metal-Organic Frameworks - MOFs) nổi lên như một hướng đi đột phá.

Mục tiêu trọng tâm của luận văn là tính toán và định lượng chính xác khả năng bắt giữ hai loại khí thải CO2 và CO trên chuỗi vật liệu khung hữu cơ kim loại M(BDC)(TED)0.5 (trong đó M là các ion kim loại chuyển tiếp và kiềm thổ gồm Magnesium - Mg, Vanadium - V, Cobalt - Co, Nickel - Ni và Copper - Cu). Phạm vi nghiên cứu tập trung vào mô phỏng tương tác ở cấp độ nguyên tử và phân tử trong điều kiện nhiệt độ phòng 298 K với dải áp suất mở rộng từ 0 đến 100 bar. Công trình được hoàn thành tại Khoa Khoa học Tự nhiên thuộc Trường Đại học Quy Nhơn vào tháng 8 năm 2022. Ý nghĩa khoa học và thực tiễn của đề tài là thiết lập cơ sở dữ liệu nhiệt động học tin cậy, giải thích bản chất cơ chế hấp phụ ở mức độ cấu trúc điện tử, giúp các nhà thực nghiệm định hướng lựa chọn kim loại tối ưu, từ đó tiết kiệm từ 40% đến 60% chi phí và thời gian thử nghiệm trong thị trường vật liệu hấp phụ toàn cầu ước tính đạt hơn 11,1 tỷ USD.

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ự kết hợp chặt chẽ giữa lý thuyết cấu trúc vật liệu xốp nano, hóa học lượng tử và nhiệt động lực học thống kê cổ điển. Hai khung lý thuyết nền tảng được áp dụng là Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) nhằm tối ưu hóa hình học ô mạng tinh thể và Phương pháp mô phỏng Monte Carlo chính tắc lớn (Grand Canonical Monte Carlo - GCMC) để mô tả trạng thái cân bằng hấp phụ pha khí - rắn.

Nghiên cứu tập trung phân tích 5 khái niệm cốt lõi:

  1. Đơn vị xây dựng cấu trúc thứ cấp (Secondary Building Units - SBUs): Cụm kim loại dạng bánh xe chèo (paddle-wheel) liên kết phối trí với các nguyên tử oxy và nitơ.
  2. Phối tử cầu nối hữu cơ hỗn hợp: Sự kết hợp giữa axit 1,4-benzenedicarboxylic (BDC) tạo mặt phẳng hai chiều và triethylenediamine (TED) đóng vai trò thanh chống không gian ba chiều.
  3. Diện tích bề mặt riêng (Specific Surface Area - SSA): Đo lường không gian tiếp xúc tính toán theo mô hình Brunauer-Emmett-Teller (BET).
  4. Thể tích lỗ xốp rỗng (Pore Volume - Vp): Không gian mao quản nano tự do bên trong cấu trúc tinh thể có kích thước dao động từ 0,8 nm đến 1,5 nm.
  5. Trường lực liên kết phi hóa trị: Mô tả thế năng tương tác Van der Waals thông qua hàm thế Lennard-Jones (LJ 12-6) kết hợp thế tĩnh điện Coulomb dựa trên phân bố điện tích nguyên tử riêng phần.

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

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã chọn lọc 5 hệ kim loại đại diện (Mg, V, Co, Ni, Cu) trong chuỗi cấu trúc M(BDC)(TED)0.5 để khảo sát tương tác với 2 phân tử khí mục tiêu là CO2 và CO, tạo thành tổng cộng 10 hệ mô phỏng nhiệt động độc lập. Mẫu mô phỏng không gian được xây dựng bằng cách nhân bản ô mạng cơ sở thành siêu ô mạng 3x3x3 (supercell) với hàng nghìn nguyên tử nhằm loại bỏ hoàn toàn sai số do hiệu ứng biên tuần hoàn.

Lý do lựa chọn phương pháp mô phỏng GCMC là vì đây là công cụ lý thuyết duy nhất cho phép duy trì cố định thế hóa học, nhiệt độ và thể tích trong khi số lượng phân tử khí hấp phụ được tự do dao động cho đến khi đạt trạng thái cân bằng. Phương pháp này vượt trội hơn mô phỏng động lực học phân tử (MD) ở khả năng tiết kiệm tài nguyên tính toán khi xử lý hệ thống vi mao quản kích thước lớn. Nguồn dữ liệu cấu trúc tinh thể ban đầu được trích xuất và chuẩn hóa từ Ngân hàng Dữ liệu Cấu trúc Cambridge (CSD).

Phân bố điện tích hạt nhân từng nguyên tử được tính toán chính xác bằng phương pháp DDEC (Density Derived Electrostatic and Chemical). Quy trình tính toán sử dụng các gói phần mềm chuyên dụng tiêu chuẩn quốc tế gồm VASP, RASPA và trực quan hóa dữ liệu qua VESTA, Avogadro, ParaView. Toàn bộ quá trình sàng lọc và khảo sát bán kính cắt tương tác thế năng được thực hiện nghiêm ngặt tại các mốc áp suất kiểm chuẩn 1 bar, 5 bar và 10 bar trong giai đoạn nghiên cứu từ năm 2021 đến năm 2022.

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

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

Quá trình mô phỏng và tính toán đã mang lại 4 phát hiện quan trọng có giá trị khoa học cao:

Thứ nhất, khả năng bắt giữ CO2 của chuỗi vật liệu thể hiện sự phụ thuộc sâu sắc vào bản chất ion kim loại trung tâm M. Cấu trúc Mg(BDC)(TED)0.5 thể hiện dung lượng hấp phụ vượt trội nhất, đạt giá trị cực đại lên tới khoảng 22,5 mmol/g ở áp suất 100 bar và nhiệt độ 298 K. Con số này cao hơn xấp xỉ 28% so với biến thể chứa Cu(BDC)(TED)0.5 và cao hơn khoảng 35% so với cấu trúc chứa Ni(BDC)(TED)0.5 trong cùng điều kiện khảo sát.

Thứ hai, đối với khí độc CO, lượng hấp phụ toàn phần bão hòa của chuỗi vật liệu đạt từ 5,8 mmol/g đến 8,2 mmol/g tại 100 bar. Cấu trúc Mg(BDC)(TED)0.5 tiếp tục dẫn đầu với dung lượng đạt khoảng 8,15 mmol/g, cao hơn gấp 7 lần so với vật liệu MIL-101(Cr) trong các công bố thực nghiệm trước đây (chỉ đạt 1,13 mmol/g ở 288 K và 1,13 bar).

Thứ ba, đặc tính diện tích bề mặt riêng tính toán của chuỗi M(BDC)(TED)0.5 dao động trong khoảng từ 1750 m2/g đến 1980 m2/g và thể tích lỗ rỗng đạt từ 0,65 cm3/g đến 0,82 cm3/g. Dữ liệu cho thấy mối tương quan tuyến tính rất chặt chẽ giữa dung lượng bắt giữ khí bão hòa với diện tích bề mặt và thể tích mao quản của tinh thể.

Thứ tư, nhiệt hấp phụ đẳng lượng (Qst) của CO2 và CO trên toàn chuỗi vật liệu nằm trong khoảng từ 22 kJ/mol đến 38 kJ/mol. Mức năng lượng này chứng minh quá trình bắt giữ khí diễn ra hoàn toàn theo cơ chế hấp phụ vật lý, giúp quá trình giải hấp và tái sinh vật liệu diễn ra dễ dàng.

Thảo luận kết quả

Nguyên nhân tạo nên sự vượt trội của cấu trúc chứa Mg là do bán kính ion Mg2+ nhỏ và đặc tính phân cực ion mạnh của liên kết Mg-O, tạo ra gradient điện trường cục bộ lớn tại các vị trí kim loại mở. Điện trường này tương tác mạnh mẽ với moment tứ cực của phân tử CO2. Khi so sánh với các vật liệu xốp nổi tiếng như MOF-177 (đạt 33,5 mmol/g ở 35 bar) hay Mg-MOF-74 (đạt 35,2 wt% ở 1 bar), chuỗi M(BDC)(TED)0.5 thể hiện độ ổn định cấu trúc cơ học vượt trội nhờ các trụ đỡ TED liên kết chéo vững chắc.

Toàn bộ dữ liệu được trực quan hóa sinh động qua hệ thống biểu đồ đường đẳng nhiệt hấp phụ thể hiện tương quan giữa dung lượng bắt giữ theo áp suất, kết hợp bảng ma trận đối sánh diện tích bề mặt riêng SSA và thể tích lỗ rỗng Vp. Các hình ảnh chụp nhanh mật độ phân tử (snapshots) chỉ ra rõ ràng: ở áp suất thấp khoảng 1 bar, các phân tử khí ưu tiên tập trung tại các vị trí gần nút kim loại mở; khi áp suất tăng lên đến 100 bar, các phân tử khí khuếch tán đồng đều và lấp đầy toàn bộ thể tích mao quản xốp. Kết quả này khẳng định tiềm năng ứng dụng thực tế to lớn của dòng vật liệu trong công nghệ phân tách khí thải công nghiệp.

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

Dựa trên các phát hiện khoa học từ mô phỏng tính toán, bốn giải pháp và khuyến nghị hành động cụ thể được đề xuất nhằm thúc đẩy ứng dụng chuỗi vật liệu vào thực tiễn:

Thứ nhất, triển khai tổng hợp thực nghiệm chuỗi vật liệu M(BDC)(TED)0.5, ưu tiên biến thể Mg và Cu bằng phương pháp cơ - hóa học xanh (mechanochemical). Các viện nghiên cứu và phòng thí nghiệm vật liệu chuyên sâu cần thực hiện trong giai đoạn 2023-2024, đặt mục tiêu giảm 70% lượng dung môi hữu cơ tiêu thụ và rút ngắn thời gian tổng hợp từ 24 giờ xuống dưới 90 phút so với phương pháp nhiệt dung môi truyền thống.

Thứ hai, nghiên cứu biến tính và chế tạo cấu trúc khung kim loại lưỡng kim (bimetallic) như Mg0.5Cu0.5(BDC)(TED)0.5 nhằm tối ưu hóa tính chất mao quản. Các nhóm nghiên cứu hóa học vật liệu tại các trường đại học chủ trì triển khai trong lộ trình 12 tháng (từ đầu năm 2024 đến đầu năm 2025), hướng tới mục tiêu gia tăng diện tích bề mặt riêng vượt mốc 2100 m2/g và nâng cao dung lượng hấp phụ CO2 thêm ít nhất 20% tại áp suất khí quyển 1 bar.

Thứ ba, thiết kế và thử nghiệm hệ module màng lọc composite MOF/polymer quy mô pilot dùng cho tháp xử lý khí thải tại các nhà máy luyện kim và nhiệt điện. Doanh nghiệp công nghệ môi trường phối hợp cùng cơ quan quản lý năng lượng triển khai trong giai đoạn 2025-2026, hướng đến chỉ tiêu kỹ thuật thu giữ hơn 85% lượng khí độc hại CO và CO2 trong dòng khí thải liên tục.

Thứ tư, thiết lập cơ sở dữ liệu số hóa mở về thông số trường lực Lennard-Jones và bộ điện tích DDEC cho hơn 50 dẫn xuất vật liệu MOFs dạng paddle-wheel. Các trung tâm tính toán hiệu năng cao đảm trách hoàn thành trong giai đoạn 2024-2025, giúp cộng đồng khoa học rút ngắn 50% thời gian sàng lọc vật liệu trước khi tiến hành thực nghiệm chế tạo.

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

Luận văn là nguồn tài liệu khoa học chuyên sâu và hữu ích cho 4 nhóm đối tượng chính:

Thứ nhất, các nhà nghiên cứu, nghiên cứu sinh và học viên cao học chuyên ngành Hóa lý thuyết, Hóa tính toán và Khoa học vật liệu: Tài liệu cung cấp quy trình chuẩn mực về phương pháp mô phỏng GCMC, kỹ thuật gán điện tích DDEC và phương pháp tối ưu hóa cấu trúc tinh thể với độ tin cậy đạt trên 95%, hỗ trợ đắc lực cho các công trình nghiên cứu và công bố quốc tế.

Thứ hai, kỹ sư công nghệ môi trường và chuyên gia vận hành hệ thống xử lý khí thải: Luận văn cung cấp bảng dữ liệu nhiệt động chi tiết về dung lượng hấp phụ CO và CO2 trong dải áp suất rộng từ 1 đến 100 bar, làm cơ sở kỹ thuật để thiết kế tháp hấp phụ rắn giúp tiết kiệm từ 30% đến 45% chi phí năng lượng vận hành so với công nghệ hấp thụ bằng amine lỏng.

Thứ ba, doanh nghiệp sản xuất và phát triển vật liệu hấp phụ công nghiệp: Tài liệu giúp các nhà quản lý R&D nắm bắt xu hướng công nghệ vật liệu nano tiên tiến trong thị trường trị giá hơn 11,1 tỷ USD, từ đó định hướng phát triển các dòng sản phẩm MOFs thương mại thay thế zeolite và than hoạt tính truyền thống.

Thứ tư, giảng viên đại học và các cơ sở đào tạo sau đại học: Luận văn đóng vai trò là tài liệu tham khảo giảng dạy chuyên đề cho các học phần Mô phỏng phân tử, Nhiệt động lực học bề mặt và Công nghệ thu giữ carbon (CCS).

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

Vật liệu M(BDC)(TED)0.5 có ưu điểm gì vượt trội so với than hoạt tính và zeolite?

Vật liệu M(BDC)(TED)0.5 sở hữu diện tích bề mặt riêng vượt trội đạt từ 1750 m2/g đến 1980 m2/g, cao gấp 2 đến 4 lần so với than hoạt tính tiêu chuẩn (thường dưới 900 m2/g). Cấu trúc vi mao quản đồng nhất cùng các nút kim loại mở tạo lực hút tĩnh điện chọn lọc mạnh mẽ với CO2 và CO ngay cả trong điều kiện độ ẩm cao.

Tại sao phương pháp mô phỏng GCMC lại được ưu tiên lựa chọn thay vì thực nghiệm thuần túy?

Mô phỏng GCMC cho phép khảo sát chính xác quá trình bắt giữ khí ở cấp độ phân tử trong dải áp suất cao từ 1 đến 100 bar mà thực nghiệm khó thực hiện an toàn. Giải pháp này định lượng nhanh dung lượng hấp phụ và nhiệt hấp phụ, giúp tiết kiệm hơn 50% chi phí hóa chất và rút ngắn thời gian nghiên cứu.

Trong chuỗi kim loại được khảo sát, kim loại nào cho hiệu quả bắt giữ CO2 tối ưu nhất?

Kim loại Magnesium (Mg) trong cấu trúc Mg(BDC)(TED)0.5 thể hiện hiệu suất bắt giữ CO2 cao nhất với dung lượng đạt khoảng 22,5 mmol/g ở 100 bar và nhiệt độ 298 K. Nguyên nhân do bán kính ion nhỏ và mật độ điện tích lớn của Mg2+ tạo ra tương tác tĩnh điện cực mạnh với phân tử CO2.

Quá trình tái sinh vật liệu MOF sau khi hấp phụ khí diễn ra như thế nào?

Nhờ nhiệt hấp phụ đẳng lượng Qst chỉ dao động trong khoảng 22 đến 38 kJ/mol (cơ chế hấp phụ vật lý thuần túy), quá trình giải phóng khí diễn ra hoàn toàn thuận nghịch. Người vận hành chỉ cần giảm áp suất nhẹ hoặc gia nhiệt ở mức 50 - 80°C là có thể giải hấp 100% lượng khí mà không làm suy giảm cấu trúc tinh thể.

Khả năng ứng dụng thực tiễn của công trình này trong xử lý khí thải công nghiệp ra sao?

Kết quả nghiên cứu cung cấp thông số nhiệt động chính xác để chế tạo màng lọc composite và tháp hấp phụ khí thải tại các nhà máy xi măng, nhiệt điện than và luyện kim. Trong thực tế, công nghệ này cho phép xử lý chọn lọc nồng độ CO và CO2 cao, đáp ứng nghiêm ngặt các quy chuẩn môi trường quốc gia.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán tính toán và dự đoán khả năng bắt giữ khí thải nhà kính của chuỗi vật liệu khung hữu cơ kim loại M(BDC)(TED)0.5. Năm đóng góp khoa học cốt lõi của công trình bao gồm:

  • Thiết lập thành công quy trình mô phỏng Monte Carlo chính tắc lớn (GCMC) chuẩn xác cho 10 hệ hấp phụ khí CO2 và CO trên chuỗi 5 kim loại Mg, V, Co, Ni, Cu.
  • Xác định Mg(BDC)(TED)0.5 là cấu trúc ưu việt nhất với khả năng bắt giữ CO2 đạt khoảng 22,5 mmol/g và CO đạt khoảng 8,15 mmol/g tại 100 bar ở nhiệt độ 298 K.
  • Khẳng định bản chất bắt giữ khí thuộc cơ chế hấp phụ vật lý thuận nghịch với nhiệt hấp phụ Qst ổn định trong khoảng 22 - 38 kJ/mol, giúp tối ưu hóa năng lượng tái sinh.
  • Chứng minh mối tương quan tuyến tính chặt chẽ giữa diện tích bề mặt riêng (1750 - 1980 m2/g) và thể tích lỗ rỗng với dung lượng lưu trữ khí thải bão hòa.
  • Xây dựng bộ thông số trường lực và phân bố điện tích DDEC đáng tin cậy phục vụ nghiên cứu và phát triển vật liệu xốp nano.

Lộ trình tiếp theo trong giai đoạn 2024-2026 tập trung vào việc chế tạo thực nghiệm màng composite pilot và kiểm thử trên dòng khí thải thực tế đa cấu tử. Quý độc giả, các viện nghiên cứu và doanh nghiệp quan tâm hãy kết nối và hợp tác cùng nhóm tác giả để ứng dụng ngay giải pháp vật liệu tiên tiến này vào các dự án thu giữ carbon và bảo vệ môi trường.