Tổng quan nghiên cứu

Vật liệu khung hữu cơ kim loại đóng vai trò quan trọng trong việc phát triển các chất xúc tác quang hóa thế hệ mới nhờ diện tích bề mặt lớn và cấu trúc lỗ xốp có thể tinh chỉnh linh hoạt. Trong số các hệ vật liệu tiên tiến, cấu trúc khung kim loại titan nổi bật với độ bền nhiệt lên tới 200 độ C và khả năng làm việc ổn định trong môi trường axit. Tinh thể MOF-901 được tổng hợp từ cụm tâm kim loại titan dạng sáu đỉnh kết hợp cùng cầu nối hữu cơ 4-aminobenzoate, sở hữu ô đơn vị lục giác chứa 336 nguyên tử với hằng số mạng a bằng 27.40 Angstrom và c bằng 8.84 Angstrom. Hệ vật liệu này có hai hệ thống lỗ xốp riêng biệt gồm lỗ xốp lớn hình tròn đường kính khoảng 12.5 Angstrom và lỗ xốp nhỏ hình elip bán kính 3.5 Angstrom, tạo điều kiện thuận lợi cho việc dung nạp các phân tử khách.

Mặc dù thực nghiệm đã chứng minh MOF-901 có độ rộng vùng cấm quang học đạt 2.65 eV và khả năng xúc tác quang vượt trội so với chất xúc tác thương mại trong phản ứng trùng hợp methyl methacrylate dưới ánh sáng khả kiến từ 340 đến 550 nm, cơ chế tương tác quỹ đạo điện tử và độ bền cơ học thực tế dưới sự hiện diện của dung môi hấp phụ vẫn là bài toán chưa có lời giải chi tiết. Đề tài tập trung khảo sát sự thay đổi cấu trúc hình học, cấu trúc điện tử và độ bền cơ học của MOF-901 khi hấp phụ phân tử methanol ở các mức nồng độ 5.40% và 10.80% khối lượng. Công trình được thực hiện tại Trường Đại học Bách Khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn từ tháng 8 năm 2019 đến tháng 2 năm 2020. Kết quả nghiên cứu mang lại ý nghĩa quan trọng trong việc làm sáng tỏ cơ chế điều biến vùng cấm và cung cấp bộ dữ liệu hằng số đàn hồi chính xác, phục vụ định hướng ứng dụng vật liệu trong công nghệ xử lý môi trường và xúc tác quang hóa.

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 cơ học lượng tử nhiều hạt, kết hợp các phép gần đúng vật lý hiện đại nhằm giải bài toán tương tác phức tạp giữa hạt nhân và điện tử trong vật rắn:

  • Phương trình Schrödinger hệ nhiều hạt và phép gần đúng Born-Oppenheimer: Tách biệt chuyển động của hạt nhân có khối lượng lớn khỏi chuyển động nhanh của các điện tử, quy bài toán về trạng thái chuyển động của điện tử trong thế tĩnh điện cố định.
  • Phương pháp trường tự hợp Hartree-Fock và định thức Slater: Mô tả hệ điện tử thông qua định thức bất đối xứng thỏa mãn nguyên lý loại trừ Pauli, làm tiền đề cho việc xấp xỉ tương tác Coulomb và hiệu ứng trao đổi.
  • Lý thuyết phiếm hàm mật độ theo hai định lý nền tảng Hohenberg-Kohn: Xác định năng lượng trạng thái cơ bản của hệ nhiều hạt hoàn toàn phụ thuộc vào mật độ điện tử không gian thay vì hàm sóng nhiều chiều.
  • Hệ phương trình Kohn-Sham và gần đúng gradient suy rộng: Áp dụng phiếm hàm Perdew-Burke-Ernzerhof để mô tả năng lượng trao đổi - tương quan không đồng nhất cục bộ, bổ sung hiệu chỉnh phân tán Grimme DFT-D3 để mô phỏng chính xác lực liên kết Van der Waals tầm xa.
  • Lý thuyết đàn hồi tinh thể học cổ điển: Dựa trên định luật Hooke tổng quát liên hệ giữa tensor ứng suất và tensor biến dạng thông qua ma trận hằng số đàn hồi 6x6, áp dụng các tiêu chuẩn ổn định cơ học Born cho hệ tinh thể đối xứng lục giác.

Các khái niệm then chốt xuyên suốt đề tài bao gồm mật độ trạng thái từng phần, độ rộng vùng cấm trực tiếp và gián tiếp, mức năng lượng Fermi, hằng số đàn hồi dao động mạng và năng lượng hấp phụ phân tử.

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

Nghiên cứu sử dụng gói phần mềm mô phỏng lượng tử VASP chuyên dụng để tính toán cấu trúc vật liệu từ nguyên lý ban đầu:

  • Mô hình tính toán và chọn mẫu: Cỡ mẫu mô phỏng bao gồm ô đơn vị chuẩn của MOF-901 với 336 nguyên tử (12 nguyên tử Ti, 12 nguyên tử N, 144 nguyên tử C, 120 nguyên tử H và 48 nguyên tử O). Mẫu nghiên cứu được khảo sát ở ba trạng thái độc lập gồm cấu trúc tinh khiết không hấp phụ, cấu trúc hấp phụ 6 phân tử methanol và cấu trúc hấp phụ 12 phân tử methanol trong một ô đơn vị.
  • Thiết lập thông số tính toán: Không gian sóng phẳng được biểu diễn bằng bộ giả thế sóng chiếu mở rộng PAW. Năng lượng cắt sóng phẳng được tối ưu hóa chuẩn xác ở mức 400 eV sau khi khảo sát dải năng lượng từ 300 eV đến 500 eV. Lưới tích phân vùng Brillouin được lấy mẫu tự động theo phương pháp Monkhorst-Pack với kích thước 1x1x4.
  • Tiêu chí hội tụ và độ tin cậy: Quá trình trường tự hợp đạt điểm dừng khi biến thiên năng lượng tổng nhỏ hơn 10 mũ trừ 5 eV và độ lớn lực tác dụng lên mỗi nguyên tử nhỏ hơn 10 mũ trừ 3 eV trên Angstrom. Lý do lựa chọn phương pháp này là vì DFT-PBE kết hợp giả thế PAW mang lại sự cân bằng tối ưu giữa độ chính xác cấu trúc hóa học và thời gian xử lý dữ liệu cho hệ tinh thể kích thước lớn trên 300 nguyên tử.
  • Thời gian thực hiện: Toàn bộ quá trình tính toán mô phỏng và phân tích dữ liệu được triển khai liên tục trong thời gian 6 tháng từ ngày 19 tháng 8 năm 2019 đến ngày 1 tháng 2 năm 2020.

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

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

Quá trình tính toán mô phỏng lượng tử đã ghi nhận 4 phát hiện khoa học mang tính đột phá về cấu trúc và tính chất của vật liệu:

  1. Sự co rút cấu trúc hình học theo trục tinh thể: Sau tối ưu hóa bằng phương pháp PBE và PBE-D3, hằng số mạng a của MOF-901 nguyên bản biến thiên rất nhỏ dưới 0.47%, trong khi hằng số mạng c giảm mạnh 7.36% so với thực nghiệm. Thể tích ô đơn vị của MOF-901 giảm 8.08% với PBE và giảm 5.28% với PBE-D3. Khi hấp phụ 6 phân tử methanol, thể tích giảm từ 1.66% đến 3.44%. Tuy nhiên, khi tăng nồng độ hấp phụ lên 12 phân tử methanol, thể tích ô mạng có xu hướng giãn nở nhẹ 4.23% ở PBE-D3 và 0.69% ở PBE do sự tăng hằng số mạng c từ 1.00% đến 1.69%.
  2. Biến dạng cục bộ của các liên kết phối trí: Độ dài liên kết Ti-O methoxy giảm từ 1.95 Angstrom ở thực nghiệm xuống 1.78 Angstrom ở mô phỏng. Đồng thời, góc nghiêng alpha giữa vector Ti-O methoxy và trục c tăng mạnh từ gần 0 độ lên 35 đến 50 độ. Ngược lại, liên kết Ti-O carboxylate bị kéo dài từ 1.98 Angstrom lên 2.11 Angstrom do sức căng cơ học từ các phân tử methanol hấp phụ bên trong lỗ xốp.
  3. Xu hướng thu hẹp độ rộng vùng cấm điện tử: Giá trị vùng cấm tính toán của MOF-901 nguyên bản là 2.07 eV với chuyển dời trực tiếp tại điểm Gamma. Khi hấp phụ 6 phân tử methanol, độ rộng vùng cấm giảm 11.11% xuống còn 1.84 eV. Khi nồng độ đạt 12 phân tử methanol, độ rộng vùng cấm tiếp tục giảm 21.41% xuống mức 1.63 eV. Mức năng lượng Fermi dịch chuyển về phía vùng năng lượng cao hơn, từ mức ban đầu lên âm 1.51 eV ở hệ 6 methanol và âm 1.10 eV ở hệ 12 methanol.
  4. Gia tăng độ bền cơ học theo nồng độ phân tử khách: Dữ liệu hằng số đàn hồi dao động mạng cho thấy giá trị C11 tăng từ 2504 kBar ở MOF-901 lên 2604 kBar ở hệ 6 methanol và đạt 2608 kBar ở hệ 12 methanol. Hằng số nén C33 tăng mạnh từ 921 kBar lên 1031 kBar (tăng 11.9%) và đạt 1152 kBar (tăng 25.1%). Năng lượng hấp phụ phân tử mang giá trị âm đạt âm 0.23 eV trên một phân tử ở phép tính PBE và âm 0.42 eV trên một phân tử khi bổ sung tương tác phân tán PBE-D3, thể hiện mức tăng liên kết lên tới 82.6%.

Thảo luận kết quả

Dữ liệu tính toán trong luận văn được mô hình hóa và trực quan hóa chi tiết thông qua các bảng đối sánh tham số mạng, đồ thị đường cong cấu trúc vùng năng lượng electron và biểu đồ phân bố mật độ trạng thái từng phần PDOS. Khi biểu diễn mật độ trạng thái, các đường cong năng lượng được bóc tách rõ ràng theo 4 thành phần cấu trúc: nguyên tử titan, oxy carboxylate, nhóm methoxy và cầu nối 4-aminobenzoate.

Biểu đồ PDOS chứng minh rằng đỉnh của vùng hóa trị và đáy của vùng dẫn trong MOF-901 chủ yếu được đóng góp bởi các phân lớp orbital 2p của cầu nối hữu cơ 4-aminobenzoate, trong khi phân lớp 3d của tâm kim loại titan chỉ đóng vai trò thứ cấp trong việc định hình biên vùng cấm. Khi các phân tử methanol đi vào lỗ xốp lớn, chúng thiết lập các liên kết phân tán Van der Waals tầm xa với vòng phenyl của phối tử. Tương tác này làm uốn cong khung hữu cơ, thay đổi độ dài liên kết Ti-O carboxylate và gián tiếp làm dịch chuyển mức năng lượng của các orbital 2p, dẫn tới sự thu hẹp độ rộng vùng cấm mà không làm mất đi trật tự tinh thể.

Độ rộng vùng cấm tính toán 2.07 eV có mức sai lệch thấp hơn 21.88% so với giá trị thực nghiệm quang học 2.65 eV. Mức chênh lệch này hoàn toàn tương đồng với các công bố lý thuyết quốc tế khi áp dụng phiếm hàm GGA cho các hệ vật liệu lai hóa tương đương như HKUST-1 hoặc UiO-66 (thường sai lệch từ 20% đến 30% do hạn chế cố hữu của việc xấp xỉ thế trao đổi). Về mặt cơ học, giá trị hằng số đàn hồi tổng C33 mang giá trị âm (âm 155 kBar và âm 157 kBar) không phản ánh sự phá hủy liên kết mà bắt nguồn từ bản chất phân lớp 2D với khoảng cách giữa các lớp là 4.42 Angstrom, nơi dao động nhiệt và tương tác trượt giữa các mặt phẳng chiếm ưu thế.

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

Dựa trên các kết quả tính toán nền tảng, luận văn đưa ra 4 khuyến nghị hành động cụ thể nhằm thúc đẩy khả năng ứng dụng thực tế của vật liệu:

  • Thực hiện tính toán nâng cao với phiếm hàm lai hóa HSE06 và gần đúng GW: Nhóm nghiên cứu vật lý tính toán cần triển khai các phép tính năng lượng kích thích để giảm sai số vùng cấm từ mức 21.88% hiện tại xuống dưới 5.0%, hoàn thành trong lộ trình 6 tháng.
  • Ứng dụng cơ chế hấp phụ dung môi để tối ưu hóa hiệu suất quang xúc tác: Các kỹ sư vật liệu tại phòng thí nghiệm cần điều chỉnh nồng độ pha tạp dung môi methanol ở mức 5.40% đến 10.80% khối lượng nhằm tinh chỉnh vùng cấm về dải 1.63 - 1.84 eV, mục tiêu nâng cao hiệu suất phản ứng trùng hợp quang hóa thêm 15% đến 20% trong vòng 9 tháng.
  • Khảo sát độ bền cơ học nhiệt động lực học phân tử AIMD: Các nhà khoa học chuyên ngành mô phỏng cần đánh giá độ ổn định của khung mạng tinh thể dưới áp suất nén thủy tĩnh từ 1 kBar đến 50 kBar và nhiệt độ làm việc từ 300 K đến 500 K trong thời gian 12 tháng.
  • Mở rộng mô hình hấp phụ đa cấu tử đối với các hợp chất khí độc hại: Trung tâm nghiên cứu vật liệu cấu trúc nano cần kiểm tra khả năng hấp phụ chọn lọc đối với CO2, CH4 và các hợp chất hữu cơ dễ bay hơi VOCs, hướng tới mục tiêu năng lượng hấp phụ âm từ âm 0.30 eV đến âm 0.60 eV trên một phân tử trong kế hoạch 18 tháng tới.

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

Nội dung luận văn mang lại giá trị học thuật và giải pháp kỹ thuật chuyên sâu cho 4 nhóm đối tượng chính:

  • Học viên cao học và nghiên cứu sinh ngành Vật lý Kỹ thuật, Vật lý Địa cầu và Khoa học Vật liệu: Nắm vững quy trình thiết lập bài toán DFT cho tinh thể kích thước lớn trên 300 nguyên tử, kỹ thuật chọn điểm k trong vùng Brillouin và phương pháp tính toán tensor đàn hồi 6 thành phần.
  • Kỹ sư nghiên cứu và phát triển vật liệu bán dẫn, quang xúc tác: Khai thác phương pháp biến tính vùng cấm phi hóa học thông qua hấp phụ phân tử để phát triển các vật liệu thu hoạch năng lượng mặt trời trong dải ánh sáng khả kiến 340 - 550 nm.
  • Chuyên gia công nghệ hóa học và môi trường: Sử dụng các thông số nhiệt động học về thể tích lỗ xốp 12.5 Angstrom và năng lượng hấp phụ âm 0.42 eV trên một phân tử để thiết kế hệ thống bẫy khí và lọc dung môi hữu cơ hiệu năng cao.
  • Giảng viên đại học chuyên ngành Vật lý Chất rắn và Hóa tính toán: Sử dụng công trình như một tài liệu tham khảo chuẩn mực về phương pháp sử dụng phần mềm VASP, phân tích mật độ trạng thái PDOS và mô hình hóa ứng suất biến dạng trong cấu trúc 2D.

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

  • Tại sao kết quả tính toán độ rộng vùng cấm bằng phương pháp PBE lại thấp hơn 21.88% so với thực nghiệm?
    • Phương pháp gần đúng gradient suy rộng GGA-PBE có xu hướng đánh giá thấp độ rộng vùng cấm của chất rắn do không tính đủ hiệu ứng tự tương tác của điện tử và bỏ qua hiệu ứng kích thích đa hạt. Mức sai lệch từ 2.65 eV thực nghiệm xuống 2.07 eV tính toán hoàn toàn nằm trong giới hạn sai số cho phép từ 20% đến 30% của lý thuyết phiếm hàm mật độ tiêu chuẩn.
  • Cơ chế nào dẫn đến việc thu hẹp độ rộng vùng cấm khi MOF-901 hấp phụ methanol?
    • Khi nồng độ methanol tăng lên 12 phân tử trong ô mạng, các phân tử này tạo tương tác phân tán tầm xa với vòng thơm của cầu nối 4-aminobenzoate. Hiện tượng này làm biến dạng hình học phối tử và kéo dài liên kết Ti-O carboxylate lên 2.11 Angstrom, làm thay đổi mật độ điện tích các orbital 2p và thu hẹp vùng cấm xuống 1.63 eV mà không cần sự xen phủ điện tử trực tiếp từ methanol.
  • Ý nghĩa vật lý của giá trị hằng số đàn hồi tổng C33 mang giá trị âm là gì?
    • Giá trị C33 tổng mang dấu âm (âm 155 kBar ở hệ 6 methanol) phản ánh đặc tính xếp lớp 2D của MOF-901 với khoảng cách giữa các lớp là 4.42 Angstrom. Liên kết giữa các lớp chủ yếu là tương tác yếu, dẫn tới sự trượt lớp và biến dạng góc liên kết Ti-O khi chịu ứng suất, trong khi hằng số đàn hồi dao động mạng C33 vẫn dương đạt 1031 kBar đảm bảo khung tinh thể ổn định.
  • Năng lượng hấp phụ âm 0.42 eV trên một phân tử thể hiện điều gì về mặt nhiệt động học?
    • Giá trị năng lượng hấp phụ âm chứng minh quá trình hấp phụ phân tử methanol vào khung mạng MOF-901 là quá trình tỏa nhiệt tự phát. Mức năng lượng âm 0.42 eV trên một phân tử nằm trong dải hấp phụ vật lý điển hình, cho phép vật liệu giữ phân tử khách ổn định ở điều kiện thường nhưng dễ dàng giải hấp để tái sinh chất xúc tác mà không phá hủy cấu trúc.
  • Cụm tâm kim loại titan đóng vai trò gì trong cấu trúc điện tử của MOF-901?
    • Cụm Ti-oxo hexameric đóng vai trò liên kết cơ học chính để định hình mạng không gian lục giác vững chắc. Tuy nhiên, phân tích PDOS chỉ ra rằng các orbital 3d của titan không đóng góp vào đỉnh vùng hóa trị và đáy vùng dẫn, cho thấy tính chất quang học của MOF-901 chủ yếu do phối tử hữu cơ kiểm soát.

Kết luận

  • Tối ưu hóa thành công mô hình cấu trúc MOF-901 nguyên bản và hệ hấp phụ methanol với quy mô 336 nguyên tử trong một ô đơn vị bằng lý thuyết DFT-PBE kết hợp hiệu chỉnh phân tán D3.
  • Xác định cơ chế co nén cấu trúc bất đẳng hướng với hằng số mạng c giảm 7.36% do sự nghiêng góc liên kết Ti-O methoxy từ 0 độ lên 35 đến 50 độ.
  • Làm sáng tỏ quy luật điều biến độ rộng vùng cấm giảm tuyến tính từ 2.07 eV xuống 1.84 eV (hệ 6 methanol) và 1.63 eV (hệ 12 methanol) dưới tác động uốn cong phối tử 4-aminobenzoate.
  • Chứng minh độ bền cơ học mạng tinh thể gia tăng từ 4% đến 26% khi nồng độ phân tử hấp phụ tăng dần, thỏa mãn đầy đủ các tiêu chuẩn ổn định Born cho hệ lục giác.
  • Xác lập năng lượng hấp phụ tự phát đạt âm 0.42 eV trên một phân tử, khẳng định tiềm năng ứng dụng linh hoạt của vật liệu trong xúc tác quang hóa.

Đóng góp lớn nhất của công trình là cung cấp cái nhìn toàn diện từ cấu trúc hình học, cơ chế điện tử đến độ bền cơ học của vật liệu MOF-901 dưới ảnh hưởng của phân tử khách. Trong giai đoạn 2020 - 2022, các hướng nghiên cứu tiếp theo cần tập trung mở rộng mô phỏng với phiếm hàm HSE06 và thử nghiệm thực tế quá trình phân hủy quang học các hợp chất hữu cơ. Quý độc giả, các nhà nghiên cứu và kỹ sư quan tâm có thể khai thác toàn bộ bộ tham số tính toán của luận văn để áp dụng vào các dự án phát triển vật liệu xúc tác quang tiên tiến.