Tổng quan luận án

Vật liệu graphene hai chiều sở hữu cấu trúc vùng năng lượng đặc biệt với các hạt tải điện là fermion không khối lượng, tán sắc tuyến tính theo phương trình Dirac tại lân cận các điểm đối xứng $K$ và $K'$ của vùng Brillouin. Nhờ đó, graphene đạt vận tốc Fermi cao ($v_F \approx 10^6 \text{ m/s}$), mở ra tiềm năng ứng dụng trong các thiết bị điện tử tốc độ cao tích hợp cùng công nghệ bán dẫn truyền thống. Tuy nhiên, việc đưa graphene vào các thiết bị điện tử số gặp hai trở ngại vật lý cốt lõi:

  1. Graphene không có vùng cấm năng lượng ($E_g = 0$), khiến việc điều khiển đóng/ngắt dòng điện trong các linh kiện logic gặp khó khăn.
  2. Graphene là một bán kim loại có mật độ trạng thái (DOS) tiệm cận không tại điểm Dirac sơ cấp (PDP), khiến mức năng lượng Fermi ($E_F$) và công thoát ($\Phi_G$) thay đổi rất nhanh theo nồng độ hạt tải điện.

Xuất phát từ các đặc tính trên, luận án của nghiên cứu sinh Vũ Hoàng Nam tập trung giải quyết hai vấn đề nghiên cứu chính được đánh số như sau:

  • Vấn đề nghiên cứu thứ nhất: Nghiên cứu cơ chế tương tác giữa siêu mạng graphene tuần hoàn tầm ngắn với bề mặt đế bán dẫn chứa các liên kết chưa bão hòa (cụ thể là đế silicon carbide - SiC) nhằm tạo thế ngoài biến thiên mạnh ở thang nguyên tử. Qua đó, xác định các cấu trúc siêu mạng có khả năng thúc đẩy sự tán xạ giữa các nón Dirac không tương đương (bẻ gãy đối xứng chiral), làm rõ nguồn gốc sự trồi lên của các nón Dirac bản sao và mô tả quy luật mở vùng cấm năng lượng tại điểm Dirac sơ cấp (PDP) cũng như các điểm Dirac thứ cấp (SDP).
  • Vấn đề nghiên cứu thứ hai: Khảo sát toàn diện tính chất tiếp giáp Schottky giữa graphene và đế bán dẫn silicon (G/Si) có và không có pha tạp dựa trên các phép tính lượng tử. Xác định các điều kiện cấu trúc để duy trì liên kết van der Waals (vdW) nhằm bảo vệ tính bán kim loại của graphene; phân tích cơ chế truyền điện tích, phân bố điện tích và các thành phần cấu thành bậc thế tiếp giáp (gồm thế truyền điện tích và thế hiệu ứng đẩy ngược); khảo sát sự hình thành vùng nghèo hạt tải trong đế và làm rõ độ nhạy của hàng rào Schottky ($S$) dưới ảnh hưởng của các trạng thái bề mặt bán dẫn.

Đối tượng nghiên cứu của luận án là cấu trúc điện tử của graphene cô lập, các cấu trúc siêu mạng graphene trên đế silicon carbide (G/SiC) và tiếp giáp dị thể giữa graphene và đế silicon (G/Si). Phạm vi nghiên cứu giới hạn ở các mô phỏng vi mô từ nguyên lý ban đầu (ab initio) bằng lý thuyết phiếm hàm mật độ (DFT), khảo sát các tính chất tĩnh, cấu trúc dải năng lượng, mật độ trạng thái và phân bố thế tĩnh điện ở thang nguyên tử.

Tổng quan tài liệu và vị trí của luận án

Luận án thực hiện điểm luận có hệ thống các hướng nghiên cứu lý thuyết và thực nghiệm trước đây về graphene và tiếp giáp kim loại/bán dẫn:

Các hướng nghiên cứu về mở vùng cấm và tái cấu trúc dải năng lượng graphene

  • Bẻ gãy đối xứng mạng phụ (đối xứng nghịch đảo): Các nghiên cứu trên hệ tiếp giáp graphene/hBN (Giovannetti và cộng sự, 2007; Hunt và cộng sự, 2013; Woods và cộng sự, 2014) chỉ ra rằng thế ngoài từ đế boron nitride lục giác (hBN) phá vỡ tính tương đương giữa hai mạng phụ carbon $A$ và $B$, từ đó mở ra một vùng cấm năng lượng tại điểm Dirac sơ cấp. Các quan sát trên siêu mạng graphene trên đế SiC (Varchon và cộng sự, 2007; Ristein và cộng sự, 2012) cũng ghi nhận hiện tượng tương tự.
  • Bẻ gãy đối xứng chiral qua siêu mạng Kekulé: Đối với siêu mạng graphene biến dạng dạng Kekulé tuần hoàn ngắn ($G-\sqrt{3}\times\sqrt{3}-R30^\circ$), các công trình của Cheianov (2009), Farjam (2009), và Giovannetti (2007) phân loại thành cấu trúc Kekulé-O (thế nhiễu loạn tại tâm vòng lục giác) và Kekulé-Y (thế nhiễu loạn tại vị trí nguyên tử carbon). Kết quả lý thuyết chỉ ra rằng chỉ có cấu trúc Kekulé-O mới tạo ra sự tán xạ ngược ngoại hồ (intervalley backscattering) giữa hai hồ $K$ và $K'$ tại điểm $\Gamma$, qua đó bẻ gãy đối xứng chiral và mở vùng cấm năng lượng.
  • Điểm Dirac thứ cấp (SDP): Các nghiên cứu của Park (2008), Bliokh (2009), và Wallbank (2013) khảo sát siêu mạng graphene dưới thế tuần hoàn biến thiên chậm ($L > 10 \text{ nm}$), ghi nhận sự xuất hiện của các nón Dirac thứ cấp tại biên vùng Brillouin siêu mạng do tán xạ nội hồ (intravalley scattering). Các công trình thực nghiệm sau đó của Yankowitz (2012) và Ponomarenko (2013) quan sát thấy SDP trên hệ graphene/hBN. Một số nghiên cứu (Song và cộng sự, 2015; Jung và cộng sự, 2015) đề xuất rằng vùng cấm tại PDP và SDP có cùng bản chất bẻ gãy đối xứng mạng phụ, trong khi Chen và cộng sự (2014) lại cho thấy sự mở vùng cấm tại hai điểm này có quy luật biến đổi khác nhau theo áp suất và góc xoay.

Các mô hình tiếp giáp kim loại/bán dẫn và tiếp giáp van der Waals hai chiều

  • Mô hình trạng thái bề mặt định xứ: Bardeen (1947), Cowley và Sze (1965) đề xuất mô hình giải thích sự sai lệch của độ cao hàng rào Schottky (SBH) so với quy tắc Schottky-Mott thông qua sự bẫy điện tích tại các trạng thái bề mặt định xứ trong vùng cấm và sự ghìm mức Fermi tại mức trung tính điện ($E_{CNL}$).
  • Mô hình trạng thái vùng cấm cảm ứng kim loại (MIGS): Heine (1965) và Tersoff (1984) chứng minh rằng đuôi hàm sóng kim loại xuyên sâu vào vùng cấm bán dẫn tạo ra các trạng thái cảm ứng ảo (VIGS), ghìm mức Fermi kim loại vào điểm quay lui ($E_B$) của bán dẫn.
  • Mô hình phân cực liên kết: Tung (2001, 2014) thiết lập lý thuyết phân cực liên kết hóa học tại vùng riêng tiếp giáp (ISR), chỉ ra rằng dipole tạo bởi các liên kết hóa học cục bộ chi phối bậc thế tiếp giáp.
  • Tiếp giáp van der Waals hai chiều: Liu và cộng sự (2018) trên tiếp giáp vdW kim loại/$\text{MoS}_2$ đạt độ nhạy $S = 0{,}96$; LaGasse và cộng sự (2019) trên tiếp giáp $\text{graphene/WSe}_2$ đạt $S = 1{,}04$. Yang và cộng sự (2012) chế tạo linh kiện barristor dựa trên tiếp giáp G/Si, điều khiển SBH thông qua điện thế cực cổng nhờ sự dịch chuyển mức Fermi $\Delta E_F$ của graphene.
  • Cấu trúc tiếp giáp G/Si: Koepke và cộng sự (2010) kết hợp thực nghiệm và DFT trên bề mặt $\text{Si}(111)-7\times7$; Li và cộng sự (2014) khảo sát bề mặt $\text{Si}(111)-4\times4$ pha tạp bổ sung nguyên tử Si.
Mô hình tiếp giáp Cơ chế truyền điện tích chính Yếu tố quy định vị trí ghìm mức Fermi Độ nhạy hàng rào ($S$)
Bardeen / Cowley-Sze Trao đổi điện tích với trạng thái bề mặt Mức trung tính điện ($E_{CNL}$) trong vùng cấm $S \to 0$ (Giới hạn Bardeen)
MIGS (Heine / Tersoff) Đuôi hàm sóng kim loại xuyên vào bán dẫn Điểm quay lui ($E_B$) của cấu trúc khối $S \ll 1$ (Phụ thuộc khối bán dẫn)
Phân cực liên kết (Tung) Dipole của liên kết hóa học tại vùng ISR Độ âm điện và cấu trúc liên kết bề mặt $S < 1$ (Do màn chắn điện môi ISR)
Tiếp giáp vdW 2D Không có liên kết hóa học, triệt tiêu MIGS Không bị ghìm (Công thoát điện cực) $S \to 1$ (Giới hạn Schottky-Mott)

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

  1. Chưa có nghiên cứu hệ thống về cơ chế bẻ gãy đối xứng chiral và tái cấu trúc dải năng lượng trong vùng nón Dirac bất đẳng hướng (biến dạng tam giác - trigonal warping) trên các siêu mạng graphene tuần hoàn tầm ngắn ($L$ cỡ vài nanomet) có chu kỳ lớn hơn mạng Kekulé ($N > 1$).
  2. Chưa phân định rõ ràng nguồn gốc vật lý của việc mở vùng cấm tại PDP và SDP khi siêu mạng chịu tác động của thế ngoài sắc nhọn thang nguyên tử.
  3. Trong tiếp giáp Schottky graphene/bán dẫn ba chiều (đại diện là G/Si), vai trò của hiệu ứng đẩy ngược (push-back effect), động thái hình thành vùng nghèo hạt tải trong đế và ảnh hưởng của các loại trạng thái bề mặt khác nhau đối với độ nhạy hàng rào Schottky chưa được mô tả định lượng đầy đủ.

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

Luận án sử dụng phương pháp tính toán lượng tử từ nguyên lý ban đầu dựa trên Lý thuyết phiếm hàm mật độ (DFT) để khảo sát các hệ vật liệu.

+-----------------------------------------------------------------------+
|                    LÝ THUYẾT PHIẾM HÀM MẬT ĐỘ (DFT)                   |
|  - Định lý Hohenberg-Kohn 1: n_0(r) xác định duy nhất thế ngoài V_ext |
|  - Định lý Hohenberg-Kohn 2: Nguyên lý biến phân năng lượng E[n(r)]   |
+-----------------------------------------------------------------------+
|                       PHƯƠNG TRÌNH KOHN-SHAM                          |
|    Giải hệ đơn hạt tương đương trong trường tự hợp hiệu dụng:         |
|    V_eff(r) = V_ext(r) + V_H(r) + V_xc(r)                             |
+-----------------------------------------------------------------------+
|                     CÔNG CỤ VÀ THAM SỐ TÍNH TOÁN                      |
|  - Phần mềm: VASP (Vienna Ab-initio Simulation Package)               |
|  - Giả thế: PAW (Projector Augmented-Wave)                            |
|  - Phiếm hàm: LDA, GGA-PBE, và optB86b-vdW (hiệu chỉnh tương tác vdW) |
|  - Kỹ thuật: Trải cấu trúc dải năng lượng (Band Unfolding)            |
+-----------------------------------------------------------------------+

Các định lý nền tảng và phương trình Kohn-Sham

  • Định lý Hohenberg-Kohn thứ nhất: Mật độ điện tích trạng thái cơ bản $n_0(\mathbf{r})$ của hệ electron tương tác chuyển động dưới tác dụng của thế ngoài $V_{ext}(\mathbf{r})$ là một hàm duy nhất xác định thế ngoài (sai khác một hằng số cộng). Hệ quả là Hamiltonian và toàn bộ các tính chất của hệ ở trạng thái cơ bản đều được xác định qua $n_0(\mathbf{r})$.
  • Định lý Hohenberg-Kohn thứ hai: Với mật độ điện tử thử $n(\mathbf{r})$ bất kỳ thỏa mãn điều kiện chuẩn hóa số hạt, phiếm hàm năng lượng toàn phần $E[n(\mathbf{r})]$ luôn đạt giá trị cực tiểu tại mật độ điện tích trạng thái cơ bản $n_0(\mathbf{r})$.
  • Phương trình Kohn-Sham: Hệ $N$ electron tương tác được quy về hệ giả định gồm $N$ electron không tương tác chuyển động trong thế hiệu dụng $V_{eff}(\mathbf{r})$: $$\left[ -\frac{\hbar^2}{2m_0}\nabla^2 + V_{eff}(\mathbf{r}) \right] \psi_i(\mathbf{r}) = \epsilon_i \psi_i(\mathbf{r})$$ Trong đó, thế hiệu dụng bao gồm thế ngoài $V_{ext}(\mathbf{r})$, thế tĩnh điện Hartree $V_H(\mathbf{r})$ và thế trao đổi - tương quan $V_{xc}(\mathbf{r}) = \frac{\delta E_{xc}[n]}{\delta n(\mathbf{r})}$. Mật độ điện tích được tính qua phép lấy tổng các orbital bị chiếm: $n(\mathbf{r}) = \sum_{i=1}^{N_{occ}} |\psi_i(\mathbf{r})|^2$.

Phép xấp xỉ trao đổi - tương quan và hiệu chỉnh van der Waals

Luận án áp dụng các phép xấp xỉ trao đổi - tương quan:

  • Xấp xỉ mật độ địa phương (LDA).
  • Xấp xỉ gradient tổng quát (GGA) với tham số hóa Perdew-Burke-Ernzerhof (PBE).
  • Phiếm hàm phi địa phương có hiệu chỉnh tương tác phân tán van der Waals: optB86b-vdW, nhằm mô tả chính xác lực hút liên kết yếu giữa lớp graphene và bề mặt đế bán dẫn.

Kỹ thuật trải cấu trúc dải năng lượng (Band Unfolding)

Khi xây dựng mô hình siêu mạng, vùng Brillouin của siêu mạng (SBZ) bị thu nhỏ và gấp lại so với vùng Brillouin nguyên tố của graphene (GBZ). Luận án áp dụng kỹ thuật trải dải năng lượng (band unfolding) để chiếu các trạng thái năng lượng từ SBZ trở lại GBZ, cho phép phân tích trực quan mức độ trồi lên của các nón Dirac bản sao và quá trình tán xạ giữa các nón Dirac.

Công cụ phần mềm và cơ sở hạ tầng tính toán

  • Phần mềm: Gói phần mềm tính toán lượng tử VASP (Vienna Ab-initio Simulation Package) sử dụng sóng phẳng và phương pháp sóng phẳng tăng cường gắn kèm máy chiếu (PAW).
  • Hạ tầng tính toán: Hệ thống tính toán hiệu năng cao tại INOMAR (ĐHQG-HCM), Phòng Tính toán Hiệu năng cao (Trường ĐH Bách Khoa - ĐHQG-HCM), Phòng thí nghiệm trọng điểm Khoa học tính toán đa tỉ lệ (Trường ĐH KHTN - ĐHQG-Hà Nội), CM Thi Lab (Trường ĐH HUTECH) và Trung tâm Khoa học Vật liệu Tính toán thuộc Viện Nghiên cứu Vật liệu (IMR), Đại học Tohoku, Nhật Bản.

Nội dung chính theo từng chương

Chương 1: Các tính chất điện tử của graphene

Chương 1 trình bày lý thuyết giải tích về cấu trúc tinh thể và cấu trúc điện tử của graphene đơn lớp cô lập.

  • Mô hình liên kết mạnh: Khảo sát mạng Bravais lục giác với hai nguyên tử carbon thuộc hai mạng phụ $A$ và $B$ trong một ô nguyên tố ($a_0 = 1{,}42 \text{ \AA}$, $a = \sqrt{3}a_0 = 2{,}46 \text{ \AA}$). Xét tương tác nhảy bậc nhất giữa các orbital $2p_z$ của các nguyên tử lân cận gần nhất với tích phân nhảy $t \approx 2{,}7 \text{ eV}$.
  • Tán sắc tuyến tính và giả spin: Khai triển năng lượng tại lân cận các điểm $K(K')$ dẫn đến phổ tán sắc tuyến tính dạng nón Dirac: $E(\mathbf{q}) = \pm \hbar v_F |\mathbf{q}|$. Hamiltonian hiệu dụng tại các hồ có dạng phương trình Dirac hai chiều $H_\kappa = \hbar v_F (\kappa \sigma_x q_x + \sigma_y q_y)$ với $\kappa = \pm 1$ và $\sigma_i$ là các ma trận Pauli. Hàm sóng mang tính chất chiral với vector giả spin hướng dọc hoặc ngược chiều vector sóng $\mathbf{q}$.
  • Mật độ trạng thái và công thoát: Mật độ trạng thái phụ thuộc tuyến tính vào năng lượng $\text{DOS}(E) = \frac{2|E|}{\pi (\hbar v_F)^2}$ và triệt tiêu tại $E = 0$. Mối liên hệ giữa năng lượng Fermi và nồng độ hạt tải điện $E_F(n) = \hbar v_F \sqrt{\pi n}$ giải thích tính nhạy cảm của mức Fermi và công thoát $\Phi_G$ trước sự pha tạp hay dịch chuyển điện tích.
  • Điều kiện mở vùng cấm: Trình bày lý thuyết về việc phá vỡ 4 trạng thái suy biến tại PDP qua hai cơ chế: (1) bẻ gãy đối xứng mạng phụ bằng cách trộn các dải $\pi(K)-\pi^*(K)$; (2) bẻ gãy đối xứng chiral bằng cách trộn các góc $K-K'$ thông qua cấu trúc siêu mạng Kekulé-O.

Chương 2: Các mô hình tiếp giáp kim loại/bán dẫn

Chương 2 phân tích các mô hình lý thuyết kinh điển và hiện đại về tiếp giáp dị thể giữa kim loại và bán dẫn (M/S).

  • Phân loại trạng thái bề mặt: Phân tích ba dạng hàm sóng trên bề mặt tinh thể bán vô hạn cắt ngang: trạng thái khối, trạng thái bề mặt định xứ (nằm trong vùng cấm) và trạng thái cộng hưởng bề mặt (nằm xen phủ trong dải khối).
  • Độ nhạy hàng rào Schottky ($S$): Xuất phát từ công thức $\Phi_{Bn} = \Phi_m - \chi_s - \Delta_{in}$, độ nhạy tiếp giáp được định nghĩa $S = \partial \Phi_{Bn}/\partial \Phi_m = 1 - \partial \Delta_{in}/\partial \Phi_m$. Khi $\partial \Delta_{in}/\partial \Phi_m = 0$, $S = 1$ (giới hạn Schottky-Mott); khi $\partial \Delta_{in}/\partial \Phi_m = 1$, $S = 0$ (giới hạn Bardeen).
  • Nguồn gốc bậc thế tiếp giáp $\Delta_{in}$:
    • Hiệu ứng đẩy ngược (push-back effect): Lực đẩy trao đổi Pauli nén đuôi đám mây điện tích tràn ra ngoài bề mặt kim loại khi tiếp xúc cơ học, làm giảm công thoát kim loại mà không kèm theo sự truyền điện tích ròng.
    • Thế truyền điện tích: Gây ra bởi sự dịch chuyển điện tích giữa kim loại với trạng thái bề mặt bán dẫn (mô hình Bardeen/Cowley-Sze), với đuôi hàm sóng kim loại xuyên sâu vào vùng cấm (mô hình MIGS của Heine/Tersoff), hoặc qua phân cực liên kết hóa học cục bộ trong vùng ISR (mô hình Tung).
  • Tiếp giáp van der Waals hai chiều: Phân tích cơ chế tiếp cận giới hạn Schottky-Mott ($S \to 1$) trong tiếp giáp vdW nhờ sự triệt tiêu các liên kết hóa học bề mặt, giảm thiểu trạng thái bẫy và tiêu tán MIGS.

Chương 3: Lý thuyết phiếm hàm mật độ

Chương 3 thiết lập khung phương pháp luận tính toán được sử dụng trong luận án:

  • Cơ sở toán học DFT: Trình bày chi tiết phép biến đổi bài toán nhiều hạt $3N$ bậc tự do về bài toán mật độ điện tích $n(\mathbf{r})$ ba chiều thông qua hệ phương trình tự hợp Kohn-Sham.
  • Phiếm hàm và hiệu chỉnh thế: Trình bày công thức toán học của các phép xấp xỉ LDA, GGA (PBE) và cơ chế tính toán năng lượng phi địa phương của phiếm hàm optB86b-vdW.
  • Phương pháp giả thế: Phân tích nguyên lý thay thế thế Coulomb kỳ dị tại vùng lõi nguyên tử bằng thế giả trơn ($V_{ps}$) và hàm sóng giả ($\tilde{\psi}$), giảm thiểu số lượng hàm sóng phẳng cơ sở.
  • Kỹ thuật tính toán trên phần mềm VASP: Chi tiết hóa quy trình tự hợp SCF, tiêu chuẩn hội tụ năng lượng và thuật toán trải dải năng lượng từ SBZ về PBZ.

Chương 4: Điều khiển sự trồi lên và tán xạ của các nón Dirac của siêu mạng graphene trên đế silicon carbide

Chương 4 giải quyết vấn đề nghiên cứu thứ nhất của luận án:

  • Mô hình tính toán: Xây dựng các mô hình tiếp giáp siêu mạng graphene tuần hoàn tầm ngắn trên bề mặt SiC chứa các liên kết chưa bão hòa:
    • $G-5\sqrt{3}\times5\sqrt{3}-R30^\circ / \text{SiC}-4\sqrt{3}\times4\sqrt{3}-R30^\circ$
    • $G-9\times9$
    • $G-\sqrt{21}\times\sqrt{21}-R10{,}893^\circ$
    • $G-7\times7$
    • $G-5\times5$
  • Cơ chế trồi lên của các nón Dirac: Kết quả trải dải năng lượng chứng minh rằng thế ngoài sắc nhọn thang nguyên tử từ các liên kết chưa bão hòa của đế SiC kích hoạt sự xuất hiện của các nón Dirac bản sao tương đương (tán xạ nội hồ) và không tương đương (tán xạ ngoại hồ).
  • Quy luật mở vùng cấm tại PDP và SDP:
    • Tại PDP: Sự mở vùng cấm xuất hiện khi tâm $\Gamma$ của SBZ gấp trùng vào điểm $K(K')$ của GBZ (các siêu mạng dạng $G-\sqrt{3N}\times\sqrt{3N}-R30^\circ$), cho phép các nón Dirac không tương đương trộn lẫn và tán xạ ngược, bẻ gãy đối xứng chiral.
    • Tại SDP: Luận án làm rõ nguồn gốc sự mở vùng cấm tại các điểm Dirac thứ cấp ở dải hóa trị và dải dẫn là do sự kết hợp giữa tán xạ nội hồ và sự bẻ gãy đối xứng mạng phụ, phân định rõ ràng sự khác biệt bản chất so với cơ chế mở vùng cấm tại PDP.
                    MÔ HÌNH TÍNH TOÁN CHƯƠNG 4 & CHƯƠNG 5
                    
  [Chương 4: Siêu mạng G/SiC]                    [Chương 5: Tiếp giáp G/Si]
  +-------------------------+                    +-------------------------+
  |    Lớp đơn Graphene     |                    |    Lớp đơn Graphene     |
  +-------------------------+                    +-------------------------+
  : Thế nguyên tử sắc nhọn  : (Tán xạ nón Dirac) :  Tương tác van der Waals: (Hiệu ứng đẩy ngược)
  +-------------------------+                    +-------------------------+
  |  Đế SiC (Dangling bonds)|                    |  Đế Si(111) đã xử lý    |
  +-------------------------+                    +-------------------------+

Chương 5: Tiếp giáp Schottky của graphene/silicon

Chương 5 giải quyết vấn đề nghiên cứu thứ hai của luận án:

  • Cấu trúc hình học tiếp giáp: Khảo sát mô hình $G-\sqrt{39}\times\sqrt{39}-R16{,}1^\circ / \text{Si}(111)-4\times4$. Kết quả tối ưu hóa hình học xác nhận rằng graphene duy trì hình thái phẳng ổn định liên kết vdW với bề mặt $\text{Si}(111)$ khi các liên kết bề mặt được bão hòa thích hợp (ví dụ bề mặt tái cấu trúc có pha tạp hoặc hydrogen hóa $\text{Si}-\text{H}$), ngăn chặn sự hình thành liên kết cộng hóa trị cục bộ làm biến dạng mạng graphene.
  • Cấu trúc điện tử và cơ chế truyền điện tích: Tính toán mật độ trạng thái hình chiếu theo lớp cho thấy nón Dirac của graphene được bảo toàn nguyên vẹn tại mức Fermi chung $E_F$. Mức Fermi của graphene bị dịch chuyển ($\Delta E_F$) do trao đổi điện tích với bề mặt silicon.
  • Phân tích bậc thế tiếp giáp $\Delta_{in}$: Luận án phân tách chính xác hai thành phần độc lập trong bậc thế tiếp giáp:
    1. Thế hiệu ứng đẩy ngược ($\Delta_{pb}$): Luôn xuất hiện do tương tác cơ học vdW giữa đám mây điện tích graphene và đế silicon, làm giảm hàm công thoát hiệu dụng.
    2. Thế truyền điện tích ($\Delta_{ct}$): Xuất phát từ sự dịch chuyển điện tích ròng giữa graphene và đế bán dẫn, phụ thuộc trực tiếp vào loại và nồng độ pha tạp của silicon.
  • Sự hình thành vùng nghèo hạt tải và độ nhạy $S$: Khác với các tiếp giáp hai chiều thuần túy, tiếp giáp G/Si hình thành vùng nghèo hạt tải điện tự do rõ rệt ở các lớp nguyên tử silicon dưới bề mặt. Độ nhạy hàng rào Schottky $S$ được cải thiện đáng kể so với tiếp giáp kim loại khối truyền thống nhưng vẫn chịu sự điều biến bởi các trạng thái bề mặt còn sót lại của silicon.

Kết quả và những đóng góp mới

  1. Về lý luận cấu trúc dải năng lượng graphene:
    • Xây dựng mô hình định lượng giải thích cơ chế trồi lên của các nón Dirac bản sao và quá trình tán xạ nội hồ/ngoại hồ trong các siêu mạng graphene tuần hoàn tầm ngắn chịu tác động của thế ngoài sắc nhọn thang nguyên tử.
    • Chứng minh bằng giải tích cấu trúc và tính toán DFT rằng việc mở vùng cấm tại điểm Dirac sơ cấp (PDP) đòi hỏi sự bẻ gãy đối xứng chiral thông qua việc gấp vùng SBZ đưa các nón không tương đương về tâm $\Gamma$, trong khi việc mở vùng cấm tại điểm Dirac thứ cấp (SDP) chủ yếu bắt nguồn từ sự giao thoa dải và bẻ gãy đối xứng mạng phụ.
  2. Về mô hình tiếp giáp Schottky graphene/bán dẫn:
    • Xác định điều kiện cấu trúc bề mặt để graphene duy trì liên kết van der Waals phẳng và bảo toàn tính chất điện tử nguyên thủy trên bề mặt silicon $\text{Si}(111)$.
    • Thiết lập mô hình bậc thế tiếp giáp tổng quát cho tiếp giáp vdW dị thể 2D/3D, tích hợp đồng thời cả thế hiệu ứng đẩy ngược ($\Delta_{pb}$) và thế truyền điện tích ($\Delta_{ct}$).
    • Làm rõ cơ chế hình thành vùng nghèo hạt tải trong đế bán dẫn 3D khi tiếp xúc với vật liệu 2D, chỉ ra mối liên hệ giữa sự dịch mức Fermi $\Delta E_F$ của graphene và độ biến thiên chiều cao hàng rào Schottky $\Delta \Phi_{Bn}$.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Các mô hình tính toán DFT trong luận án được thực hiện ở nhiệt độ 0 Kelvin ($T = 0 \text{ K}$) và giả định các cấu trúc bề mặt hoàn hảo theo chu kỳ siêu mạng. Chưa xét đến ảnh hưởng của dao động mạng (phonon) ở nhiệt độ phòng, cũng như các khuyết tật mạng ngẫu nhiên trong điều kiện chế tạo thực nghiệm.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng tính toán DFT có xét đến hiệu ứng tán xạ electron-phonon ở nhiệt độ phòng trên tiếp giáp G/Si và G/SiC.
    • Khảo sát động học truyền hạt tải điện qua hàng rào Schottky dưới tác dụng của điện trường ngoài phụ thuộc thời gian.
    • Nghiên cứu ứng dụng các cấu trúc siêu mạng tuần hoàn ngắn trong việc thiết kế các linh kiện thung lũng điện tử (valleytronics).

Giá trị tham khảo

Luận án là tài liệu tham khảo cho các đối tượng sau:

  • Nghiên cứu sinh và học viên cao học ngành Khoa học Vật liệu, Vật lý Chất rắn, Vật lý Tính toán: Có thể tham khảo phương pháp luận DFT, kỹ thuật trải dải năng lượng (band unfolding) trong VASP và phương pháp phân tích điện tích bề mặt.
  • Các nhóm nghiên cứu vật liệu bán dẫn và linh kiện nano: Cung cấp thông số định lượng về cấu trúc tiếp giáp, khoảng cách liên kết, thế tiếp giáp và cơ chế điều khiển hàng rào Schottky cho các thiết bị logic dựa trên graphene tích hợp trên silicon.

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

1. Tại sao graphene cô lập không thể ứng dụng trực tiếp trong các linh kiện logic số?

Do cấu trúc vùng năng lượng của graphene có vùng dẫn và vùng hóa trị chạm nhau tại các điểm Dirac $K$ và $K'$, tạo nên vùng cấm năng lượng bằng không ($E_g = 0$). Điều này khiến thiết bị không thể chuyển hoàn toàn sang trạng thái ngắt dòng (OFF state), dẫn đến tỷ số dòng ON/OFF thấp.

2. Sự khác nhau giữa cơ chế bẻ gãy đối xứng mạng phụ và bẻ gãy đối xứng chiral trong việc mở vùng cấm graphene là gì?

  • Bẻ gãy đối xứng mạng phụ (nghịch đảo): Làm mất tính tương đương giữa hai mạng con carbon $A$ và $B$ (ví dụ: do thế ngoài tác dụng khác nhau lên hai vị trí), gây trộn dải $\pi(K)$ và $\pi^*(K)$ trong cùng một hồ.
  • Bẻ gãy đối xứng chiral: Đòi hỏi sự tán xạ ngược ngoại hồ kết hợp hai hồ không tương đương $K$ và $K'$ vào cùng một điểm (tâm $\Gamma$ của SBZ), chỉ xảy ra trên các cấu trúc siêu mạng có chu kỳ thích hợp (như cấu trúc kiểu Kekulé-O).

3. Hiệu ứng đẩy ngược (push-back effect) ảnh hưởng thế nào đến bậc thế tiếp giáp?

Khi graphene tiếp xúc vật lý vdW với bề mặt đế, lực đẩy trao đổi Pauli nén phần đuôi hàm sóng điện tử tràn ra ngoài bề mặt trở lại phía trong. Sự tái phân bố điện tích này tạo ra một dipole bề mặt gọi là thế đẩy ngược ($\Delta_{pb}$), làm giảm công thoát hiệu dụng của tiếp giáp mà không cần có sự truyền điện tích ròng giữa hai vật liệu.

4. Tại sao độ cao hàng rào Schottky trong tiếp giáp graphene/silicon có thể điều chỉnh được bằng điện trường ngoài?

Mật độ trạng thái của graphene rất nhỏ tại lân cận điểm Dirac ($\text{DOS} \to 0$). Khi có sự trao đổi điện tích hoặc tác động của điện áp cổng, mức Fermi của graphene bị dịch chuyển một lượng $\Delta E_F$ đáng kể ($E_F \propto \sqrt{n}$), dẫn đến sự thay đổi trực tiếp của công thoát $\Phi_G$ và điều biến độ cao hàng rào Schottky $\Delta \Phi_{Bn} \approx \Delta E_F$.

Kết luận

Luận án của nghiên cứu sinh Vũ Hoàng Nam đã thực hiện nghiên cứu chuyên sâu về tương tác giữa graphene và các đế bán dẫn (SiC và Si) bằng phương pháp lý thuyết phiếm hàm mật độ (DFT). Công trình đã làm sáng tỏ cơ chế điều khiển sự trồi lên của các nón Dirac và quy luật mở vùng cấm năng lượng tại điểm Dirac sơ cấp cũng như thứ cấp trên các siêu mạng graphene tuần hoàn tầm ngắn. Đồng thời, luận án đã xây dựng mô hình tiếp giáp Schottky G/Si hoàn chỉnh, phân tích định lượng bậc thế tiếp giáp gồm thế đẩy ngược và thế truyền điện tích, cung cấp cơ sở vật lý vi mô cho việc tích hợp graphene vào công nghệ linh kiện bán dẫn.