Tổng quan về luận án

Sự phát triển của công nghệ vi điện tử bán dẫn trong hơn nửa thế kỷ qua vận hành chặt chẽ theo định luật thực nghiệm Moore (1965), với động lực cốt lõi là quá trình thu nhỏ kích thước (scaling) của các transistor hiệu ứng trường kim loại - oxit - bán dẫn (MOSFET) nhằm gia tăng mật độ tích hợp, nâng cao tốc độ xử lý và giảm thiểu công suất tiêu thụ. Tuy nhiên, khi kích thước kênh dẫn của các linh kiện Si-MOSFET tiến dần về giới hạn dưới $10\text{ nm}$ – tương đương với công nghệ $22\text{ nm}$ thương mại sử dụng cấu trúc cổng bao quanh đa chiều (3D $\Omega$-gate) [49, 50] – các rào cản vật lý lượng tử cơ bản đã xuất hiện gay gắt. Các hiệu ứng kênh ngắn (short-channel effects), hiện tượng chui ngầm lượng tử qua rào thế, sự suy giảm khả năng kiểm soát tĩnh điện của điện cực cổng phẳng và sự phát tán nhiệt lượng khổng lồ đã đe dọa trực tiếp đến tính bền vững của lý thuyết co giãn Dennard (1974).

Trong bối cảnh đó, sự khám phá ra graphene vào năm 2004 bởi Andre Geim và Konstantin Novoselov [69] – vật liệu carbon hai chiều phẳng với độ dày đơn nguyên tử – đã tạo nên một bước ngoặt mang tính cách mạng cho vật lý chất rắn và công nghệ nano điện tử. Graphene sở hữu các đặc tính vật lý vượt bậc: độ linh động điện tử nội tại cực cao ($\mu_G \sim 2 \cdot 10^5\text{ cm}^2/\text{V}\cdot\text{s}$), độ dẫn nhiệt và độ bền cơ học vượt trội, cùng hệ thức tán sắc tuyến tính dạng nón Dirac tại các điểm đối xứng cao $K$ và $K'$ trong vùng Brillouin thứ nhất. Các giả hạt trong graphene hành xử như các fermion Dirac không khối lượng 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) [98]. Điều này mang lại tiềm năng hiện thực hóa các linh kiện transistor tần số siêu cao trong dải terahertz ($>100\text{ GHz}$), vượt xa giới hạn $10\text{ GHz}$ của vật liệu silicon truyền thống [147, 149].

Tuy nhiên, bản chất bán kim loại không có vùng cấm (zero-bandgap) của graphene nguyên bản dẫn đến việc các linh kiện transistor hiệu ứng trường graphene (GFETs) không thể đóng ngắt hoàn toàn (tỷ số dòng $I_{\text{on}}/I_{\text{off}}$ thấp), gây rò rỉ dòng điện nghiêm trọng trong mạch logic kỹ thuật số. Đồng thời, các hiệu ứng tiếp xúc tại ranh giới dị thể kim loại - graphene (M-G interface) và cơ chế phản ứng quang - điện tử dưới tác động của thế tuần hoàn ngoài vẫn còn tồn tại nhiều khoảng trống lý thuyết vi mô chưa được làm sáng tỏ. Luận án tiến sĩ với đề tài "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ử" được thực hiện nhằm giải quyết trực tiếp các điểm nghẽn học thuật này thông qua các câu hỏi nghiên cứu và giả thuyết khoa học xác định:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế vi mô nào chi phối sự biến đổi cấu trúc vùng năng lượng, sự định xứ trạng thái điện tử $p_z$ và đặc trưng hấp thụ quang học của graphene khi chịu tác động của thế vô hướng tuần hoàn một chiều trong cấu trúc siêu mạng graphene (GSLs)?
  • Câu hỏi nghiên cứu 2 (RQ2): Bản chất tương tác cơ học lượng tử tại mặt phân giới tiếp xúc kim loại - graphene (M-G) ảnh hưởng như thế nào đến sự phân bố mật độ trạng thái $p_z$-DOS, cấu trúc dịch mức Fermi và điện trở tiếp xúc nội tại?
  • Câu hỏi nghiên cứu 3 (RQ3): Sự tích hợp giữa mô hình hàm tự năng tiếp xúc $\Sigma_{G-M}$ và hình thức luận hàm Green không cân bằng (NEGF) có thể mô tả chính xác các đặc trưng truyền dẫn tĩnh điện và đặc tuyến $I_{DS}-V_{GS}$, $I_{DS}-V_{DS}$ của linh kiện GFET kênh ngắn ($L_c = 40 - 100\text{ nm}$) có điện cực cổng dây nano GaN như thế nào?

Khung lý thuyết của luận án được xây dựng vững chắc dựa trên phương pháp gần đúng liên kết chặt (Tight-Binding Approximation - TBA), phương trình sóng tựa Dirac tương đối tính, hình thức luận đáp ứng tuyến tính Kubo cho tính toán độ dẫn quang, và phương pháp hàm Green không cân bằng (NEGF) tự phối với phương trình Poisson tĩnh điện.

Phạm vi nghiên cứu bao quát từ cấu trúc tinh thể nguyên tử, dải nano graphene (GNRs) biên zigzag và armchair, cấu trúc siêu mạng (A-GSLs, Z-GSLs) với rào thế $U_b$, tiếp xúc của graphene với mạng lập phương tâm mặt FCC (111) và lục giác xếp chặt HCP (0001) của hàng loạt kim loại ($\text{Cu, Au, Pt, Pd, Ti, Ag, Al, Ir, Co, Cd, Ru}$), đến việc mô phỏng linh kiện GFET thực tế với chiều dài kênh $L_c = 40 - 100\text{ nm}$. Các đột phá của luận án đã đóng góp trực tiếp vào gói phần mềm mô phỏng lượng tử OPEDEVS (Opto-Electronic Devices Simulation), đồng thời tạo ra 3 công bố quốc tế uy tín trên các tạp chí chuyên ngành hàng đầu như Applied Physics Letters và Journal of Physics: Condensed Matter.

Literature Review và Positioning

Lý thuyết cấu trúc vùng năng lượng của mạng tinh thể lục giác graphene bắt đầu từ công trình kinh điển của P. R. Wallace (1946) [98] khi xem xét graphene như một mô hình cơ sở cho than chì. Trong suốt hơn năm thập kỷ, các định lý nhiệt động học của Peierls và Landau [76, 105] đã khiến giới vật lý tin rằng màng tinh thể hai chiều không thể tồn tại ổn định ở trạng thái tự do do các thăng giáng nhiệt sẽ phá hủy trật tự tinh thể tầm xa. Định kiến này bị phá vỡ hoàn toàn vào năm 2004 khi nhóm Andre Geim và Konstantin Novoselov tại Đại học Manchester phân lập thành công đơn lớp graphene bằng phương pháp bóc tách cơ học [69]. Khám phá này đã kích hoạt làn sóng nghiên cứu bùng nổ trên toàn cầu [150, 152].

Wallace (1946)               Geim & Novoselov (2004)            IBM / Liao et al. (2010)
Cấu trúc dải năng lượng   ──> Bóc tách đơn lớp graphene   ───>   GFET siêu cao tần (>100 GHz)
dạng nón Dirac lý thuyết     và hiệu ứng trường thực nghiệm       Cổng tự sắp xếp GaN nanowire
         │                                │                                   │
         ▼                                ▼                                   ▼
Novoselov et al. (2005)      Tworzydło et al. (2006)             Luận án này (TS. Đỗ Vân Nam lab)
Lý thuyết Klein tunneling    Độ dẫn cực tiểu 4e²/(πh)            Mô hình TBA-NEGF đóng gói OPEDEVS,
và Hall dị thường            thực nghiệm vs lý thuyết           giải mã GSLs, tiếp xúc M-G & GFET

Trong dòng chảy lý thuyết truyền dẫn lượng tử, hiện tượng chui ngầm Klein (Klein tunneling) đối với các fermion Dirac tương đối tính đã được Katsnelson, Geim và Novoselov (2006) [85] tiên đoán và sau đó được Young và Kim (2009) [11] kiểm chứng thực nghiệm, xác nhận rằng hạt tải trong graphene có thể xuyên qua các rào thế tĩnh điện cao và rộng với xác suất gần bằng 1 ở góc tới trực diện.

Một tranh luận học thuật sâu sắc khác liên quan đến giới hạn độ dẫn điện cực tiểu của graphene: trong khi các phép đo thực nghiệm của Novoselov et al. (2005) [70] ghi nhận giá trị điện trở lượng tử cực đại xấp xỉ giá trị $R_Q = h/4e^2$ (tương đương độ dẫn $\sigma_{\min} \approx 4e^2/h$), công trình lý thuyết của Tworzydło et al. (2006) [57] dựa trên công thức Landauer lại chứng minh rằng độ dẫn suất giới hạn lượng tử thuần túy trong giới hạn không mất trật tự là $\sigma_{\min} = 4e^2/(\pi h)$, được hình thành bởi sự truyền dẫn qua các trạng thái mờ (evanescent modes) với vector sóng ảo $ik$.

Về phương diện điều biến cấu trúc năng lượng, các nghiên cứu siêu mạng graphene (GSLs) của Park et al. (2008) và Brey & Fertig (2009) [19-24, 62] đã chỉ ra khả năng làm biến dạng nón Dirac và tạo ra các điểm Dirac phụ dưới thế tuần hoàn. Tuy nhiên, cơ chế định xứ vi mô của các hàm sóng trạng thái $p_z$ và ảnh hưởng trực tiếp của chúng lên độ dẫn quang học phụ thuộc phân cực photon chưa được mô hình hóa một cách tường minh và thống nhất.

Đối với vấn đề tiếp xúc kim loại - graphene, các công trình thực nghiệm của Huard et al. (2008) [120] và tính toán nguyên lý đầu ab-initio của Giovannetti et al. (2008) đã chỉ ra hiện tượng dịch chuyển mức Fermi và sự hình thành các chuyển tiếp $p-n$ hoặc $p-p$ dưới chân điện cực. Mặc dù vậy, các tính toán mật độ hàm mật độ (DFT) từ nguyên lý đầu thường đòi hỏi tài nguyên tính toán khổng lồ và không thể mở rộng để mô phỏng sự vận chuyển lượng tử của toàn bộ linh kiện GFET kích thước hàng chục đến hàng trăm nanomet.

Ở hướng phát triển linh kiện thực nghiệm, nhóm nghiên cứu tại IBM (2010) [147, 149] đã đạt kỷ lục tần số cắt $100\text{ GHz}$ cho transistor graphene kích thước wafer $2 \times 2\text{ cm}^2$, và nhóm Lei Liao et al. (2010) [82] đã công bố linh kiện GFET hiệu năng cao sử dụng dây nano GaN ($\kappa \approx 10$) làm điện cực cổng tự sắp xếp để giảm thiểu điện trở ký sinh.

Luận án định vị chính xác tại điểm giao cắt giữa vật lý chất rắn lý thuyết và mô phỏng linh kiện ứng dụng: kế thừa các nền tảng giải tích của Wallace, Tworzydło và Katsnelson, đồng thời phát triển một khung mô hình hiệu dụng Tight-Binding kết hợp hàm tự năng tiếp xúc $\Sigma_{G-M}$ và hình thức luận NEGF. Cách tiếp cận này vượt qua giới hạn quy mô của các tính toán DFT nguyên lý đầu (như VASP hay SIESTA), khắc phục điểm yếu của các mô hình bán cổ điển (vốn bỏ qua tính chất pha lượng tử), và tái lập hoàn hảo các dữ liệu đo đạc thực nghiệm từ các công trình quốc tế của Lei Liao et al. và Huard et al.

Đó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à hoàn thiện lý thuyết cấu trúc điện tử và động học lượng tử của vật liệu carbon hai chiều thông qua các đóng góp cụ thể:

  1. Mở rộng mô hình Tight-Binding cho cấu trúc siêu mạng graphene (GSLs): Thiết lập ma trận Hamiltonian lượng tử hóa lần hai tổng quát cho các siêu mạng A-GSLs (hướng armchair) và Z-GSLs (hướng zigzag) dưới tác động của điện thế tuần hoàn $U(x)$ biên độ $U_b$. Luận án chứng minh sự tái cấu trúc của các mặt đẳng năng lượng, sự hình thành các nón Dirac bất đối xứng mới và hiện tượng ghim mặt năng lượng (energy surface pinning) dọc theo các phương vector sóng $k_x, k_y$.
  2. Khám phá hiện tượng định xứ dị thường (exotic localization): Xác định bằng chứng giải tích và số học về sự định xứ không gian kỳ lạ của hàm sóng điện tử $p_z$ tại các vị trí cực trị của hàm thế trong một chu kỳ siêu mạng, phá vỡ tính chất giải định xứ truyền thống của mạng graphene phẳng.
  3. Lý thuyết hóa cơ chế suy giảm độ dẫn quang và hiệu ứng khóa Pauli: Phân tích vi mô các phần tử ma trận chuyển mức quang lượng tử $M_{vc}(\mathbf{k}) = \langle \psi_v | \hat{v} | \psi_c \rangle$ dựa trên hình thức luận Kubo. Luận án chứng minh sự sụt giảm mạnh mẽ của độ dẫn quang $\sigma(\omega)$ trong miền năng lượng photon $(0, U_b)$, đồng thời giải thích tường minh cơ chế ngăn chặn chuyển dời điện tử liên dải do hiện tượng khóa Pauli (Pauli blocking) cùng tính bất đối xứng theo trạng thái phân cực photon dọc/ngang.
  4. Mô hình hóa tương tác vi mô tại mặt phân giới Kim loại - Graphene (M-G): Xây dựng mô hình tán xạ lượng tử dựa trên việc ghép nối ma trận Hamiltonian riêng phần giữa màng graphene và các bề mặt kim loại FCC (111) ($\text{Cu, Au, Pt, Ag, Al, Ir}$) và HCP (0001) ($\text{Co, Cd, Ru, Ti}$). Mô hình hóa thành công tương tác lai hóa orbitan $p_z - d$ và $p_z - s$ thông qua ma trận tự năng tiếp xúc $\Sigma_{G-M} = \Lambda - i\Gamma$, cho phép lượng hóa chính xác sự dịch mức Fermi và điện trở tiếp xúc nội tại mà không cần giải toàn bộ phương trình Kohn-Sham đa hạt.
                  ┌─────────────────────────────────────────┐
                  │   Khung Phân Tích Tích Hợp Đa Cấp Độ    │
                  └────────────────────┬────────────────────┘
                                       │
         ┌─────────────────────────────┼─────────────────────────────┐
         ▼                             ▼                             ▼
┌──────────────────┐          ┌──────────────────┐          ┌──────────────────┐
│   Cấp Vi Mô      │          │  Cấp Trung Gian  │          │   Cấp Linh Kiện  │
│  (Atomic Scale)  │          │  (Mesoscopic)    │          │  (Device Level)  │
├──────────────────┤          ├──────────────────┤          ├──────────────────┤
│• DFT / VASP      │          │• TBA Hamiltonian │          │• Package OPEDEVS │
│• Lai hóa p_z - d │ ───────> │• Hình thức Kubo  │ ───────> │• NEGF + Poisson  │
│• p_z-DOS,        │          │• Khóa Pauli      │          │• Đặc tuyến       │
│  dịch Fermi      │          │• Tự năng Σ_{G-M} │          │  I_DS - V_GS, V_DS│
└──────────────────┘          └──────────────────┘          └──────────────────┘

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 chặt chẽ:

$$\begin{aligned} \text{Hamiltonian Liên kết chặt (TBA)} &\longrightarrow \text{Hình thức luận Đáp ứng Tuyến tính Kubo (Quang học)} \ &\searrow \text{Hình thức luận Hàm Green phi cân bằng NEGF (Truyền dẫn GFET)} \end{aligned}$$

Hệ thống phương trình vi phân và tích phân lượng tử được thiết lập với các điều kiện biên tường minh:

  • Đối với dải nano GNRs: Sử dụng điều kiện biên triệt tiêu hàm sóng tại các mép dải biên zigzag ($M_{\text{zline}}$) và armchair ($M_{\text{aline}}$), kết hợp ma trận biến đổi pha Peierls Landau Gauge $\vec{A} = (-By, 0, 0)$ để lượng hóa hiệu ứng từ trường đều $B$ tác động trực giao.
  • Đối với tiếp xúc M-G: Thiết lập phương trình truyền qua Landauer-Büttiker đa kênh với hàm tự năng tiếp xúc $\Sigma_{S,D}(E)$ mô tả chính xác hiệu ứng mở rộng mức năng lượng $\Gamma = \text{Im}[\Sigma]$ và dịch chuyển mức năng lượng tĩnh điện $\Lambda = \text{Re}[\Sigma]$.
  • Điều kiện biên linh kiện GFET: Giải phương trình Poisson 2D/3D tự phối (self-consistent loop) với điện thế đặt tại cổng $V_{GS}$, nguồn-máng $V_{DS}$, phản ánh hằng số điện môi cao $\kappa \sim 10$ của dây nano GaN.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng định lượng (Positivist-Rationalist Epistemology), kết hợp các mô hình giải tích chính xác của vật lý lượng tử với các thuật toán mô phỏng số hiệu năng cao. Thiết kế nghiên cứu phân tầng đa quy mô (multi-scale simulation design) được triển khai qua ba cấp độ liên hoàn:

  1. Cấp độ nguyên tử (Atomic Scale): Khảo sát cấu trúc mạng tinh thể tổ ong gồm hai mạng con tam giác đối xứng $A$ và $B$ với khoảng cách liên kết $a_{CC} \approx 1.42\text{ \AA}$, diện tích ô cơ sở $S_0 = \frac{3\sqrt{3}}{2}a_{CC}^2$, tính toán lai hóa $sp^2$ (liên kết $\sigma$) và orbitan $2p_z$ (liên kết $\pi$, năng lượng nhảy lân cận gần nhất $t_{CC} \approx -2.67\text{ eV}$, tích phân xen phủ $E_{pp\pi} \approx -2.7\text{ eV}$).
  2. Cấp độ trung gian (Mesoscopic Scale): Xây dựng ma trận Hamiltonian kích thước lớn cho siêu mạng GSLs và cấu trúc kim loại - graphene - kim loại (M-G-M), tính toán tích phân ma trận chuyển mức quang và mật độ trạng thái điện tử (DOS).
  3. Cấp độ linh kiện (Device Scale): Tích hợp các tham số vi mô vào module GFET của gói phần mềm OPEDEVS để mô phỏng quá trình tiêm hạt tải và đặc trưng Von-Ampe của transistor thực tế.
                         THUẬT TOÁN TỰ PHỐI NEGF - POISSON
                               (Package OPEDEVS)
                                       │
                                       ▼
                       ┌───────────────────────────────┐
                       │  Giả thiết thế ban đầu V(r)   │
                       └───────────────┬───────────────┘
                                       │
                                       ▼
                       ┌───────────────────────────────┐
                       │  Tính tự năng Σ_S(E), Σ_D(E)  │
                       │     và tự năng tiếp xúc Σ_GM  │
                       └───────────────┬───────────────┘
                                       │
                                       ▼
                       ┌───────────────────────────────┐
                       │  Tính hàm Green phi cân bằng  │
                       │  G(E) = [EI - H - Σ_S - Σ_D]⁻¹│
                       └───────────────┬───────────────┘
                                       │
                                       ▼
                       ┌───────────────────────────────┐
                       │  Tính mật độ hạt tải n(r)     │
                       │     và nồng độ lỗ trống p(r)  │
                       └───────────────┬───────────────┘
                                       │
                                       ▼
                       ┌───────────────────────────────┐
                       │   Giải phương trình Poisson   │
                       │     ∇·[ε(r)∇V(r)] = -ρ(r)     │
                       │     để tìm thế mới V_new(r)   │
                       └───────────────┬───────────────┘
                                       │
                                       ▼
                              /─────────────────\
                             < |V_new - V| < 10⁻⁵? >
                              \─────────────────/
                                  │           │
                                 No          Yes
                                  │           │
                                  └───┐       ▼
                                      │  ┌───────────────────────────────┐
                                      └──│  Tính dòng I_DS (Landauer)    │
                                         │  Xuất đặc tuyến I-V & độ dẫn G│
                                         └───────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình tính toán và mô phỏng được chuẩn hóa với độ tin cậy và tính chặt chẽ học thuật cao:

  • Tính toán cấu trúc dải và phổ quang: Áp dụng biến đổi Fourier từ không gian thực sang không gian mạng đảo $(\mathbf{k}$-space) trong vùng Brillouin thứ nhất hình lục giác với diện tích $\Omega_{BZ} = \frac{8\pi^2}{3\sqrt{3}a_{CC}^2}$. Khai triển Hamiltonian tựa Dirac quanh hai điểm đối xứng cao $K, K' = \left(\frac{2\pi}{3a_{CC}}, \pm \frac{2\pi}{3\sqrt{3}a_{CC}}\right)$:

    $$\hat{H}_K = v_F \vec{\sigma} \cdot \vec{p} = v_F (\hat{p}_x \sigma_x + \hat{p}_y \sigma_y)$$

    trong đó $\sigma_{x,y,z}$ là các ma trận Pauli. Hàm mật độ trạng thái điện tử (DOS) ở lân cận điểm Dirac được chuẩn hóa giải tích theo biểu thức:

    $$\rho(E) = \frac{4|E|}{\pi (\hbar v_F)^2}$$

  • Mô hình hóa từ trường bằng pha Peierls: Khi đặt từ trường đều $B$ vuông góc với mặt dải nano, năng lượng nhảy $t_{CC}$ được điều chỉnh thông qua tích phân Peierls:

    $$T_{ij}(B) = t_{CC} \exp\left( i \frac{e}{\hbar} \int_{\mathbf{R}_i}^{\mathbf{R}_j} \vec{A} \cdot d\vec{r} \right)$$

    với hàm thế vector Landau $\vec{A} = (-By, 0, 0)$, cho phép phân tích sự hình thành các mức Landau và sự tách dải năng lượng trong dải nano biên zigzag và armchair.

  • Thuật toán tự phối NEGF - Poisson trong OPEDEVS:

    1. Khởi tạo thế tĩnh điện ban đầu $V^{(0)}(\vec{r})$.

    2. Xây dựng ma trận Hamiltonian của kênh dẫn graphene $H$ kết hợp ma trận tự năng của các điện cực tiếp xúc kim loại $\Sigma_S(E), \Sigma_D(E)$ và tự năng tương tác $\Sigma_{G-M}$.

    3. Tính toán hàm Green trễ (retarded Green's function):

      $$G(E) = \left[ (E + i\eta)I - H - \Sigma_S(E) - \Sigma_D(E) - \Sigma_{G-M} \right]^{-1}$$

    4. Xác định hàm mật độ phổ hạt tải $A(E) = i[G(E) - G^\dagger(E)]$ và tích phân để tính mật độ điện tử $n(\vec{r})$, nồng độ lỗ trống $p(\vec{r})$.

    5. Giải phương trình Poisson không gian 3 chiều: $\nabla \cdot [\varepsilon(\vec{r}) \nabla V(\vec{r})] = -\rho(\vec{r}) = -e[p(\vec{r}) - n(\vec{r}) + N_D^+ - N_A^-]$ để tìm thế nghiệm mới $V^{(k+1)}(\vec{r})$.

    6. Lặp chu trình tự phối cho đến khi đạt tiêu chuẩn hội tụ sai số thế tĩnh điện $\max |V^{(k+1)} - V^{(k)}| < 10^{-5}\text{ eV}$.

    7. Tính toán dòng điện truyền dẫn qua kênh bằng công thức tích phân Landauer:

      $$I_{DS} = \frac{2e}{h} \int_{-\infty}^{+\infty} T(E) [f_S(E - \mu_S) - f_D(E - \mu_D)] dE$$

      với hàm truyền qua $T(E) = \text{Tr}[\Gamma_S G \Gamma_D G^\dagger]$.

Data và phân tích

  • Mẫu cấu trúc hình học: Chiều dài kênh dẫn $L_c = 40\text{ nm}, 60\text{ nm}, 100\text{ nm}$; chu kỳ siêu mạng $N = 2N_1 = 30$ (đối với A-GSLs) và $N = 2N_1 = 40$ (đối với Z-GSLs); điện áp thế tuần hoàn $U_b$ biến thiên từ $0$ đến $6U_0$.
  • Tham số tiếp xúc kim loại: Khảo sát các tổ hợp $\text{Cu-G-Cu}$, $\text{Au-G-Au}$, $\text{Pt-G-Pt}$, $\text{Pd-G-Pd}$, $\text{Ti-G-Ti}$ trên các ô mạng tương thích chứa 2 nguyên tử (như Cu-FCC, Co-HCP) và 8 nguyên tử (như Ag, Al, Ir, Pt, Au-FCC; Cd, Ru, Ti-HCP). Sử dụng mã VASP 4 để trích xuất hàm mật độ trạng thái chiếu $p_z$-DOS và các tích phân nhảy năng lượng $t_{pz-s}$, $t_{pz-d}$.
  • Kiểm tra tính vững (Robustness & Validation): Kết quả tính toán dòng điện cực tiểu $I_{\min}$, dòng cực đại $I_{\max}$, độ dẫn $G(V_{GS})$ và đặc tuyến Von-Ampe $I_{DS}-V_{DS}$ từ module GFET của gói OPEDEVS được đối chiếu trực tiếp với các số liệu thực nghiệm đo đạc từ nhóm nghiên cứu quốc tế của Lei Liao et al. [82], cho thấy sự trùng khớp tuyệt đối về xu hướng vật lý và độ lớn bậc dòng điện.

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

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

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        4 PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                 │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. ĐỊNH XỨ KỲ LẠ TRONG SIÊU MẠNG (GSLs)                                               │
│    Hàm sóng p_z bị giam giữ dị thường tại các vùng cực trị của thế tuần hoàn U(x);   │
│    Tái tạo các nón Dirac bất đối xứng và ghim mặt năng lượng tại biên vùng Brillouin. │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 2. DẬP TẮT ĐỘ DẪN QUANG DO KHÓA PAULI                                                  │
│    Độ dẫn quang σ(ω) triệt tiêu hoàn toàn trong dải photon (0, U_b) do hiệu ứng      │
│    khóa Pauli (Pauli blocking); Độ dẫn phụ thuộc phân cực ánh sáng rõ rệt.             │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 3. PHÂN HÓA BẢN CHẤT TIẾP XÚC KIM LOẠI - GRAPHENE (M-G)                                │
│    • Nhóm tương tác yếu (Cu, Au): Nón Dirac bảo toàn, điện trở tiếp xúc cao.           │
│    • Nhóm lai hóa mạnh (Ti, Pd, Pt): Mất tính chất nón Dirac, dịch Fermi sâu,           │
│      tạo kênh tiêm hạt tải ohmic dẫn điện lý tưởng.                                   │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 4. LÀM RÕ CƠ CHẾ DÒNG CỰC TIỂU VÀ ĐẶC TRƯNG GFET KÊNH NGẮN                              │
│    • Re[Σ_{G-M}] quyết định sự dịch điểm Dirac và điện áp trung hòa điện tích V_Dirac.│
│    • Im[Σ_{G-M}] quy định độ mở rộng mức năng lượng và giá trị dòng rò cực tiểu I_min.│
└────────────────────────────────────────────────────────────────────────────────────────┘
  1. Hiện tượng định xứ kỳ lạ (exotic localization) của trạng thái điện tử trong siêu mạng GSLs: Khác với tính chất chuyển động tự do của điện tử trong graphene nguyên bản, sự áp đặt một trường thế vô hướng tuần hoàn một chiều đã bẻ gãy tính đối xứng chuẩn, làm xuất hiện sự định xứ mạnh mẽ của các trạng thái điện tử $p_z$ tại các vị trí nhất định trong chu kỳ thế tĩnh điện. Mật độ xác suất tìm thấy điện tử $P_n(k_y, x)$ (với $k_x = 0, n = 1, 2, 3, 4$) cho thấy các điện tử bị giam giữ dị thường trong các giếng thế, tạo ra các đỉnh cộng hưởng sắc nét trên phổ mật độ trạng thái (DOS) trong phạm vi năng lượng cỡ $1\text{ eV}$. Cấu trúc dải năng lượng của A-GSLs và Z-GSLs xuất hiện sự ghim lại (pinning) của các mặt tán sắc dọc theo phương $k_x, k_y$, đồng thời hình thành nên các nón Dirac mới với độ dốc vận tốc nhóm dị hướng.
  2. Sự dập tắt độ dẫn quang trong miền năng lượng photon $(0, U_b)$ do hiệu ứng khóa Pauli: Tính toán độ dẫn quang động học $\sigma(\omega)$ thông qua hình thức luận Kubo đã phát hiện sự suy giảm nghiêm trọng của khả năng hấp thụ quang học trong cửa sổ năng lượng photon $\hbar\omega \in (0, U_b)$. Phân tích các phần tử ma trận chuyển mức quang lượng tử $M_{vc}$ khẳng định cơ chế ngăn chặn chuyển dời điện tử liên dải từ dải hóa trị $\pi$ lên dải dẫn $\pi^*$ do trạng thái đích đã bị chiếm đóng bởi các điện tử bị điều biến thế năng tĩnh điện (cơ chế khóa Pauli - Pauli blocking tương tự như trường hợp graphene bị pha tạp với ngưỡng $2E_F$). Đặc biệt, độ dẫn quang thể hiện sự phụ thuộc chọn lọc mạnh mẽ vào trạng thái phân cực của photon (phân cực song song hoặc vuông góc với phương điều biến thế).
  3. Quy luật phân hóa cơ chế tiếp xúc Kim loại - Graphene (M-G): Phân tích cấu trúc dải điện tử và $p_z$-DOS của các tổ hợp tiếp xúc dị thể M-G cho thấy sự phân chia thành hai nhóm kim loại rõ rệt:
    • Nhóm tương tác yếu (như Cu, Au): Các orbitan $d$ của kim loại nằm sâu dưới mức Fermi, liên kết giữa graphene và kim loại chủ yếu là tương tác van der Waals yếu kết hợp dịch chuyển điện tích nhỏ; cấu trúc nón Dirac của graphene vẫn được bảo toàn gần như nguyên vẹn, dẫn đến rào cản thế năng tiếp xúc và điện trở tiếp xúc nội tại tương đối lớn.
    • Nhóm tương tác hóa học / lai hóa mạnh (như Ti, Pd, Pt): Có sự xen phủ và lai hóa mạnh mẽ giữa các orbitan $d$ chưa lấp đầy của kim loại với orbitan $p_z$ của graphene; cấu trúc nón Dirac bị phá hủy hoàn toàn tại vùng tiếp xúc, mức Fermi bị dịch sâu vào vùng dẫn hoặc vùng hóa trị (tạo pha tạp kiểu $n$ hoặc $p$ cực mạnh), mở ra mật độ trạng thái dẫn điện lớn và tạo điều kiện tiêm điện tử hiệu quả với điện trở suất tiếp xúc cực thấp.
  4. Tác động định lượng của tự năng tiếp xúc $\Sigma_{G-M}$ lên đặc tuyến linh kiện GFET: Mô phỏng tự phối NEGF - Poisson trên cấu trúc GFET kênh ngắn ($L_c = 40\text{ nm}, 60\text{ nm}$) chứng minh rằng:
    • Phần thực của tự năng $\text{Re}[\Sigma_{G-M}]$ (chẳng hạn khi khảo sát tại giá trị $\text{Re}[\Sigma_{G-M}] = -0.1\text{ eV}$) quyết định mức độ uốn cong thế năng tĩnh điện và sự dịch chuyển vị trí điểm trung hòa điện tích (điểm Dirac) trên trục điện áp cổng $V_{GS}$.
    • Phần ảo $\text{Im}[\Sigma_{G-M}] = -\Gamma$ (mô tả thời gian sống và tốc độ tán xạ/tiêm hạt tải từ điện cực) chi phối trực tiếp độ rộng mở mức năng lượng, làm tăng giá trị dòng cực tiểu $I_{\min}$ tại điểm Dirac và làm mềm độ dốc chuyển mạch trong đặc tuyến truyền dẫn $I_{DS}-V_{GS}$ và đặc tuyến ra $I_{DS}-V_{DS}$.

Implications đa chiều

  • Ý nghĩa học thuật và lý thuyết: Cung cấp bức tranh toàn diện và sâu sắc về hành vi cơ học lượng tử của fermion Dirac hai chiều dưới tác động đồng thời của trường thế ngoài, từ trường trực giao và ranh giới tiếp xúc kim loại. Đặt nền móng lý thuyết vững chắc cho ngành quang điện tử học graphene (Graphene Optoelectronics).
  • Đổi mới phương pháp luận tính toán: Module GFET trong gói phần mềm OPEDEVS do nhóm phát triển đã thiết lập một công cụ mô phỏng lượng tử mạnh mẽ, kết hợp hiệu quả giữa ma trận thưa, thuật toán đệ quy NEGF và giải Poisson tự phối, cho phép tính toán chính xác các linh kiện nano kích thước hàng chục nanomet với chi phí thời gian và bộ nhớ giảm hàng trăm lần so với các phần mềm DFT ab-initio truyền thống.
  • Ứng dụng thực tiễn trong công nghiệp bán dẫn: Cung cấp hướng dẫn công nghệ mang tính quyết định cho các kỹ sư thiết kế linh kiện: khuyến nghị lựa chọn các kim loại tiếp xúc nhóm lai hóa mạnh ($\text{Ti, Pd}$) để chế tạo điện cực nguồn/máng nhằm tối thiểu hóa điện trở tiếp xúc ký sinh, hoặc sử dụng nhóm kim loại tiếp xúc yếu ($\text{Au, Cu}$) khi cần bảo tồn độ linh động Dirac; định hướng tích hợp điện môi cổng hằng số cao (như dây nano GaN $\kappa \approx 10$) để tăng cường khả năng kiểm soát tĩnh điện kênh dẫn.

Limitations và Future Research

  • Giới hạn về mô hình tán xạ không đàn hồi: Mô hình NEGF hiện tại chủ yếu khảo sát quá trình truyền dẫn đạn đạo (ballistic transport) và tán xạ đàn hồi do thế tĩnh điện/điện cực, chưa tính toán đầy đủ tương tác tán xạ electron - phonon ở nhiệt độ phòng ($T = 300\text{ K}$) cũng như tán xạ do tạp chất tích điện ngẫu nhiên trong lớp oxit đế.
  • Giới hạn về cấu trúc biên thực tế: Các tính toán dải nano GNRs và siêu mạng GSLs giả định cấu trúc biên lý tưởng hoàn hảo (zigzag hoặc armchair thuần túy). Trong thực tế chế tạo bằng phương pháp khắc lithography hóa học, sự mất trật tự biên (edge disorder) và các khuyết tật mạng (vacancies, Stone-Wales) có thể làm biến tính phổ năng lượng.
  • Định hướng nghiên cứu tương lai:
    1. Tích hợp các cơ chế tán xạ không đàn hồi electron - phonon và tán xạ quang học phonon vào thuật toán tự phối của gói OPEDEVS để mô phỏng chính xác hiện tượng bão hòa vận tốc hạt tải ở điện trường cao.
    2. Mở rộng khung lý thuyết để khảo sát các hiệu ứng nhiệt điện lượng tử (Seebeck effect) và hiệu ứng spin-orbit coupling trong các cấu trúc dị thể Van der Waals hai chiều thế hệ mới (như Graphene/h-BN hoặc Graphene/TMDs).
    3. Hoàn thiện giao diện đồ họa người dùng (GUI) và tối ưu hóa tính toán song song đa luồng (MPI/OpenMP) cho module GFET trong OPEDEVS, hướng tới thương mại hóa công cụ mô phỏng linh kiện nano.

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

  • Tác động học thuật quốc tế: Các kết quả cốt lõi của luận án đã được công bố trên các tạp chí vật lý chuyên ngành hàng đầu thế giới của Viện Vật lý Hoa Kỳ (AIP - Applied Physics Letters) và Viện Vật lý Vương quốc Anh (IOP - Journal of Physics: Condensed Matter). Các công trình này đã nhận được sự quan tâm, trích dẫn và sử dụng trực tiếp bởi nhiều nhóm nghiên cứu lý thuyết và thực nghiệm quốc tế, khẳng định trình độ hội nhập học thuật cao của nghiên cứu vật lý chất rắn tại Việt Nam.
  • Chuyển đổi công nghệ và công nghiệp bán dẫn: Cung cấp cơ sở khoa học vi mô cho các tập đoàn công nghệ (như Intel, IBM, TSMC) trong việc đánh giá tính khả thi của vật liệu graphene trong lộ trình thay thế silicon (Beyond-CMOS roadmap). Các kết quả về độ dẫn quang và siêu mạng GSLs mở ra triển vọng chế tạo các bộ tách sóng quang siêu nhạy, bộ điều biến quang terahertz và pin quang điện thế hệ mới.
  • Đóng góp phát triển phần mềm nội địa: Khẳng định năng lực tự chủ công nghệ tính toán thông qua việc xây dựng và hoàn thiện gói phần mềm mã nguồn chuyên dụng OPEDEVS, giảm thiểu sự phụ thuộc vào các bộ phần mềm thương mại đắt đỏ của nước ngoài.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận hệ thống phương pháp luận hoàn chỉnh từ lý thuyết lượng tử hóa lần hai, biến đổi Peierls, hình thức luận Kubo đến giải thuật NEGF phức tạp; sử dụng trực tiếp các công thức giải tích và ma trận Hamiltonian đã được chuẩn hóa.
  • Kỹ sư R&D Linh kiện Bán dẫn & Vi điện tử: Nắm bắt quy luật tương tác mặt phân giới kim loại - graphene để lựa chọn kim loại chế tạo lớp tiếp xúc ohmic hoặc Schottky tối ưu; hiểu rõ tác động của điện dung cổng và hằng số điện môi dây nano GaN lên tần số cắt của GFET.
  • Nhà phát triển Phần mềm Mô phỏng Lượng tử: Khai thác thuật toán tự phối NEGF - Poisson, kỹ thuật ma trận thưa và module GFET trong package OPEDEVS để ứng dụng cho các vật liệu hai chiều tương tự (như Silicene, Germanene, Phosphorene, MoS₂).

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 thiết lập mô hình vi mô giải thích sự định xứ dị thường của điện tử $p_z$ và sự dập tắt độ dẫn quang trong siêu mạng graphene (GSLs) dựa trên hình thức luận Kubo và hiệu ứng khóa Pauli. Luận án đã mở rộng trực tiếp lý thuyết trường Dirac tương đối tính của Wallace (1946) và phương trình tựa Dirac hai thành phần spinor cho hệ tuần hoàn ngoài một chiều, chứng minh rằng điện thế vô hướng tuần hoàn $U(x)$ không chỉ làm biến dạng vận tốc Fermi mà còn tạo ra các điểm Dirac phụ bất đối xứng và ngăn chặn quá trình chuyển dời quang học liên dải trong cửa sổ photon $(0, U_b)$.

2. Đổi mới phương pháp luận nổi bật so với các nghiên cứu quốc tế trước đây là gì?

So với các nghiên cứu DFT thuần túy (như Giovannetti et al., 2008 - vốn bị giới hạn kích thước tính toán dưới vài trăm nguyên tử) và các mô hình bán cổ điển (vốn bỏ qua tính chất pha lượng tử), luận án đã sáng tạo phương pháp tiếp cận đa quy mô (multi-scale): trích xuất tương tác orbitan vi mô từ DFT/VASP để xây dựng ma trận tự năng hiệu dụng $\Sigma_{G-M}$, sau đó tích hợp vào thuật toán NEGF tự phối Poisson trong package OPEDEVS. Phương pháp này cho phép mô phỏng chính xác đặc tính truyền dẫn lượng tử của toàn bộ cấu trúc linh kiện GFET kích thước thực tế ($L_c = 40 - 100\text{ nm}$) với tốc độ tính toán vượt trội.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi số liệu nào?

Phát hiện bất ngờ nhất là sự phân hóa đối nghịch sâu sắc giữa hai nhóm kim loại tiếp xúc: Kim loại quý như Cu và Au dù có độ dẫn điện khối tuyệt vời nhưng lại tạo tiếp xúc kém hiệu quả với graphene do tương tác yếu, giữ nguyên nón Dirac và gây điện trở tiếp xúc lớn; ngược lại, các kim loại hoạt tính như Ti, Pd dù có cấu trúc dải phức tạp nhưng nhờ sự lai hóa mạnh giữa orbitan $d$ và $p_z$ đã phá hủy nón Dirac tại vùng tiếp xúc, tạo mật độ trạng thái $p_z$-DOS khổng lồ tại mức Fermi, biến vùng tiếp xúc thành kênh tiêm hạt tải ohmic hoàn hảo.

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

Có. Luận án cung cấp chi tiết toàn bộ cơ sở ma trận Hamiltonian của mạng graphene, dải nano zigzag/armchair, siêu mạng GSLs, bảng ma trận vận tốc chuyển mức quang, các giá trị tích phân nhảy năng lượng $t_{pz-s}, t_{pz-d}$ cho từng cặp kim loại (Bảng 3.1) và hướng dẫn chi tiết quy trình cài đặt, thiết lập tham số đầu vào cho module GFET trong gói phần mềm OPEDEVS (Phần Phụ lục).

5. Lộ trình nghiên cứu 10 năm được phác thảo như thế nào?

Lộ trình hướng tới ba mục tiêu chiến lược: (1) Nâng cấp package OPEDEVS tích hợp đầy đủ tán xạ electron - phonon nhiệt độ phòng và các hiệu ứng tán xạ quang học; (2) Mở rộng mô phỏng sang các linh kiện spin-FET và van điều khiển thung lũng (valleytronics) dựa trên siêu mạng dị thể graphene xoắn góc ma thuật (twisted bilayer graphene); (3) Hợp tác với các phòng thí nghiệm thực nghiệm vi chế tạo để hiện thực hóa linh kiện GFET cổng dây nano GaN siêu cao tần hoạt động trong dải sóng terahertz ($>100\text{ GHz}$).

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 kết luận khoa học cốt lõi:

  1. Hệ thống hóa toàn diện và phát triển các mô hình giải tích cấu trúc vùng năng lượng, hệ thức tán sắc nón Dirac, mật độ trạng thái (DOS) và hiệu ứng từ trường Peierls trên các màng đơn lớp graphene và dải nano GNRs biên zigzag/armchair.
  2. Khám phá bản chất hiện tượng định xứ dị thường của điện tử $p_z$ và sự ghim mặt năng lượng trong siêu mạng graphene (GSLs) chịu tác động của thế tuần hoàn một chiều.
  3. Lý giải tường minh cơ chế dập tắt độ dẫn quang động học $\sigma(\omega)$ trong miền photon $(0, U_b)$ thông qua hình thức luận Kubo và hiệu ứng khóa Pauli, xác lập quy luật chọn lọc phân cực photon.
  4. Xây dựng thành công mô hình ma trận tự năng $\Sigma_{G-M}$ mô tả chính xác tương tác lai hóa orbitan $p_z - d/s$, phân loại rõ nét cơ chế tiếp xúc ohmic/barrier giữa graphene và các bề mặt kim loại FCC (111), HCP (0001).
  5. Phát triển hoàn thiện module GFET trong gói phần mềm mô phỏng lượng tử OPEDEVS, thực hiện thành công việc mô phỏng tự phối NEGF - Poisson cho linh kiện GFET kênh ngắn cổng GaN nanowire, làm rõ vai trò quyết định của tự năng tiếp xúc $\Sigma_{G-M}$ lên đặc tuyến Von-Ampe và kiểm chứng chuẩn xác với các dữ liệu thực nghiệm quốc tế.

Nghiên cứu đã khẳng định vị thế tiên phong trong việc làm chủ lý thuyết truyền dẫn lượng tử tiên tiến, đóng góp quan trọng vào kho tàng tri thức vật lý chất rắn tính toán và mở ra những triển vọng to lớn cho sự phát triển của công nghệ vi điện tử - quang điện tử thế hệ mới.