Tổng quan nghiên cứu

Theo thống kê từ Tổ chức Y tế Thế giới (WHO), ô nhiễm không khí gây ra khoảng 4,2 triệu ca tử vong sớm mỗi năm và chiếm tỷ lệ 7,6% tổng số ca tử vong trên toàn cầu. Các loại khí thải độc hại như carbon dioxide (CO2), carbon monoxide (CO), nitrogen monoxide (NO), nitrogen dioxide (NO2) và ammonia (NH3) phát sinh mạnh mẽ từ hoạt động công nghiệp cùng phương tiện giao thông, gây tổn thương nghiêm trọng đến hệ hô hấp và tim mạch của con người. Để kiểm soát và cảnh báo sớm, việc phát triển các cảm biến khí có độ nhạy cao và khả năng chọn lọc tốt là yêu cầu cấp thiết.

Graphene, vật liệu hai chiều (2D) với cấu trúc một lớp nguyên tử carbon liên kết sp2, sở hữu độ linh động điện tử vượt trội nhưng lại thiếu độ mở vùng cấm (zero bandgap) và độ chọn lọc khí tự nhiên còn hạn chế. Nhằm khắc phục nhược điểm này, luận văn tập trung nghiên cứu việc kết hợp graphene với các chất nền điện môi phù hợp gồm hexagonal boron nitride (h-BN) và silicon dioxide (α-SiO2) để tạo nên cấu trúc dị thể (heterostructure). Mục tiêu chính của đề tài là khảo sát cấu hình tối ưu, đánh giá năng lượng hấp phụ, khoảng cách tương tác và phân tích mật độ trạng thái điện tử (DOS) của 5 loại khí độc cùng hơi nước (H2O). Nghiên cứu được thực hiện trong chương trình Thạc sĩ Công nghệ Nano (mã ngành 8440140.11QTD) tại Trường Đại học Việt Nhật thuộc Đại học Quốc gia Hà Nội trong giai đoạn năm 2018 đến năm 2020. Kết quả nghiên cứu góp phần mở rộng vùng cấm của graphene lên mức 0,053 eV, nâng cao độ chọn lọc phát hiện khí độc và thúc đẩy ứng dụng cơ chế hấp phụ vật lý (physisorption) nhằm giảm thiểu rác thải điện tử từ các loại cảm biến dùng một lần.

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 thuyết phiếm hàm mật độ (Density Functional Theory - DFT) của Hohenberg - Kohn và phương pháp tính toán Kohn - Sham. Để giải quyết tương tác giữa các hệ điện tử nhiều hạt, phép gần đúng mật độ cục bộ (LDA) kết hợp phép gần đúng gradient tổng quát (GGA-PBE) được sử dụng để tối ưu hóa hàm thế trao đổi - tương quan. Đồng thời, hiệu chỉnh lực phân tán van der Waals theo mô hình revPBE-vdW được tích hợp nhằm mô tả chuẩn xác các tương tác tầm xa giữa phân tử khí và bề mặt vật liệu 2D.

Mô hình nghiên cứu bao gồm hai hệ cấu trúc dị thể chính: graphene trên đế hexagonal boron nitride (G/h-BN) và graphene trên đế silicon dioxide (G/α-SiO2) tinh thể α-quartz với mặt phẳng định hướng (0001). Ba khái niệm then chốt chi phối cơ chế tương tác trong mô hình gồm: năng lượng xếp chồng (stacking stability), năng lượng hấp phụ (adsorption energy) và mật độ trạng thái điện tử (DOS/PDOS). Các thông số về orbital phân tử bị chiếm cao nhất (HOMO) và orbital phân tử trống thấp nhất (LUMO) được phân tích để làm sáng tỏ cơ chế dịch chuyển điện tích.

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

Nguồn dữ liệu tính toán lượng tử từ nguyên lý ban đầu (ab initio) được tạo lập thông qua gói phần mềm mô phỏng chuyên dụng Vienna Ab initio Simulation Package (VASP), kết hợp cùng công cụ trực quan hóa VESTA. Nghiên cứu sử dụng siêu ô mạng (supercell) 3×3 chứa 36 nguyên tử cho hệ G/h-BN với hằng số mạng a = b = 7,50 Å và siêu ô mạng 2×2 cho hệ G/α-SiO2 với hằng số mạng a = b = 9,82 Å. Việc lựa chọn cỡ mẫu siêu ô mạng này nhằm đảm bảo không gian phân bố cho các phân tử khí, triệt tiêu tương tác giả giữa các phân tử lân cận trong điều kiện biên tuần hoàn. Vùng chân không được thiết lập trên 20 Å theo phương trục z (cụ thể là 21,47 Å cho G/h-BN và 18 Å cho G/α-SiO2) để loại bỏ hoàn toàn tương tác giữa các phiến tuần hoàn lặp lại.

Phương pháp lấy mẫu không gian khảo sát 4 vị trí hình học hấp phụ gồm: vị trí đỉnh carbon 1 (Top 1), vị trí đỉnh carbon 2 (Top 2), vị trí tâm vòng lục giác (Hollow) và vị trí giữa liên kết carbon - carbon (Bridge) dưới 2 hướng định hướng thẳng đứng (Upright) và nằm ngang (Horizontal). Năng lượng cắt sóng phẳng (cut-off energy) được chọn ở mức 520 eV cho hệ G/h-BN và 600 eV cho hệ G/α-SiO2. Lưới tích phân K-points trong vùng Brillouin được thiết lập tương ứng là 8×8×1 và 6×6×1. Quy trình tối ưu hóa cấu hình được kiểm soát chặt chẽ với tiêu chuẩn hội tụ năng lượng đạt 10^-5 eV và lực Hellman - Feynman tác dụng lên mỗi nguyên tử nhỏ hơn 0,01 eV/Å. Timeline tính toán mô phỏng được triển khai xuyên suốt 24 tháng từ năm 2018 đến năm 2020.

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

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

Thứ nhất, phân tích độ lệch mạng tinh thể (lattice mismatch) cho thấy hệ G/h-BN có mức bất tương xứng rất nhỏ chỉ khoảng 1,70%, trong khi hệ G/α-SiO2 đạt khoảng 0,48%. Cấu hình xếp chồng kiểu AB với nguyên tử C nằm trực tiếp trên nguyên tử B (C on top B) đạt độ bền vững tối ưu với năng lượng xếp chồng Estack = -0,025 eV/nguyên tử, thấp hơn mức -0,023 eV/nguyên tử của cấu hình C trên N. Cấu trúc G/h-BN tối ưu đã mở thành công vùng cấm của graphene lên 0,053 eV với mức năng lượng Fermi đạt -1,692 eV.

Thứ hai, đối với nhóm khí gốc carbon, phân tử CO2 hấp phụ bền nhất tại vị trí Bridge theo hướng nằm ngang với năng lượng hấp phụ Eads = -0,224 eV và khoảng cách từ tâm phân tử đến bề mặt graphene là d = 3,222 Å. Phân tử CO cũng đạt trạng thái cân bằng tại vị trí Bridge theo hướng nằm ngang với năng lượng hấp phụ Eads = -0,142 eV và khoảng cách d = 3,415 Å.

Thứ ba, đối với nhóm khí gốc nitơ, phân tử thuận từ NO hấp phụ thuận lợi nhất tại vị trí Top 2 theo phương ngang với Eads = -0,185 eV và khoảng cách d = 3,204 Å. Đáng chú ý, phân tử NO2 thể hiện tương tác mạnh nhất tại vị trí Top 1 với năng lượng hấp phụ đạt -0,365 eV và khoảng cách tương tác ngắn nhất là 3,021 Å, vượt trội hơn khoảng 63% về năng lượng liên kết so với CO2 và cao hơn 157% so với CO.

Thảo luận kết quả

Toàn bộ dữ liệu tính toán có thể được hệ thống hóa trực quan thông qua các bảng tổng hợp năng lượng hấp phụ tương ứng với từng tọa độ không gian và các biểu đồ mật độ trạng thái (DOS, PDOS). Biểu đồ DOS chỉ ra rằng vùng cấm của các phân tử khí không bị phá vỡ sau khi hấp phụ, đồng thời mức năng lượng Fermi của toàn hệ dịch chuyển từ -1,692 eV xuống -2,044 eV khi có mặt phân tử khí. Đối với NO2, orbital LUMO (3π*) nằm dưới điểm Dirac tạo điều kiện cho quá trình chuyển dịch một phần điện tích từ nguyên tử oxy của NO2 sang các nguyên tử carbon của graphene.

Khoảng cách hấp phụ trong tất cả các trường hợp đều duy trì ở mức trên 3,0 Å, khẳng định bản chất tương tác là hấp phụ vật lý chi phối bởi lực van der Waals. So với graphene tự do trong các nghiên cứu trước đây vốn có độ chọn lọc thấp, sự xuất hiện của chất nền h-BN và SiO2 giúp phân biệt rõ rệt mức biến thiên trở kháng điện tử giữa các loại khí. Đặc biệt, năng lượng hấp phụ vừa phải của NO2 (-0,365 eV) và CO (-0,142 eV) là điều kiện lý tưởng để cảm biến đạt độ nhạy cao trong môi trường thực tế nhưng vẫn dễ dàng nhả khí để tái sử dụng ở nhiệt độ phòng.

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

  • Triển khai chế tạo thực nghiệm màng dị thể G/h-BN bằng kỹ thuật lắng đọng pha hơi hóa học (CVD) hoặc bóc tách cơ học: Các phòng thí nghiệm vật liệu nano cần thực hiện quy trình chế tạo thử nghiệm với mục tiêu đạt độ nhạy phát hiện khí NO2 ở nồng độ dưới 1 ppm trong khung thời gian 12 tháng.
  • Phát triển mạch đo trở kháng vi sai tích hợp cho cảm biến khí nano: Các kỹ sư điện tử và vi cơ điện tử (MEMS) cần thiết kế module đo tín hiệu trở kháng đa kênh với mục tiêu giảm thời gian hồi phục tín hiệu cảm biến xuống dưới 15 giây sau mỗi chu kỳ nhả khí, hoàn thành thử nghiệm trong vòng 6 đến 9 tháng.
  • Mở rộng mô hình tính toán DFT với các chất nền bán dẫn 2D thế hệ mới: Các nhóm nghiên cứu vật liệu tính toán cần ứng dụng phương pháp revPBE-vdW để khảo sát thêm các cấu trúc như MoS2, WS2 hoặc phosphorene nhằm nâng cao năng lượng hấp phụ chọn lọc cho khí NH3 và CO thêm khoảng 20% đến 30% trong lộ trình 18 tháng.
  • Ban hành quy chuẩn kỹ thuật cho thiết bị quan trắc khí thải tự động sử dụng cảm biến vật liệu mới: Bộ Tài nguyên và Môi trường phối hợp cùng Bộ Khoa học và Công nghệ xây dựng bộ tiêu chuẩn kiểm định cảm biến nano, hướng tới mục tiêu giảm sai số đo đạc nồng độ khí độc công nghiệp xuống dưới 5% trong thời hạn 24 tháng.

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

  • Nhà nghiên cứu và học viên cao học chuyên ngành Công nghệ Nano, Khoa học Vật liệu: Tiếp cận quy trình tính toán DFT hoàn chỉnh từ khâu thiết lập siêu ô mạng, tối ưu hóa cấu hình dị thể 2D đến kỹ thuật phân tích mật độ trạng thái điện tử (DOS/PDOS) trên phần mềm VASP.
  • Kỹ sư phát triển cảm biến khí và linh kiện bán dẫn: Khai thác các thông số định lượng về vị trí hấp phụ tối ưu, năng lượng liên kết và khoảng cách tương tác để thiết kế cấu trúc màng cảm biến có khả năng tái sử dụng cao.
  • Chuyên gia phân tích và quản lý trạm quan trắc môi trường: Nắm bắt cơ chế chọn lọc tín hiệu điện tử của vật liệu nano đối với các khí nhà kính và khí độc hại, từ đó lựa chọn giải pháp thiết bị quan trắc phù hợp.
  • Doanh nghiệp sản xuất thiết bị đo lường thông minh và IoT: Định hướng ứng dụng công nghệ vật liệu 2D vào sản xuất cảm biến thương mại giúp tiết kiệm chi phí nguyên liệu và giảm thiểu rác thải linh kiện điện tử.

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

Tại sao cần đưa chất nền h-BN hoặc SiO2 vào dưới lớp graphene trong nghiên cứu cảm biến khí?

Graphene tự do có cấu trúc dải năng lượng không có vùng cấm và khả năng chọn lọc khí kém. Việc đặt graphene lên chất nền h-BN giúp mở vùng cấm khoảng 0,053 eV nhờ phá vỡ tính đối xứng mạng con, đồng thời giữ bề mặt nhẵn không có bẫy điện tích, giúp tăng cường độ ổn định và độ chọn lọc tín hiệu trở kháng khi phát hiện khí độc.

Cơ chế hấp phụ vật lý mang lại ưu thế gì cho việc chế tạo cảm biến khí?

Hấp phụ vật lý có năng lượng liên kết trung bình từ -0,142 eV đến -0,365 eV và khoảng cách tương tác lớn trên 3,0 Å. Đặc điểm này giúp phân tử khí bám vào bề mặt mà không làm đứt gãy liên kết hóa học, cho phép cảm biến nhanh chóng giải hấp phụ để tái sử dụng nhiều lần mà không bị suy hao hiệu năng.

Vì sao cấu hình xếp chồng AB với carbon nằm trên boron lại ổn định hơn cấu hình carbon trên nitơ?

Nguyên tử nitơ có độ âm điện cao hơn boron (3,04 so với 2,04) nên tích điện âm trong khi boron mang điện tích dương nhẹ. Ở cấu hình C on top B, các anion N nằm đúng tâm vòng lục giác của graphene, làm giảm tối đa lực đẩy tĩnh điện với đám mây điện tử C-C và tạo năng lượng xếp chồng ổn định hơn mức -0,025 eV/nguyên tử.

Khí độc nào có tương tác mạnh nhất trên bề mặt dị thể G/h-BN?

Khí NO2 thể hiện tương tác mạnh nhất với năng lượng hấp phụ đạt -0,365 eV tại vị trí Top 1 và khoảng cách hấp phụ là 3,021 Å. Cấu trúc thuận từ cùng orbital LUMO (3π*) nằm dưới mức Fermi giúp quá trình nhận điện tử từ mạng graphene diễn ra mạnh mẽ hơn so với các khí CO2, CO hay NO.

Phương pháp hiệu chỉnh revPBE-vdW đóng vai trò gì trong độ chính xác của nghiên cứu?

Các phiếm hàm DFT truyền thống thường bỏ sót hoặc đánh giá không đầy đủ lực tán xạ van der Waals tầm xa. Việc bổ sung hiệu chỉnh revPBE-vdW trong gói VASP giúp tính toán chính xác năng lượng tương tác giữa phân tử khí và bề mặt 2D, đảm bảo độ lệch khoảng cách hấp phụ được kiểm soát chặt chẽ dưới 0,05 Å.

Kết luận

  • Luận văn xác định thành công cấu hình dị thể G/h-BN ổn định nhất theo kiểu xếp chồng AB (C on top B) với năng lượng xếp chồng Estack = -0,025 eV/nguyên tử và độ lệch mạng thấp chỉ 1,70%.
  • Làm sáng tỏ cơ chế mở vùng cấm của graphene đạt mức 0,053 eV trên nền chất nền h-BN, tạo tiền đề quan trọng cho việc ứng dụng vật liệu 2D vào các thiết bị bán dẫn nano.
  • Khẳng định bản chất hấp phụ vật lý thuận nghịch của 5 loại khí độc (CO2, CO, NO, NO2, NH3) và hơi nước H2O với khoảng cách tương tác đều duy trì trên 3,0 Å.
  • Chứng minh tính chọn lọc vượt trội của hệ vật liệu đối với khí NO2 với năng lượng hấp phụ đạt -0,365 eV thông qua phân tích cơ chế dịch chuyển điện tích trên biểu đồ DOS và PDOS.
  • Cung cấp cơ sở dữ liệu tính toán lượng tử chuẩn xác, định hướng cho các nghiên cứu thực nghiệm chế tạo cảm biến khí nano tái sử dụng trong vòng 12 đến 24 tháng tới; các đơn vị nghiên cứu và doanh nghiệp công nghệ nên đẩy mạnh ứng dụng mô hình dị thể này vào thực tiễn sản xuất thiết bị quan trắc môi trường.