Tổng quan nghiên cứu

Kỷ nguyên vật liệu nano hai chiều chứng kiến sự phát triển bùng nổ kể từ khi graphene được khám phá với khả năng hấp thụ 2,3% quang phổ nhìn thấy và độ linh động điện tử vượt trội. Tuy nhiên, rào cản vùng cấm bằng 0 eV của graphene cũng như các đơn lớp nhóm IVA như germanene đã hạn chế đáng kể khả năng ứng dụng trong các linh kiện bán dẫn logic. Để khắc phục điểm nghẽn này, việc tạo ra các dị cấu trúc phẳng hai chiều là giải pháp đột phá nhằm kiểm soát cấu trúc điện tử và kích hoạt từ tính nội tại. Luận văn thạc sĩ chuyên ngành Vật lý chất rắn với mã số 8440104 được thực hiện tại Trường Đại học Sư phạm thuộc Đại học Thái Nguyên, hoàn thành và nghiệm thu vào tháng 11 năm 2023 với chỉ số trùng lặp Turnitin được kiểm duyệt nghiêm ngặt ở mức 12%.

Nghiên cứu tập trung giải quyết bài toán dự đoán cấu trúc tinh thể, kiểm tra độ bền động học và phân tích chi tiết tính chất điện tử cùng từ tính của các dị cấu trúc ngang hai chiều giữa đơn lớp germanene và các đơn lớp pnictogen nhóm VA gồm arsenene, antimonene và AsSb. Phạm vi nghiên cứu bao quát cả hai cấu hình tiếp xúc biên dạng armchair và zigzag với 7 tỷ lệ hóa học biến thiên. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc chứng minh khả năng mở dải năng lượng cấm từ mức 0 eV lên cực đại 0,87 eV, đồng thời tạo ra mômen từ tự phát đạt tới 0,50 mu_B dọc theo đường biên zigzag mà không cần pha tạp ngoại lai. Với độ tương thích mạng lý tưởng dưới 5,78%, nghiên cứu mở ra triển vọng to lớn trong việc phát triển các transistor hiệu ứng trường thế hệ mới và các thiết bị điện tử học spin tiêu thụ năng lượng cực thấp.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng cơ học lượng tử hệ nhiều hạt và lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT). Quá trình giải bài toán vi mô áp dụng gần đúng Born-Oppenheimer nhằm tách biệt chuyển động của electron khỏi hệ hạt nhân nguyên tử nặng hơn hàng nghìn lần, kết hợp nguyên lý biến phân để xác định trạng thái cơ bản có năng lượng cực tiểu. Nghiên cứu vận dụng chặt chẽ hai định lý nền tảng Hohenberg-Kohn nhằm chứng minh mật độ electron xác định duy nhất thế ngoài, cùng hệ phương trình tự hợp Kohn-Sham để đưa hệ electron tương tác phức tạp về hệ các hạt độc lập chuyển động trong thế hiệu dụng. Để xử lý năng lượng tương quan trao đổi, mô hình gần đúng gradient suy rộng (GGA) với tham số hóa Perdew-Burke-Ernzerhof (PBE) được áp dụng nhờ tính chuẩn xác vượt trội so với mô hình Thomas-Fermi cổ điển.

Các khái niệm then chốt được phân tích chuyên sâu gồm: mô hình sóng giả tăng cường xạ ảnh PAW (Projector Augmented Wave) tái lập tương tác electron hóa trị và ion lõi; hiện tượng lai hóa quỹ đạo sp2 xen kẽ sp3 tạo nên độ mấp mô cấu trúc dao động từ 0,69 Å đến 1,64 Å; cấu trúc dải năng lượng với các điểm đối xứng cao Gamma, K, X trong vùng Brillouin; nón Dirac đặc trưng của trạng thái bán kim loại; sự phân cực spin tự phát và hiệu ứng phân tách Rashba xuất phát từ sự phá vỡ tính đối xứng không gian tại bề mặt tiếp xúc.

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

Toàn bộ các phép tính toán mô phỏng nguyên lý ban đầu được thực thi trên phần mềm VASP (Vienna Ab initio Simulation Package), kết hợp các công cụ chuyên dụng gồm p4vasp, VESTA và Origin để xử lý số liệu, phân tích mật độ trạng thái và trực quan hóa tinh thể.

Về cỡ mẫu và phương pháp chọn mẫu mô hình: Nghiên cứu thiết lập hệ siêu ô cơ sở (supercell) một chiều tuần hoàn gồm 32 nguyên tử, tương ứng với 8 đơn vị chu kỳ cấu trúc thành phần. Phương pháp chọn mẫu cấu hình được tiến hành tuần tự theo chuỗi tỷ lệ hóa học Gem-X8-m với m biến thiên từ 1 đến 7 (trong đó X đại diện cho As, Sb hoặc AsSb), tương ứng với nồng độ germanene trải dài từ 12,5% đến 87,5%. Tổng cộng có 42 mô hình dị cấu trúc độc lập được tối ưu hóa toàn phần cho cả hai hướng tiếp xúc armchair và zigzag. Lý do lựa chọn phương pháp DFT với phiếm hàm GGA-PBE là khả năng cân bằng tối ưu giữa độ chính xác vật lý của cấu trúc tinh thể với sai số hằng số mạng dưới 1,5% so với thực nghiệm và thời gian tính toán trên hệ thống máy tính hiệu năng cao.

Các thông số kỹ thuật được kiểm soát chặt chẽ: năng lượng ngưỡng sóng phẳng (cut-off energy) đạt mức 500 eV; lưới vector sóng k trong không gian nghịch đảo lấy theo phương pháp Monkhorst-Pack với kích thước 2x12x1; tiêu chuẩn hội tụ lực Hellmann-Feynman tác dụng lên mỗi nguyên tử đạt giá trị dưới 0,01 eV/Å; sai số năng lượng trong các vòng lặp tự hợp dừng lại ở mức dưới 10^-6 eV. Toàn bộ quy trình mô phỏng và phân tích dữ liệu được triển khai hoàn chỉnh trong giai đoạn 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

Thứ nhất, về độ tương thích mạng và năng lượng tạo thành: Hằng số mạng của đơn lớp germanene đạt 4,06 Å, tạo ra độ sai lệch mạng cực nhỏ chỉ 0,73% khi ghép nối với antimonene (4,09 Å), 2,53% với AsSb và 5,78% với arsenene (3,60 Å). Phân tích nhiệt động học chỉ ra rằng hệ germanene-antimonene có năng lượng hình thành thấp hơn từ 25% đến 40% so với hệ germanene-arsenene, đồng thời cấu trúc định hướng zigzag có độ ổn định năng lượng cao hơn đáng kể so với hướng armchair.

Thứ hai, về khả năng mở vùng cấm linh hoạt ở hướng armchair: Toàn bộ các dị cấu trúc armchair đều chuyển đổi thành công từ tính chất bán kim loại của germanene (vùng cấm 0 eV) sang trạng thái bán dẫn hoàn chỉnh. Độ rộng vùng cấm đạt giá trị cực đại tại nồng độ 25% germanene (m = 2), lần lượt đạt 0,74 eV ở hệ Ge-As, 0,76 eV ở hệ Ge-Sb và đạt đỉnh 0,87 eV ở hệ Ge-AsSb. Khi tỷ lệ germanene tăng lên 87,5% (m = 7), vùng cấm giảm tuyến tính xuống mức 0,03 eV đến 0,05 eV. Nghiên cứu cũng xác nhận sự chuyển dịch linh hoạt giữa khe cấm trực tiếp (tại điểm Gamma hoặc X) và khe cấm gián tiếp tùy thuộc vào tỷ lệ m.

Thứ ba, về sự xuất hiện của từ tính tự phát dọc theo biên zigzag: Trái ngược với cấu hình armchair không phân cực spin, hướng zigzag kích hoạt từ tính nội tại mạnh mẽ. Hệ Ge-Sb duy trì mômen từ bão hòa cực kỳ ổn định trong khoảng từ 0,49 đến 0,50 mu_B trên toàn bộ dải nồng độ từ 25% đến 75% germanene (m từ 2 đến 6). Hệ Ge-As bắt đầu xuất hiện từ tính từ m = 3 (0,13 mu_B) và đạt đỉnh 0,50 mu_B ở m = 6. Hệ Ge-AsSb ghi nhận mômen từ ổn định từ 0,46 đến 0,49 mu_B trong dải m từ 2 đến 6 trước khi giảm về 0,01 mu_B ở m = 7.

Thứ tư, về các trạng thái điện tử tôpô đặc biệt: Ở cấu hình zigzag nồng độ cao (m = 7), hai hệ Ge7-As1 và Ge7-Sb1 thể hiện đặc tính của chất bán dẫn Weyl đơn cực với độ rộng vùng cấm siêu hẹp lần lượt là 22,5 meV và 19,0 meV, đi kèm hiệu ứng phân tách Rashba rõ rệt tại vùng hóa trị.

Thảo luận kết quả

Cơ chế vật lý dẫn đến sự mở vùng cấm xuất phát từ việc phá vỡ tính đối xứng mạng tổ ong nguyên bản và sự lai hóa orbital sp3 mạnh mẽ giữa các trạng thái Ge-pz với các orbital p của arsen và antimon tại mặt phân giới. Kết quả phân tích mật độ trạng thái (DOS) cho thấy các mức năng lượng cận mức Fermi hoàn toàn do các orbital nguyên tử Ge chi phối, chứng minh khe cấm được mở ra trực tiếp ngay trong lớp germanene. Về nguồn gốc từ tính biên zigzag, bản đồ mật độ spin khẳng định từ tính bắt nguồn chủ yếu từ các nguyên tử Ge ở đường biên phân cách một chiều do sở hữu liên kết pz chưa bão hòa, tạo nên sự phân tách năng lượng rõ rệt giữa hai kênh spin-up và spin-down.

So với các dị cấu trúc dọc liên kết van der Waals yếu hoặc các dị cấu trúc ngang kinh điển như graphene-stanene vốn không có từ tính nội tại, công trình này lần đầu tiên chứng minh khả năng tự phát sinh từ tính mà không cần cấy ion từ tính bên ngoài. Về phương diện trực quan hóa, toàn bộ các quy luật trên được thể hiện sinh động qua bảng đối sánh hằng số mạng, đồ thị đường biểu diễn sự suy giảm tuyến tính của độ rộng vùng cấm theo tỷ lệ phần trăm germanene, giản đồ cấu trúc dải năng lượng E-k minh họa vị trí dời điểm VBM/CBM, và bản đồ 3D phân bố mật độ spin khu trú sắc nét dọc theo đường biên một chiều.

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

Thứ nhất, tối ưu hóa quy trình chế tạo bằng phương pháp epitaxy chùm phân tử (MBE): Các phòng thí nghiệm trọng điểm và viện nghiên cứu vật liệu nano cần chuẩn hóa quy trình lắng đọng hai giai đoạn với sai số nhiệt độ đế dưới 1,5 độ C, nhằm kiểm soát độ sắc nét của mặt phân giới dị cấu trúc zigzag đạt độ hoàn hảo trên 95% trong lộ trình 2024 - 2025.

Thứ hai, thiết kế chế tạo transistor hiệu ứng trường kênh đôi (FET) thế hệ mới: Doanh nghiệp bán dẫn và các nhóm phát triển vi mạch cần ứng dụng dị cấu trúc armchair Ge-AsSb với độ rộng vùng cấm 0,87 eV làm vật liệu kênh dẫn, hướng tới mục tiêu nâng cao tỷ số dòng on/off vượt ngưỡng 10^5 và hạ thấp điện thế chuyển mạch xuống dưới 0,5 V trước năm 2026.

Thứ ba, phát triển các phần tử van spin và bộ nhớ từ trường không tiêu tán: Các chuyên gia công nghệ điện tử spin cần khai thác cấu hình zigzag Ge-Sb có mômen từ ổn định 0,50 mu_B để chế tạo các thiết bị tiếp xúc từ trở chui hầm (TMR), đặt mục tiêu nâng tỷ số từ trở lên trên 150% ở nhiệt độ phòng và mở rộng mật độ lưu trữ đạt trên 10 Terabit trên mỗi inch vuông trong giai đoạn 2025 - 2027.

Thứ tư, mở rộng khảo sát ảnh hưởng của ứng suất cơ học và điện trường ngoài: Các nhà vật lý tính toán tại các trường đại học cần tiếp tục mô phỏng hành vi của vật liệu dưới tác động của biến dạng cơ học hai trục từ -5% đến +5% và điện trường ngoài vuông góc lên tới 0,5 V/Å, nhằm thiết lập cẩm nang điều khiển pha lượng tử hoàn tất trong vòng 12 tháng tới.

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

Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý chất rắn, Khoa học vật liệu: Tiếp cận quy trình mô phỏng lượng tử chuẩn mực trên phần mềm VASP, phương pháp thiết lập siêu ô 32 nguyên tử và giải thuật phân tích orbital chuyên sâu.

Kỹ sư R&D trong ngành vi mạch bán dẫn và quang điện tử: Khai thác kho dữ liệu vùng cấm biến thiên linh hoạt từ 0,03 eV đến 0,87 eV để ứng dụng vào thiết kế cảm biến quang học, pin mặt trời nano và đi-ốt phát quang dải sóng ngắn.

Các nhà khoa học phát triển công nghệ điện tử học spin (Spintronics): Vận dụng các phát hiện về từ tính tự phát biên một chiều đạt 0,50 mu_B và hiệu ứng phân tách Rashba để kiến tạo các cổng logic lượng tử tiêu thụ điện năng tối thiểu.

Giảng viên đại học và chuyên gia công nghệ vật liệu hai chiều: Sử dụng các mô hình tinh thể, cơ chế lai hóa sp2-sp3 và phổ phonon làm học liệu giảng dạy tiên tiến cho các học phần vật lý chất rắn và vật lý nano chuyên sâu.

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

Tại sao dị cấu trúc ngang hai chiều được đánh giá cao hơn dị cấu trúc dọc?

Dị cấu trúc dọc chủ yếu hình thành nhờ lực van der Waals yếu, trong khi dị cấu trúc ngang tạo ra các liên kết cộng hóa trị vững chắc trong cùng một mặt phẳng nguyên tử. Nghiên cứu thực nghiệm năm 2020 của nhóm tác giả Ogikubo đã chứng minh liên kết cộng hóa trị này mang lại độ bền cơ học cao, đồng thời cho phép tích hợp liền mạch trên bề mặt phẳng của các linh kiện vi mạch nano với kích thước siêu mỏng dưới 1 nm.

Bằng cách nào dị cấu trúc armchair mở được vùng cấm của germanene lên đến 0,87 eV?

Đơn lớp germanene vốn là chất bán kim loại với vùng cấm 0 eV do tính đối xứng mạng tổ ong phẳng. Khi kết hợp với các đơn lớp nhóm VA theo hướng armchair, liên kết biên phá vỡ tính đối xứng này và kích hoạt lai hóa sp3 mạnh mẽ. Dữ liệu tính toán chứng minh ở nồng độ 25% germanene, độ rộng vùng cấm được mở rộng tối đa đạt 0,87 eV ở hệ Ge-AsSb và 0,76 eV ở hệ Ge-Sb, đáp ứng xuất sắc tiêu chuẩn linh kiện bán dẫn.

Từ tính tự phát trong dị cấu trúc hướng zigzag bắt nguồn từ đâu?

Phân tích mật độ trạng thái và mật độ spin xác nhận từ tính không xuất phát từ ion kim loại chuyển tiếp mà được sinh ra tại ranh giới tiếp xúc một chiều. Các nguyên tử Ge tại đường biên zigzag sở hữu orbital pz chưa bão hòa, dẫn đến sự phân cực spin mãnh liệt quanh mức Fermi và tạo ra mômen từ bão hòa ổn định từ 0,46 đến 0,50 mu_B trên dải nồng độ từ 25% đến 75% germanene.

Độ không phối hợp hằng số mạng tác động thế nào đến khả năng chế tạo thực nghiệm?

Độ sai lệch mạng càng nhỏ thì ứng suất cơ học nội tại càng thấp, ngăn ngừa khuyết tật mạng khi tổng hợp. Với độ không phối hợp chỉ 0,73% giữa germanene (4,06 Å) và antimonene (4,09 Å), hệ Ge-Sb sở hữu năng lượng hình thành thấp nhất, chứng tỏ độ thuận lợi nhiệt động lực học vượt trội để lắng đọng thực nghiệm thông qua phương pháp chùm phân tử epitaxy.

Luận văn đã áp dụng các tiêu chí hội tụ nào để đảm bảo tính xác thực của phép tính DFT?

Công trình thiết lập năng lượng cắt 500 eV, mạng lưới điểm k Monkhorst-Pack 2x12x1 và phiếm hàm GGA-PBE trên phần mềm VASP. Tiêu chuẩn dừng tự hợp năng lượng đạt mức dưới 10^-6 eV và lực tác dụng Hellmann-Feynman trên mỗi nguyên tử nhỏ hơn 0,01 eV/Å, đảm bảo sai số hằng số mạng dưới 1,5% so với các dữ liệu thực nghiệm đã công bố.

Kết luận

  • Thiết lập thành công mô hình cấu trúc lượng tử cho 42 biến thể dị cấu trúc ngang phẳng giữa germanene và các đơn lớp pnictogen nhóm VA gồm arsenene, antimonene và AsSb.
  • Chứng minh tính tương thích mạng tối ưu giữa germanene và antimonene với độ sai lệch chỉ 0,73%, tạo điều kiện hoàn hảo cho quá trình tổng hợp thực nghiệm.
  • Kiểm soát linh hoạt độ rộng vùng cấm bán dẫn từ 0,03 eV lên tới đỉnh 0,87 eV ở cấu hình armchair thông qua biến thiên nồng độ germanene từ 12,5% đến 87,5%.
  • Phát hiện từ tính tự phát bền vững đạt 0,50 mu_B dọc theo đường biên zigzag nhờ sự đóng góp chủ đạo của các orbital Ge-pz chưa bão hòa.
  • Ghi nhận trạng thái bán dẫn Weyl đơn cực với vùng cấm siêu hẹp 19 meV đến 22,5 meV và hiệu ứng phân tách Rashba mở đường cho các ứng dụng tôpô lượng tử.

Trong giai đoạn 2024 - 2026, các nghiên cứu tiếp nối cần tập trung vào việc hiện thực hóa quy trình chế tạo epitaxy chùm phân tử và kiểm chứng các tính chất quang điện tử trong phòng thí nghiệm. Quý độc giả và các nhà khoa học quan tâm hãy kết nối, trích dẫn và tải toàn văn tài liệu luận văn để cùng hợp tác thúc đẩy sự phát triển của công nghệ vật liệu bán dẫn tiên tiến!