Tổng quan về luận án

Nghiên cứu về vật liệu hai chiều (2D) khởi sắc mạnh mẽ từ khi graphene được cô lập thực nghiệm, mở ra kỷ nguyên của các hạt tải Dirac fermion phi khối lượng với vận tốc Fermi đạt xấp xỉ $v_F \approx 10^6\text{ m/s}$ và tán sắc tuyến tính dạng nón $E(\mathbf{q}) = \pm \hbar v_F |\mathbf{q}|$. Mặc dù sở hữu độ linh động điện tử vượt trội, graphene nguyên bản thiếu một vùng cấm năng lượng ($E_g = 0\text{ eV}$), cản trở trực tiếp khả năng ứng dụng trong các linh kiện logic số (digital electronics) đòi hỏi tỷ số đóng ngắt ($I_{on}/I_{off}$) cao. Đồng thời, do mật độ trạng thái (DOS) triệt tiêu tuyến tính tiệm cận 0 tại điểm trung hòa điện ($\text{DOS}(E) \propto |E|$), graphene có vị trí mức Fermi $E_F(n) = \hbar v_F \sqrt{\pi n}$ vô cùng nhạy bén với nồng độ hạt tải $n$, cho phép biến điệu công thoát thuận lợi hơn kim loại 3D truyền thống. Luận án tiến sĩ chuyên ngành Khoa học Vật liệu (Mã số: 62 44 01 22) của tác giả Vũ Hoàng Nam tại Viện Công nghệ Nano - Trường Đại học Khoa học Tự nhiên, ĐHQG-TP.HCM, dưới sự hướng dẫn khoa học của PGS.TS Cao Minh Thì và TS. Lê Minh Hưng, đã giải quyết thấu đáo hai nghịch lý cốt lõi: "(i) chế tác graphene thành bán dẫn bằng cách tái cấu trúc vùng năng lượng của nó để mở một vùng cấm năng lượng đáng kể tại điểm Dirac; và (ii) ngược lại, bảo vệ tính chất bán kim loại và khả năng thay đổi công thoát của graphene để tăng sự nhạy của hàng rào tiếp giáp (hàng rào Schottky) giữa graphene và bán dẫn."

                       ┌────────────────────────────────────────────────────────┐
                       │  TƯƠNG TÁC GRAPHENE & ĐẾ BÁN DẪN (TÍNH TOÁN DFT)       │
                       └───────────────────┬────────────────────────────────────┘
                                           │
         ┌─────────────────────────────────┴─────────────────────────────────┐
         ▼                                                                   ▼
┌─────────────────────────────────┐                         ┌─────────────────────────────────┐
│ VẤN ĐỀ 1: KỸ THUẬT VÙNG NĂNG LƯỢNG│                         │ VẤN ĐỀ 2: TIẾP GIÁP SCHOTTKY    │
│  - Bẻ gãy đối xứng Sublattice   │                         │  - Tiếp giáp Graphene / Si(111) │
│  - Bẻ gãy đối xứng Chiral       │                         │  - Khử ghim Fermi (vdW vs MIGS) │
│  - Tán xạ ngoại hồ (K-K')       │                         │  - Dipole tiếp giáp & Push-back │
│  - Siêu mạng Kekulé & SiC       │                         │  - Độ nhạy hàng rào S -> 1      │
└─────────────────────────────────┘                         └─────────────────────────────────┘

Khoảng trống học thuật (research gap) được định vị cụ thể: các nghiên cứu trước đây (Giovannetti et al., 2007; Hunt et al., 2013; Ponomarenko et al., 2011) chủ yếu tập trung vào việc bẻ gãy đối xứng mạng phụ trên siêu mạng tuần hoàn tầm xa ($L > 10\text{ nm}$) dưới trường thế trơn mịn hoặc tiếp xúc van der Waals yếu với boron nitride lục giác (hBN), nhưng chưa giải thích được bản chất vật lý của sự trồi lên và tán xạ giữa các nón Dirac không tương đương ở siêu mạng tuần hoàn tầm ngắn ($L \sim 1 - 3\text{ nm}$) chịu thế ngoài sắc nhọn ở thang nguyên tử. Đồng thời, cơ chế truyền điện tích và vi mô dipole tại mặt phân giới tiếp giáp Schottky 2D/3D (graphene/silicon) vẫn bị quy giản cơ học theo các mô hình tiếp giáp kim loại/bán dẫn 3D cổ điển mà bỏ qua biến điệu mức Fermi nội tại của graphene.

Các câu hỏi nghiên cứu và giả thuyết tương ứng được thiết lập chặt chẽ:

  1. Câu hỏi 1 (Q1): Cơ chế lượng tử nào chi phối sự trồi lên, trộn lẫn và tán xạ giữa các nón Dirac sơ cấp (PDP) và thứ cấp (SDP) trong các siêu mạng graphene tầm ngắn trên đế bán dẫn silicon carbide (SiC)?
    • Giả thuyết 1 (H1): Một trường thế ngoài sắc nhọn từ các liên kết chưa bão hòa của đế SiC trên các siêu mạng đối xứng chọn lọc ($G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$) sẽ cung cấp vector truyền xung lượng $|\mathbf{Q}| \sim |\mathbf{K}| = 4\pi / (3\sqrt{3}a_0)$, bẻ gãy đối xứng chiral thông qua tán xạ ngược ngoại hồ (intervalley scattering) để mở vùng cấm năng lượng thực sự tại PDP.
  2. Câu hỏi 2 (Q2): Bản chất của bậc thế tiếp giáp và sự truyền điện tích qua mặt phân giới graphene/silicon (G/Si) tuân theo mô hình trạng thái bề mặt, MIGS, phân cực liên kết hay tiếp xúc van der Waals thuần túy?
    • Giả thuyết 2 (H2): Tiếp giáp G/Si đạt được liên kết van der Waals ổn định khi bề mặt silicon được bão hòa hóa trị hoặc tái cấu trúc thích hợp, qua đó triệt tiêu các trạng thái vùng cấm cảm ứng kim loại (MIGS) và hiện tượng ghim mức Fermi, đưa độ nhạy hàng rào Schottky $S = \partial \Phi_{Bn} / \partial \phi_m$ tiệm cận giới hạn Schottky-Mott lý tưởng ($S \to 1$).

Khung lý thuyết của luận án tích hợp: Lý thuyết phiếm hàm mật độ (Hohenberg-Kohn, 1964; Kohn-Sham, 1965), Lý thuyết Hamiltonian liên kết mạnh (Tight-Binding Model) với ma trận nhảy bậc nhất $t \approx 2.7\text{ eV}$, Mô hình tiếp giáp Schottky-Mott, Mô hình trạng thái bề mặt Bardeen - Cowley - Sze, Lý thuyết trạng thái vùng cấm cảm ứng kim loại (MIGS - Heine, 1965; Tersoff, 1984), và Lý thuyết phân cực liên kết (Bond Polarization Theory - Tung, 1992, 2001). Phạm vi nghiên cứu bao quát các siêu mạng graphene tuần hoàn tầm ngắn ($G-\sqrt{3}\times\sqrt{3}-R30^\circ$, $G-5\sqrt{3}\times 5\sqrt{3}-R30^\circ$, $G-9\times 9$, $G-\sqrt{21}\times\sqrt{21}-R10.893^\circ$, $G-7\times 7$) trên đế SiC và các cấu trúc tiếp giáp $G-\sqrt{39}\times\sqrt{39}-R16.1^\circ/\text{Si}(111)-4\times 4$ với các kiểu pha tạp bề mặt đa dạng, cung cấp bức tranh toàn cảnh chính xác về mặt lượng tử cho vật lý tiếp xúc 2D/3D.

Literature Review và Positioning

Lịch sử nghiên cứu tiếp giáp kim loại/bán dẫn và điều khiển dải năng lượng graphene ghi nhận những dòng nghiên cứu chính với các tranh luận học thuật sâu sắc:

                            TIẾN TRÌNH LÝ THUYẾT TIẾP GIÁP & VÙNG NĂNG LƯỢNG
 1947 - 1965                    1984 - 2001                    2007 - 2013                    2018 - Nay
 ┌──────────────────────────┐   ┌──────────────────────────┐   ┌──────────────────────────┐   ┌──────────────────────────┐
 │ MÔ HÌNH 3D CỔ ĐIỂN       │   │ MÔ HÌNH MIGS & PHÂN CỰC  │   │ GRAPHENE & ĐẾ 2D         │   │ TIẾP GIÁP vdW 2D/3D      │
 │ - Schottky-Mott (S=1)    │──►│ - Tersoff (MIGS & CNL)   │──►│ - Giovannetti (G/hBN gap)│──►│ - Liu et al. (MoS2 S=0.96│
 │ - Bardeen (Trạng thái BM)│   │ - Tung (Bond Polarization│   │ - Ponomarenko (SDP)      │   │ - Luận án: Khử ghim Fermi│
 │ - Cowley-Sze (Dipole)    │   │   và ISR Dipole)         │   │ - Hunt et al. (Sublattice│   │   & bẻ gãy Chiral trên Si│
 └──────────────────────────┘   └──────────────────────────┘   └──────────────────────────┘   └──────────────────────────┘
  1. Dòng nghiên cứu mở vùng cấm graphene: Dựa trên phương trình Dirac tuyến tính cho 2D graphene $H = \hbar v_F (\tau_z \sigma_x k_x + \sigma_y k_y)$, vùng cấm có thể hình thành qua hai cơ chế tách biệt: bẻ gãy đối xứng nghịch đảo mạng phụ (trộn $\pi(K)-\pi^*(K)$) hoặc bẻ gãy đối xứng chiral (trộn góc $K-K'$). Giovannetti et al. (2007) và Hunt et al. (2013) đã chứng minh trên tiếp giáp graphene/hBN rằng độ lệch mạng $1.8%$ tạo ra thế tuần hoàn làm bất đối xứng hai nguyên tử Carbon A và B, mở ra khoảng trống năng lượng tại PDP. Tuy nhiên, Wallbank et al. (2013) và Mucha-Kruczyński et al. (2013) tranh luận rằng các thế tuần hoàn tầm xa chỉ tạo ra tán xạ nội hồ (intravalley), hình thành các điểm Dirac thứ cấp (SDP) mà không thể kích hoạt tán xạ ngược ngoại hồ (intervalley backscattering). Đối với mạng Kekulé, Farjam & Rafii-Tabar (2009) cùng Cheianov et al. (2009) phân biệt hai cấu trúc: Kekulé-O (thế tác động tại tâm vòng hexagonal làm biến dạng 6 liên kết C-C) và Kekulé-Y (thế tác động trực tiếp lên nguyên tử C). Chỉ có Kekulé-O mới tạo ra tán xạ ngoại hồ bẻ gãy đối xứng chiral tại PDP. Vấn đề học thuật bỏ ngỏ là chưa có nghiên cứu nào làm rõ hành vi của các siêu mạng lớn hơn Kekulé ($G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$ với $N > 1$) khi chịu thế ngoài sắc nhọn.
  2. Dòng nghiên cứu tiếp giáp Schottky và hiện tượng ghim mức Fermi (Fermi-level pinning): Trong các tiếp giáp M/S kim loại 3D truyền thống, độ nhạy của hàng rào Schottky $S = \partial \Phi_{Bn} / \partial \phi_m$ luôn suy giảm nghiêm trọng từ giới hạn Schottky-Mott ($S = 1$) về gần giới hạn Bardeen ($S \approx 0$). Cowley & Sze (1965) gán hiện tượng này cho mật độ bẫy trạng thái bề mặt cao ($D_{it}$), trong khi Heine (1965) và Tersoff (1984) phát triển lý thuyết MIGS, chứng minh phần đuôi hàm sóng kim loại suy giảm theo dạng hàm mũ phức $\kappa = -iq$ thâm nhập vào vùng cấm bán dẫn, ghim mức Fermi vào điểm quay lui (branching point $E_B$ hay Charge Neutrality Level $\phi_{CNL}$): $$E(q) = V_0 + E_0 + \left(\frac{\hbar^2 q^2}{2m_0}\right) \pm \left[V_1^2 - 4E_0\left(\frac{\hbar^2 q^2}{2m_0}\right)\right]^{1/2}$$ Tung (2001) đưa ra lý thuyết phân cực liên kết, khẳng định dipole hình thành từ các liên kết hóa học trong vùng tiếp giáp riêng (ISR) quyết định $S = [1 + e^2 N_{ad} d_{in} / (\epsilon_{in}(E_g + \kappa))]^{-1}$.

Vị thế học thuật của luận án được xác lập rõ ràng khi đối chiếu với 2 nghiên cứu quốc tế mang tính đột phá:

  • So với công trình của Liu et al. (Nature, 2018) về tiếp xúc van der Waals giữa kim loại 2D và $\text{MoS}_2$ đạt độ nhạy $S = 0.96$ (vượt trội hoàn toàn so với tiếp xúc hóa học có $S = 0.09$), luận án mở rộng biên độ từ tiếp giáp 2D/2D sang tiếp giáp 2D/3D (graphene/silicon), giải mã cách thức bảo toàn độ nhạy $S \to 1$ ngay cả khi xuất hiện sự truyền điện tích thực.
  • So với nghiên cứu của LaGasse et al. (2019) trên tiếp giáp graphene/$\text{WSe}_2$ đạt $S = 1.04$ nhờ dịch dải năng lượng bằng điện trường ngoài cực cổng, luận án của NCS. Vũ Hoàng Nam làm sáng tỏ bản chất vi mô: biến thiên hàng rào $\Delta \Phi_B$ thực chất được dẫn dắt bởi sự dịch mức Fermi nội tại $\Delta E_F$ của graphene do mật độ trạng thái tiệm cận 0, độc lập với việc hình thành dipole tiếp giáp kiểu 3D cổ điển.

Đóng góp lý thuyết và khung phân tích

                   CƠ CHẾ LƯỢNG TỬ ĐIỀU BIẾN TÍNH CHẤT ĐIỆN TỬ GRAPHENE
┌───────────────────────────────────────────────────────────────────────────────────────┐
│ 1. BẺ GÃY ĐỐI XỨNG MẠNG PHỤ (Sublattice / Inversion Symmetry Breaking)                │
│    - Tương tác bất đối xứng giữa mạng con A và B: H_AA != H_BB                        │
│    - Trộn trạng thái pi(K) và pi*(K) -> Mở vùng cấm trực tiếp Eg tại PDP             │
├───────────────────────────────────────────────────────────────────────────────────────┤
│ 2. BẺ GÃY ĐỐI XỨNG CHIRAL (Chiral Symmetry Breaking via Intervalley Scattering)       │
│    - Vector xung lượng |Q| ~ |K| từ thế ngoài sắc nhọn đế bán dẫn                     │
│    - Gấp vùng Brillouin SBZ -> Trộn lẫn 2 nón Dirac không tương đương K và K' tại gamma │
│    - Triệt tiêu xuyên ngầm Klein -> Mở vùng cấm tại PDP và SDP                        │
├───────────────────────────────────────────────────────────────────────────────────────┤
│ 3. KHỬ GHIM FERMI TRONG TIẾP GIÁP vdW 2D/3D (De-pinning Schottky Barrier)             │
│    - Loại bỏ liên kết hóa học cục bộ & trạng thái bề mặt bẫy điện tích (D_it -> 0)    │
│    - Triệt tiêu trạng thái MIGS trong vùng cấm bán dẫn 3D                             │
│    - Điều hưởng hàng rào tuyến tính: Delta Phi_B ~ Delta E_F (Đạt giới hạn S -> 1)    │
└───────────────────────────────────────────────────────────────────────────────────────┘

Đóng góp cho lý thuyết

Luận án đã mở rộng và tái định hình sâu sắc các lý thuyết vật lý chất rắn nền tảng:

  1. Mở rộng Lý thuyết Tán sắc Dirac và Đối xứng Chiral: Luận án xây dựng mô hình toán học giải thích hiện tượng tán xạ ngoại hồ trong các siêu mạng graphene tuần hoàn tầm ngắn. Khi siêu mạng thỏa mãn điều kiện đối xứng hình học $G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$, tâm của vùng Brillouin siêu mạng (SBZ - điểm $\gamma$) trùng khít hoàn toàn với các góc $K$ và $K'$ của vùng Brillouin nguyên tố (GBZ). Thế ngoài chu kỳ thang nguyên tử đóng vai trò như một toán tử nhiễu loạn cung cấp vector chuyển xung lượng $|\mathbf{Q}| \approx 4\pi / (3\sqrt{3}a_0)$, kích hoạt sự giao thoa lượng tử giữa hai nón Dirac có giả spin đối song ($\mathbf{s}K$ và $\mathbf{s}{K'}$). Hiện tượng này bẻ gãy đối xứng chiral, loại bỏ hoàn toàn nghịch lý xuyên ngầm Klein (Klein tunneling), mở ra vùng cấm năng lượng $E_g$ xác định tại PDP.
  2. Xác lập Mô hình Tiếp giáp Schottky 2D/3D Phi Cổ Điển: Thách thức trực tiếp mô hình Cowley - Sze và lý thuyết MIGS của Tersoff khi áp dụng cho vật liệu có mật độ trạng thái thấp. Luận án chứng minh rằng trong tiếp giáp van der Waals graphene/bán dẫn, công thức hàng rào Schottky: $$\Phi_{Bn} = \phi_m - \chi_s - \Delta_{in}$$ không bị chi phối bởi dipole bề mặt cố định $\Delta_{in} = -(e Q_m d_{in})/\epsilon_{in}$. Thay vào đó, do DOS của graphene triệt tiêu tại lân cận điểm Dirac, mức Fermi bị dịch chuyển mạnh một khoảng $\Delta E_F = \hbar v_F \sqrt{\pi n}$. Khi đó, sự biến điệu của hàng rào Schottky tuân theo hệ thức đột phá: $$\Delta \Phi_{Bn} \approx \Delta E_F$$ chứng minh trên phương diện lượng tử rằng tiếp giáp vdW 2D/3D hoàn toàn có thể đạt tới giới hạn Schottky-Mott lý tưởng ($S = 1$) mà không bị ghim mức Fermi.

Khung phân tích độc đáo

Khung phân tích của công trình được kiến tạo từ sự tích hợp liên ngành giữa 3 trụ cột lý thuyết: Cơ học lượng tử cấu trúc dải (Quantum Band Theory), Điện động lực học bề mặt mặt phân giới (Interface Electrostatics), và Nhiệt động học tính toán phiếm hàm mật độ.

Các đóng góp khái niệm mới được định nghĩa tường minh:

  • Thế cảm ứng phân cực đảo (Induced Reverse Polarization Potential): Sự phân bố lại mật độ điện tích phẳng $\Delta \rho(z)$ tại mặt phân giới do hiệu ứng đẩy ngược Pauli (push-back / pillow effect) nén phần đuôi hàm sóng điện tử, tạo ra bước nhảy thế tĩnh điện đối kháng với thế truyền điện tích.
  • Hiện tượng trồi nón Dirac (Dirac Cone Emergence & Unfolding): Quá trình các dải năng lượng bị gấp (folded bands) trong vùng SBZ được chiếu và trải ngược về vùng GBZ thông qua toán tử chiếu hàm sóng Bloch, làm hiển lộ nguồn gốc vật lý thực của các trạng thái phân bố tại các điểm Dirac thứ cấp (SDP).
  • Điều kiện biên (Boundary Conditions): Cơ chế bảo toàn tính chất bán kim loại chỉ duy trì nghiêm ngặt khi khoảng cách tiếp giáp $d_{vdW} \ge 3.0\text{ \AA}$ và bề mặt chất nền silicon được thụ động hóa hoàn toàn các liên kết treo (dangling bonds) bằng hydro hoặc tái cấu trúc bề mặt có trật tự cao.

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

                        QUY TRÌNH TÍNH TOÁN DFT & XỬ LÝ DỮ LIỆU LƯỢNG TỬ
┌─────────────────────────┐     ┌─────────────────────────┐     ┌─────────────────────────┐
│ KHỞI TẠO CẤU TRÚC       │     │ TỐI ƯU HÓA HÌNH HỌC     │     │ TÍNH TOÁN ĐIỆN TỬ SCF   │
│ - Mô hình phiếm hàm vdW │────►│ - Phần mềm VASP         │────►│ - Năng lượng cắt Encut  │
│ - Ghép siêu mạng G/SiC  │     │ - Lực ion < 0.01 eV/Å   │     │ - Lưới k-points dày đặc │
│ - Siêu mạng G/Si(111)   │     │ - Tiêu chuẩn E < 10^-6eV│     │ - Hiệu chỉnh optB86b-vdW│
└─────────────────────────┘     └─────────────────────────┘     └─────────────────────────┘
                                                                             │
                                ┌────────────────────────────────────────────┘
                                ▼
┌─────────────────────────────────────────────────────────────────────────────────────────┐
│ KỸ THUẬT PHÂN TÍCH CHUYÊN SÂU                                                           │
│  - Kỹ thuật trải dải năng lượng (Band Unfolding) từ SBZ về GBZ nguyên tố                │
│  - Phân tích hiệu mật độ điện tích trung bình phẳng: Delta rho(z)                       │
│  - Chiết xuất thế tĩnh điện cảm ứng Hartre & Dipole tiếp giáp Delta_in                  │
│  - Tính toán độ cao hàng rào Schottky (SBH) và tham số độ nhạy S = d(Phi_B)/d(phi_m)    │
└─────────────────────────────────────────────────────────────────────────────────────────┘

Thiết kế nghiên cứu

Nghiên cứu được thiết kế trên nền tảng triết lý khoa học thực chứng lượng tử (Quantum Positivism) và hiện thực phản tư (Critical Realism). Phép toán được thực hiện thông qua mô hình đa quy mô (multi-level design) kết hợp giữa cấu trúc mạng nguyên tử vi mô và giản đồ dải năng lượng vĩ mô. Không gian siêu mạng được mô hình hóa bằng kỹ thuật siêu ô tuần hoàn (supercell slab model) với lớp chân không dày hơn $15\text{ \AA}$ nhằm loại bỏ hoàn toàn tương tác giả giữa các lớp định kỳ theo phương trục $z$.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng lượng tử tuân thủ nghiêm ngặt các tiêu chuẩn tính toán chất rắn:

  1. Thuật toán giải phương trình Kohn-Sham tự hợp (Self-Consistent Field - SCF): Sử dụng gói phần mềm chuyên dụng VASP (Vienna Ab-initio Simulation Package). Tương tác trao đổi - tương quan điện tử được xử lý bằng phép xấp xỉ Gradient tổng quát (GGA-PBE), phép xấp xỉ mật độ địa phương (LDA), và đặc biệt là phiếm hàm hiệu chỉnh lực phân tán van der Waals không cục bộ optB86b-vdW nhằm mô tả chính xác tương tác liên lớp yếu.
  2. Phương pháp giả thế (Pseudopotential Method): Sử dụng sóng phẳng tăng cường xạ ảnh (Projector Augmented Wave - PAW). Năng lượng cắt sóng phẳng (cutoff energy) được thiết lập tối ưu tại mức $E_{cut} = 400 - 500\text{ eV}$.
  3. Lấy mẫu không gian mạng đảo (k-point sampling): Lưới k-points theo sơ đồ Monkhorst-Pack được phân bố siêu mịn (từ $9\times 9\times 1$ đến $15\times 15\times 1$ cho các cấu trúc siêu mạng lớn). Tiêu chí hội tụ năng lượng điện tử đạt mức khắt khe $10^{-6}\text{ eV}$, và lực tác dụng lên từng nguyên tử sau tối ưu hóa hình học nhỏ hơn $0.01\text{ eV/\AA}$.
  4. Kỹ thuật trải dải năng lượng (Band Unfolding Technique): Sử dụng thuật toán chiếu véc-tơ hàm sóng riêng từ vùng Brillouin siêu mạng (SBZ) về vùng Brillouin nguyên tố (GBZ), cho phép bóc tách cấu trúc dải phức tạp, nhận diện trực quan vị trí trồi lên của các nón Dirac bản sao và quá trình mở vùng cấm năng lượng.

Data và phân tích

Luận án tiến hành phân tích định lượng trên hệ thống dữ liệu đồ sộ:

  • Khảo sát các cấu trúc siêu mạng tiếp giáp $G/\text{SiC}$ với các thông số đối xứng: $G-5\sqrt{3}\times 5\sqrt{3}-R30^\circ / \text{SiC}-4\sqrt{3}\times 4\sqrt{3}-R30^\circ$, $G-\sqrt{21}\times\sqrt{21}-R10.893^\circ$, $G-7\times 7$, $G-9\times 9$.
  • Hệ tiếp giáp $G-\sqrt{39}\times\sqrt{39}-R16.1^\circ/\text{Si}(111)-4\times 4$ với độ lệch mạng (lattice mismatch) cực thấp ($< 1%$).
  • Phân tích hiệu mật độ điện tích tích hợp phẳng: $$\Delta \rho(z) = \rho_{G/Si}(z) - \rho_G(z) - \rho_{Si}(z)$$ và thế cảm ứng tĩnh điện $V_H(z)$ giải từ phương trình Poisson một chiều: $$\frac{d^2 V_H(z)}{dz^2} = -\frac{\Delta \rho(z)}{\epsilon_0}$$
  • Kiểm tra tính vững chắc (Robustness checks): So sánh đối chuẩn kết quả giữa ba phiếm hàm optB86b-vdW, PBE và LDA, chứng minh rằng các tương tác phân tán tầm xa đóng vai trò sống còn trong việc duy trì cấu trúc dải tuyến tính của graphene mà các phép xấp xỉ truyền thống không thể mô tả chuẩn xác.

Phát hiện đột phá và implications

                   CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN
┌──────────────────────────────────────────────────────────────────────────────┐
│ PHÁT HIỆN 1: ĐIỀU KIỆN HÌNH HỌC MỞ VÙNG CẤM TẠI ĐIỂM DIRAC SƠ CẤP (PDP)      │
│ - Chỉ các siêu mạng dạng G-sqrt(3N)xsqrt(3N)-R phi bẻ gãy đối xứng chiral    │
│ - Tán xạ ngược ngoại hồ (K-K') mở vùng cấm thực sự Eg                        │
├──────────────────────────────────────────────────────────────────────────────┤
│ PHÁT HIỆN 2: BẢN CHẤT VẬT LÝ CỦA ĐIỂM DIRAC THỨ CẤP (SDP)                    │
│ - SDP hình thành do thế tuần hoàn trơn mịn; PDP mở vùng cấm do thế sắc nhọn  │
│ - Nguồn gốc vật lý phân kỳ hoàn toàn giữa SDP và PDP                         │
├──────────────────────────────────────────────────────────────────────────────┤
│ PHÁT HIỆN 3: BẢO TOÀN LIÊN KẾT vdW TRÊN ĐẾ Si(111)                           │
│ - Trái ngược dự đoán cũ về liên kết cộng hóa trị gây gợn sóng graphene       │
│ - Đế Si(111)-4x4 tái cấu trúc duy trì graphene phẳng tuyệt đối (d ~ 3.3 Å)   │
├──────────────────────────────────────────────────────────────────────────────┤
│ PHÁT HIỆN 4: ĐẠT GIỚI HẠN SCHOTTKY-MOTT LÝ TƯỞNG (S -> 1)                    │
│ - Triệt tiêu trạng thái MIGS và bẫy điện tích bề mặt                         │
│ - Biến thiên hàng rào tỷ lệ thuận trực tiếp với độ dịch Fermi: Delta Phi ~ EF │
└──────────────────────────────────────────────────────────────────────────────┘

Những phát hiện then chốt

Luận án mang lại 4 phát hiện khoa học mang tính bước ngoặt:

  1. Quy luật hình học điều khiển tán xạ ngoại hồ tại PDP: Tính toán DFT chứng minh rằng không phải mọi siêu mạng graphene đều có khả năng mở vùng cấm khi tương tác với đế. Chỉ những cấu trúc siêu mạng có tính chất đối xứng $G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$ (như $G-5\sqrt{3}\times 5\sqrt{3}-R30^\circ$) mới làm điểm $\gamma$ của SBZ gấp đúng vào điểm $K(K')$ của GBZ. Dưới tác động của trường thế sắc nhọn từ các liên kết treo bề mặt đế $\text{SiC}$, các nón Dirac không tương đương bị trộn lẫn, kích hoạt tán xạ ngược ngoại hồ ($K-K'$) và bẻ gãy đối xứng chiral, hình thành vùng cấm năng lượng rõ rệt. Ngược lại, các siêu mạng như $G-7\times 7$ hay $G-9\times 9$ chủ yếu duy trì tán xạ nội hồ.
  2. Giải mã nguồn gốc vật lý phân kỳ giữa PDP và SDP: Trái ngược với các công bố quốc tế trước đây (ví dụ: Amet et al., 2013) cho rằng việc mở vùng cấm tại PDP và SDP có cùng chung bản chất bẻ gãy đối xứng mạng phụ, luận án chứng minh: vùng cấm tại PDP bắt nguồn từ sự bẻ gãy đối xứng chiral do tán xạ ngoại hồ thang nguyên tử, trong khi các điểm Dirac thứ cấp (SDP) xuất hiện thuần túy do sự gấp dải và tán xạ nội hồ dưới trường thế tuần hoàn dài, hoàn toàn có thể bị triệt tiêu vùng cấm bởi hiệu ứng xuyên ngầm Klein.
  3. Cơ chế liên kết van der Waals phẳng trên đế $\text{Si}(111)$: Bác bỏ quan điểm của các nghiên cứu DFT sơ khai (vốn cho rằng graphene trên silicon luôn hình thành các liên kết cộng hóa trị cục bộ $C-Si$ làm gợn sóng và phá hủy nón Dirac). Luận án khẳng định trên bề mặt $\text{Si}(111)-4\times 4$ tái cấu trúc hoặc bề mặt $\text{Si}(111)$ thụ động hóa Hydro ($\text{Si-H}$), khoảng cách tiếp giáp ổn định ở mức van der Waals ($d \approx 3.25 - 3.40\text{ \AA}$), bảo toàn nguyên vẹn cấu trúc điện tử Dirac và độ dốc tán sắc tuyến tính của graphene.
  4. Hiện tượng khử ghim mức Fermi và tiếp cận giới hạn Schottky-Mott ($S \approx 1$): Trong tiếp giáp $G/\text{Si}(111)$, do sự vắng mặt của các trạng thái vùng cấm cảm ứng kim loại (MIGS) và trạng thái bề mặt định xứ trong vùng cấm silicon, mức Fermi của graphene không bị ghim. Bậc thế tiếp giáp thực chất được phân tách thành hai thành phần đối kháng: thế hiệu ứng đẩy ngược ($\Delta_{pb}$) và thế truyền điện tích ($\Delta_{ct}$). Khi thay đổi công thoát đế thông qua pha tạp, sự dịch chuyển mức Fermi $\Delta E_F$ trực tiếp biến điệu chiều cao hàng rào Schottky với độ nhạy $S = \partial \Phi_{Bn} / \partial \phi_m \approx 0.95 - 1.00$, hiện thực hóa giới hạn Schottky-Mott lý tưởng trên tiếp giáp 2D/3D.

Implications đa chiều

  • Về mặt lý thuyết: Thống nhất lý thuyết tiếp giáp bề mặt bằng cách tích hợp trạng thái Dirac 2D vào hệ phương trình tĩnh điện Schottky-Mott, mở ra hướng tiếp cận chuẩn xác cho các hệ dị thể 2D/3D (2D/3D heterostructures).
  • Về mặt phương pháp luận: Kỹ thuật trải dải năng lượng (band unfolding) kết hợp phiếm hàm hiệu chỉnh phân tán optB86b-vdW thiết lập chuẩn mực tính toán mới cho các hệ siêu mạng phức hợp có kích thước nguyên tử lớn.
  • Về mặt ứng dụng thực tiễn: Cung cấp nguyên lý thiết kế cho linh kiện Transistor rào chắn graphene (Graphene Barristor), cảm biến quang điện tử siêu nhạy và pin mặt trời dị thể graphene/silicon với hiệu suất thu gom hạt tải tối ưu nhờ kiểm soát hoàn hảo hàng rào Schottky.
  • Khuyến nghị chính sách R&D: Định hướng cho các phòng thí nghiệm quốc gia và doanh nghiệp bán dẫn đầu tư vào công nghệ epitaxy graphene trực tiếp trên đế silicon và SiC quy mô wafer, đón đầu công nghệ hậu silicon (Beyond-CMOS).

Limitations và Future Research

  1. Giới hạn quy mô tính toán nhiệt độ: Các tính toán DFT trong luận án được thực hiện tại trạng thái cơ bản ($T = 0\text{ K}$), chưa tính toán trực tiếp ảnh hưởng của dao động mạng (phonon-electron scattering) ở nhiệt độ phòng ($T = 300\text{ K}$) đến độ nhạy $S$ và độ mở vùng cấm.
  2. Giới hạn về khuyết tật ngẫu nhiên: Mô hình siêu mạng giả định cấu trúc tinh thể hoàn hảo không tì vết, trong khi thực tế màng graphene chế tạo bằng phương pháp CVD luôn tồn tại biên hạt (grain boundaries), khuyết tật điểm (Stone-Wales defects) và tạp chất môi trường.
  3. Chương trình nghiên cứu 5-10 năm tới:
    • Mở rộng tính toán DFT kết hợp phương trình GW và Bethe-Salpeter (GW-BSE) để đánh giá chuẩn xác hiệu ứng kích thích quang và exciton tại mặt phân giới tiếp giáp.
    • Nghiên cứu động lực học phân tử lượng tử (ab-initio Molecular Dynamics - AIMD) khảo sát độ bền nhiệt động học của tiếp giáp $G/\text{Si}$ ở điều kiện vận hành thực tế ($T = 300 - 500\text{ K}$).
    • Mở rộng khảo sát tiếp giáp van der Waals giữa graphene với các bán dẫn oxit năng lượng rộng ($\text{Ga}_2\text{O}_3$, $\text{TiO}_2$) và bán dẫn 2D thế hệ mới (MXenes, Phosphorene).

Tác động và ảnh hưởng

                                HỆ THỐNG TÁC ĐỘNG ĐA LĨNH VỰC
┌─────────────────────────┐     ┌─────────────────────────┐     ┌─────────────────────────┐
│ TÁC ĐỘNG HỌC THUẬT      │     │ CÔNG NGHIỆP BÁN DẪN     │     │ LỢI ÍCH XÃ HỘI & MÔI TRƯỜNG│
│ - Chuẩn hóa mô hình     │     │ - Tích hợp Graphene vào │     │ - Nền tảng chip tiêu thụ│
│   tiếp giáp vdW 2D/3D   │────►│   dây chuyền Si-Fab     │────►│   năng lượng siêu thấp  │
│ - Nâng cao chỉ số trích │     │ - Thương mại hóa dòng   │     │ - Thiết bị quang điện   │
│   dẫn ngành Vật liệu    │     │   chip Barristor logic  │     │   hiệu suất cao bền vững│
└─────────────────────────┘     └─────────────────────────┘     └─────────────────────────┘
  • Tác động học thuật: Định hình lại các mô hình giáo khoa kinh điển về tiếp xúc kim loại/bán dẫn, cung cấp cơ sở lý thuyết chuẩn xác cho hàng trăm công trình nghiên cứu tiếp nối về cấu trúc dải của vật liệu Dirac 2D.
  • Chuyển đổi công nghiệp: Giải tỏa nút thắt công nghệ trong việc tích hợp graphene vào dây chuyền sản xuất vi mạch silicon hiện hữu, mở đường cho việc chế tạo các linh kiện chuyển mạch tốc độ Terahertz ($10^{12}\text{ Hz}$) với điện áp vận hành cực thấp.
  • Lợi ích kinh tế - xã hội: Giảm thiểu tiêu hao năng lượng trong các trung tâm dữ liệu (data centers) toàn cầu nhờ các bộ vi xử lý dựa trên dị thể graphene/bán dẫn siêu mát, đóng góp trực tiếp vào mục tiêu phát triển công nghệ xanh và điện tử bền vững.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà nghiên cứu sau tiến sĩ (Postdocs): Thụ hưởng phương pháp luận tính toán dải năng lượng chi tiết, quy trình giải mã ma trận hàm sóng và kỹ thuật trải dải (band unfolding) trong VASP để áp dụng cho các hệ vật liệu mới.
  • Giáo sư & Chuyên gia Vật lý chất rắn: Sở hữu hệ thống bằng chứng lượng tử chuẩn xác để hoàn thiện lý thuyết tiếp xúc mặt phân giới 2D/3D và cơ chế tán xạ Dirac fermion.
  • Kỹ sư R&D trong công nghiệp bán dẫn: Nắm bắt các thông số kỹ thuật tối ưu (khoảng cách lớp vdW, cấu trúc tái cấu trúc bề mặt đế Si) để thiết kế quy trình lắng đọng màng mỏng và chế tạo linh kiện Barristor hiệu năng cao.
  • Nhà hoạch định chiến lược Khoa học & Công nghệ: Có luận cứ khoa học vững chắc để xây dựng các chương trình tài trợ trọng điểm quốc gia về công nghệ vật liệu bán dẫn tiên tiến và vi điện tử thế hệ mới.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
    • Trả lời: Đó là việc chứng minh sự bẻ gãy đối xứng chiral và mở vùng cấm tại điểm Dirac sơ cấp (PDP) của graphene thông qua cơ chế tán xạ ngược ngoại hồ ($K-K'$) dưới tác động của trường thế ngoài thang nguyên tử sắc nhọn trên các siêu mạng đối xứng chọn lọc $G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$, đồng thời phát hiện quy luật biến điệu hàng rào Schottky 2D/3D $\Delta \Phi_{Bn} \approx \Delta E_F$ đạt giới hạn Schottky-Mott lý tưởng ($S \to 1$).
  2. Đột phá về mặt phương pháp luận so với các nghiên cứu trước đây thể hiện ở điểm nào?
    • Trả lời: Khác với các nghiên cứu dùng mô hình tiếp xúc 3D của Cowley - Sze hay phép xấp xỉ thế trơn mịn (Mucha-Kruczyński et al.), luận án kết hợp phiếm hàm phân tán không cục bộ optB86b-vdW với thuật toán trải dải (band unfolding) từ SBZ về GBZ trên các siêu ô mạng nguyên tử lớn tiếp giáp SiC và $\text{Si}(111)$, cho phép quan sát vi mô trực tiếp quá trình phân kỳ và tái cấu trúc các nón Dirac bản sao.
  3. Phát hiện nào mang tính bất ngờ (counter-intuitive) nhất?
    • Trả lời: Sự phân kỳ về nguồn gốc vật lý giữa điểm Dirac sơ cấp (PDP) và điểm Dirac thứ cấp (SDP). Trong khi PDP đòi hỏi thế ngoài sắc nhọn để bẻ gãy đối xứng chiral, thì SDP lại sinh ra thuần túy do thế tuần hoàn trơn và chịu sự chi phối của hiệu ứng xuyên ngầm Klein. Đồng thời, sự truyền điện tích qua tiếp giáp vdW G/Si không hề làm giảm độ nhạy $S$ của hàng rào như lý thuyết tiếp xúc kim loại 3D truyền thống dự đoán.
  4. Quy trình nghiên cứu có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?
    • Trả lời: Hoàn toàn đầy đủ. Toàn bộ tham số mạng ban đầu, bán kính cắt sóng phẳng ($E_{cut} = 400 - 500\text{ eV}$), lưới tích phân không gian k, cấu hình giả thế PAW, tiêu chuẩn hội tụ lực ion ($< 0.01\text{ eV/\AA}$) và các bước thiết lập phiếm hàm optB86b-vdW trong VASP đều được tường minh hóa chi tiết.
  5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?
    • Trả lời: Mở rộng tính toán sang các hiệu ứng tương quan điện tử mạnh (Strongly Correlated Electron Systems), vật lý góc ma thuật (Magic-Angle Twisted Bilayer Graphene), mô phỏng truyền hạt tải lượng tử phụ thuộc thời gian (TD-NEGF) và tích hợp graphene trên các đế bán dẫn siêu rộng nhằm hiện thực hóa linh kiện công suất Terahertz.

Kết luận

Luận án tiến sĩ của NCS. Vũ Hoàng Nam đại diện cho một công trình khoa học mẫu mực, chuẩn xác và công phu thuộc chuyên ngành Khoa học Vật liệu, mang lại 6 đóng góp cốt lõi:

  1. Xác lập bản đồ điều khiển đối xứng chiral và cơ chế tán xạ ngoại hồ ($K-K'$) trên các siêu mạng graphene tuần hoàn tầm ngắn $G-\sqrt{3N}\times\sqrt{3N}-R\phi^\circ$.
  2. Chứng minh điều kiện hình học và vật lý để mở vùng cấm năng lượng thực chất $E_g$ tại điểm Dirac sơ cấp (PDP) bằng trường thế nguyên tử sắc nhọn từ đế SiC.
  3. Giải mã sự khác biệt bản chất giữa cơ chế mở vùng cấm tại PDP và sự xuất hiện của các điểm Dirac thứ cấp (SDP).
  4. Bác bỏ các định kiến về liên kết cộng hóa trị phá hủy màng, chứng minh trạng thái tiếp xúc van der Waals phẳng hoàn hảo của graphene trên bề mặt $\text{Si}(111)-4\times 4$.
  5. Khám phá cơ chế phân tách dipole tiếp giáp (thế đẩy ngược Pauli $\Delta_{pb}$ và thế truyền điện tích $\Delta_{ct}$), chứng minh hiện tượng giải phóng ghim mức Fermi trên tiếp xúc 2D/3D.
  6. Thiết lập mô hình giải tích cho độ nhạy hàng rào Schottky $\Delta \Phi_{Bn} \approx \Delta E_F$, đưa tiếp giáp graphene/silicon tiệm cận giới hạn Schottky-Mott lý tưởng ($S \to 1$).

Công trình không chỉ mở ra 3 hướng nghiên cứu mới về điện tử học thung lũng (Valleytronics), vật lý tiếp xúc dị thể 2D/3D và linh kiện logic vượt giới hạn CMOS, mà còn định hình chuẩn mực tính toán lượng tử chính xác cao cho ngành khoa học vật liệu nano trong nước và quốc tế.