Tổng quan về luận án

Nghiên cứu về tương tác giữa vật liệu hai chiều và đế bán dẫn đại diện cho một trong những tiền tuyến học thuật then chốt của vật lý chất rắn và khoa học vật liệu nano hiện đại. 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 nghiên cứu sinh Vũ Hoàng Nam, thực hiện tại Viện Công nghệ Nano - Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia TP. Hồ Chí Minh dưới sự hướng dẫn của PGS.TS Cao Minh Thì và TS. Lê Minh Hưng, mang tựa đề: "Nghiên cứu tương tác giữa graphene và đế bán dẫn bằng lý thuyết phiếm hàm mật độ" (2023). Công trình tập trung giải quyết đồng thời hai nghịch lý vật lý cơ bản mang tính nền tảng: thứ nhất, sự thiếu vắng vùng cấm năng lượng ($E_g = 0$) tại điểm Dirac sơ cấp (Primitive Dirac Point - PDP) do cấu trúc đối xứng dải năng lượng ngăn cản graphene hoạt động như một kênh dẫn bán dẫn trong linh kiện số logic; thứ hai, sự mâu thuẫn giữa hành vi dịch chuyển mức Fermi ($\Delta E_F$) đặc thù của bán kim loại với các mô hình tiếp giáp Schottky kim loại/bán dẫn (Metal/Semiconductor - M/S) truyền thống vốn bị chi phối bởi hiện tượng ghìm mức Fermi (Fermi-level pinning).

Khoảng trống nghiên cứu (research gap) được xác định rõ nét từ thực trạng y văn: các nghiên cứu trước đây chủ yếu khảo sát thế tuần hoàn trơn mịn hoặc siêu mạng tầm xa ($L > 10\text{ nm}$), chưa làm rõ cơ chế tán xạ ngoại hồ (intervalley scattering) đối với siêu mạng tuần hoàn tầm ngắn chịu tác động của trường thế ngoài sắc nhọn quy mô nguyên tử; đồng thời chưa có mô hình lượng tử giải thích thấu đáo cơ chế truyền điện tích, vai trò của thế đẩy ngược (push-back effect) và động học của vùng nghèo hạt tải trong tiếp giáp van der Waals (vdW) giữa graphene và đế bán dẫn silicon (G/Si).

Để giải quyết triệt để vấn đề, luận án xây dựng hai hệ thống câu hỏi và giả thuyết nghiên cứu cụ thể:

  1. Hệ câu hỏi $Q_1$ & Giả thuyết $H_1$: Thế ngoài sắc nhọn từ liên kết chưa bão hòa của đế silicon carbide (SiC) tác động như thế nào lên các siêu mạng graphene tuần hoàn tầm ngắn? Giả thuyết $H_1$ khẳng định việc lựa chọn cấu hình siêu mạng có tâm vùng Brillouin siêu mạng (SBZ) gấp vào góc $K(K')$ của vùng Brillouin nguyên tố (GBZ) sẽ bẻ gãy đối xứng chiral, kích hoạt tán xạ ngược ngoại hồ và mở vùng cấm năng lượng hiệu dụng tại cả điểm Dirac sơ cấp lẫn thứ cấp (Secondary Dirac Point - SDP).
  2. Hệ câu hỏi $Q_2$ & Giả thuyết $H_2$: Cơ chế vận chuyển điện tích và sự hình thành dipole tiếp giáp tại mặt phân giới G/Si diễn ra theo mô hình nào? Giả thuyết $H_2$ cho rằng tiếp giáp vdW G/Si bảo toàn tính bán kim loại của graphene, cho phép điều biến công thoát thông qua dịch chuyển mức Fermi $\Delta E_F$, đư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$) thay vì bị giới hạn Bardeen ($S \to 0$).

Khung lý thuyết của luận án được thiết lập vững chắc trên nền tảng cơ học lượng tử của phương trình Dirac cho fermion tương đối tính không khối lượng, kết hợp lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) của Hohenberg-Kohn và Kohn-Sham. Phạm vi nghiên cứu bao quát các cấu hình giao diện nguyên tử phức tạp: từ các siêu mạng $G-5\sqrt{3}\times 5\sqrt{3}\text{-R}30^\circ/\text{SiC}$, $G-9\times 9$, $G-\sqrt{21}\times\sqrt{21}\text{-R}10.893^\circ$, $G-7\times 7$ cho đến tiếp giáp dị thể $G-\sqrt{39}\times\sqrt{39}\text{-R}16.1^\circ/\text{Si}(111)\text{-}4\times 4$. Ý nghĩa đột phá của luận án thể hiện ở việc định lượng chính xác độ mở vùng cấm, giải mã bản chất các trạng thái phân cực giao diện và cung cấp nguyên lý thiết kế cho thế hệ linh kiện barristor siêu nhanh vượt qua giới hạn bán dẫn truyền thống.

Literature Review và Positioning

Lịch sử nghiên cứu cấu trúc điện tử của graphene khởi nguồn từ mô hình liên kết mạnh (Tight-Binding Model) của Wallace (1947) và Slonczewski - Weiss (1958), chứng minh rằng mạng lục giác hai chiều với hai nguyên tử carbon $A$ và $B$ trong một ô cơ sở tạo ra sự suy biến tuyến tính tại hai góc $K$ và $K'$ của GBZ với biểu thức năng lượng: $$E(k) = \pm t \sqrt{1 + 4\cos\left(\frac{3a_0}{2}k_x\right)\cos\left(\frac{\sqrt{3}a_0}{2}k_y\right) + 4\cos^2\left(\frac{\sqrt{3}a_0}{2}k_y\right)}$$ trong đó năng lượng nhảy lân cận $t \approx 2.7\text{ eV}$ và khoảng cách liên kết $a_0 \approx 1.42\text{ \AA}$. Nghiên cứu của Novoselov và Geim (2004) đã hiện thực hóa thực nghiệm graphene, xác nhận vận tốc Fermi $v_F = 3ta_0/(2\hbar) \approx 10^6\text{ m/s}$ và sự hiện diện của các fermion Dirac không khối lượng.

Tuy nhiên, cuộc tranh luận học thuật về cơ chế mở vùng cấm và điều khiển dải năng lượng nhanh chóng phân hóa thành các luồng quan điểm đối lập:

  • Trường phái bẻ gãy đối xứng mạng phụ (Sublattice Symmetry Breaking): Giovannetti cùng các cộng sự (2007) và Hunt cùng các cộng sự (2013) lập luận rằng việc đặt graphene lên đế boron nitride lục giác (hBN) với độ lệch mạng $1.8%$ tạo ra trường thế bất đối xứng tác động lên mạng phụ $A$ và $B$, dẫn đến sự trộn lẫn dải $\pi(K)-\pi^*(K)$ và mở vùng cấm tại PDP.
  • Trường phái bẻ gãy đối xứng chiral (Chiral Symmetry Breaking): Ngược lại, các nghiên cứu lý thuyết của Cheianov cùng các cộng sự (2006) và Bena (2008) chỉ ra rằng việc bẻ gãy đối xứng nghịch đảo không đủ để kiểm soát tán sắc phi tuyến ở mức năng lượng cao; cần phải có sự trộn lẫn góc $K$ và $K'$ (trộn $K-K'$) thông qua các siêu mạng biến dạng Kekulé (Kekulé-O và Kekulé-Y). Song, các công trình này tồn tại mâu thuẫn lớn khi khảo sát hành vi tại điểm Dirac thứ cấp (SDP): trong khi Park cùng các cộng sự (2008) cho rằng độ mở vùng cấm tại PDP và SDP có chung nguồn gốc từ thế mạng phụ, các quan sát thực nghiệm của Yankowitz cùng các cộng sự (2012) lại phát hiện sự bất đối xứng rõ rệt dưới tác động của áp suất và góc xoay.
                      ┌──────────────────────────────────────────┐
                      │    Tranh luận học thuật tiếp giáp M/S   │
                      └────────────────────┬─────────────────────┘
                                           │
             ┌─────────────────────────────┴─────────────────────────────┐
             ▼                                                           ▼
┌─────────────────────────┐                                 ┌─────────────────────────┐
│ Mô hình Cổ điển (3D)    │                                 │ Mô hình vdW Mới (2D)    │
├─────────────────────────┤                                 ├─────────────────────────┤
│ - Bardeen (1947)        │                                 │ - Liu et al. (2018)     │
│ - Cowley & Sze (1965)   │                                 │   S = 0.96 (vdW MoS2)   │
│ - Heine/Tersoff (MIGS)  │                                 │ - LaGasse et al. (2019) │
│ - Tung (Bond Polariz.)  │                                 │   S = 1.04 (vdW WSe2)   │
│ => Ghìm Fermi (S -> 0)  │                                 │ => Schottky-Mott (S->1) │
└─────────────────────────┘                                 └─────────────────────────┘

Trong lĩnh vực tiếp giáp kim loại/bán dẫn, các mô hình cổ điển của Bardeen (1947), Cowley & Sze (1965), lý thuyết trạng thái vùng cấm cảm ứng kim loại (MIGS) của Heine (1965) và Tersoff (1984), cùng lý thuyết phân cực liên kết của Tung (2000, 2001) đều thống nhất rằng mật độ trạng thái bề mặt cao ($D_{it}$) hoặc sự liên tục của đuôi hàm sóng kim loại vào vùng cấm bán dẫn sẽ ghìm chặt mức Fermi ($E_F$) tại mức trung tính điện ($E_{CNL}$) hoặc điểm quay lui ($E_B$), khiến độ nhạy $S \ll 1$. Khi chuyển sang vật liệu hai chiều, Liu cùng các cộng sự (2018) trên tiếp giáp vdW kim loại/$\text{MoS}2$ đạt được $S = 0.96$ (so với $S = 0.09$ của tiếp giáp hóa học); LaGasse cùng các cộng sự (2019) trên tiếp giáp graphene/$\text{WSe}2$ ghi nhận $S = 1.04$. Tuy nhiên, với tiếp giáp graphene/silicon ba chiều (G/Si), Yang cùng các cộng sự (2012) khi chế tạo linh kiện barristor lại quan sát thấy sự biến thiên phức tạp của độ cao hàng rào Schottky ($\Phi{Bn}$) phụ thuộc nồng độ pha tạp, mâu thuẫn với giá trị lý thuyết Schottky-Mott thuần túy $\Phi{Bn} = \Phi_M - \chi_S = 4.6\text{ eV} - 4.05\text{ eV} = 0.55\text{ eV}$.

Luận án của NCS. Vũ Hoàng Nam định vị chính xác vào điểm giao thoa của các tranh luận trên: thiết lập cơ chế vi mô giải thích sự trồi lên và tán xạ của các nón Dirac bản sao trong siêu mạng tầm ngắn, đồng thời giải bài toán động học truyền điện tích và thế đẩy ngược tại tiếp giáp G/Si nhằm dung hòa sự khác biệt giữa vật lý tiếp giáp 2D vdW và bán dẫn khối 3D.

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

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

Luận án đã mở rộng và tái cấu trúc các lý thuyết vật lý chất rắn nền tảng thông qua các luận điểm khoa học chuẩn xác:

  1. Mở rộng lý thuyết tán xạ Dirac đối với siêu mạng graphene tuần hoàn tầm ngắn: Công trình chứng minh rằng điều kiện tiên quyết để bẻ gãy đối xứng chiral tại PDP không chỉ giới hạn ở cấu trúc Kekulé-O nhỏ nhất ($G-\sqrt{3}\times\sqrt{3}\text{-R}30^\circ$) như các công bố kinh điển của Cheianov hay Pereira, mà hoàn toàn có thể hiện thực hóa trên các siêu mạng mở rộng $G-\sqrt{3}N\times\sqrt{3}N\text{-R}\phi^\circ$ ($N > 1$). Luận án chứng minh vector xung lượng kích thước $|Q| \sim 1/L$ từ thế ngoài sắc nhọn nguyên tử cho phép vận chuyển các nón Dirac bản sao không tương đương về tâm $\gamma$ của SBZ, giải tỏa sự suy biến và tạo lập vùng cấm năng lượng thực sự.
  2. Xác lập cơ chế phân kỳ nguồn gốc vật lý giữa PDP và SDP: Luận án bác bỏ giả định phổ biến rằng vùng cấm tại PDP và SDP có chung bản chất bẻ gãy đối xứng mạng phụ. Bằng chứng định lượng chỉ ra rằng: vùng cấm tại PDP bắt nguồn từ sự tán xạ ngược ngoại hồ (intervalley backscattering) bẻ gãy đối xứng chiral, trong khi cấu trúc dải tại SDP hình thành từ sự cạnh tranh giữa tán xạ nội hồ (intravalley scattering) mở dải và hiệu ứng xuyên hầm Klein đóng dải.
  3. Phát triển khung lý thuyết tiếp giáp Schottky bán kim loại/bán dẫn: Luận án điều chỉnh mô hình Bardeen - Cowley - Sze và lý thuyết phân cực liên kết của Tung cho hệ graphene/bán dẫn. Công thức độ cao hàng rào Schottky được tổng quát hóa thành: $$\Phi_{Bn} = S(\Phi_{G0} + \Delta E_F - \chi_S) + (1 - S)(E_g - \phi_{0}) - \Delta_{push-back}$$ chứng minh rằng mật độ trạng thái tiệm cận không tại PDP ($DOS \sim |E|$) cho phép mức Fermi dịch chuyển tự do $\Delta E_F = \hbar v_F \sqrt{\pi n}$, đóng vai trò quyết định độ nhạy tiếp giáp thay vì bị ghìm cứng như kim loại khối.
                      ┌──────────────────────────────────────────┐
                      │    Khung Tích Hợp Đa Lý Thuyết DFT       │
                      └────────────────────┬─────────────────────┘
                                           │
         ┌─────────────────────────────────┼─────────────────────────────────┐
         ▼                                 ▼                                 ▼
┌──────────────────┐             ┌──────────────────┐             ┌──────────────────┐
│ Dirac-Weyl & TB  │             │ Hohenberg-Kohn & │             │ Schottky-Mott &  │
│ (Fermion Chiral) │             │ Kohn-Sham (DFT)  │             │ Bardeen-Tung (MS)│
└────────┬─────────┘             └────────┬─────────┘             └────────┬─────────┘
         │                                │                                │
         └────────────────────────────────┼────────────────────────────────┘
                                          │
                                          ▼
                      ┌──────────────────────────────────────────┐
                      │ Mô hình Tái cấu trúc Dải & Tiếp giáp vdW │
                      │ - Trải dải SBZ -> GBZ (Band Unfolding)   │
                      │ - Lớp nghèo & Thế phân cực Pauli         │
                      └──────────────────────────────────────────┘

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

Luận án thiết lập một khung phân tích tích hợp ba trục lý thuyết: Lý thuyết tán sắc Dirac-Weyl hạt tải tương đối tính, Lý thuyết phiếm hàm mật độ lượng tử đa hạt, và Cơ học tiếp giáp dị thể bán dẫn.

Cách tiếp cận phân tích mới nổi bật ở việc kết hợp kỹ thuật chiếu hàm sóng giả thế (Projector Augmented-Wave) với giải thuật trải dải năng lượng (Band Unfolding Technique). Phương pháp này cho phép ánh xạ ma trận mật độ xác suất và vector sóng từ vùng Brillouin siêu mạng (SBZ) bị gấp phức tạp trở lại vùng Brillouin nguyên tố (GBZ), cô lập hoàn toàn tương tác giữa các nón Dirac bản sao.

Điều kiện biên xác định rõ: Khung phân tích áp dụng nghiêm ngặt cho tương tác vật lý vdW có khoảng cách cân bằng $d > 3.0\text{ \AA}$, thế ngoài có gradient không gian sắc nhọn ở thang nguyên tử, và nồng độ pha tạp nền của đế bán dẫn $N_D < 10^{19}\text{ cm}^{-3}$.

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

Thiết kế nghiên cứu

Nghiên cứu được định hình theo triết học duy thực phản biện (Critical Realism) và thực chứng lượng tử (Quantum Positivism), xuất phát từ nguyên lý ban đầu (ab initio) không phụ thuộc vào tham số thực nghiệm ngoại suy. Thiết kế nghiên cứu đa tầng (multi-level design) được cấu trúc chặt chẽ:

  • Tầng nguyên tử/hình học: Tối ưu hóa cấu trúc mạng dị thể, tìm trạng thái hồi phục cấu trúc năng lượng cực tiểu thông qua thuật toán Quasi-Newton (RMM-DIIS).
  • Tầng điện tử/dải năng lượng: Giải hệ phương trình tự hợp Kohn-Sham để thu nhận phân bố mật độ điện tích $n(\mathbf{r})$, cấu trúc dải năng lượng $E(\mathbf{k})$ và mật độ trạng thái (DOS).
  • Tầng tiếp giáp/thế tĩnh điện: Tính toán hàm thế tĩnh điện trung bình phẳng $\overline{V}(z)$, hàm hiệu mật độ điện tích $\Delta \rho(z) = \rho_{G/Si}(z) - \rho_G(z) - \rho_{Si}(z)$, và dịch chuyển thế chân không để xác định chính xác bậc thế tiếp giáp $\Delta_{in}$.
┌────────────────────────────────────────────────────────────────────────────┐
│                    Quy trình Tính toán Lượng tử DFT-VASP                  │
└─────────────────────────────────────┬──────────────────────────────────────┘
                                      │
                                      ▼
┌────────────────────────────────────────────────────────────────────────────┐
│ 1. Xây dựng Cấu hình Siêu mạng & Tối ưu hóa Hình học                      │
│    - G/SiC: G-5√3x5√3-R30°, G-9x9, G-√21x√21-R10.893°, G-7x7              │
│    - G/Si: G-√39x√39-R16.1° / Si(111)-4x4 (Mismatch < 1.5%)                │
└─────────────────────────────────────┬──────────────────────────────────────┘
                                      │
                                      ▼
┌────────────────────────────────────────────────────────────────────────────┐
│ 2. Vòng lặp Tự hợp Kohn-Sham (SCF) với Phiếm hàm optB86b-vdW               │
│    - Năng lượng hội tụ: ΔE < 10^-6 eV | Lực nguyên tử: F < 0.01 eV/Å       │
│    - Vacuum layer > 15 Å | Dipole correction dọc trục z                    │
└─────────────────────────────────────┬──────────────────────────────────────┘
                                      │
                                      ▼
┌────────────────────────────────────────────────────────────────────────────┐
│ 3. Phân tích Dải năng lượng & Điện động học Giao diện                      │
│    - Trải dải SBZ -> GBZ (Band Unfolding) để quan sát phân kỳ nón Dirac    │
│    - Tích phân thế Hartree V(z), mật độ điện tích phẳng Δρ(z)              │
│    - Tách bậc thế đẩy ngược Pauli Δpush-back và thế truyền điện tích ΔCT   │
└────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình tính toán tuân thủ các chuẩn mực quốc tế cao nhất trong vật lý tính toán:

  • Công cụ thực thi: Bộ phần mềm VASP (Vienna Ab initio Simulation Package), tích hợp các thuật toán giải phương trình Kohn-Sham: $$\left[ -\frac{\hbar^2}{2m}\nabla^2 + V_{ext}(\mathbf{r}) + V_H(\mathbf{r}) + V_{xc}(\mathbf{r}) \right] \psi_i(\mathbf{r}) = \varepsilon_i \psi_i(\mathbf{r})$$
  • Hiệu chỉnh tương tác phân tán: Phiếm hàm phi định xứ optB86b-vdW được sử dụng làm trung tâm để mô tả chính xác lực phân tán van der Waals, khắc phục triệt để nhược điểm của các phiếm hàm chuẩn LDA (thường đánh giá quá cao lực liên kết) và GGA-PBE (đánh giá quá thấp liên kết vdW).
  • Mô hình phiến (Slab Model): Thiết lập phiến bán dẫn có độ dày tối thiểu 6-12 lớp nguyên tử, các liên kết treo ở mặt đáy được thụ động hóa hoàn toàn bằng nguyên tử hydro ($-\text{H}$). Khoảng chân không (vacuum layer) được thiết lập $> 15\text{ \AA}$ để triệt tiêu hoàn toàn tương tác giả giữa các phiến tuần hoàn tuần hoàn theo trục $z$.
  • Độ tin cậy và chuẩn hóa: Năng lượng cắt sóng phẳng ($E_{cut}$) được kiểm tra hội tụ nghiêm ngặt ở mức $400 - 500\text{ eV}$. Lưới không gian $k$ (Monkhorst-Pack grid) được phân bố dày đặc (lên tới $12\times 12\times 1$ cho ô siêu mạng), tiêu chuẩn hội tụ năng lượng điện tử đạt $\Delta E < 10^{-6}\text{ eV}$, và tiêu chuẩn hội tụ lực tác dụng lên mỗi nguyên tử đạt $|F| < 0.01\text{ eV/\AA}$.

Data và phân tích

Phân tích dữ liệu lượng tử được thực hiện với các kiểm tra độ vững (robustness checks) toàn diện:

  • Kiểm chứng đa phiếm hàm: So sánh chéo kết quả giữa optB86b-vdW, LDA và PBE trên cùng các cấu hình $G/\text{Si}$ để đánh giá độ nhạy của khoảng cách giao diện và năng lượng liên kết.
  • Khử sai số giả tạo: Áp dụng thuật toán hiệu chỉnh lưỡng cực (dipole correction) theo phương vuông góc với mặt phân giới nhằm loại bỏ sai số thế tĩnh điện do sự tái phân bố điện tích bất đối xứng gây ra.
  • Trải dải phổ năng lượng: Phép chiếu hàm sóng Bloch từ không gian siêu mạng $E_{SBZ}(\mathbf{K})$ sang không gian nguyên tố $E_{GBZ}(\mathbf{k})$ được định lượng hóa bằng trọng số phổ $P_{K}(\mathbf{k}) = \sum_m |\langle \psi_{m\mathbf{K}} | \phi_{n\mathbf{k}} \rangle|^2$, cho phép hiển thị trực quan và sắc nét sự biến đổi của nón Dirac mà không bị nhiễu bởi các dải bị gấp.

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

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

  1. Khám phá quy luật điều khiển nón Dirac bản sao trong siêu mạng graphene tuần hoàn tầm ngắn: Luận án chứng minh rằng trên bề mặt $\text{SiC}(0001)$, thế ngoài sắc nhọn từ các liên kết chưa bão hòa kích hoạt sự trồi lên của các nón Dirac bản sao tại các điểm đối xứng cao. Đối với siêu mạng $G-5\sqrt{3}\times 5\sqrt{3}\text{-R}30^\circ/\text{SiC}-4\sqrt{3}\times 4\sqrt{3}\text{-R}30^\circ$, tâm $\gamma$ của SBZ gấp chính xác vào điểm $K(K')$ của GBZ, tạo ra sự giao thoa và tán xạ ngược ngoại hồ cực mạnh giữa các nón Dirac có tính đối xứng chiral ngược nhau ($s_K$ và $s_{K'}$), bẻ gãy đối xứng chiral và mở một vùng cấm năng lượng trực tiếp tại PDP đạt giá trị danh định $E_g \approx 0.12 - 0.35\text{ eV}$ tùy thuộc vào mức độ tái cấu trúc.
  2. Giải mã cơ chế vật lý độc lập của điểm Dirac thứ cấp (SDP): Khác biệt với PDP, tại các siêu mạng $G-9\times 9$ và $G-7\times 7$, sự tán xạ nội hồ chiếm ưu thế. Luận án phát hiện rằng các nón Dirac thứ cấp xuất hiện tại biên của SBZ do sự giao thoa dải năng lượng không bị triệt tiêu bởi hiệu ứng xuyên hầm Klein. Sự mở vùng cấm tại SDP phụ thuộc vào tensor biến dạng mạng cục bộ chứ không đồng nhất với cơ chế mở dải tại PDP, giải quyết triệt để bất đồng y văn tồn tại suốt một thập kỷ giữa các nhóm nghiên cứu quốc tế.
  3. Xác lập trạng thái tiếp giáp vdW thuần túy trên bề mặt $\text{Si}(111)\text{-}4\times 4$: Tính toán DFT với phiếm hàm optB86b-vdW khẳng định rằng việc tái cấu trúc bề mặt $\text{Si}(111)\text{-}4\times 4$ với các nguyên tử Si bổ sung tạo ra bề mặt phẳng tối ưu, duy trì khoảng cách vdW ổn định $d_{eq} \approx 3.28\text{ \AA}$ và độ lệch mạng $< 1.5%$. Cấu hình này bảo toàn nguyên vẹn tính chất bán kim loại và độ tán sắc tuyến tính của graphene, bác bỏ các dự đoán sai lệch trước đó về sự hình thành liên kết cộng hóa trị phá hủy mạng $\text{sp}^2$.
  4. Tách biệt định lượng hiệu ứng đẩy ngược Pauli và thế truyền điện tích: Luận án lần đầu tiên định lượng độc lập hai thành phần tạo nên bậc thế tiếp giáp $\Delta_{in}$ tại giao diện G/Si:
    • Thế đẩy ngược Pauli ($\Delta_{push-back} \approx 0.21\text{ eV}$): sinh ra do sự nén đuôi hàm sóng điện tử tự do khi hai đám mây điện tích tiếp cận nhau.
    • Thế truyền điện tích ($\Delta_{CT}$): phụ thuộc tuyến tính vào mật độ điện tích dịch chuyển qua mặt phân giới.
  5. Tiệm cận giới hạn Schottky-Mott ($S \approx 0.92 - 0.98$) trong tiếp giáp G/Si có điều biến: Khi bề mặt Si được thụ động hóa hoặc kiểm soát trạng thái bề mặt định xứ, mật độ bẫy điện tích $D_{it} \to 0$, mức Fermi của graphene dịch chuyển tự do $\Delta E_F$ theo nồng độ hạt tải cảm ứng, đưa tham số độ nhạy hàng rào Schottky đạt $S = 0.95$, tiệm cận sát giới hạn Schottky-Mott lý tưởng ($S=1$) và vượt trội hoàn toàn so với tiếp giáp kim loại khối/Si truyền thống ($S \approx 0.27$).
┌────────────────────────────────────────────────────────────────────────────┐
│           Bản chất Phân kỳ Năng lượng tại Giao diện Siêu mạng             │
├─────────────────────────────────────┬──────────────────────────────────────┤
│ Điểm Dirac Sơ cấp (PDP)             │ Điểm Dirac Thứ cấp (SDP)             │
├─────────────────────────────────────┼──────────────────────────────────────┤
│ - Vị trí: Góc K(K') của GBZ         │ - Vị trí: Cạnh/biên của SBZ          │
│ - Cơ chế: Bẻ gãy đối xứng Chiral    │ - Cơ chế: Cạnh tranh tán xạ nội hồ   │
│ - Tán xạ: Ngoại hồ (Intervalley)    │   và xuyên hầm Klein                 │
│ - Nguồn gốc: Giao thoa nón phi tương│ - Nguồn gốc: Tái cấu trúc dải tuần   │
│   đương tại điểm γ của SBZ          │   hoàn & tensor biến dạng cục bộ     │
└─────────────────────────────────────┴──────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình lượng tử hoàn chỉnh thống nhất cơ học tương đối tính Dirac trong vật liệu 2D với vật lý tiếp giáp bán dẫn 3D, làm sáng tỏ vai trò của đối xứng chiral và tương tác vdW.
  • Về mặt phương pháp: Chuẩn hóa quy trình tính toán DFT kết hợp trải dải năng lượng (Band Unfolding) và hiệu chỉnh vdW phi định xứ cho các hệ dị thể diện tích lớn ($> 100$ nguyên tử), có thể chuyển giao trực tiếp để nghiên cứu các họ vật liệu 2D khác như TMDs, MXenes, phosphorene.
  • Về mặt thực tiễn và công nghệ: Mở ra lộ trình công nghệ khả thi để chế tạo các transistor barristor graphene/silicon có tốc độ đóng ngắt cực cao, tổn hao năng lượng thấp, và pin quang điện graphene/silicon hiệu suất cao tương thích hoàn hảo với dây chuyền chế tạo CMOS hiện nay.
  • Về mặt chính sách và đầu tư: Cung cấp cơ sở khoa học định lượng để các trung tâm công nghệ cao quốc gia và các doanh nghiệp bán dẫn định hướng đầu tư nghiên cứu vật liệu mới tích hợp trên nền tảng công nghệ silicon hiện hữu.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn nội tại mang tính học thuật:

  1. Giới hạn xấp xỉ của phiếm hàm DFT: Các phiếm hàm trao đổi - tương quan tiêu chuẩn (kể cả optB86b-vdW) có xu hướng đánh giá thấp (underestimate) độ rộng vùng cấm năng lượng của chất bán dẫn so với thực nghiệm do thiếu số hạng tự tương tác (self-interaction error).
  2. Kích thước siêu mạng và tài nguyên tính toán: Do giới hạn năng lực siêu máy tính, các siêu mạng khảo sát dừng lại ở các cấu hình tuần hoàn có chu kỳ ngắn ($L < 3\text{ nm}$); chưa thể mô phỏng trực tiếp các góc xoay moiré cực nhỏ ($\theta < 1.1^\circ$) vốn đòi hỏi hàng chục nghìn nguyên tử trong một ô đơn vị.
  3. Ảnh hưởng của nhiệt độ và khuyết tật thực tế: Các tính toán được thực hiện tại nhiệt độ tuyệt đối $T = 0\text{ K}$ trong môi trường chân không lý tưởng, chưa bao quát đầy đủ sự tán xạ phonon nhiệt độ phòng và các sai hỏng mạng ngẫu nhiên (chất pha tạp cục bộ, tạp chất môi trường).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Nâng cấp tính toán dải năng lượng bằng phương pháp lý thuyết nhiễu loạn nhiều hạt $GW$ kết hợp phương trình Bethe-Salpeter (BSE) để thu nhận chính xác hiệu ứng kích thích quang học và exciton.
  • Hướng 2: Mở rộng mô hình tiếp giáp vdW sang các cấu hình xoay góc ngẫu nhiên (Twisted Bilayer / Moiré Heterostructures) trên đế bán dẫn oxit phân cực cao ($\text{HfO}_2$, $\beta\text{-Ga}_2\text{O}_3$).
  • Hướng 3: Khảo sát động học vận chuyển hạt tải phụ thuộc thời gian (Time-Dependent DFT) dưới tác động của xung laser siêu ngắn (femtosecond laser pulse).
  • Hướng 4: Nghiên cứu thực nghiệm đối chứng bằng kỹ thuật hiển vi quét đầu dò xuyên hầm quét (STM/STS) và quang phổ phát xạ quang góc phân giải (ARPES) trên các mẫu chế tạo bằng epitaxy chùm phân tử (MBE).

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

  • Ảnh hưởng học thuật: Các kết quả của luận án đóng góp trực tiếp vào kho tàng tri thức của vật lý chất rắn và khoa học vật liệu nano, dự kiến tạo ra chỉ số trích dẫn cao trong các tạp chí chuyên ngành quốc tế uy tín (hệ ISI/Scopus Q1 như Physical Review B, Applied Surface Science, Carbon).
  • Chuyển dịch công nghiệp bán dẫn: Cung cấp giải pháp kỹ thuật cốt lõi giúp ngành công nghiệp vi điện tử vượt qua "điểm nghẽn" vật lý của định luật Moore, mở đường cho việc tích hợp điện cực trong suốt, siêu mỏng và linh hoạt bằng graphene vào công nghệ chip bán dẫn silicon.
  • Lợi ích xã hội và môi trường: Các thiết bị bán dẫn và pin mặt trời dị thể G/Si hiệu suất cao phát triển từ nguyên lý của luận án sẽ cắt giảm tiêu thụ năng lượng toàn cầu, giảm phát thải nhiệt trong các trung tâm dữ liệu siêu lớn (Data Centers).
  • Vị thế quốc tế: Khẳng định năng lực nghiên cứu vật lý tính toán đỉnh cao của các nhà khoa học Việt Nam, hội nhập bình đẳng với các nhóm nghiên cứu hàng đầu tại Nhật Bản (Viện IMR, Đại học Tohoku) và thế giới.

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

  • Nghiên cứu sinh và Giới học giả Trẻ: Tiếp cận quy trình phương pháp luận mẫu mực về mô phỏng lượng tử DFT, kỹ thuật trải dải năng lượng và cách xử lý tương tác bề mặt vdW phức tạp.
  • Các Nhà Khoa học Vật liệu & Vật lý Chất rắn: Khai thác các dữ liệu định lượng về dải năng lượng, độ cao hàng rào Schottky và cơ chế tán xạ Dirac để phát triển các mô hình lý thuyết mở rộng.
  • Kỹ sư R&D trong Công nghiệp Bán dẫn: Ứng dụng các thông số tiếp xúc giao diện, độ dày vdW và điều kiện pha tạp để tối ưu hóa quy trình lắng đọng hóa hơi hóa học (CVD) và chế tạo transistor barristor.
  • Nhà hoạch định Chính sách Khoa học & Công nghệ: Có luận cứ xác đáng để phê duyệt các dự án đầu tư phòng thí nghiệm tính toán hiệu năng cao (HPC) phục vụ thiết kế vật liệu chiến lược quốc gia.

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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Đối xứng Chiral và Tán xạ Dirac sang các siêu mạng tuần hoàn tầm ngắn $G-\sqrt{3}N\times\sqrt{3}N\text{-R}\phi^\circ$ ($N > 1$) dưới tác động của thế ngoài nguyên tử sắc nhọn. Luận án đã chứng minh cơ chế bẻ gãy đối xứng chiral thông qua việc gấp điểm $\gamma$ của SBZ vào điểm $K(K')$ của GBZ, tạo ra sự tán xạ ngược ngoại hồ điều khiển được, mở ra vùng cấm năng lượng tại PDP mà không phá hủy độ linh động của hạt tải.

2. Đột phá phương pháp luận của luận án khi so sánh với ít nhất hai nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Giovannetti et al. (2007) (chỉ khảo sát ô nguyên tố trên đế hBN bằng phiếm hàm LDA tiêu chuẩn gây sai số liên kết) và công trình của Park et al. (2008) (sử dụng thế nhiễu loạn trơn mịn liên tục 1D không phản ánh đúng cấu trúc nguyên tử thực), luận án của NCS. Vũ Hoàng Nam đã tạo ra bước tiến vượt bậc:

  • Sử dụng phiếm hàm phi định xứ tiên tiến optB86b-vdW giải quyết chính xác năng lượng tương tác phân tán.
  • Triển khai kỹ thuật Band Unfolding trên siêu mạng 2D thực tế trên bề mặt bán dẫn 3D ($\text{SiC}$, $\text{Si}$), cho phép phân tách rành mạch phổ năng lượng của dải chính và dải bản sao.

3. Phát hiện nào trong dữ liệu tính toán mang lại sự bất ngờ lớn nhất?

Phát hiện bất ngờ nhất là sự phân ly hoàn toàn về nguồn gốc vật lý giữa độ mở vùng cấm tại PDP và SDP. Trái ngược với quan niệm phổ biến cho rằng hai vùng cấm này đồng biến và cùng do thế mạng phụ quy định, dữ liệu DFT chỉ ra rằng vùng cấm PDP chịu sự chi phối của đối xứng chiral và tán xạ ngoại hồ, trong khi SDP được điều khiển bởi tính chất bất đẳng hướng của nón Dirac (trigonal warping) và sự cạnh tranh tán xạ nội hồ dưới hiệu ứng xuyên hầm Klein.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết:

  • Cấu hình mạng tinh thể, vector mạng thực ($\mathbf{a}_1, \mathbf{a}_2$) và mạng đảo ($\mathbf{b}_1, \mathbf{b}_2$).
  • Năng lượng cắt sóng phẳng ($E_{cut} = 400 - 500\text{ eV}$), lưới mẫu điểm $k$, tiêu chuẩn hội tụ lực và năng lượng.
  • Các bước tạo file đầu vào (POSCAR, INCAR, KPOINTS, POTCAR) cho phần mềm VASP, cách thiết lập thế giả PAW và quy chuẩn trải dải năng lượng.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1-3 năm): Tích hợp hiệu ứng kích thích nhiều hạt $GW\text{-BSE}$ và khảo sát tiếp giáp vdW với các bán dẫn oxit dải rộng.
  • Giai đoạn 2 (4-6 năm): Mô phỏng cấu trúc góc xoay ma thuật (Moiré Valleytronics) và hiện tượng siêu dẫn tương quan ở giao diện 2D/3D.
  • Giai đoạn 3 (7-10 năm): Chuyển giao quy trình thiết kế vật liệu sang chế tạo thực nghiệm và thử nghiệm tiền công nghiệp chip logic barristor G/Si tích hợp mật độ cao.

Kết luận

  1. Luận án giải quyết thành công bài toán mở vùng cấm năng lượng tại điểm Dirac sơ cấp (PDP) của graphene bằng cơ chế bẻ gãy đối xứng chiral trên các siêu mạng tuần hoàn tầm ngắn $G-\sqrt{3}N\times\sqrt{3}N\text{-R}\phi^\circ$ nhờ thế ngoài sắc nhọn từ đế bán dẫn SiC.
  2. Thiết lập cơ sở vật lý khẳng định tính độc lập về bản chất dải năng lượng giữa PDP (tán xạ ngoại hồ bẻ gãy đối xứng chiral) và SDP (tán xạ nội hồ cạnh tranh với xuyên hầm Klein).
  3. Xác định cấu trúc tối ưu của tiếp giáp vdW phẳng $G-\sqrt{39}\times\sqrt{39}\text{-R}16.1^\circ/\text{Si}(111)\text{-}4\times 4$, bảo toàn tính bán kim loại và cấu trúc dải tuyến tính của graphene.
  4. Tách biệt và định lượng chính xác hai thành phần của bậc thế giao diện: thế đẩy ngược Pauli ($\Delta_{push-back}$) và thế truyền điện tích ($\Delta_{CT}$).
  5. Chứng minh độ nhạy hàng rào Schottky của tiếp giáp G/Si đạt $S \approx 0.95$, tiệm cận giới hạn Schottky-Mott lý tưởng nhờ khả năng dịch chuyển mức Fermi $\Delta E_F$ linh hoạt của graphene.
  6. Mở ra ba hướng nghiên cứu đột phá: Điện tử học thung lũng (Valleytronics), Linh kiện logic Barristor thế hệ mới, và Vật liệu dị thể Moiré 2D/3D hiệu năng cao.