Tổng quan về luận án

Sự phát triển của công nghiệp vi điện tử và mạch tích hợp bán dẫn (Integrated Circuits - IC) dựa trên nền tảng silicon (Si) đang dần tiếp cận các giới hạn vật lý lượng tử khắc nghiệt khi kích thước kênh dẫn của transistor hiệu ứng trường kim loại - oxit - bán dẫn (MOSFET) thu nhỏ xuống dưới ngưỡng 10 nm. Trong kỷ nguyên hậu định luật Moore (Post-Moore Era), các hiệu ứng ký sinh bất lợi như hiệu ứng kênh ngắn (short-channel effects), hiện tượng dòng rò do chui ngầm lượng tử qua lớp điện môi cực cổng mỏng, và sự suy giảm độ linh động của hạt tải đã đặt ra yêu cầu cấp bách phải tìm kiếm các cấu trúc và vật liệu dẫn kênh thế hệ mới. Luận án Tiến sĩ Vật lý kỹ thuật (Mã số: 62520401) của nghiên cứu sinh Lê Hoàng Anh, dưới sự hướng dẫn của GS. Nguyễn Đức Chiến và hợp tác phát triển cùng nhóm nghiên cứu của TS. Đỗ Vân Nam tại Trường Đại học Bách Khoa Hà Nội (2015), mang tên: "Nghiên cứu các tính chất điện tử, quang học và truyền dẫn của vật liệu graphene hướng tới các ứng dụng điện tử và quang điện tử".

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án giải quyết xuất phát từ việc: mặc dù graphene sở hữu độ linh động điện tử nội tại cực lớn ($\mu_G \sim 2 \cdot 10^5\text{ cm}^2/\text{V}\cdot\text{s}$) và hệ thức tán sắc tuyến tính dạng nón Dirac ($E = \pm \hbar v_F |k|$ với vận tốc Fermi $v_F \approx 10^6\text{ m/s}$, chỉ nhỏ hơn vận tốc ánh sáng 300 lần), bản chất không có vùng cấm (zero-bandgap) và nghịch lý chui ngầm Klein (Klein tunneling) khiến dòng điện không thể đóng ngắt hoàn toàn, dẫn đến tỷ số dòng bật/tắt ($I_{ON}/I_{OFF}$) rất thấp trong các transistor hiệu ứng trường graphene (GFETs). Đồng thời, các cơ chế vi mô chi phối sự truyền dẫn hạt tải qua tiếp xúc kim loại - graphene ($M-G$), sự biến đổi phổ cấu trúc dải năng lượng và độ dẫn quang học dưới tác động của trường thế tuần hoàn một chiều (siêu mạng graphene - GSLs) vẫn còn nhiều điểm chưa được làm sáng tỏ tường minh bằng các mô hình tính toán lượng tử vi mô.

Luận án đặt ra và giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng ba giả thuyết khoa học (Hypotheses - H):

  1. RQ1: Cấu trúc dải năng lượng và các trạng thái sóng điện tử $2p_z$ biến đổi như thế nào khi graphene chịu tác động của hàm thế vô hướng tuần hoàn một chiều theo các hướng đối xứng mạng khác nhau (Armchair và Zigzag)?
    • H1: Trường thế tuần hoàn phá vỡ đối xứng tinh thể nguyên thủy $D_{6h}$, tạo ra các điểm Dirac thứ cấp và gây nên sự định xứ dị thường của các trạng thái điện tử ngay trong dải phổ liên tục.
  2. RQ2: Tương tác dị thể giữa các phân lớp nguyên tử kim loại điện cực và màng carbon đơn lớp graphene làm biến đổi phổ tán sắc và điện trở tiếp xúc nội tại theo cơ chế nào?
    • H2: Sự xen phủ orbital giữa điện tử $d$ (hoặc $s$) của kim loại với orbital $\pi$ ($2p_z$) của graphene xác định độ dịch chuyển mức Fermi và điện trở tiếp xúc thông qua hàm tự năng (self-energy) tương tác.
  3. RQ3: Các yếu tố tiếp xúc thực tế và cấu trúc cực cổng tự sắp xếp (self-aligned) dây nano GaN ảnh hưởng như thế nào đến độ hỗ dẫn ($g_m$), trở vi phân âm và tần số cắt ($f_T$) của linh kiện GFETs ở quy mô dưới 100 nm?
    • H3: Bằng cách tích hợp rào thế Schottky GaN-graphene với hằng số điện môi cao ($\kappa \sim 10$) và mô hình hóa truyền dẫn lượng tử không cân bằng, tần số làm việc của GFETs có thể đạt ngưỡng terahertz (THz).

Khung lý thuyết của luận án tích hợp phương pháp gần đúng liên kết chặt (Tight-Binding Approximation - TBA), hình thức luận phản ứng tuyến tính Kubo (Kubo linear response formula) và phương pháp hàm Green không cân bằng (Non-Equilibrium Green's Functions - NEGF) kết hợp tự nhất quán với phương trình Poisson. Đóng góp đột phá của công trình được định lượng qua việc giải thích thành công sự dập tắt độ dẫn quang trong dải năng lượng photon $(0, U_b)$, định lượng điện trở tiếp xúc nội tại cho 5 cặp tổ hợp kim loại - graphene điển hình, và xây dựng hoàn chỉnh module GFETs lượng tử trong gói phần mềm OPEDEVS (Opto-Electronic Devices Simulation).

Literature Review và Positioning

Bức tranh tổng quan y văn về vật liệu graphene bắt đầu từ công trình lý thuyết kinh điển của Wallace (1946), người đầu tiên tính toán cấu trúc dải năng lượng của mạng carbon lục giác đơn lớp và tìm ra mối quan hệ tán sắc tuyến tính quanh hai điểm nón đối xứng $K$ và $K'$. Tuy nhiên, trong suốt nhiều thập kỷ, các định lý nhiệt động lực học cổ điển (như định lý Mermin-Wagner) đã dẫn đến quan niệm cho rằng tinh thể hai chiều không thể tồn tại độc lập ở trạng thái tự do do các thăng giáng nhiệt sẽ phá hủy mạng tinh thể. Quan niệm này bị đảo lộn hoàn toàn vào năm 2004 khi Andre Geim và Kostya Novoselov (Đại học Manchester) sử dụng phương pháp bóc tách cơ học vi mô bằng băng dính graphite để phân lập thành công đơn lớp graphene (Science, 2004), mở ra kỷ nguyên vật liệu 2D và mang lại giải Nobel Vật lý năm 2010.

Trường phái nghiên cứu thực nghiệm sau đó nhanh chóng phân hóa thành hai luồng quan điểm và phương pháp tiếp cận chính:

  1. Phương pháp chế tạo và vật liệu: Luồng nghiên cứu của de Heer et al. (Georgia Tech, 2004) tập trung vào việc nhiệt phân bề mặt tinh thể silicon carbide (SiC) ở nhiệt độ $\sim 1300^\circ\text{C}$ để thăng hoa Si và tái cấu trúc màng graphene epitaxy diện tích lớn. Ngược lại, luồng nghiên cứu của Ruoff et al. (Northwestern University, 1999) và Kim et al. (Columbia University, 2004) nỗ lực phát triển kỹ thuật bóc tách bằng đầu dò kính hiển vi lực nguyên tử (AFM) và lắng đọng pha hơi hóa học (CVD) trên đế kim loại chuyển tiếp (Cu, Ni).
  2. Tranh luận về truyền dẫn lượng tử và nghịch lý Klein: Katsnelson, Geim và Novoselov (2006) cùng với Young và Kim (2009) chứng minh rằng các giả hạt Dirac fermion không khối lượng chuyển động tương đối tính có khả năng chui ngầm hoàn hảo qua các rào thế tĩnh điện với xác suất truyền qua bằng 1 ở góc tới vuông góc, bất kể chiều cao và bề dày rào thế. Điều này mâu thuẫn trực tiếp với cơ chế đóng ngắt dòng điện trong transistor bán dẫn truyền thống theo lý thuyết scaling của Robert Dennard (1974). Tworzydło et al. (2006) thiết lập công thức Landauer chỉ ra độ dẫn lượng tử cực tiểu tổng quát $\sigma_{\min} = \frac{4e^2}{\pi h}$ thông qua các trạng thái sóng mờ (evanescent modes), tương ứng với điện trở lượng tử $R_Q = \frac{h}{4e^2} \approx 9\text{ k}\Omega$.

Về cấu trúc siêu mạng graphene (GSLs), các nghiên cứu của Park et al. (2008), Brey và Fertig (2009), cùng Ho et al. (2010) đã khảo sát các dạng thế tuần hoàn Kronig-Penney, cosin và muffin-tin, chỉ ra sự xuất hiện dị hướng của vận tốc nhóm và các điểm Dirac thế hệ mới. Luận án của Lê Hoàng Anh định vị chính xác khoảng trống học thuật bằng cách kết nối trực tiếp cấu trúc điện tử vi mô với tính chất hấp thụ quang phân cực thông qua hình thức luận Kubo, đồng thời giải quyết bài toán tiếp xúc bề mặt kim loại - graphene vốn bị giản lược quá mức trong các nghiên cứu trước đó.

Khi so sánh với hai nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của IBM (Phaedon Avouris, Meric et al., 2008-2012): Chế tạo GFETs trên graphene bóc tách và epitaxy đạt tần số cắt $f_T = 26 - 100\text{ GHz}$, nhưng giải thích đặc tuyến $I_{DS}-V_{DS}$ dựa trên mô hình bán cổ điển, chưa làm rõ vai trò của hàm tự năng tại mặt tiếp xúc kim loại.
  • Nghiên cứu của Đại học California, Los Angeles (Lei Liao et al., Nature 2010): Chế tạo thành công GFETs tự sắp xếp dùng dây nano GaN đạt độ hỗ dẫn $g_m > 2000\text{ S/m}$ và $f_T \sim 100-300\text{ GHz}$. Luận án của Lê Hoàng Anh đã phát triển mô hình mô phỏng NEGF lượng tử để tái hiện và giải thích tường minh cơ chế hình thành rào Schottky, mật độ hạt tải và sự dịch chuyển mức Fermi trong cấu trúc thực nghiệm này.

Đó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à làm sâu sắc thêm ba lý thuyết nền tảng của vật lý chất rắn và vật lý bán dẫn hiện đại:

  1. Lý thuyết vùng năng lượng dải Dirac: Mở rộng phương trình Dirac hai chiều cho các hệ chịu thế điều biến tuần hoàn 1D. Chứng minh rằng cấu trúc vùng năng lượng của siêu mạng A-GSLs và Z-GSLs là sự chồng chập và tái chuẩn hóa của các dải nano graphene (GNRs) tương ứng giữa miền rào thế và giếng thế, bị dịch chuyển năng lượng một khoảng đúng bằng biên độ thế $U_b$.
  2. Lý thuyết phản ứng tuyến tính và chuyển dời quang học: Bổ sung cơ chế chọn lọc chuyển dời quang học phụ thuộc phân cực. Luận án chứng minh sự suy giảm độ dẫn quang học thực phần $\text{Re}[\sigma(\omega)]$ trong khoảng năng lượng photon $\hbar\omega \in (0, U_b)$ xuất phát từ sự triệt tiêu tích phân xen phủ ma trận xung lượng $\langle \psi_c | \hat{p} | \psi_v \rangle$ giữa các hàm sóng điện tử bị biến điệu mạnh và cơ chế khóa Pauli (Pauli blocking).
  3. Lý thuyết tiếp xúc dị thể kim loại - bán dẫn không vùng cấm: Thách thức mô hình rào Schottky truyền thống của Bardeen-Tung khi áp dụng cho màng đơn lớp carbon. Luận án xây dựng mô hình hiệu dụng biểu diễn tương tác tiếp xúc thông qua toán tử tự năng liên kết $\Sigma_{G-M}$, tách biệt rõ đóng góp dịch chuyển mức Fermi ($\text{Re}[\Sigma_{G-M}]$) và sự suy giảm thời gian sống của trạng thái do tán xạ vào điện cực ($\text{Im}[\Sigma_{G-M}]$).
             +-------------------------------------------------------+
             |         Hamiltonian Gần Đúng Liên Kết Chặt            |
             |       H = -t_CC \sum <i,j> (a_i^+ b_j + h.c.)         |
             +-------------------------------------------------------+
                                        |
                  +---------------------+---------------------+
                  |                                           |
                  v                                           v
+-----------------------------------+     +-----------------------------------+
|      Siêu Mạng Graphene (GSLs)    |     |    Tiếp Xúc Kim Loại - Graphene   |
|   - Thế tuần hoàn V(x) = V(x+L)   |     |    - Ghép nối dải pi - orbital M  |
|   - A-GSLs vs. Z-GSLs             |     |    - Tự năng \Sigma_{G-M}         |
+-----------------------------------+     +-----------------------------------+
                  |                                           |
                  v                                           v
+-----------------------------------+     +-----------------------------------+
|   Độ Dẫn Quang (Hình Thức Kubo)   |     |      Mô Phỏng Lượng Tử NEGF       |
|   - Suy giảm độ dẫn vùng (0, U_b) |     |      - Module GFET / OPEDEVS      |
|   - Bất đối xứng phân cực quang   |     |      - Tự nhất quán NEGF-Poisson  |
+-----------------------------------+     +-----------------------------------+

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

Khung phân tích của luận án kết hợp đa phương pháp trên nhiều thang đo:

  • Phương pháp liên kết chặt lượng tử hóa lần thứ hai: Sử dụng hệ cơ sở orbital $2p_z$ định xứ trên hai mạng con tam giác $A$ và $B$ với tích phân nhảy $t_{CC} \approx 2.67\text{ eV}$ và khoảng cách nguyên tử $a_{CC} \approx 1.42\text{ \AA}$. Biểu thức ma trận Hamiltonian trong không gian xung lượng: $$H(\mathbf{k}) = \begin{pmatrix} 0 & h(\mathbf{k}) \ h^*(\mathbf{k}) & 0 \end{pmatrix}, \quad h(\mathbf{k}) = -t_{CC} \sum_{j=1}^3 e^{i \mathbf{k} \cdot \mathbf{e}_j}$$ cho phép giải chính xác cấu trúc dải năng lượng mà không cần tham số hóa phức tạp như các phương pháp bán kinh nghiệm khác.
  • Hình thức luận hàm Green không cân bằng (NEGF): Phương trình hàm Green trễ của kênh dẫn dưới tác động của phân cực nguồn - máng và điện cực cổng: $$G(E) = \left[ (E + i\eta)I - H_{channel} - \Sigma_S(E) - \Sigma_D(E) - U_{scf} \right]^{-1}$$ trong đó $\Sigma_S, \Sigma_D$ là các ma trận tự năng của điện cực nguồn/máng và $U_{scf}$ là thế tĩnh điện tự nhất quán thu được từ phương trình Poisson 3D giải trên lưới vi phân hữu hạn.
  • Điều kiện biên (Boundary Conditions): Luận án khảo sát tường minh hai điều kiện biên hình học nghiêm ngặt: biên ngoằn ngoèo (zigzag) bảo toàn trạng thái biên định xứ tại mức Fermi, và biên tay vịn (armchair) dẫn đến việc mở vùng cấm năng lượng phụ thuộc vào số lượng chuỗi nguyên tử $N$.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng lượng tử (Quantum Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp mô phỏng số lượng tử chính xác (computational quantum mechanics) làm công cụ thực nghiệm số. Thiết kế nghiên cứu bao gồm ba cấp độ tích hợp:

  1. Cấp độ nguyên tử/orbital: Tính toán cấu trúc vi mô bằng phương pháp phiếm hàm mật độ từ nguyên lý đầu (Ab initio DFT via VASP code) cho các siêu ô mạng tiếp xúc kim loại (FCC: Cu, Ag, Au, Pt, Ir và HCP: Co, Cd, Ru, Ti).
  2. Cấp độ dải năng lượng và quang học: Mô hình hóa liên kết chặt kết hợp hình thức luận Kubo giải tích và số trị để tính tensor độ dẫn quang $\sigma_{\alpha\beta}(\omega)$.
  3. Cấp độ linh kiện nano: Mô phỏng tự nhất quán NEGF-Poisson giải quyết bài toán truyền dẫn đạn đạo (ballistic transport) và bán đạn đạo qua kênh dẫn GFETs.

Quy trình nghiên cứu rigorous

Quy trình tính toán và kiểm chứng được chuẩn hóa qua 4 bước:

[BƯỚC 1: Cấu Trúc Vi Mô]
Khởi tạo cấu trúc tinh thể -> DFT / VASP -> Trích xuất thông số nhảy tích phân t_pz-s, t_pz-d

[BƯỚC 2: Siêu Mạng GSLs]
Xây dựng siêu ô mạng A-GSL / Z-GSL -> Thiết lập thế tuần hoàn V(x) -> Chéo hóa ma trận Hamiltonian -> Phổ E(k)

[BƯỚC 3: Quang Học Kubo]
Tính ma trận vận tốc v_x, v_y -> Tính tích phân chuyển dời quang học -> Tensor độ dẫn quang \sigma(\omega)

[BƯỚC 4: Linh Kiện GFETs]
Module GFET / OPEDEVS -> Lặp tự nhất quán Poisson-NEGF -> Hội tụ điện tích/thế năng -> Trích xuất I_DS - V_DS, I_DS - V_GS, g_m, f_T

Độ tin cậy và tính hợp thức của nghiên cứu được kiểm soát chặt chẽ thông qua việc đối chuẩn (benchmarking) trực tiếp với các số liệu thực nghiệm:

  • Giá trị vận tốc Fermi $v_F = \frac{3}{2} \frac{t_{CC} a_{CC}}{\hbar} = 0.98 \times 10^6\text{ m/s}$ hoàn toàn khớp với thực nghiệm tán xạ quang điện tử phân giải góc (ARPES).
  • Độ dẫn cực tiểu $\sigma_{\min} = 4e^2/\pi h$ và hằng số hấp thụ quang học $\pi \alpha \approx 2.3%$ (với $\alpha = e^2/\hbar c \approx 1/137$) khớp chính xác với đo đạc quang phổ của Nair et al. (Science, 2008).

Data và phân tích

Các tham số mô phỏng và mẫu tính toán được thiết lập chi tiết:

  • Kích thước siêu mạng GSLs: Chu kỳ thế $L$ biến thiên từ $3\text{ nm}$ đến $20\text{ nm}$, độ rộng ô đơn vị $N = 2N_1 = 30$ đến $40$ chuỗi nguyên tử; biên độ thế tĩnh điện $U_b$ khảo sát từ $0\text{ eV}$ đến $6U_0$ ($U_0 \approx 0.1\text{ eV}$).
  • Thông số tiếp xúc kim loại - graphene: Cấu trúc mạng lập phương tâm mặt mặt (111) của Cu (ô 2 nguyên tử), Ag, Al, Ir, Pt, Au (ô 8 nguyên tử); mạng lục giác xếp chặt (0001) của Co (ô 2 nguyên tử), Cd, Ru, Ti (ô 8 nguyên tử).
  • Thông số linh kiện GFETs: Chiều dài kênh dẫn $L_c = 20\text{ nm}, 40\text{ nm}, 60\text{ nm}, 100\text{ nm}$; lớp điện môi dây nano GaN có hằng số điện môi $\kappa = 10$, bán kính tương đương với độ dày điện môi $t_{ins} = 16.4\text{ nm}$, điện dung cực cổng $C_{ins} = 486\text{ nF/cm}^2$, nồng độ pha tạp donor của dây GaN $n \approx 2 \times 10^{19}\text{ cm}^{-3}$.
  • Công cụ phần mềm: Tự xây dựng và phát triển các module C/C++, Fortran trong gói OPEDEVS, thuật toán đại số ma trận thưa đệ quy (Recursive Green's Function - RGF) và song song hóa MPI/OpenMP.

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

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 đột phá:

  1. Sự định xứ kỳ lạ (anomalous localization) của trạng thái điện tử trong GSLs: Phân tích mật độ xác suất $|P_n(k_y, x)|^2$ của điện tử $p_z$ trong một chu kỳ thế chứng minh rằng một số trạng thái điện tử bị dập tắt biên độ hoàn toàn ở miền rào thế hoặc giếng thế. Khác với cơ chế giam cầm cổ điển và cơ chế Anderson, hiện tượng định xứ này:

    • Xuất hiện ngay cả ở các mức năng lượng cao hơn đỉnh rào thế ($E > U_b$).
    • Không bị cô lập trong vùng cấm mà nhúng sâu (embedded) bên trong dải phổ liên tục của các trạng thái truyền qua (extended states).
    • Đóng vai trò như các trạng thái cộng hưởng Fabry-Perot lượng tử điều khiển sự chui ngầm chọn lọc qua cấu trúc tuần hoàn.

  2. Quy luật suy giảm độ dẫn quang học và tính bất đối xứng phân cực: Tính toán độ dẫn quang $\sigma(\omega)$ dựa trên hình thức luận Kubo chỉ ra rằng trong miền năng lượng photon $\hbar\omega \in (0, U_b)$, độ dẫn quang của GSLs bị suy giảm mạnh so với giá trị phổ quát $\sigma_0 = e^2/4\hbar$ của graphene tự do. Nguyên nhân là do sự biến điệu của hàm sóng làm giảm xác suất chuyển dời giữa dải hóa trị $\pi$ và dải dẫn $\pi^*$, kết hợp với hiệu ứng khóa Pauli. Đồng thời, độ dẫn quang thể hiện tính dị hướng rõ rệt đối với ánh sáng phân cực dọc theo hướng biến thiên của thế so với hướng vuông góc.

  3. Bản chất truyền dẫn và điện trở nội tại tại tiếp xúc kim loại - graphene: Khảo sát đặc tuyến $I-V$ và xác suất truyền qua $T(E, k)$ cho thấy tương tác $M-G$ phân hóa rõ thành hai nhóm:

    • Nhóm tương tác yếu (Cu, Au, Pt): Giữ nguyên hình nón Dirac nhưng làm dịch chuyển mức Fermi (tạo tiếp xúc kiểu $p-p$ hoặc $p-n$). Điện trở tiếp xúc thấp, đặc tuyến $I-V$ có tính tuyến tính cao.
    • Nhóm tương tác mạnh (Ti, Pd): Phá hủy cấu trúc nón Dirac nguyên thủy, lai hóa mạnh orbital $d-\pi$, tạo ra mật độ trạng thái phụ tại mức Fermi nhưng làm tăng điện trở tán xạ tiếp xúc.

  4. Hiệu năng siêu cao tần của GFETs cấu trúc tự sắp xếp: Mô phỏng tự nhất quán NEGF cho linh kiện GFETs với điện cực cổng dây nano GaN và cực nguồn/máng tự sắp xếp chứng minh:

    • Giảm thiểu triệt để trở kháng ký sinh lối vào nhờ liên kết tự sắp xếp hoàn hảo giữa màng Pt và dây GaN.
    • Độ hỗ dẫn đạt giá trị kỷ lục $g_m > 2000\text{ S/m}$, mật độ dòng bão hòa $I_{ON} \sim 5\text{ mA/mm}$.
    • Tần số cắt nội tại $f_T$ tăng tuyến tính khi thu nhỏ kênh dẫn dưới $40\text{ nm}$ và đạt ngưỡng $1000\text{ GHz}$ ($1\text{ THz}$) ở kênh dẫn $20\text{ nm}$.
    • Giải thích hiện tượng uốn cong của đặc tuyến $I_{DS}-V_{DS}$ tại điện áp cổng nhỏ là do sự chuyển đổi loại hạt tải từ $p$ sang $n$ và sự xuất hiện của vùng trở vi phân âm (negative differential resistance) do hiệu ứng chui ngầm Klein.
Thông số / Đặc tính MOSFET Silicon truyền thống (Công nghệ 22-32nm) GFETs kênh dẫn thông thường (Meric et al. 2008) GFETs tự sắp xếp dây GaN (Luận án Lê Hoàng Anh / Liao 2010)
Vật liệu kênh dẫn Silicon khối (3D Si) Graphene bóc tách / CVD Graphene đơn lớp tự sắp xếp
Độ linh động hạt tải ($\mu$) $\sim 500 - 1.400\text{ cm}^2/\text{V}\cdot\text{s}$ $\sim 3.000 - 10.000\text{ cm}^2/\text{V}\cdot\text{s}$ $> 20.000 - 200.000\text{ cm}^2/\text{V}\cdot\text{s}$
Độ hỗ dẫn ($g_m$) $\sim 800 - 1.200\text{ S/m}$ $\sim 150 - 500\text{ S/m}$ $> 2.000\text{ S/m}$
Tần số cắt nội tại ($f_T$) $\sim 53\text{ GHz}$ ($L_g = 550\text{ nm}$) $\sim 26 - 100\text{ GHz}$ ($L_g = 240\text{ nm}$) $\sim 300 - 1.000\text{ GHz}$ ($L_g = 20 - 100\text{ nm}$)
Cơ chế điện cực cổng Planar Gate / Omega-Gate 3D Top-gate phẳng thông thường Dây nano GaN tự sắp xếp ($\kappa \sim 10$)
Tỷ số bật/tắt ($I_{ON}/I_{OFF}$) $> 10^5 - 10^6$ (Logic số xuất sắc) $\sim 2 - 10$ (Rất thấp) $\sim 5 - 20$ (Tối ưu cho Analog/RF)

Implications đa chiều

  • Ý nghĩa học thuật và phương pháp luận: Đặt nền móng lý thuyết vững chắc cho ngành quang điện tử 2D; tích hợp thành công module GFET vào package OPEDEVS, tạo công cụ mô phỏng lượng tử mạnh mẽ cho cộng đồng vật lý tính toán tại Việt Nam và quốc tế.
  • Ý nghĩa thực tiễn công nghệ: Chỉ dẫn trực tiếp cho kỹ thuật chế tạo mạch tích hợp cao tần (RF/Microwave ICs), cảm biến quang học hồng ngoại và photodetector graphene siêu nhạy hoạt động dải rộng từ ánh sáng nhìn thấy đến terahertz.
  • Khuyến nghị chính sách công nghệ: Cần định hướng đầu tư phát triển công nghệ bán dẫn quốc gia vào phân khúc chuyên biệt: linh kiện siêu cao tần (RF) và quang điện tử quân sự/viễn thông sử dụng vật liệu 2D, thay vì cạnh tranh trực tiếp ở mảng logic số silicon thương mại quy mô lớn.

Limitations và Future Research

Luận án thừa nhận 4 giới hạn nghiên cứu xuất phát từ các xấp xỉ vật lý và điều kiện tính toán:

  1. Bỏ qua hiệu ứng tán xạ không đàn hồi và tán xạ phonon: Mô hình NEGF áp dụng trong chương 4 chủ yếu dựa trên giả định truyền dẫn đạn đạo thuần túy (ballistic transport). Trong thực tế ở nhiệt độ phòng ($T = 300\text{ K}$), tán xạ với phonon âm thanh và phonon quang học của lớp đế SiO2/GaN sẽ làm suy giảm vận tốc trôi bão hòa và độ hỗ dẫn thực tế của linh kiện.
  2. Xấp xỉ một orbital $2p_z$: Mô hình liên kết chặt chưa tính đến sự lai hóa của các orbital $\sigma$ ($sp^2$) khi màng graphene bị biến dạng uốn cong (rippling/wrinkling) hoặc chịu ứng suất cơ học mạnh.
  3. Mô hình tiếp xúc kim loại phẳng lý tưởng: Chưa tính toán tường minh ảnh hưởng của các sai hỏng cấu trúc nguyên tử (khuyết tật nút khuyết, biên hạt, tạp chất hấp phụ ngẫu nhiên) tại mặt phân giới $M-G$.
  4. Giới hạn quy mô chiều dài kênh dẫn: Chưa mở rộng mô phỏng cho các kênh dẫn $> 500\text{ nm}$ do giới hạn về chi phí tính toán của thuật toán ma trận nghịch đảo trong NEGF.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm 4 hướng cụ thể:

  • Tích hợp tán xạ electron-phonon và tán xạ tạp chất tích điện vào phương trình hàm Green thông qua các ma trận tự năng tán xạ không đàn hồi $\Sigma_{scatt}(E)$.
  • Khảo sát các cấu trúc Moiré xoay góc (Twisted Bilayer Graphene - TBG) ở góc ma thuật (magic angle $\theta \approx 1.1^\circ$) để nghiên cứu trạng thái siêu dẫn nhiệt độ cao và cách điện tương quan.
  • Phát triển mô phỏng linh kiện phototransistor kết hợp graphene với các vật liệu bán dẫn phân lớp chuyển tiếp (TMDCs như $\text{MoS}_2, \text{WSe}_2$).
  • Hoàn thiện gói phần mềm OPEDEVS thành nền tảng mã nguồn mở phục vụ đào tạo và nghiên cứu bán dẫn.

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

Các công trình công bố từ luận án trên các tạp chí quốc tế chuyên ngành thuộc Viện Vật lý Hoa Kỳ (AIP) và Viện Vật lý Vương quốc Anh (IOP) như Applied Physics Letters [Appl. Phys. Lett. 101, 161605 (2012)] và Journal of Physics: Condensed Matter đã khẳng định chất lượng học thuật đạt chuẩn quốc tế:

  • Tác động học thuật: Các bài báo đã được trích dẫn và sử dụng bởi nhiều nhóm nghiên cứu hàng đầu thế giới về vật lý chất rắn và linh kiện graphene, khẳng định vị thế của nghiên cứu vật lý bán dẫn Việt Nam trên bản đồ khoa học thế giới.
  • Chuyển đổi công nghiệp bán dẫn: Cung cấp cơ sở dữ liệu tính toán vật liệu cho các tập đoàn vi điện tử trong việc lựa chọn kim loại tiếp xúc tối ưu (Ti, Pt, Pd, Au) và thiết kế cấu trúc transistor RF hoạt động ở băng tần sóng milimet và terahertz.
  • Lợi ích xã hội và an ninh quốc phòng: Công nghệ GFETs và cảm biến quang điện tử graphene mở đường cho việc chế tạo các hệ thống radar mảng pha quét điện tử tốc độ cao, hệ thống viễn thông 5G/6G, và thiết bị nhìn đêm thế hệ mới tiêu thụ năng lượng cực thấp.

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Vật lý/Vật liệu: Tiếp cận phương pháp luận chuẩn mực về mô hình hóa liên kết chặt, hình thức luận Kubo và công cụ NEGF lượng tử.
  • Các nhà khoa học và chuyên gia vật lý tính toán: Thừa hưởng các module thuật toán mô phỏng trong package OPEDEVS để mở rộng cho các hệ vật liệu hai chiều mới (Silicene, Phosphorene, TMDCs).
  • Kỹ sư R&D công nghệ bán dẫn: Nhận được chỉ dẫn định lượng về thông số điện trở tiếp xúc, quy trình chế tạo tự sắp xếp và giới hạn tần số hoạt động của linh kiện transistor nano.
  • Nhà hoạch định chính sách phát triển công nghệ cao: Có cơ sở khoa học tin cậy để định hình chiến lược nghiên cứu vi mạch quốc gia, ưu tiên nguồn lực cho các công nghệ bán dẫn chuyên dụng đón đầu xu thế hậu silicon.

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 lý thuyết độc đáo nhất là việc phát hiện và chứng minh bản chất của sự định xứ dị thường của các trạng thái điện tử $p_z$ trong siêu mạng graphene (GSLs). Luận án đã mở rộng lý thuyết Dirac-Weyl hai chiều và lý thuyết giam cầm lượng tử cổ điển, chứng minh rằng các trạng thái định xứ có thể tồn tại ở mức năng lượng cao hơn bờ rào thế ($E > U_b$) và nhúng trực tiếp trong dải phổ liên tục của các trạng thái mở rộng mà không cần sự tồn tại của một vùng cấm năng lượng hữu hạn.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?

So với các mô hình bán cổ điển của Meric et al. (2008) hay mô hình tiếp xúc Schottky thuần túy của Thiele et al. (2010), luận án đã:

  • Tích hợp thành công mô hình tự năng tương tác vi mô $\Sigma_{G-M}$ rút ra từ tính toán nguyên lý đầu DFT vào hệ phương trình truyền dẫn lượng tử NEGF.
  • Xây dựng thuật toán giải tự nhất quán 3D giữa phương trình Poisson và phương trình hàm Green không cân bằng trong gói phần mềm OPEDEVS, cho phép mô phỏng chính xác cấu trúc thực nghiệm GFETs dây nano GaN tự sắp xếp với chi phí tính toán tối ưu hơn nhiều so với các phần mềm thương mại.

3. Phát hiện thực nghiệm/mô phỏng nào bất ngờ nhất và bằng chứng dữ liệu đi kèm?

Phát hiện bất ngờ nhất là sự suy giảm mạnh của độ dẫn quang học $\text{Re}[\sigma(\omega)]$ trong dải photon $(0, U_b)$. Về mặt trực giác cổ điển, sự có mặt của các rào thế tuần hoàn sẽ làm tăng mật độ trạng thái chuyển dời. Tuy nhiên, dữ liệu tính toán giải tích từ hình thức luận Kubo chứng minh rằng sự biến điệu không gian của hàm sóng điện tử làm triệt tiêu các yếu tố ma trận chuyển dời lưỡng cực quang học giữa dải $\pi$ và $\pi^*$, kết hợp với cơ chế khóa Pauli, dẫn đến độ dẫn quang bị dập tắt gần như hoàn toàn ở năng lượng thấp trước khi phục hồi về giá trị $\sigma_0 = e^2/4\hbar$ khi $\hbar\omega > U_b$.

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

Có. Luận án cung cấp chi tiết toàn bộ hệ cơ sở toán học tại phần Phụ lục:

  • Bảng ma trận Hamiltonian và ma trận vận tốc cho siêu mạng A-GSLs và Z-GSLs.
  • Bảng thông số năng lượng nhảy $t_{CC} = 2.67\text{ eV}$, $a_{CC} = 1.42\text{ \AA}$, các tích phân nhảy $t_{pz-s}$ và $t_{pz-d}$ cho 10 cấu hình tiếp xúc kim loại - graphene.
  • Quy trình thuật toán lặp tự nhất quán Poisson-NEGF và hướng dẫn sử dụng chi tiết module GFET trong package OPEDEVS.

5. Kế hoạch nghiên cứu 10 năm được phác thảo như thế nào?

Kế hoạch nghiên cứu 10 năm tập trung vào 3 trục chính:

  1. Nâng cấp bộ công cụ OPEDEVS hỗ trợ tương tác spin-quỹ đạo (Spin-Orbit Coupling) và các hiệu ứng nhiệt - điện tử (thermoelectric transport).
  2. Thiết kế linh kiện logic tiết kiệm năng lượng dựa trên cấu trúc siêu mạng Moiré và dị thể Van der Waals (Graphene/h-BN/TMDC).
  3. Hợp tác với các phòng thí nghiệm chế tạo bán dẫn quốc tế để hiện thực hóa các mạch tích hợp RF Terahertz và màng photodetector phân cực siêu nhạy ứng dụng trong viễn thông 6G.

Kết luận

  1. Thiết lập hoàn chỉnh bức tranh cấu trúc vùng năng lượng của siêu mạng graphene (A-GSLs và Z-GSLs) chịu thế tuần hoàn 1D bằng phương pháp gần đúng liên kết chặt, làm sáng tỏ sự hình thành các nón Dirac thế hệ mới và cơ chế tán sắc dị hướng.
  2. Phát hiện và giải thích bản chất vật lý của hiện tượng định xứ dị thường của các trạng thái điện tử $2p_z$ nhúng trong dải phổ liên tục, đóng vai trò trạng thái cộng hưởng Fabry-Perot lượng tử.
  3. Giải thích tường minh cơ chế suy giảm độ dẫn quang học trong dải photon $(0, U_b)$ và sự phụ thuộc trạng thái phân cực ánh sáng thông qua hình thức luận phản ứng tuyến tính Kubo.
  4. Đề xuất mô hình hiệu dụng mô tả liên kết điện tử tại mặt tiếp xúc giữa graphene và 10 cấu hình kim loại điển hình (Cu, Ag, Au, Pt, Ir, Co, Cd, Ru, Ti), cho phép xác định chính xác điện trở tiếp xúc nội tại.
  5. Phát triển thành công module GFETs trong gói phần mềm mô phỏng lượng tử OPEDEVS, thực hiện mô phỏng tự nhất quán NEGF-Poisson cho transistor kênh dẫn nano tự sắp xếp dùng dây GaN, chứng minh khả năng hoạt động ở tần số terahertz ($f_T \sim 1000\text{ GHz}$).
  6. Công trình tạo tiền đề vững chắc cho việc mở ra 3 nhánh nghiên cứu mũi nhọn: quang điện tử 2D phân cực, mạch tích hợp nano siêu cao tần (RF Terahertz), và phần mềm mô phỏng lượng tử chuyên dụng phục vụ chiến lược khoa học công nghệ quốc gia.