Tổng quan về luận án

Sự phát triển vượt bậc của công nghệ vi điện tử bán dẫn trong hơn nửa thế kỷ qua dựa trên quy luật thu nhỏ kích thước (scaling) các transistor hiệu ứng trường kim loại - oxit - bán dẫn (MOSFET) theo định luật Moore (1965) và nguyên lý Dennard (1974). Tuy nhiên, khi tiến vào giới hạn dưới 10 nm, các linh kiện silicon truyền thống đối mặt với hàng loạt rào cản vật lý nghiêm trọng: hiệu ứng kênh ngắn (short-channel effects), dòng rò lượng tử qua lớp cách điện cổng, sự suy giảm độ linh động hạt tải do tán xạ pha tạp và mật độ tỏa nhiệt Joule vượt ngưỡng kiểm soát. Trong bối cảnh đó, graphene – mạng tinh thể carbon đơn lớp hai chiều (2D) dạng tổ ong với liên kết lai hóa $sp^2$ – nổi lên như một ứng viên đột phá nhờ độ linh động điện tử nội tại cực cao ($\mu \sim 2 \times 10^5\text{ cm}^2/\text{V}\cdot\text{s}$ ở nhiệt độ phòng) và vận tốc Fermi $v_F \approx 10^6\text{ m/s}$, mở ra triển vọng hiện thực hóa các mạch tích hợp hoạt động ở dải tần số Terahertz (THz).

Mặc dù sở hữu các đặc tính vật lý vượt trội, việc ứng dụng graphene vào các linh kiện điện tử và quang điện tử nano gặp phải ba khoảng trống nghiên cứu (research gaps) mang tính cốt lõi:

  1. Thiếu hụt mô hình lý thuyết lượng tử toàn diện giải thích vi mô sự điều biến cấu trúc vùng năng lượng và động học tán sắc của các giả hạt Dirac khi chịu tác động bởi trường thế tĩnh điện biến thiên tuần hoàn một chiều (Graphene Superlattices - GSLs).
  2. Bản chất vi mô của sự truyền dẫn điện tử và điện trở tiếp xúc tại bề mặt phân giới kim loại - graphene (M-G interface) trong cấu trúc linh kiện thực tế chưa được mô hình hóa chính xác qua hình thức luận tự năng lượng lượng tử.
  3. Thiếu công cụ mô phỏng lượng tử tự nhất quán (self-consistent quantum transport simulation) kết hợp phương pháp hàm Green không cân bằng (NEGF) để đánh giá linh kiện transistor hiệu ứng trường graphene (GFETs) có cấu trúc cổng tự sắp xếp (self-aligned) kích thước sub-100 nm.

Để giải quyết các rào cản trên, luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Trường thế tĩnh điện tuần hoàn 1D làm biến đổi phổ tán sắc Dirac, tính đối xứng tinh thể và các trạng thái định xứ điện tử trong GSLs như thế nào?
    • H1: Trường thế tuần hoàn phá vỡ tính đẳng hướng của vận tốc Fermi, dẫn đến sự định xứ dị thường của các điện tử $2p_z$ ngay cả ở các mức năng lượng cao hơn bờ rào thế tĩnh điện.
  • RQ2: Cơ chế kích thích quang học và độ dẫn quang động của GSLs thay đổi ra sao dưới bức xạ điện từ phân cực trong dải năng lượng photon thấp?
    • H2: Sự suy giảm độ dẫn quang xuất hiện trong miền năng lượng photon $(0, U_b)$ do hiệu ứng dập tắt ma trận chuyển mức quang học và quy tắc chọn lọc phân cực.
  • RQ3: Tương tác liên kết điện tử và tự năng tiếp xúc $\Sigma_{G-M}$ tại các điện cực kim loại ảnh hưởng thế nào đến đặc trưng $I_{DS}-V_{GS}$ và $I_{DS}-V_{DS}$ của linh kiện GFETs kênh ngắn?
    • H3: Sự kết hợp giữa thế tự nhất quán Poisson-Schrödinger và tự năng tiếp xúc kim loại phản ánh trung thực điểm trung hòa điện tích (Dirac point shift) và độ hỗ dẫn $g_m$ của linh kiện GFETs.

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), phương trình sóng Dirac hai chiều, hình thức luận đáp ứng tuyến tính Kubo (Kubo Formalism) và phương pháp hàm Green không cân bằng (Non-Equilibrium Green's Functions - NEGF). Phạm vi nghiên cứu bao quát các cấu hình dải nano graphene biên armchair (A-GSLs) và biên zigzag (Z-GSLs), các tiếp xúc tiếp mặt với kim loại họ FCC (Cu, Au, Pt, Pd) và HPC (Ti), cùng mô hình linh kiện GFETs kênh dẫn chiều dài từ $40\text{ nm}$ đến $100\text{ nm}$ sử dụng điện cực cổng dây nano GaN có hằng số điện môi cao $\kappa \approx 10$.

Literature Review và Positioning

Cấu trúc vùng năng lượng của điện tử trong mạng tinh thể graphene lần đầu tiên được Wallace khảo cứu vào năm 1946 bằng phương pháp liên kết chặt như một mô hình toán học nền tảng cho than chì graphite. Tuy nhiên, suốt hơn nửa thế kỷ, quan điểm nhiệt động học của Landau và Peierls cho rằng màng tinh thể hai chiều không thể tồn tại độc lập ở điều kiện thường do thăng giáng nhiệt phá hủy trật tự tinh thể tầm xa. Bước ngoặt lịch sử diễn ra vào năm 2004 khi Andre Geim và Kostya Novoselov (Đại học Manchester) cô lập thành công đơn lớp graphene bằng phương pháp bóc tách cơ học, quan sát thấy hiệu ứng trường điện từ và hiệu ứng Hall lượng tử dị thường với bậc nhảy $\Delta\sigma_{xy} = 4e^2/h$ ở nhiệt độ phòng.

Trong dòng chảy lý thuyết truyền dẫn, một tranh luận học thuật lớn xoay quanh hiện tượng chui ngầm Klein (Klein Tunneling) được Katsnelson, Geim và Novoselov (2006) đặt ra: các fermion Dirac không khối lượng có khả năng xuyên qua các rào thế tĩnh điện cao vô hạn với xác suất truyền qua bằng 1 ở góc tới vuông góc, khiến việc giam cầm điện tử bằng cổng tĩnh điện thông thường trở nên bất khả thi. Nhóm Tworzydło et al. (2006) đã chứng minh độ dẫn cực tiểu phổ quát đạt giá trị $\sigma_{min} = \frac{4e^2}{\pi h}$ thông qua các trạng thái mờ (evanescent modes). Đối với siêu mạng graphene, Park et al. (2008, 2009) và Brey & Fertig (2009) phát hiện sự tái cấu trúc nón Dirac và xuất hiện các điểm Dirac phụ (satellite Dirac points) dưới thế tuần hoàn 1D. Tuy nhiên, các nghiên cứu này chủ yếu dùng phép gần đúng Dirac liên tục, bỏ qua các hiệu ứng gián đoạn mạng nguyên tử và không làm rõ được bức tranh chuyển mức quang học.

So với các nghiên cứu thực nghiệm quốc tế tiên phong như công trình của Lei Liao et al. (Nature, 2010) về GFETs tự sắp xếp dùng dây nano GaN đạt độ hỗ dẫn $g_m > 2000\text{ S/m}$ và mô phỏng đạn đạo của Jyotsna Chauhan et al. (IEEE TED, 2011), luận án định vị sự đóng góp độc đáo ở ba điểm:

  1. Xây dựng ma trận Hamiltonian kích thước lớn trong không gian vector sóng cho siêu mạng A-GSL và Z-GSL, phát hiện bản chất định xứ dị thường của điện tử $p_z$.
  2. Ứng dụng hình thức luận Kubo giải thích sự triệt tiêu độ dẫn quang trong dải $(0, U_b)$, làm sáng tỏ thực nghiệm của Li et al. (2008) và Phaedon Avouris et al. (IBM, 2012).
  3. Đề xuất mô hình tự năng hiệu dụng $\Sigma_{G-M}$ gắn kết tính chất vi mô bề mặt kim loại vào hệ phương trình vận chuyển NEGF trong gói phần mềm mô phỏng thiết bị quang - điện tử OPEDEVS.

Đó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 thuyết liên kết chặt của Wallace và hình thức luận trường Dirac lượng tử hóa lần hai cho hệ thấp chiều chịu thế tuần hoàn, mang lại các bước tiến lý thuyết cụ thể:

  • Mở rộng thuyết tán sắc Dirac trong thế siêu mạng: Thiết lập hệ ma trận Hamiltonian tuần hoàn $H_{total} = H + V$, chứng minh rằng thế vô hướng một chiều $U(x)$ điều biến mạnh mẽ hàm sóng điện tử, bẻ cong các mặt đẳng năng lượng và ghim (pinning) các trạng thái năng lượng tại biên vùng Brillouin mới.
  • Phát hiện hiện tượng định xứ dị thường (Anomalous Electronic Localization): Khác với cơ chế giam cầm lượng tử cổ điển (chỉ xảy ra trong giếng thế), luận án chứng minh sự tồn tại của các trạng thái định xứ của điện tử $2p_z$ nằm ở mức năng lượng bên ngoài khoảng biến thiên của hàm thế ($E > U_b$). Các trạng thái này nhúng trực tiếp vào phổ liên tục của các trạng thái mở rộng mà không làm mất tính kết hợp pha, đóng vai trò như các mode cộng hưởng Fabry-Pérot nội tại.
  • Lý thuyết chuyển dời quang học trong GSLs: Mở rộng công thức độ dẫn quang động học Kubo, chỉ ra rằng việc phá vỡ đối xứng tinh thể $D_{6h}$ sang đối xứng định hướng $1\text{D}$ làm suy giảm xen phủ hàm sóng giữa các dải $\pi$ và $\pi^*$. Điều này dẫn tới hiện tượng "khóa chuyển mức" (optical transition blocking) tương tự khóa Pauli trong miền photon $\hbar\omega \in (0, U_b)$.

$$\sigma_{\alpha\beta}(\omega) = \frac{ie^2\hbar}{\Omega} \sum_{n \neq m} \frac{f(E_n) - f(E_m)}{E_n - E_m} \frac{\langle n | v_\alpha | m \rangle \langle m | v_\beta | n \rangle}{E_n - E_m + \hbar\omega + i\eta}$$

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết: Mô hình Tight-Binding orbital $\pi$ kết cặp pha Peierls, Hình thức luận động học Kubo và Thuyết hàm Green không cân bằng (NEGF) Keldysh.

Điều kiện biên và phạm vi áp dụng của khung phân tích được xác định rõ:

  1. Giới hạn nhiệt độ thấp đến nhiệt độ phòng ($T = 0\text{ K} - 300\text{ K}$) trong chế độ truyền dẫn đạn đạo (ballistic transport).
  2. Biên thế tĩnh điện mượt được mô tả qua hàm phân bố Fermi-Dirac $f(x) = [1 + \exp(\alpha x)]^{-1}$.
  3. Bề mặt tiếp xúc kim loại - graphene được giải ma trận tự năng lượng $\Sigma_{G-M}$ riêng phần dựa trên dữ liệu cấu trúc vùng nguyên lý đầu (DFT-VASP).

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

Thiết kế nghiên cứu

Luận án tuân thủ lập trường nhận thức luận thực chứng (Positivism) kết hợp phương pháp mô phỏng vật lý tính toán lượng tử chính xác cao (Computational Quantum Physics). Thiết kế nghiên cứu đa tầng (multi-scale design) đi từ cấp độ orbital nguyên tử ($10^{-10}\text{ m}$) đến cấu trúc siêu mạng ($10^{-8}\text{ m}$) và cấu trúc linh kiện hoàn chỉnh ($10^{-7}\text{ m}$).

Cấp độ phân tích Đối tượng vật lý Công cụ / Phương pháp Thông số đầu vào then chốt
Nguyên tử / Dải nano Đơn lớp Graphene, GNRs (Armchair, Zigzag) Tight-Binding, Lượng tử hóa lần 2 $a_{CC} = 1.42\text{ \AA}$, $t_{CC} = 2.67\text{ eV}$
Siêu mạng (GSLs) A-GSLs, Z-GSLs chịu thế $U(x)$ Hamiltonian $N \times N$, Hình thức luận Kubo Chu kỳ $T = 3N a_{CC}$, Biên độ $U_b = 0 - 6U_0$
Tiếp xúc phân giới Tiếp xúc M-G (Cu, Au, Pt, Pd, Ti) DFT-VASP 4.6 kết hợp mô hình tiếp xúc Ô mạng FCC (111), HPC (0001), Tự năng $\Sigma_{G-M}$
Linh kiện nano GFETs cổng dây nano GaN NEGF tự nhất quán (Gói OPEDEVS) $L_c = 40-100\text{ nm}$, $\kappa = 10$, $V_{DS}, V_{GS}$

Quy trình nghiên cứu rigorous

Quy trình tính toán tuân thủ quy chuẩn toán lý chặt chẽ:

  1. Rời rạc hóa không gian thực và biến đổi Fourier: Chuyển đổi toán tử sinh hủy $a_i, b_j$ sang không gian sóng đảo $k$, thiết lập ma trận Hamiltonian dải năng lượng: $$H(k) = \begin{pmatrix} 0 & h(k) \ h^*(k) & 0 \end{pmatrix}$$ với $h(k) = -t_{CC}\sum_{j=1}^3 e^{i \mathbf{k}\cdot\mathbf{e}_j}$, giải tích trị riêng $E(k) = \pm |h(k)|$.
  2. Thuật toán ma trận thưa và đệ quy: Trong cấu hình GSLs, ô cơ sở mở rộng chứa $N$ nguyên tử dẫn tới ma trận kích thước lớn. Luận án triển khai thuật toán đường chéo hóa ma trận khối thưa (sparse block-tridiagonal matrix solver) để tối ưu thời gian tính toán.
  3. Chu trình tự nhất quán NEGF-Poisson: Hàm Green chậm $G^R(E)$ và hàm mật độ điện tích $\rho(\mathbf{r})$ được lặp tự nhất quán qua phương trình Poisson tĩnh điện: $$G^R(E) = \left[ (E + i\eta)I - H_{TB} - \Sigma_S(E) - \Sigma_D(E) - U_{sc}(\mathbf{r}) \right]^{-1}$$ $$\nabla \cdot [\varepsilon(\mathbf{r}) \nabla U_{sc}(\mathbf{r})] = -\rho(\mathbf{r}) = -e \int \frac{dE}{2\pi} G^<(E; \mathbf{r}, \mathbf{r})$$
                                < |U_new - U_old| < ε? >

Data và phân tích

Dữ liệu mô phỏng được kiểm chứng chéo (triangulation) giữa ba phương pháp: giải tích giải phương trình Dirac, tính toán số trị Tight-Binding trên gói phần mềm tự phát triển, và so chuẩn với kết quả nguyên lý đầu từ VASP 4.6. Độ hội tụ tự nhất quán của thế tĩnh điện đạt chuẩn sai số $\Delta U < 10^{-5}\text{ eV}$.

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

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

  1. Sự hình thành và dị hướng nón Dirac phụ trong GSLs: Khi tăng biên độ thế $U_b$ từ $0\text{ eV}$ lên $6U_0$, các nón Dirac mới xuất hiện tại các điểm đối xứng của vùng Brillouin siêu mạng. Vận tốc Fermi theo phương song song với rào thế $v_y$ bị suy giảm mạnh về gần 0, trong khi vận tốc trực giao $v_x$ duy trì giá trị không đổi $v_F \approx 10^6\text{ m/s}$, tạo ra tính dẫn hướng dòng điện tử dị hướng cực hạn.
  2. Bằng chứng định lượng về trạng thái định xứ dị thường: Phân tích mật độ xác suất $|\psi_n(k_y, x)|^2$ của điện tử $p_z$ chứng minh sự cô lập mật độ điện tích tại các miền cục bộ dọc theo chu kỳ thế. Đặc biệt, các mức năng lượng $E > U_b$ vẫn thể hiện hàm sóng bị thắt nút hoàn toàn, chứng minh sự tồn tại của các trạng thái giả liên kết (quasi-bound states) trong continuum.
  3. Hiện tượng triệt tiêu độ dẫn quang trong miền $(0, U_b)$: Độ dẫn quang động học $\sigma(\omega)$ của graphene đồng nhất giữ giá trị hằng số phổ quát $\sigma_0 = \frac{e^2}{4\hbar}$. Tuy nhiên, trong GSLs, $\sigma(\omega)$ giảm sâu tiệm cận 0 trong khoảng photon năng lượng nhỏ $\hbar\omega < U_b$ do triệt tiêu yếu tố ma trận vận tốc $\langle n | v | m \rangle$, kèm theo tính phụ thuộc nhạy bén vào góc phân cực ánh sáng tới.
  4. Bản chất tiếp xúc kim loại - graphene: Khảo sát 5 tổ hợp tiếp xúc kim loại (Cu, Au, Pt, Pd, Ti) chỉ ra rằng tương tác $d-\pi$ làm dịch chuyển mức Fermi và mở rộng phần ảo của tự năng $\Gamma = \text{Im}(\Sigma_{G-M})$. Kim loại Pd và Pt thể hiện khả năng tiêm hạt tải tối ưu với điện trở tiếp xúc nội tại nhỏ nhất, duy trì cấu trúc nón Dirac không bị biến dạng nghiêm trọng.
  5. Mô phỏng GFETs kênh ngắn sub-100 nm: Mô phỏng thành công linh kiện GFETs kênh $L_c = 40-60\text{ nm}$ sử dụng điện cực cổng GaN ($\kappa \approx 10$). Kết quả tính toán giải thích chính xác đường đặc trưng $I_{DS}-V_{GS}$ dạng chữ V bất đối xứng và hiện tượng bão hòa dòng yếu do hiệu ứng chui ngầm Klein ở trạng thái khóa kênh.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp khung động học lượng tử chuẩn xác cho các hệ graphene dị thể và siêu mạng Moiré, đặt tiền đề giải thích các pha tương quan điện tử mạnh trong cấu trúc góc xoay (twisted bilayer graphene).
  • Về mặt phương pháp luận: Khẳng định tính ưu việt của việc đóng gói module GFETs trong nền tảng OPEDEVS, cho phép mô phỏng linh kiện 2D phức tạp với chi phí tính toán thấp hơn hàng chục lần so với DFT nguyên lý đầu thuần túy.
  • Về mặt ứng dụng kỹ thuật: Cung cấp quy tắc thiết kế (design guidelines) then chốt cho các kỹ sư bán dẫn: lựa chọn kim loại tiếp xúc tiếp xúc tiếp xúc có công thoát phù hợp (Pd, Pt) và cấu hình cổng điện môi cao để tối đa hóa độ hỗ dẫn $g_m > 2000\text{ S/m}$ trong các bộ tiền khuếch đại RF siêu cao tần.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu có tính biên:

  1. Mô hình truyền dẫn lượng tử NEGF chủ yếu xét trong chế độ đạn đạo (ballistic), chưa tích hợp đầy đủ các kênh tán xạ phonon nội tại (acoustic & optical phonons) và tán xạ tạp chất tích điện ở nhiệt độ cao ($T > 300\text{ K}$).
  2. Giản đồ thế tĩnh điện siêu mạng $U(x)$ được giả định là biến thiên 1D lý tưởng, chưa bao quát các thăng giáng thế ngẫu nhiên 2D (charge puddles) sinh ra từ lớp đế điện môi $\text{SiO}_2$.
  3. Tiếp xúc kim loại - graphene được mô hình hóa qua ma trận tự năng lượng hiệu dụng một mức, chưa xét đến biến dạng cơ học cục bộ (strain) tại rìa tiếp xúc kim loại.

Định hướng nghiên cứu tương lai:

  • Tích hợp cơ chế tán xạ không đàn hồi electron-phonon vào thuật toán tự nhất quán NEGF của gói OPEDEVS.
  • Khảo sát các cấu trúc siêu mạng Moiré hai chiều (2D GSLs) trên nền vật liệu cách điện tinh thể boron nitride lục giác (h-BN).
  • Thiết kế và mô phỏng các bộ điều biến quang tử Terahertz tích hợp dựa trên hiệu ứng triệt tiêu độ dẫn quang trong GSLs.

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

Các kết quả nghiên cứu trong luận án đã được công bố trên các tạp chí vật lý quốc tế danh tiếng thuộc hệ thống ISI:

  • Bài báo về tính chất điện tử và sự định xứ dị thường trong GSLs công bố trên Applied Physics Letters (AIP).
  • Bài báo về đặc trưng hấp thụ và độ dẫn quang học của GSLs công bố trên Journal of Physics: Condensed Matter (IOP).
  • Công trình về mô hình liên kết điện tử và điện trở tiếp xúc kim loại - graphene đăng trên Applied Physics Letters.

Về mặt công nghệ, nghiên cứu đóng góp nền tảng tính toán hỗ trợ các viện nghiên cứu và doanh nghiệp bán dẫn trong việc định hình kiến trúc bóng bán dẫn hậu silicon (Beyond-Silicon Roadmap), thúc đẩy chế tạo các cảm biến quang học nano và mạch tích hợp tần số siêu cao phục vụ viễn thông thế hệ mới (5G/6G).

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Nắm bắt phương pháp luận toán lý chặt chẽ từ giải tích Dirac đến kỹ thuật tính ma trận thưa trong mô hình liên kết chặt và NEGF.
  • Chuyên gia Vật lý lý thuyết & Vật lý chất rắn: Khai thác các phát hiện về trạng thái định xứ dị thường và cơ chế điều biến quang trong siêu mạng 2D để mở rộng cho các vật liệu đơn lớp mới như TMDcs (MoS2, WS2) hay Phosphorene.
  • Kỹ sư R&D Bán dẫn & Thiết bị vi điện tử: Sử dụng trực tiếp module GFETs trong gói OPEDEVS để thiết kế tối ưu hóa độ dày lớp cách điện cổng, lựa chọn vật liệu kim loại nguồn/máng và kiểm soát dòng rò linh kiện.

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à phát hiện hiện tượng định xứ dị thường của các điện tử $2p_z$ trong siêu mạng graphene. Luận án đã mở rộng thuyết liên kết chặt của Wallace (1946) và phương trình Dirac cho vật liệu 2D, chứng minh rằng sự điều biến của trường thế tuần hoàn 1D có thể tạo ra các trạng thái giả liên kết (quasi-bound states) có hàm sóng bị dập tắt biên độ hoàn toàn tại một vùng không gian, tồn tại ngay ở mức năng lượng cao hơn bờ rào thế ($E > U_b$) và nhúng trong phổ liên tục của các trạng thái mở rộng.

2. Phương pháp nghiên cứu có điểm gì cải tiến so với các công trình quốc tế tiền nhiệm?

So với các nghiên cứu của Park et al. (2008) chỉ dùng gần đúng Dirac dải liên tục hoặc Huard et al. (2008) tiếp cận bán cổ điển, luận án đã:

  1. Xây dựng ma trận Hamiltonian nguyên tử đầy đủ có xét đến pha Peierls cho ô đơn vị siêu mạng kích thước lớn ($N$ nguyên tử).
  2. Tích hợp ma trận tự năng lượng tiếp xúc $\Sigma_{G-M}$ rút ra từ cấu trúc vi mô DFT vào hệ phương trình NEGF tự nhất quán Poisson-Schrödinger trong gói OPEDEVS, mô tả chính xác cả hiệu ứng lượng tử đạn đạo lẫn ảnh hưởng của điện cực kim loại thực tế.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng số liệu minh chứng?

Phát hiện bất ngờ nhất là sự suy giảm mạnh độ dẫn quang học $\sigma(\omega)$ về gần 0 trong dải photon $(0, U_b)$. Trong khi graphene tự do luôn có độ dẫn quang hằng số $\sigma_0 = \frac{e^2}{4\hbar} \approx 6.08 \times 10^{-5}\text{ \Omega}^{-1}$ bất kể tần số, thì trong GSLs, sự biến điệu của hàm sóng làm triệt tiêu các tích phân xen phủ ma trận chuyển mức dipole giữa dải hóa trị $\pi$ và dải dẫn $\pi^*$, tạo nên một "cửa sổ quang học trong suốt nhân tạo" phụ thuộc vào biên độ thế $U_b$.

4. Giao thức mô phỏng có khả năng tái lặp (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ biểu thức giải tích ma trận Hamiltonian cho A-GSLs và Z-GSLs trong phụ lục, thông số liên kết chuẩn hóa ($t_{CC} = 2.67\text{ eV}$, $a_{CC} = 1.42\text{ \AA}$), thuật toán lặp tự nhất quán Poisson-NEGF và hướng dẫn vận hành Module GFETs trong packages OPEDEVS. Mọi kết quả phổ năng lượng và đặc tuyến truyền dẫn đều có thể tái lập chính xác.

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

Lộ trình 10 năm tập trung vào:

  • Mở rộng mô hình sang các siêu mạng Moiré 2D tạo bởi góc xoay ma thuật (Magic-Angle Twisted Bilayer Graphene).
  • Tích hợp tán xạ đa hạt (electron-electron, electron-phonon) trong NEGF.
  • Thiết kế mạch tích hợp lai GFET-CMOS và bộ nhớ lượng tử dựa trên trạng thái định xứ topo.

Kết luận

  1. Hệ thống hóa toàn diện cơ sở vật lý lượng tử: Xác lập bức tranh giải tích và số trị hoàn chỉnh về cấu trúc vùng năng lượng dải $2p_z$ của graphene đơn lớp, dải nano biên zigzag/armchair và các cấu trúc siêu mạng 1D.
  2. Khám phá cơ chế định xứ dị thường: Chứng minh sự tồn tại của các trạng thái điện tử định xứ trên bờ rào thế tĩnh điện, mở rộng hiểu biết kinh điển về nghịch lý chui ngầm Klein trong các vật liệu Dirac hai chiều.
  3. Làm sáng tỏ động học quang học GSLs: Giải thích hiện tượng triệt tiêu độ dẫn quang trong miền photon $(0, U_b)$ và tính dị hướng phân cực quang học bằng hình thức luận Kubo.
  4. Mô hình hóa thành công tiếp xúc kim loại - graphene: Đề xuất phương pháp tự năng $\Sigma_{G-M}$ lượng hóa điện trở tiếp xúc và dịch chuyển mức Fermi cho 5 họ kim loại chuyển tiếp chủ đạo.
  5. Phát triển nền tảng mô phỏng linh kiện OPEDEVS: Đóng gói module GFETs lượng tử tự nhất quán, mô phỏng thành công đặc tính tĩnh và động của transistor graphene kênh ngắn sub-100 nm có cổng tự sắp xếp dây nano GaN.
  6. Mở ra hướng phát triển công nghệ linh kiện Terahertz: Đặt nền móng lý thuyết và công cụ thiết kế tin cậy cho thế hệ linh kiện điện tử và quang điện tử nano tốc độ siêu cao hậu silicon.