Tổng quan nghiên cứu

Ô nhiễm không khí do các hợp chất hữu cơ dễ bay hơi (VOC) đang trở thành một thách thức môi trường và y tế nghiêm trọng trên phạm vi toàn cầu. Theo các thống kê y khoa, trong khí thở ra của con người có chứa hơn 200 hợp chất VOC khác nhau, trong đó nhiều chất như propenal, acetone, ethylamine hay isopentane đóng vai trò là dấu ấn sinh học đặc trưng cho các bệnh lý hiểm nghèo, đặc biệt là ung thư phổi giai đoạn đầu. Việc phát hiện sớm các phân tử khí này ở nồng độ siêu vết thông qua hơi thở giúp gia tăng tỷ lệ điều trị thành công lên trên 70%, đồng thời giảm thiểu đáng kể chi phí điều trị xâm lấn so với các kỹ thuật chụp cắt lớp hay cộng hưởng từ truyền thống.

Nghiên cứu được thực hiện tại Phòng Thí nghiệm Hóa học tính toán và mô phỏng thuộc Trường Đại học Quy Nhơn, hoàn thành vào tháng 10 năm 2023. Đề tài tập trung vào việc khảo sát khả năng hấp phụ của 5 phân tử VOC tiêu biểu trên bề mặt vật liệu hai chiều dichalcogenide kim loại chuyển tiếp MoX2 (với X là lưu huỳnh S và selen Se) ở cả hai trạng thái tinh khiết và pha tạp sắt (Fe). Vật liệu MoS2 dạng đơn lớp sở hữu năng lượng vùng cấm trực tiếp khoảng 1,9 eV (so với 1,3 eV ở dạng khối), mang lại diện tích tiếp xúc riêng lớn và độ nhạy dẫn truyền điện tử vượt trội.

Mục tiêu cốt lõi của công trình là xác định cấu hình hấp phụ bền vững nhất, tính toán các thông số nhiệt động học bề mặt và làm sáng tỏ bản chất các tương tác liên phân tử ở cấp độ lượng tử. Kết quả nghiên cứu không chỉ bổ sung cơ sở dữ liệu học thuật quan trọng cho ngành Hóa lý thuyết và Hóa lý mà còn mở ra định hướng thiết kế các cảm biến khí y sinh có độ nhạy cao, phục vụ chẩn đoán bệnh không xâm lấn trong thực tiễn y khoa.

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 xây dựng trên nền tảng cơ học lượng tử hiện đại, kết hợp đồng bộ bốn trụ cột lý thuyết chính:

  1. Thuyết phiếm hàm mật độ (DFT): Dựa trên hai định lý nền tảng Hohenberg-Kohn và hệ phương trình tự hợp Kohn-Sham, DFT biểu diễn năng lượng toàn phần của hệ nhiều electron dưới dạng một phiếm hàm của mật độ electron. Phương pháp này cho phép giải quyết chính xác cấu trúc điện tử của các hệ vật liệu tuần hoàn phức tạp với chi phí tính toán tối ưu.
  2. Thuyết nguyên tử trong phân tử (AIM) của Richard Bader: Lý thuyết sử dụng ma trận Hessian mật độ electron để xác định các điểm tới hạn liên kết (BCP), điểm tới hạn vòng (RCP) và điểm tới hạn lồng (CCP). Dấu và độ lớn của toán tử Laplace mật độ electron $\nabla^2\rho(r)$, cùng mật độ năng lượng electron khu trú $H(r)$, giúp phân loại chính xác bản chất liên kết là cộng hóa trị, liên kết hydro hay tương tác van der Waals.
  3. Phân tích orbital liên kết thích hợp (NBO): Phương pháp chuyển đổi hàm sóng thành các orbital nguyên tử tự nhiên (NAO) và orbital liên kết tự nhiên (NBO). Năng lượng bền hóa bậc hai $E^{(2)}$ định lượng mức độ xen phủ và truyền điện tích giữa các orbital phân tử bị chiếm (cho electron) và orbital phản liên kết còn trống (nhận electron).
  4. Giản đồ thế tĩnh điện phân tử (MEP): Khảo sát sự phân bố điện tích bề mặt với các giá trị thế năng tĩnh điện đẳng trị khoảng 0,01 au, cho phép dự đoán chính xác các vị trí giàu electron (tâm nucleophile) và nghèo electron (tâm electrophile) trên từng phân tử VOC.

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

Nghiên cứu tiến hành khảo sát trên bộ mẫu gồm 5 phân tử VOC có cấu trúc đa dạng (acetone, 1-propanol, propenal, isopentane, ethylamine) tương tác với 4 mô hình bề mặt vật liệu: MoS2 tinh khiết, MoSe2 tinh khiết, MoS2 pha tạp Fe và MoSe2 pha tạp Fe. Phương pháp chọn mẫu cấu trúc tuân theo nguyên tắc cực tiểu hóa năng lượng tự do bề mặt, thiết lập hơn 20 hướng tiếp cận không gian khác nhau giữa các nhóm chức hoạt động của VOC và các tâm hấp phụ trên vật liệu.

Lý do lựa chọn phương pháp phân tích:

  • Tối ưu hóa hình học hệ tuần hoàn: Sử dụng gói phần mềm mô phỏng VASP với phiếm hàm vdW-DF2 nhằm hiệu chỉnh chính xác tương tác phân tán van der Waals tầm xa, vốn chiếm ưu thế trong các hệ vật liệu hai chiều xếp lớp. Năng lượng cắt sóng phẳng được thiết lập ở mức 400 eV đến 500 eV nhằm đảm bảo độ hội tụ năng lượng tuyệt đối.
  • Tính toán lượng tử phân tử và cụm bề mặt: Sử dụng phần mềm Gaussian 09 kết hợp bộ hàm cơ sở mở rộng LANL2DZ cho các nguyên tử kim loại chuyển tiếp (Mo, Fe) và bộ cơ sở Pople chuẩn cho các nguyên tử phi kim (C, H, O, N, S, Se).
  • Định lượng năng lượng tương tác: Năng lượng hấp phụ ($E_A$) và năng lượng tương tác ($E_I$) được tính toán cùng với năng lượng biến dạng cấu trúc ($E_D$) của cả phân tử và bề mặt, loại trừ hoàn toàn sai số dịch chuyển cơ sở. Toàn bộ chuỗi tính toán và phân tích topo được thực hiện liên tục trong khoảng thời gian từ năm 2022 đến cuối năm 2023.

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 lượng tử đã thu được những kết quả đột phá về khả năng lưu giữ và biến tính điện tử của vật liệu:

  • Hấp phụ trên bề mặt MoX2 tinh khiết mang bản chất vật lý thuần túy: Năng lượng hấp phụ $E_A$ của cả 5 phân tử VOC trên MoS2 và MoSe2 nguyên bản dao động ở mức yếu đến trung bình, trong khoảng từ -0,15 eV đến -0,48 eV. Khoảng cách cân bằng giữa các nguyên tử của phân tử VOC và lớp nguyên tử chalcogen (S hoặc Se) trên bề mặt đều lớn hơn 2,85 Å, chứng minh hệ liên kết chủ yếu được duy trì bởi lực phân tán van der Waals.
  • Pha tạp Fe làm tăng đột biến độ bền hấp phụ: Khi thay thế một nguyên tử Mo trên bề mặt bằng nguyên tử Fe, năng lượng hấp phụ $E_A$ tăng mạnh từ 150% đến trên 320%, đạt giá trị âm sâu từ -0,85 eV đến -1,82 eV. Khoảng cách tương tác giữa tâm dị tố (O, N) của phân tử VOC với nguyên tử Fe bị rút ngắn đáng kể xuống dưới 2,15 Å, chuyển hóa tương tác sang cơ chế hấp phụ hóa học mạnh mẽ.
  • Tính chọn lọc cao đối với các phân tử chứa dị tố: Ethylamine (chứa nhóm amine mang nguyên tử N) và propenal, acetone (chứa nhóm carbonyl mang nguyên tử O) thể hiện năng lượng liên kết cao hơn từ 45% đến 65% so với hydrocarbon mạch nhánh bão hòa như isopentane. Độ phân cực lớn của các liên kết C=O và C-N tạo điều kiện cho quá trình ghép cặp orbital diễn ra thuận lợi.

Thảo luận kết quả

Các phát hiện lượng tử được minh chứng chi tiết thông qua các biểu đồ phân bố mật độ electron và bảng tổng hợp thông số topo học:

  • Cơ chế chuyển dịch điện tích NBO: Khi phân tử VOC tiếp cận bề mặt pha tạp Fe, xảy ra sự dịch chuyển mật độ electron (EDT) rõ rệt từ các orbital cặp electron tự do LP của O hoặc N sang orbital d còn trống của kim loại Fe. Năng lượng ổn định hóa bậc hai $E^{(2)}$ của tương tác cho nhận này đạt giá trị trên 18,5 kcal/mol, giải thích tại sao năng lượng biến dạng $E_D$ của phân tử tăng lên khoảng 0,12 eV đến 0,25 eV do liên kết nội phân tử bị kéo dãn nhẹ.
  • Đặc trưng topo AIM: Tại các điểm tới hạn liên kết BCP hình thành giữa Fe và nguyên tử O/N, giá trị mật độ electron $\rho(r)$ đạt từ 0,045 au đến 0,082 au, trong khi toán tử Laplace $\nabla^2\rho(r) > 0$ kết hợp với tổng mật độ năng lượng $H(r) < 0$. Dữ liệu này khẳng định sự xuất hiện của liên kết phối trí mang một phần tính chất cộng hóa trị phân cực, vượt trội hoàn toàn so với tương tác tĩnh điện thông thường.
  • So sánh với các nghiên cứu tiền nhiệm: Khả năng hấp phụ khí của MoSe2 pha tạp Fe cho thấy hiệu năng tương đương hoặc vượt trội khoảng 15% so với các hệ SnS2 pha tạp Ru hay MoSe2 pha tạp Al từng được công bố trong các báo cáo quốc tế gần đây. Dữ liệu khi trình bày dưới dạng đồ thị tương quan năng lượng - khoảng cách và biểu đồ cột so sánh năng lượng vùng cấm thể hiện rõ tính ưu việt của vật liệu pha tạp Fe trong việc thu bắt các phân tử khí độc.

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

Dựa trên các kết quả tính toán nhiệt động học và cấu trúc điện tử, luận văn đưa ra 4 khuyến nghị thực tiễn nhằm phát triển ứng dụng vật liệu trong công nghệ cảm biến và xử lý môi trường:

  1. Ứng dụng công nghệ màng mỏng chế tạo cảm biến hơi thở y tế: Các viện nghiên cứu công nghệ nano và doanh nghiệp y sinh cần đẩy mạnh tích hợp đơn lớp MoS2 pha tạp Fe vào các vi chip cảm biến khí quang điện trở. Mục tiêu là chế tạo thiết bị phân tích hơi thở cầm tay có khả năng phát hiện các phân tử chỉ điểm ung thư phổi như propenal và acetone ở ngưỡng nồng độ siêu vết dưới 50 ppb, với thời gian đáp ứng tín hiệu nhanh dưới 5 giây trong lộ trình 18 tháng tới.
  2. Tối ưu hóa quy trình tổng hợp vật liệu pha tạp: Các phòng thí nghiệm hóa học vật liệu nên áp dụng phương pháp lắng đọng hơi hóa học (CVD) có kiểm soát nồng độ tiền chất để điều chỉnh tỷ lệ pha tạp Fe trên bề mặt MoX2 ở mức tối ưu từ 1,5% đến 3,0% nguyên tử. Tỷ lệ này đảm bảo tối đa hóa mật độ tâm hoạt tính mà không làm phá vỡ cấu trúc mạng tinh thể 2D hay làm suy giảm độ linh động của hạt mang điện.
  3. Phát triển màng lọc xử lý khí thải công nghiệp: Bộ Tài nguyên và Môi trường cùng các ban quản lý khu công nghiệp cần thí điểm ứng dụng vật liệu composite chứa MoS2 biến tính Fe để thu hồi các dung môi hữu cơ dễ bay hơi độc hại (như ethylamine, acetone) tại các nhà máy sản xuất sơn, in ấn và dược phẩm, hướng tới mục tiêu giảm thiểu trên 85% lượng khí VOC phát thải ra môi trường trước năm 2026.
  4. Mở rộng mô hình sàng lọc tính toán tự động: Các nhóm nghiên cứu Hóa học tính toán nên kết hợp phương pháp DFT với các thuật toán học máy (Machine Learning) để mở rộng cơ sở dữ liệu tính toán cho hơn 30 kim loại chuyển tiếp khác (như Co, Ni, Cu, Ru, Pt), rút ngắn 60% thời gian thử nghiệm thực tế trong việc tìm kiếm các vật liệu xúc tác quang và cảm biến thế hệ mới.

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

Công trình nghiên cứu mang giá trị học thuật chuyên sâu và khả năng ứng dụng thực tế cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Các nhà nghiên cứu Hóa học lượng tử và Khoa học vật liệu tính toán: Cung cấp quy trình chuẩn mực từ việc xây dựng mô hình bề mặt định kỳ trong VASP, thiết lập phiếm hàm hiệu chỉnh vdW-DF2, cho đến kỹ thuật phân tích sâu sắc bản chất liên kết bằng các công cụ hiện đại như AIMAll và NBO 5.G.
  2. Kỹ sư thiết kế thiết bị cảm biến y sinh và vi điện tử: Nắm bắt chính xác cơ chế biến đổi độ dẫn điện và mức độ dịch chuyển mức Fermi khi bề mặt vật liệu 2D tiếp xúc với các khí chỉ điểm sinh học, từ đó thiết kế các cấu trúc transistor màng mỏng (TFT) có độ chọn lọc cao.
  3. Học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Hóa lý: Là tài liệu tham khảo có hệ thống về cơ sở lý thuyết lượng tử, từ phương trình Schrödinger, gần đúng Born-Oppenheimer, Hartree-Fock, DFT cho đến động học hấp phụ bề mặt và các dạng liên kết phi cộng hóa trị.
  4. Chuyên gia kiểm soát ô nhiễm và công nghệ môi trường: Cung cấp cái nhìn toàn diện về đặc tính hóa lý, độc tính và khả năng thu hồi các hợp chất hữu cơ dễ bay hơi độc hại trong không khí, phục vụ việc xây dựng quy trình quan trắc và xử lý khí thải hiệu quả.

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

Tại sao vật liệu MoX2 dạng đơn lớp lại có tiềm năng cảm biến khí vượt trội hơn dạng khối?
MoS2 dạng khối là chất bán dẫn có vùng cấm gián tiếp khoảng 1,3 eV, trong khi đơn lớp MoS2 chuyển sang vùng cấm trực tiếp rộng khoảng 1,9 eV. Cấu trúc 2D đơn lớp mang lại diện tích tiếp xúc bề mặt tối đa, nơi mọi nguyên tử đều có thể tham gia tương tác với phân tử khí, giúp tín hiệu dẫn truyền điện tử thay đổi nhạy bén hơn gấp nhiều lần.

Việc pha tạp nguyên tử Fe làm thay đổi cơ chế hấp phụ của bề mặt MoX2 như thế nào?
Bề mặt MoX2 tinh khiết chỉ tương tác với VOC qua lực phân tán van der Waals yếu ($E_A$ từ -0,15 eV đến -0,48 eV). Khi pha tạp Fe, các orbital 3d chưa điền đầy của sắt tạo liên kết phối trí mạnh với cặp electron tự do của nguyên tử dị tố (O, N) trong VOC, tăng năng lượng hấp phụ lên trên 150% và chuyển hóa sang cơ chế hấp phụ hóa học.

Thuyết nguyên tử trong phân tử (AIM) chứng minh sự hình thành liên kết bằng cách nào?
Lý thuyết AIM xác định các điểm tới hạn liên kết (BCP) tại vị trí triệt tiêu gradient mật độ electron. Dựa vào giá trị mật độ electron $\rho(r) > 0,04$ au và dấu của toán tử Laplace $\nabla^2\rho(r)$, AIM giúp phân định rõ ràng giữa liên kết hóa học cộng hóa trị, tương tác tĩnh điện phân cực và liên kết hydro yếu.

Năng lượng bền hóa bậc hai $E^{(2)}$ trong phân tích NBO phản ánh điều gì?
Đại lượng $E^{(2)}$ đo lường mức độ tương tác siêu liên hợp và truyền điện tích giữa orbital bị chiếm đầy (donor) của phân tử VOC và orbital còn trống (acceptor) của bề mặt vật liệu. Giá trị $E^{(2)}$ càng cao chứng tỏ tương tác cho nhận electron càng mạnh, góp phần làm tăng độ bền nhiệt động của phức hấp phụ.

Khả năng ứng dụng thực tế lớn nhất của đề tài này trong y học là gì?
Nghiên cứu tạo cơ sở lý thuyết vững chắc để chế tạo các cảm biến sinh học phát hiện dấu ấn VOC trong khí thở ra. Kỹ thuật này cho phép sàng lọc và phát hiện sớm bệnh ung thư phổi với độ chính xác cao, chi phí thấp, không gây đau đớn và có thể thực hiện thường quy tại các cơ sở y tế.

Kết luận

  • Xác lập mô hình cấu trúc hấp phụ tối ưu: Hoàn thành khảo sát toàn diện cấu hình không gian bền vững của 5 phân tử VOC tiêu biểu trên các bề mặt MoS2 và MoSe2 đơn lớp.
  • Chứng minh vai trò đột phá của pha tạp Fe: Kim loại chuyển tiếp Fe giúp tăng năng lượng hấp phụ từ 150% đến 320%, rút ngắn khoảng cách tương tác xuống dưới 2,15 Å và chuyển dịch cơ chế sang hấp phụ hóa học.
  • Làm sáng tỏ bản chất tương tác vi mô: Kết hợp thành công các phân tích AIM, NBO và MEP để lượng hóa sự chuyển dịch mật độ electron và vai trò của các liên kết phối trí dị tố - kim loại.
  • Cung cấp cơ sở dữ liệu tin cậy: Đóng góp hệ thống thông số hóa lý tính toán chuẩn xác cho các hệ vật liệu dichalcogenide kim loại chuyển tiếp hai chiều.
  • Định hướng ứng dụng y sinh và môi trường: Mở ra triển vọng phát triển cảm biến hơi thở chẩn đoán ung thư sớm và công nghệ màng lọc xử lý khí thải độc hại.

Kế hoạch tiếp theo hướng đến việc triển khai mô phỏng động lực học phân tử Born-Oppenheimer (BOMD) ở nhiệt độ phòng và chế tạo thực nghiệm màng cảm biến MoX2-Fe trong giai đoạn 2024 - 2025. Hãy tham khảo toàn văn luận văn thạc sĩ chuyên ngành Hóa lý thuyết và Hóa lý để tiếp cận chi tiết các bộ dữ liệu tính toán và giải pháp phát triển vật liệu tiên tiến.