Tổng quan nghiên cứu

Kể từ khi được phân lập thành công vào năm 2004, Graphene đã trở thành tâm điểm của ngành vật lý vật liệu và công nghệ bán dẫn toàn cầu. Với cấu trúc phẳng hai chiều tạo bởi các nguyên tử Carbon lai hóa sp2, Graphene sở hữu những đặc tính vật lý vượt trội như độ dẫn nhiệt xấp xỉ 5000 W/m.K (cao gấp 12 lần đồng), suất Young đạt 1.0 TPa, độ bền kéo 125 GPa (gấp hơn 100 lần thép cứng nhất) và độ linh động hạt tải lý thuyết lên tới 200.000 cm2/Vs ở mật độ hạt tải 10^12 cm-2. Tuy nhiên, rào cản lớn nhất ngăn Graphene thay thế vật liệu Silicon truyền thống trong các linh kiện bán dẫn thế hệ mới là khe năng lượng nguyên thủy bằng 0 eV. Khi không có vùng cấm năng lượng, transistor làm từ Graphene không thể đóng ngắt dòng điện hiệu quả, khiến tỷ số dòng bật trên tắt chỉ dao động ở mức 100 đến 400, không đủ đáp ứng tiêu chuẩn hoạt động của vi mạch logic ở nhiệt độ phòng.

Nhằm vượt qua hạn chế này, luận văn thạc sĩ chuyên ngành Vật lý (mã số 60.04) của học viên Nguyễn Thị Len, dưới sự hướng dẫn khoa học của Tiến sĩ Nguyễn Tiến Cường tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội năm 2017, phối hợp cùng nhóm nghiên cứu tại Đại học Osaka (Nhật Bản), đã tập trung giải quyết bài toán cốt lõi: Thiết kế, mô phỏng và điều khiển khe năng lượng của Graphene thông qua cấu trúc lai Armchair – Zigzag. Mục tiêu cụ thể là khảo sát cấu trúc dải năng lượng, mật độ trạng thái điện tử và phổ truyền electron trên các cấu trúc hình học tiên tiến như kênh dẫn góc vuông 90 độ, dải gấp khúc, dải chữ U và cấu trúc đục lỗ, đồng thời đánh giá tác động của biến dạng cơ học lên các thông số điện tử. Nghiên cứu mang ý nghĩa then chốt khi chứng minh giải pháp mở rộng khe năng lượng từ 0.1 eV đến hơn 0.8 eV mà không cần thu hẹp kích thước kênh dẫn xuống dưới 5 nm hay sử dụng kỹ thuật pha tạp hóa học vốn làm suy giảm độ linh động của hạt tải điện.

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

Khung lý thuyết áp dụng

Khung lý thuyết của nghiên cứu được xây dựng trên hai trụ cột vững chắc của cơ học lượng tử tính toán: Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) và phương pháp hàm Green không cân bằng (Non-Equilibrium Green Function - NEGF).

Lý thuyết phiếm hàm mật độ xuất phát từ việc giải bài toán hệ nhiều hạt dựa trên gần đúng Born-Oppenheimer nhằm tách biệt chuyển động của electron và hạt nhân. Trọng tâm của DFT là hai định lý nền tảng của Hohenberg-Kohn (1964), khẳng định mật độ electron trạng thái cơ bản xác định duy nhất thế năng ngoài và toàn bộ tính chất năng lượng của hệ. Dựa trên đó, hệ phương trình Kohn-Sham (1965) đã chuyển đổi bài toán tương tác phức tạp của hệ electron sang hệ các hạt không tương tác chuyển động trong một thế hiệu dụng. Năng lượng tương quan trao đổi được mô tả chính xác thông qua gần đúng mật độ địa phương (LDA) và gần đúng gradient suy rộng (GGA).

Để mô tả hiện tượng truyền dẫn lượng tử qua cấu trúc mở, phương pháp hàm Green không cân bằng kết hợp kỹ thuật đường bao Keldysh được áp dụng. Phương pháp này liên kết chặt chẽ với DFT để xử lý bài toán vận chuyển electron qua hệ bán tuần hoàn gồm ba vùng: điện cực trái (Left), vùng tán xạ trung tâm (Center) và điện cực phải (Right). Các khái niệm cốt lõi như cấu trúc vùng năng lượng, mật độ trạng thái điện tử (DOS), hàm truyền lượng tử và orbital pz tạo liên kết pi phi định xứ với khoảng cách liên kết C-C chuẩn 0.142 nm và diện tích ô đơn vị 0.052 nm2 đóng vai trò then chốt trong việc phân tích các tính chất điện tử của mạng Graphene.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng lượng tử từ nguyên lý ban đầu (ab-initio) thông qua gói phần mềm mã nguồn mở OpenMX (Open source package for Material explorer) phiên bản tối ưu hóa cho máy tính song song MPI/OpenMP, kết hợp chương trình TranMain để tính toán động học truyền electron.

Cỡ mẫu mô phỏng bao gồm hơn 30 cấu hình dải Graphene nano đại diện với các thông số hình học đa dạng: Dải Zigzag (8-ZGNRs), chuỗi dải Armchair có độ rộng từ N = 6 đến N = 14 hàng nguyên tử, cấu trúc góc 90 độ, kênh dẫn gấp khúc và 3 dạng đục lỗ (lỗ tròn, lỗ vuông, lỗ tam giác) được khảo sát dưới các mức biến dạng kéo căng từ 1% đến 10%. Phương pháp chọn mẫu cấu hình được xây dựng có chủ đích nhằm phản ánh đầy đủ hai kiểu biên cơ bản (Armchair và Zigzag) cùng các dạng tiếp giáp lai hóa phổ biến trong thực nghiệm chế tạo vi cơ.

Lý do lựa chọn phương pháp DFT kết hợp NEGF là vì đây là chuẩn mực cao nhất hiện nay để dự đoán chính xác cấu trúc vi mô và tính chất truyền hạt tải ở quy mô dưới nanomet với sai số năng lượng dưới 0.05 eV mà không phụ thuộc vào các tham số giả định từ thực nghiệm. Quy trình phân tích tuân thủ nghiêm ngặt 3 bước: Tối ưu hóa hệ số ma trận cơ sở (basis set), xác định năng lượng cắt SCF trên 150 Rydberg và bán kính cutoff từ 5.0 đến 7.0 a.u.; tính toán cấu trúc dải cho hai kênh dẫn; giải hệ phương trình NEGF tự hội tụ trên vùng trung tâm để xuất phổ truyền electron. Toàn bộ quá trình tính toán và kiểm chuẩn mô hình được thực hiện liên tục trong khoảng thời gian 12 tháng.

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: Dải Graphene biên Zigzag thuần túy (8-ZGNRs) luôn duy trì tính chất kim loại với khe năng lượng bằng 0 eV và mật độ trạng thái đạt cực đại tại mức Fermi. Ngược lại, dải biên Armchair (N-AGNRs) thể hiện sự phân hóa rõ rệt: Các dải có số hàng nguyên tử N = 3p hoặc N = 3p+1 mang tính bán dẫn với khe năng lượng mở rộng từ 0.35 eV đến 1.20 eV tùy theo bề rộng dải, trong khi dải N = 3p+2 thể hiện tính chất bán kim loại với khe năng lượng xấp xỉ 0 eV.

Phát hiện thứ hai: Khi tích hợp hai cấu trúc biên vào dạng lai hóa góc vuông 90 độ và dải chữ U, tính đối xứng của mạng tinh thể bị phá vỡ, tạo ra rào thế tán xạ lượng tử mạnh mẽ tại nút giao tiếp. Cấu trúc lai này đã mở thành công khe năng lượng dao động từ 0.25 eV đến 0.78 eV mà không cần thu hẹp kích thước dải xuống dưới ngưỡng giới hạn công nghệ 5 nm, cho phép nâng tỷ số chuyển mạch dòng điện bật trên tắt lên hơn 1000 lần so với dải phẳng thông thường.

Phát hiện thứ ba: Cấu trúc đục lỗ nano (nanoperforations) làm thay đổi căn bản phổ truyền electron của dải Graphene. Với cùng diện tích lỗ bị khuyết tật, cấu trúc đục lỗ hình tam giác và hình vuông tạo ra sự giam hãm lượng tử cục bộ mạnh hơn, giúp mở rộng vùng cấm năng lượng cao hơn từ 40% đến 65% so với dạng lỗ tròn.

Phát hiện thứ tư: Biến dạng cơ học đơn trục tác động trực tiếp và có tính quy luật lên cấu trúc dải năng lượng. Khi đặt lực kéo căng từ 2% đến 8%, độ rộng khe năng lượng của hệ 7-AGNRs và hệ góc vuông biến thiên tuyến tính với tốc độ điều biến đạt khoảng 0.06 eV trên mỗi 1% độ biến dạng.

Thảo luận kết quả

Nguyên nhân vật lý sâu xa của hiện tượng mở khe năng lượng trong cấu trúc lai Armchair – Zigzag là do sự gián đoạn của liên kết pi giải tỏa và sự thay đổi góc liên kết giữa các orbital lai hóa sp2 tại khu vực chuyển tiếp. Khi dòng electron di chuyển qua khúc quanh 90 độ hoặc vùng biên lai hóa, sự lệch pha của hàm sóng phản xạ và truyền qua triệt tiêu các trạng thái điện tử ở lân cận mức năng lượng Fermi, dẫn đến sự hình thành vùng cấm năng lượng thực sự.

So với phương pháp pha tạp hóa học truyền thống (như gắn các nhóm H, F, Cl hoặc phân tử BN) thường làm biến dạng mạng tinh thể và khiến độ linh động hạt tải sụt giảm hơn 60%, giải pháp tạo cấu trúc lai thuần Carbon giúp bảo toàn liên kết C-C bền vững và duy trì độ linh động ở mức cao trên 10.000 cm2/Vs. Kết quả mô phỏng cấu trúc chữ U hoàn toàn tương thích với các quan sát thực nghiệm chế tạo bằng chùm ion hội tụ (FIB), trong đó transistor đạt tỷ số bật tắt lên tới 10^5, vượt trội gấp nhiều lần mức 400 của transistor Graphene hai lớp phân cực điện trường ngoài.

Dữ liệu nghiên cứu được biểu diễn trực quan và khoa học thông qua các biểu đồ phân tán năng lượng E(k), đồ thị mật độ trạng thái DOS, phổ hàm truyền T(E) dọc theo trục năng lượng và bảng ma trận tham số tối ưu hóa năng lượng cutoff, giúp làm sáng tỏ toàn diện cơ chế điều khiển dòng electron trong kênh dẫn nano.

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

Thứ nhất: Triển khai ứng dụng công nghệ quang khắc chùm ion hội tụ (Focused Ion Beam - FIB) để chế tạo thử nghiệm các kênh dẫn lai Armchair - Zigzag dạng chữ U và góc 90 độ với độ phân giải biên dưới 2 nm. Mục tiêu là chế tạo thành công transistor hiệu ứng trường Graphene đạt tỷ số dòng điện bật trên tắt vượt mốc 10^5 trước năm 2028. Chủ thể thực hiện là các phòng thí nghiệm quốc gia về công nghệ nano phối hợp cùng các viện nghiên cứu vật lý chất rắn.

Thứ hai: Tích hợp kỹ thuật kỹ thuật biến dạng cơ học (strain engineering) lên các đế vi cơ dẻo như PDMS hoặc PET. Áp dụng mức biến dạng kéo căng định hướng từ 3% đến 5% nhằm tinh chỉnh động khe năng lượng trong khoảng 0.3 eV đến 0.6 eV, hoàn thiện module cảm biến áp điện siêu nhạy trong giai đoạn 2027 - 2029. Chủ thể chịu trách nhiệm là các trung tâm nghiên cứu và phát triển linh kiện MEMS/NEMS.

Thứ ba: Nâng cấp nền tảng tính toán lượng tử bằng cách tích hợp các thuật toán học máy (Machine Learning) với gói phần mềm OpenMX để tự động hóa quá trình sàng lọc hơn 10.000 cấu hình lai hình học phức tạp. Giải pháp này giúp rút ngắn 70% thời gian tính toán mô phỏng trong vòng 18 tháng, giao cho các nhóm nghiên cứu Vật lý tính toán và Tin học lượng tử đảm nhiệm.

Thứ tư: Chuẩn hóa quy trình xử lý biên dải Graphene bằng phương pháp thụ động hóa Hydro (H-passivation) ở cấp độ nguyên tử nhằm loại bỏ triệt để các bẫy điện tích biên, giảm suy hao tán xạ điện tử xuống dưới mức 5% trước năm 2028. Nhiệm vụ này do các kỹ sư công nghệ vật liệu bề mặt và bán dẫn thực hiện.

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

Nhóm thứ nhất: Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý chất rắn, Khoa học vật liệu và Vật lý tính toán. Luận văn cung cấp tài liệu tham khảo chuẩn mực về phương pháp kết hợp DFT và NEGF, quy trình thiết lập tham số trong OpenMX và phương pháp xử lý dữ liệu truyền dẫn lượng tử ab-initio.

Nhóm thứ hai: Kỹ sư thiết kế linh kiện bán dẫn và mạch tích hợp nano. Tài liệu là cẩm nang hữu ích để hiểu rõ cơ chế mở khe năng lượng không cần pha tạp, hỗ trợ nghiên cứu phát triển các thế hệ transistor hiệu ứng trường FET thay thế Silicon ở các tiến trình dưới 2 nm.

Nhóm thứ ba: Chuyên gia nghiên cứu cảm biến nano và thiết bị quang điện tử dẻo. Những phát hiện về sự tương tác giữa cấu trúc đục lỗ, biến dạng cơ học và độ nhạy bề mặt 2630 m2/g cung cấp nền tảng vững chắc để thiết kế các cảm biến khí, cảm biến sinh học và màn hình OLED uốn dẻo thế hệ mới.

Nhóm thứ tư: Nhà quản lý công nghệ, giám đốc R&D tại các tập đoàn vi điện tử và các quỹ đầu tư khoa học công nghệ. Luận văn cung cấp bức tranh toàn cảnh về tính khả thi, lộ trình kỹ thuật và tiềm năng thương mại hóa của vật liệu Graphene 2D trong ngành công nghiệp vi mạch giai đoạn 5 đến 10 năm tới.

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

Câu hỏi 1: Tại sao Graphene nguyên thủy có độ dẫn điện cực cao nhưng không thể dùng trực tiếp làm transistor logic? Trả lời: Graphene nguyên thủy là vật liệu bán kim loại có khe năng lượng bằng 0 eV, khiến dòng electron liên tục di chuyển qua vật liệu mà không thể ngắt hoàn toàn. Khi ứng dụng vào transistor logic, linh kiện không thể đạt được trạng thái tắt hoàn chỉnh, dẫn đến dòng rò lớn và tỷ số bật trên tắt rất thấp, gây tiêu hao năng lượng nghiêm trọng.

Câu hỏi 2: Cấu trúc lai Armchair – Zigzag mở khe năng lượng cho Graphene bằng cơ chế nào? Trả lời: Cấu trúc lai kết hợp hai dạng biên khác nhau (như khúc quanh 90 độ hoặc dải chữ U) làm phá vỡ tính đối xứng tịnh tiến của mạng tổ ong. Sự gián đoạn liên kết pi tạo ra hiện tượng giao thoa sóng lượng tử và bẫy thế cục bộ, giúp hình thành vùng cấm năng lượng từ 0.25 eV đến 0.78 eV mà vẫn giữ nguyên độ linh động cao.

Câu hỏi 3: Sự kết hợp giữa DFT và hàm Green không cân bằng (NEGF) trong OpenMX mang lại lợi ích gì? Trả lời: DFT đóng vai trò tính toán chính xác thế hiệu dụng và cấu trúc điện tử trạng thái cơ bản từ nguyên lý ban đầu, trong khi NEGF mô tả thống kê lượng tử không cân bằng của dòng electron qua các điện cực bán vô hạn. Sự kết hợp này mô phỏng chân thực quá trình truyền dẫn electron qua kênh dẫn nano với độ chính xác cao dưới 0.05 eV.

Câu hỏi 4: Biến dạng cơ học ảnh hưởng ra sao đến tính chất điện tử của dải Graphene nano? Trả lời: Biến dạng cơ học làm thay đổi khoảng cách liên kết C-C chuẩn 0.142 nm và độ xen phủ của các orbital pz. Khi chịu lực kéo căng từ 2% đến 8%, độ rộng khe năng lượng thay đổi tuyến tính khoảng 0.06 eV trên mỗi 1% biến dạng, tạo cơ sở chế tạo các linh kiện cảm biến áp điện có thể điều chỉnh dải hoạt động.

Câu hỏi 5: So với phương pháp pha tạp hóa học, giải pháp tạo cấu trúc lai có ưu thế vượt trội nào? Trả lời: Phương pháp pha tạp hóa học đưa nguyên tử lạ vào mạng Carbon thường phá hủy cấu trúc tinh thể nguyên thủy và làm giảm hơn 60% độ linh động của hạt tải. Trong khi đó, cấu trúc lai thuần túy sử dụng liên kết C-C đồng nhất, giúp mở khe năng lượng hiệu quả đồng thời duy trì độ bền cơ học 125 GPa và độ linh động hạt tải trên 10.000 cm2/Vs.

Kết luận

  • Luận văn đã chứng minh thành công cơ chế mở và điều khiển khe năng lượng của Graphene trong dải từ 0.1 eV đến hơn 0.8 eV thông qua việc thiết kế các cấu trúc lai Armchair – Zigzag như dạng góc 90 độ, dạng chữ U và đục lỗ nano.
  • Khẳng định tính ưu việt của phương pháp tiếp cận lượng tử từ nguyên lý ban đầu thông qua sự kết hợp giữa lý thuyết phiếm hàm mật độ (DFT) và hàm Green không cân bằng (NEGF) trên phần mềm OpenMX.
  • Làm sáng tỏ quy luật biến đổi cấu trúc dải năng lượng dưới tác động của biến dạng cơ học đơn trục từ 1% đến 10%, mở ra hướng tinh chỉnh khe năng lượng chủ động.
  • Cung cấp cơ sở khoa học vững chắc giải thích tỷ số đóng ngắt dòng điện vượt trội lên tới 10^5 của kênh dẫn Graphene chữ U trong các thử nghiệm thực nghiệm chế tạo bằng chùm ion hội tụ.
  • Thiết lập hệ thống thông số mô phỏng chuẩn xác bao gồm bán kính cutoff, năng lượng cắt SCF và cấu hình hình học tối ưu cho các nghiên cứu vật liệu nano 2D tiếp theo.

Đóng góp lớn nhất của luận văn là đã giải quyết trọn vẹn điểm nghẽn vật lý về khe năng lượng của Graphene mà không cần đánh đổi độ linh động hạt tải hay độ bền cơ học đặc trưng. Lộ trình nghiên cứu tiếp theo cần tập trung hoàn thiện quy trình chế tạo mẫu thực nghiệm bằng công nghệ FIB và tích hợp mạch tích hợp trong vòng 24 tháng tới. Các viện nghiên cứu, trường đại học và doanh nghiệp công nghệ vi mạch bán dẫn nên đẩy mạnh ứng dụng mô hình tính toán này để sớm hiện thực hóa các dòng chip Graphene tốc độ cao thế hệ mới.