Tổng quan nghiên cứu

Sự gia tăng không ngừng của nồng độ khí nhà kính trong khí quyển đang là một trong những thách thức môi trường nghiêm trọng nhất của nhân loại. Dữ liệu quan trắ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ỹ ghi nhận nồng độ carbon dioxide trong khí quyển đã tăng liên tục từ mức 280 ppm vào đầu thế kỷ XIX lên 408,65 ppm vào tháng 7/2019 và chạm đỉnh 419 ppm vào tháng 5/2021. Ủy ban Liên chính phủ về Biến đổi Khí hậu cảnh báo nồng độ này có thể đạt tới 950 ppm vào năm 2100 nếu các hoạt động phát thải không được kiểm soát. Bên cạnh đó, carbon monoxide cũng là khí thải độc hại nguy hiểm phát sinh từ quá trình đốt cháy không hoàn toàn nhiên liệu hóa thạch, gây tổn thương hệ thần kinh nghiêm trọng khi hít phải.

Các công nghệ thu giữ truyền thống như hấp thụ hóa học bằng dung dịch amine lỏng (monoethanolamine, diethanolamine) có thể loại bỏ tới 98% carbon dioxide nhưng bộc lộ nhiều hạn chế 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 nguy hại. Trước bối cảnh thị trường vật liệu hấp phụ toàn cầu đạt khoảng 11,1 tỷ USD vào năm 2019 (tăng 38% so với năm 2013), việc phát triển vật liệu rắn xốp nano có hiệu suất cao trở thành hướng đi tất yếu.

Nghiên cứu này tập trung tính toán định lượng khả năng bắt giữ hai khí thải carbon dioxide và carbon monoxide của chuỗi khung hữu cơ kim loại M(BDC)(TED)0.5 (với M là các kim loại chuyển tiếp và kiềm thổ gồm Magnesium, Vanadium, Nickel, Cobalt và Copper). Nghiên cứu được 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 trong giai đoạn 2021–2022 ở dải áp suất từ 0 đến 100 bar tại nhiệt độ phòng 298 K. Kết quả đề tài mang ý nghĩa khoa học sâu sắc, làm sáng tỏ cơ chế tương tác ở cấp độ nguyên tử và định hướng tối ưu hóa việc lựa chọn tâm kim loại để phát triển vật liệu xử lý khí thải hiệu quả.

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 thống kê và hóa học lượng tử thông qua hai trụ cột lý thuyết chính:

  1. Lý thuyết mô phỏng Monte Carlo chính tắc lớn: Mô hình toán học ngẫu nhiên xử lý hệ thống vi mô trong điều kiện thế hóa học, thể tích và nhiệt độ không đổi, cho phép số lượng hạt dao động tự do nhằm mô tả chân thực quá trình hấp phụ khí trong vật liệu xốp.
  2. Lý thuyết phiếm hàm mật độ kết hợp phương pháp xác định điện tích riêng phần dựa trên mật độ hạt: Công cụ hóa học tính toán chính xác để tối ưu hóa hình học ô mạng tinh thể và phân bổ trường tĩnh điện nguyên tử.
  3. Mô hình thế tương tác Lennard-Jones 12-6 và lực tĩnh điện Coulomb: Mô tả chính xác thế năng Van der Waals và tương tác điện tích giữa phân tử khí với thành mao quản vật liệu.

Các khái niệm chuyên ngành cốt lõi được áp dụng xuyên suốt bao gồm: Khung hữu cơ kim loại (vật liệu tinh thể nano xốp tạo bởi các ion kim loại và cầu nối hữu cơ), Đơn vị xây dựng cấu trúc thứ cấp (khối đa diện phối trí kim loại - oxy), Diện tích bề mặt riêng tính theo phương pháp Brunauer-Emmett-Teller, Thể tích lỗ rỗng, và Nhiệt hấp phụ đẳng tích.

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

Nguồn dữ liệu của đề tài được thiết lập từ 5 mô hình ô mạng cơ sở của chuỗi M(BDC)(TED)0.5 tương ứng với 5 kim loại Mg, V, Ni, Co, Cu. Cỡ mẫu mô phỏng được mở rộng qua phép lặp không gian ba chiều 3x3x3 ô mạng cơ sở, tạo ra một siêu ô mạng chứa hàng nghìn tâm tương tác nguyên tử. Phương pháp chọn mẫu áp dụng thuật toán lấy mẫu xác suất Metropolis với hơn 100.000 chu kỳ chạy thử nghiệm cân bằng và lấy mẫu.

Lý do lựa chọn phương pháp Monte Carlo chính tắc lớn thay vì mô phỏng động học phân tử là nhờ khả năng cho phép số lượng phân tử khí thay đổi linh hoạt theo ba bước chuyển trạng thái (chèn hạt, xóa hạt, dịch chuyển/quay hạt) mà không phải giải các phương trình đạo hàm theo thời gian, giúp tiết kiệm hơn 60% chi phí tính toán đối với hệ mao quản lớn. Toàn bộ tính toán mô phỏng được thực hiện trên hệ thống phần mềm tính toán chuyên dụng RASPA và VASP, kết hợp các công cụ trực quan hóa cấu trúc VESTA, Avogadro và GNUPLOT trong khoảng thời gian từ tháng 9/2021 đến tháng 8/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 đã mang lại 4 phát hiện quan trọng về khả năng bắt giữ khí thải của chuỗi vật liệu:

  1. Khả năng hấp phụ carbon dioxide vượt trội so với carbon monoxide: Trên toàn bộ 5 cấu trúc khảo sát ở nhiệt độ 298 K, dung lượng hấp phụ carbon dioxide cao hơn từ 250% đến 400% so với carbon monoxide tại dải áp suất cao 50–100 bar. Nguyên nhân xuất phát từ momen tứ cực mạnh của phân tử carbon dioxide so với momen lưỡng cực yếu của carbon monoxide.
  2. Ảnh hưởng rõ rệt của tâm kim loại phối trí: Biến thể Mg(BDC)(TED)0.5 và Ni(BDC)(TED)0.5 thể hiện dung lượng bắt giữ carbon dioxide tại vùng áp suất thấp 1–5 bar cao hơn khoảng 30% so với các biến thể chứa Đồng (Cu) và Cobalt (Co), nhờ ái lực liên kết ion mạnh của tâm Magnesium và Nickel đối với các nguyên tử oxy của khí thải.
  3. Bản chất hấp phụ vật lý thuận nghịch: Nhiệt hấp phụ đẳng tích của carbon dioxide dao động trong khoảng từ 25 kJ/mol đến 45 kJ/mol, nằm hoàn toàn trong thang năng lượng hấp phụ vật lý. Khi dung lượng hấp phụ bề mặt tăng lên, nhiệt hấp phụ có xu hướng giảm dần do các vị trí tương tác mạnh nhất tại nút kim loại và cầu nối hữu cơ được lấp đầy trước.
  4. Tương quan tuyến tính giữa dung lượng hấp phụ bão hòa với đặc trưng rỗng: Tại áp suất bão hòa 100 bar, dung lượng hấp phụ toàn phần tương quan thuận chặt chẽ với diện tích bề mặt riêng (dao động từ 1500 m2/g đến 1952 m2/g) và thể tích lỗ rỗng của cấu trúc khung.

Thảo luận kết quả

Cơ chế hấp phụ được giải thích thông qua cấu trúc hình học paddle-wheel đặc trưng của chuỗi M(BDC)(TED)0.5. Sự kết hợp giữa phối tử benzen-dicarboxylate và triethylenediamine tạo nên mạng lưới lỗ xốp 3D đồng nhất. Khi phân tử khí khuếch tán vào mao quản, trường tĩnh điện do phương pháp DDEC tính toán cho thấy lực hút tập trung mạnh nhất tại lân cận các nút kim loại mở và các vòng thơm liên kết.

Khi so sánh với các vật liệu xốp khác, dung lượng bắt giữ khí của chuỗi M(BDC)(TED)0.5 tuy thấp hơn các kỷ lục lý thuyết như DUT-60 (diện tích bề mặt 7839 m2/g) hay MOF-210 (diện tích bề mặt 10450 m2/g, hấp phụ 74,2 wt% carbon dioxide ở 50 bar), nhưng lại thể hiện độ ổn định cấu trúc và độ bền cơ hóa học vượt trội hơn các dòng zeolite hoặc than hoạt tính thông thường.

Dữ liệu mô phỏng được trình bày rõ nét thông qua biểu đồ đường đẳng nhiệt hấp phụ thể hiện tương quan giữa áp suất từ 0 đến 100 bar với dung lượng hấp phụ bề mặt và hấp phụ toàn phần, cùng bảng tổng hợp các thông số tối ưu hóa ô mạng, diện tích bề mặt riêng và thể tích lỗ rỗng. Không gian phân bố khí được mô phỏng trực quan bằng ảnh chụp 3D mật độ phân tử trong ô mạng tinh thể.

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

Dựa trên các kết quả tính toán mô phỏng, đề tài đưa ra 4 nhóm giải pháp mang tính ứng dụng thực tiễn cao:

  1. Tối ưu hóa quy trình tổng hợp thực nghiệm bằng phương pháp cơ - hóa học: Các nhóm nghiên cứu thực nghiệm cần ứng dụng kỹ thuật nghiền bi tự động không dùng dung môi lỏng để tổng hợp chuỗi tinh thể M(BDC)(TED)0.5 trong thời gian từ 30 đến 120 phút, hướng đến mục tiêu đạt diện tích bề mặt riêng trên 1800 m2/g và giảm thiểu 50% chi phí hóa chất độc hại trong 6–12 tháng tới.
  2. Thiết kế và chế tạo hệ vật liệu khung kim loại lưỡng kim: Phòng thí nghiệm vật liệu nano cần triển khai tổng hợp hệ vật liệu lai như Cu-Ni(BDC)(TED)0.5 trong khung thời gian 12–18 tháng, nhằm nâng cao dung lượng lưu trữ và tách khí thêm 20% đến 25% so với cấu trúc đơn kim truyền thống.
  3. Phát triển màng hỗn hợp ma trận polymer tích hợp tinh thể MOF quy mô pilot: Kỹ sư công nghệ hóa chất và các trung tâm nghiên cứu môi trường cần phối hợp thử nghiệm chế tạo màng lọc hỗn hợp trong 18–24 tháng, đặt mục tiêu đạt hiệu suất tách carbon dioxide/nitơ trên 85% và cắt giảm 35% năng lượng tái sinh so với tháp hấp thụ amine.
  4. Mở rộng cơ sở dữ liệu mô phỏng sàng lọc thông lượng cao: Các nhà khoa học dữ liệu vật liệu tính toán cần tiếp tục mở rộng mô hình GCMC cho hơn 50 cấu trúc MOF biến thể với các hỗn hợp khí đa cấu tử phức tạp chứa hơi nước và lưu huỳnh dioxit trong vòng 6 tháng tới.

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

Nội dung và phương pháp trong công trình này mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Hóa lý thuyết, Vật lý tính toán và Khoa học vật liệu: Tài liệu cung cấp quy trình thiết lập trường lực, phương pháp tính điện tích DDEC và thuật toán mô phỏng Monte Carlo chính tắc lớn hoàn chỉnh để ứng dụng vào nghiên cứu tương tác phân tử.
  2. Kỹ sư công nghệ môi trường và chuyên gia xử lý khí thải công nghiệp: Công trình cung cấp các thông số định lượng chuẩn xác về dung lượng bắt giữ carbon dioxide và carbon monoxide ở các dải áp suất khác nhau, hỗ trợ tính toán thiết kế tháp hấp phụ xử lý khí thải.
  3. Doanh nghiệp sản xuất vật liệu tiên tiến và các dự án bẫy carbon: Luận văn là cơ sở khoa học tin cậy giúp các đơn vị R&D đánh giá hiệu quả kinh tế - kỹ thuật của từng loại ion kim loại trước khi đầu tư kinh phí tổng hợp quy mô lớn.
  4. Giảng viên và nhà nghiên cứu tại các trường đại học: Công trình là tài liệu tham khảo giá trị cho các bài giảng chuyên đề về vật liệu nano xốp, hóa học phối trí và ứng dụng công nghệ thông tin trong mô phỏng vật liệu.

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

Phương pháp mô phỏng Monte Carlo chính tắc lớn có ưu thế gì so với thực nghiệm truyền thống? Mô phỏng GCMC cho phép khảo sát chính xác quá trình hấp phụ ở cấp độ phân tử trong dải áp suất rộng từ 0 đến 100 bar mà không tốn kém hóa chất thí nghiệm. Phương pháp giúp tiết kiệm hơn 70% thời gian nghiên cứu và sàng lọc cấu trúc tối ưu trước khi chế tạo thực tế.

Tại sao chuỗi vật liệu M(BDC)(TED)0.5 lại được lựa chọn để nghiên cứu bắt giữ khí thải? Chuỗi vật liệu này sở hữu cấu trúc mạng xốp 3 chiều vững chắc với các đơn vị xây dựng SBU dạng paddle-wheel, diện tích bề mặt riêng đạt tới 1952 m2/g và độ ổn định nhiệt cao, tạo điều kiện lý tưởng cho việc khuếch tán và lưu trữ các phân tử khí nhỏ.

Kim loại nào trong chuỗi khảo sát thể hiện hiệu năng hấp phụ carbon dioxide triển vọng nhất? Magnesium và Nickel là hai kim loại thể hiện hiệu năng bắt giữ carbon dioxide vượt trội ở vùng áp suất thấp nhờ mật độ điện tích cục bộ cao, giúp tạo lực hút tĩnh điện mạnh với phân tử khí mà vẫn đảm bảo tính thuận nghịch hoàn toàn.

Vật liệu khung hữu cơ kim loại này có khả năng tái sinh và sử dụng lại nhiều lần không? Vật liệu có khả năng tái sinh tuyệt vời vì nhiệt hấp phụ nằm trong khoảng 25–45 kJ/mol thuộc dạng hấp phụ vật lý. Quá trình giải hấp diễn ra nhanh chóng khi giảm áp suất về 1 bar ở nhiệt độ phòng mà không gây phá hủy mạng tinh thể.

Làm thế nào để áp dụng kết quả tính toán mô phỏng vào dây chuyền lọc khí công nghiệp? Các kỹ sư có thể kết hợp hạt nano MOF vào chất nền polymer để tạo màng lọc composite hoặc nén thành các viên xốp định hình đặt trong tháp hấp phụ dao động áp suất, giúp xử lý trực tiếp dòng khí thải ở quy mô nhà máy.

Kết luận

  • Công trình đã thiết lập thành công mô hình tính toán định lượng khả năng bắt giữ carbon dioxide và carbon monoxide của chuỗi 5 vật liệu M(BDC)(TED)0.5 bằng phương pháp mô phỏng Monte Carlo chính tắc lớn kết hợp lý thuyết phiếm hàm mật độ.
  • Nghiên cứu xác định dung lượng hấp phụ carbon dioxide vượt trội hơn carbon monoxide từ 250% đến 400% ở áp suất cao, đồng thời chứng minh biến thể chứa Magnesium và Nickel có lực bắt giữ khí mạnh nhất ở áp suất thấp.
  • Giá trị nhiệt hấp phụ đẳng tích từ 25 kJ/mol đến 45 kJ/mol khẳng định cơ chế hấp phụ vật lý bền vững, cho phép tái sinh chất hấp phụ với mức tiêu hao năng lượng thấp.
  • Đề xuất lộ trình 6–24 tháng tiếp theo tập trung vào việc tổng hợp xanh bằng phương pháp cơ - hóa học, chế tạo màng lọc hỗn hợp và thử nghiệm trong môi trường khí thải thực tế.
  • Các nhà nghiên cứu và doanh nghiệp công nghệ môi trường hãy tận dụng ngay cơ sở dữ liệu tính toán từ luận văn để tối ưu hóa việc lựa chọn vật liệu và thúc đẩy ứng dụng công nghệ bẫy carbon hiện đại.