Tổng quan nghiên cứu

Kể từ khi graphene được khám phá và phân tách thành công vào năm 2004 với độ linh động hạt tải vượt trội trên 15.000 cm2/(V·s) cùng mô đun Young đạt mức xấp xỉ 1.000 GPa, làn sóng nghiên cứu vật liệu hai chiều (2D) đã bùng nổ mạnh mẽ trên toàn cầu. Tuy nhiên, rào cản lớn nhất của graphene chính là vùng cấm bằng 0 eV và tương tác spin - quỹ đạo cực kỳ yếu, khiến việc ứng dụng vào các linh kiện bán dẫn như transistor hiệu ứng trường gặp nhiều trở ngại do không thể tắt mở dòng điện một cách tối ưu. Nhằm khắc phục hạn chế này, hướng nghiên cứu chuyển dịch mạnh mẽ sang họ vật liệu bán dẫn hai chiều đơn lớp, tiêu biểu là các hợp chất monochalcogenide nhóm III (MX với M = Ga, In; X = S, Se, Te) và các cấu trúc biến thể bất đối xứng Janus (M2XY và GaInXO).

Mục tiêu cốt lõi của nghiên cứu là ứng dụng lý thuyết phiếm hàm mật độ (DFT) để giải mã có hệ thống cấu trúc hình học, độ bền nhiệt động, phổ dao động mạng phonon, cấu trúc vùng năng lượng điện tử và cơ chế điều khiển tính chất điện tử dưới tác động của biến dạng cơ học cũng như điện trường ngoài. Phạm vi khảo sát tập trung vào 6 cấu trúc đơn lớp đối xứng MX, các cấu trúc Janus đơn lớp M2XY cùng 6 cấu hình dị thể GaInXO trong giai đoạn thực hiện 2019-2022. Ý nghĩa của công trình được khẳng định qua việc công bố thành công 08 bài báo khoa học chất lượng cao, bao gồm 06 công trình trên các tạp chí quốc tế thuộc danh mục SCIE uy tín. Kết quả đạt được cung cấp cơ sở dữ liệu vật lý toàn diện với hằng số mạng biến thiên từ 3,59 Å đến 4,38 Å, xác lập tiền đề khoa học vững chắc cho việc thiết kế linh kiện quang điện tử, cảm biến khí và công nghệ xúc tác quang hóa thế thái mới.

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 phiếm hàm mật độ (DFT) thông qua hai định lý cơ bản Hohenberg - Kohn (1964) và hệ phương trình tự nhất quán Kohn - Sham (1965). Thay vì giải hàm sóng đa hạt N-electron phức tạp với 3N biến số không gian, DFT chuyển đổi bài toán sang việc mô tả mật độ electron phụ thuộc duy nhất vào 3 biến tọa độ không gian, giúp tiết kiệm chi phí tính toán đồng thời giữ nguyên tính chính xác vật lý. Năng lượng trao đổi - tương quan được xử lý thông qua phép gần đúng gradient suy rộng (GGA) với phiếm hàm Perdew - Burke - Ernzerhof (PBE) và phiếm hàm lai Heyd - Scuseria - Ernzerhof (HSE06) với thông số sàng lọc tương tác tầm ngắn khoảng 0,15 Å-1 nhằm khắc phục hiện tượng đánh giá thấp độ rộng vùng cấm của GGA.

Bên cạnh đó, hiệu ứng tương tác spin - quỹ đạo (SOC) được tích hợp trong phương trình Kohn - Sham thông qua phương pháp PBE+SOC để phân tích chính xác sự phân tách mức năng lượng tại các điểm đối xứng cao. Tương tác van der Waals yếu giữa các lớp nguyên tử được hiệu chỉnh bằng phương pháp bán thực nghiệm DFT-D2 của Grimme. Các khái niệm cốt lõi như vùng Brillouin với các điểm đối xứng đặc trưng (Gamma, K, M), cực đại vùng hóa trị (VBM), cực tiểu vùng dẫn (CBM), mật độ trạng thái riêng phần (PDOS) và công thoát điện tử tạo nên khung lý thuyết chuẩn xác để đánh giá toàn diện hành vi lượng tử của hệ vật liệu.

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

Toàn bộ quá trình tính toán số trị được tiến hành tự động hóa trên gói phần mềm mô phỏng mã nguồn mở Quantum ESPRESSO. Không gian k trong vùng Brillouin được lấy mẫu kỹ lưỡng bằng phương pháp lưới Monkhorst - Pack với kích thước (15 x 15 x 1) k-mesh. Năng lượng ngưỡng cắt cho hệ sóng phẳng được thiết lập ở mức 500 eV, đảm bảo tiêu chuẩn hội tụ năng lượng tuyệt đối đạt mức 10^-6 eV giữa các bước lặp và ngưỡng lực tác dụng lên mỗi nguyên tử nhỏ hơn 10^-3 eV/Å. Nhằm triệt tiêu hoàn toàn tương tác giả định giữa các lớp lân cận do điều kiện biên tuần hoàn, một khoảng chân không dày 20 Å được thiết lập theo phương vuông góc với mặt phẳng đơn lớp (trục z).

Để khảo sát độ ổn định động học, phổ tán sắc phonon được tính toán thông qua lý thuyết nhiễu loạn phiếm hàm mật độ (DFPT) trên siêu ô cơ sở kích thước (6 x 6 x 1) gồm 144 nguyên tử. Độ bền nhiệt động được kiểm chứng qua mô phỏng động lực học phân tử lượng tử (AIMD) sử dụng thuật toán nhiệt động Nosé - Hoover ở nhiệt độ phòng 300 K với bước thời gian 1 femtigiây kéo dài liên tục trong 5 picogiây. Các hệ số đàn hồi độc lập Cij được xác định bằng cách áp dụng biến dạng kéo nén đơn trục và song trục trong biên độ nhỏ từ -0,015 đến +0,015 với bước nhảy 0,005 (gồm 7 điểm dữ liệu) nhằm đối chiếu với tiêu chuẩn cơ học Born - Huang. Toàn bộ tiến trình nghiên cứu kéo dài 3 năm (2019-2022) dưới sự tài trợ từ Quỹ Đổi mới sáng tạo Vingroup qua đề tài VINIF.2020.TS.100.

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

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

Phát hiện thứ nhất: Về mặt hình học và độ bền cơ - nhiệt, 6 cấu trúc đơn lớp MX (GaS, GaSe, GaTe, InS, InSe, InTe) kết tinh ở nhóm đối xứng không gian D3h dạng lục giác với hằng số mạng a tăng dần từ 3,59 Å (đối với GaS) lên 4,38 Å (đối với InTe) do sự gia tăng bán kính nguyên tử từ lưu huỳnh đến tellur. Năng lượng liên kết kết tụ Ecoh dao động trong khoảng từ -3,5 eV/nguyên tử đến -5,2 eV/nguyên tử, khẳng định cấu trúc có độ bền liên kết nội tại cao. Phổ phonon DFPT của toàn bộ hệ vật liệu không xuất hiện tần số ảo âm trong toàn bộ vùng Brillouin, kết hợp cùng mô phỏng AIMD tại 300 K chứng minh mạng tinh thể duy trì cấu trúc nguyên vẹn mà không bị phá vỡ liên kết.

Phát hiện thứ hai: Về đặc trưng vùng năng lượng, các đơn lớp MX là chất bán dẫn có vùng cấm gián tiếp với độ rộng dải năng lượng tính theo PBE đạt từ 1,25 eV đến 2,60 eV, và tăng lên mức 1,85 eV đến 3,45 eV khi sử dụng phiếm hàm lai HSE06. Khi chuyển sang cấu trúc bất đối xứng Janus M2XY (như Ga2SSe, Ga2STe, In2SeTe) và GaInXO (như SGaInO, SeGaInO, OGaInS), sự phá vỡ tính đối xứng gương qua mặt phẳng làm biến đổi căn bản cấu trúc dải năng lượng. Một số cấu hình Janus chuyển hóa thành công từ vùng cấm gián tiếp sang vùng cấm trực tiếp tại điểm Gamma, mang lại ưu thế phát xạ quang tử vượt trội.

Phát hiện thứ ba: Sự xuất hiện của lưỡng cực điện nội tại tạo ra sự chênh lệch thế tĩnh điện giữa hai bề mặt của đơn lớp Janus GaInXO, dẫn đến độ chênh lệch công thoát điện tử Delta Phi đạt từ 1,50 eV đến 2,85 eV. Độ linh động hạt tải của điện tử dọc theo các trục truyền dẫn đạt giá trị ấn tượng, vượt ngưỡng 1.200 cm2/(V·s) trong các cấu hình tối ưu hóa. Khi áp dụng biến dạng cơ học song trục từ -10% đến +10%, độ rộng vùng cấm suy giảm tuyến tính hoặc phi tuyến, tạo ra bước chuyển pha bán dẫn - kim loại rõ rệt ở mức biến dạng nén trên 8%.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự thay đổi tính chất điện tử trong cấu trúc Janus bắt nguồn từ sự khác biệt độ âm điện giữa các nguyên tử chalcogen (S, Se, Te) và oxy, dẫn đến sự phân bố bất đối xứng của mật độ điện tích không gian theo phương vuông góc. Sự phân cực điện trường nội tại này hoạt động như một điện trường định hướng tự phát, uốn cong các dải năng lượng biên CBM và VBM.

Khi so sánh với các nghiên cứu tương tự trên vật liệu dichalcogenide kim loại chuyển tiếp (TMDs) như MoS2 hay WS2, cấu trúc Janus monochalcogenide nhóm III thể hiện khả năng phản ứng nhạy bén hơn đáng kể dưới tác động của ứng suất cơ học ngoài. Trong khi TMDs đòi hỏi mức biến dạng kéo lớn trên 10% để chuyển đổi loại vùng cấm, hệ đơn lớp GaInXO chỉ cần mức biến dạng điều biến khoảng 4% đến 6% để đạt được sự chuyển đổi tương đương.

Dữ liệu nghiên cứu này có thể được trình diễn một cách trực quan thông qua các đồ thị phân tán dải năng lượng E(k) kết hợp biểu đồ mật độ trạng thái riêng phần (PDOS), thể hiện rõ sự đóng góp của các orbital nguyên tử Ga-4p, In-5p và chalcogen-p tại đỉnh vùng hóa trị. Đồng thời, các bảng so sánh hằng số đàn hồi C11, C12 cùng bản đồ đường mức thế tĩnh điện 2D là minh chứng rõ nét cho sự phân cực không gian của hệ cấu trúc.

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

Thứ nhất, triển khai thử nghiệm tổng hợp thực nghiệm các cấu trúc đơn lớp Janus GaInXO và M2XY thông qua kỹ thuật lắng đọng pha hơi hóa học (CVD) hoặc epitaxy chùm phân tử (MBE). Mục tiêu định lượng cần hướng tới là đạt được màng tinh thể đơn lớp đồng nhất có kích thước diện tích trên 10 micromet vuông với độ dày dưới 1,0 nm trong khung thời gian 12 đến 18 tháng, do các phòng thí nghiệm trọng điểm quốc gia về công nghệ nano chủ trì thực hiện.

Thứ hai, thiết kế và chế tạo thử nghiệm linh kiện transistor hiệu ứng trường (FET) thế hệ mới trên cơ sở kênh dẫn đơn lớp Janus monochalcogenide. Chỉ số mục tiêu cần đạt là tỷ số đóng mở dòng điện (Ion/Ioff) vượt ngưỡng 10^6 và độ linh động hạt tải thực tế duy trì trên 800 cm2/(V·s) tại nhiệt độ phòng trong vòng 24 tháng, do các nhóm nghiên cứu kỹ thuật vi điện tử bán dẫn đảm nhận.

Thứ ba, ứng dụng đơn lớp Janus GaInXO làm vật liệu quang xúc tác phân tách nước sinh khí hydro sạch. Cần tận dụng điện trường nội tại để tối ưu hóa sự phân tách cặp electron - lỗ trống, điều chỉnh vị trí mép vùng CBM và VBM bao bọc thế oxy hóa khử chuẩn của nước (-4,44 eV và -5,67 eV so với mức chân không), hướng đến mục tiêu đạt hiệu suất chuyển đổi năng lượng mặt trời thành hydro (STH) trên 15% trong vòng 18 tháng tại các viện nghiên cứu năng lượng tái tạo.

Thứ tư, mở rộng mô hình tính toán đa quy mô bằng cách kết hợp DFT với thế tương tác máy học (Machine Learning Interatomic Potentials). Giải pháp này cho phép tăng tốc độ sàng lọc cấu trúc lên gấp 100 lần, hoàn thành việc khảo sát hơn 50 cấu trúc dị lớp van der Waals trong vòng 12 tháng, do các nhóm nghiên cứu vật lý tính toán và khoa học dữ liệu phối hợp triển khai.

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

Nhóm nghiên cứu sinh và học viên cao học chuyên ngành Vật lý lý thuyết, Vật lý chất rắn và Khoa học vật liệu: Tài liệu cung cấp quy trình chuẩn hóa từ việc xây dựng mô hình nguyên tử, thiết lập tham số tính toán DFT chuyên sâu trên Quantum ESPRESSO đến phân tích phổ phonon DFPT và mô phỏng AIMD.

Nhóm kỹ sư thiết kế linh kiện bán dẫn vi điện tử và công nghệ nano: Luận văn cung cấp tập dữ liệu toàn diện về cấu trúc vùng năng lượng, công thoát điện tử, khối lượng hiệu dụng và độ linh động hạt tải, hỗ trợ trực tiếp cho quá trình mô phỏng và phát triển transistor FET thế hệ tiếp theo.

Nhóm nhà khoa học lĩnh vực năng lượng tái tạo và hóa học quang xúc tác: Luận văn là tài liệu tham khảo giá trị về cơ chế phân cực bề mặt và thế tĩnh điện nội tại trong vật liệu Janus, giúp định hướng thiết kế các chất xúc tác quang hóa phân tách nước tạo hydro hiệu quả cao.

Nhóm giảng viên đại học giảng dạy các học phần Vật lý chất rắn nâng cao và Vật lý tính toán: Công trình đóng vai trò như một giáo trình chuyên khảo thực nghiệm số, cung cấp nhiều ví dụ minh họa sinh động về mối liên hệ giữa cấu trúc vi mô và tính chất vĩ mô của vật liệu thấp chiều.

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

Vật liệu monochalcogenide nhóm III có ưu điểm gì nổi bật so với graphene truyền thống? Khác với graphene có vùng cấm bằng 0 eV gây khó khăn cho việc đóng ngắt dòng điện trong transistor, đơn lớp monochalcogenide nhóm III sở hữu độ rộng vùng cấm nội tại phù hợp từ 1,25 eV đến 3,45 eV. Điều này cho phép vật liệu duy trì tỷ số đóng ngắt dòng điện cao trong khi vẫn giữ được độ linh động hạt tải lớn trên 1.000 cm2/(V·s).

Tại sao cấu trúc bất đối xứng Janus lại tạo ra các tính chất vật lý vượt trội? Cấu trúc Janus được tạo thành bằng cách phá vỡ tính đối xứng gương của màng đơn lớp, thay thế một mặt nguyên tử chalcogen bằng nguyên tử khác. Sự bất đối xứng này sinh ra điện trường nội tại vuông góc với mặt phẳng, tạo ra chênh lệch công thoát điện tử từ 1,5 eV đến 2,8 eV, giúp phân tách các hạt mang điện và nâng cao hiệu suất quang điện.

Biến dạng cơ học điều biến cấu trúc dải năng lượng của vật liệu như thế nào? Dưới tác động của biến dạng cơ học song trục từ -10% đến +10%, khoảng cách giữa các nguyên tử thay đổi dẫn đến sự tái phân bố mật độ điện tích xen phủ. Kết quả tính toán cho thấy biến dạng kéo làm dịch chuyển dải dẫn xuống thấp và đẩy dải hóa trị lên cao, thu hẹp vùng cấm và kích thích chuyển pha bán dẫn sang kim loại khi biến dạng vượt 8%.

Phiếm hàm lai HSE06 mang lại ưu thế gì so với phép gần đúng chuẩn GGA-PBE? Phép gần đúng GGA-PBE thường đánh giá thấp độ rộng vùng cấm của chất bán dẫn từ 30% đến 50% do sai số tự tương tác điện tử. Phiếm hàm lai HSE06 trộn thêm 25% thành phần trao đổi Hartree - Fock tầm ngắn, giúp hiệu chỉnh chính xác thế trao đổi và mang lại kết quả vùng cấm tiệm cận sát với số liệu đo đạc thực nghiệm.

Khả năng ứng dụng thực tế rõ nét nhất của cấu trúc Janus GaInXO là gì? Cấu trúc Janus GaInXO có tiềm năng ứng dụng lớn nhất trong lĩnh vực quang xúc tác tách nước sinh hydro nhờ vị trí dải năng lượng CBM cao hơn thế khử H+/H2 và VBM thấp hơn thế oxy hóa O2/H2O, kết hợp cùng điện trường nội tại giúp giảm thiểu hiện tượng tái hợp quang sinh giữa electron và lỗ trống.

Kết luận

  • Luận văn đã thiết lập bức tranh khoa học toàn diện về cấu trúc tinh thể, phổ phonon và các tham số đàn hồi của 6 đơn lớp monochalcogenide nhóm III và các biến thể Janus.
  • Khẳng định 100% các cấu trúc khảo sát đạt độ ổn định động học và nhiệt động vững chắc ở nhiệt độ phòng 300 K thông qua phân tích DFPT và AIMD.
  • Phát hiện cơ chế chuyển đổi vùng cấm từ gián tiếp sang trực tiếp và chuyển pha kim loại dưới tác động của biến dạng cơ học song trục biên độ từ -10% đến +10%.
  • Xác định độ chênh lệch thế tĩnh điện bề mặt đạt mức 1,50 eV đến 2,85 eV trong cấu trúc Janus GaInXO, mở ra hướng ứng dụng đột phá cho pin quang điện và tách nước quang hóa.
  • Thành tựu nghiên cứu được bảo chứng qua 08 công trình khoa học công bố trên các tạp chí chuyên ngành uy tín, trong đó có 06 bài báo thuộc danh mục SCIE.

Lộ trình nghiên cứu tiếp theo cần tập trung vào việc mô phỏng dị lớp tiếp xúc bán dẫn van der Waals trong 12 tháng tới và tiến hành thử nghiệm chế tạo thực nghiệm màng đơn lớp trong 24 tháng. Hãy kết nối ngay với các nhóm nghiên cứu vật lý tính toán để đón đầu công nghệ vật liệu hai chiều tương lai!