Tổng quan luận án

Luận án tiến sĩ của nghiên cứu sinh Lê Hoàng Anh, thực hiện dưới sự hướng dẫn khoa học của PGS. TS. Nguyễn Đức Chiến tại Đại học Bách Khoa Hà Nội (năm 2015), thuộc chuyên ngành Vật lý kỹ thuật (mã số: 62520401), tập trung giải quyết các bài toán vật lý lý thuyết và mô phỏng số lượng tử liên quan đến trạng thái điện tử, tính chất quang học và cơ chế truyền dẫn qua tiếp xúc của vật liệu graphene trong các cấu trúc vi mô và linh kiện bán dẫn nano.

Tính cấp thiết và khoảng trống nghiên cứu

Ngành công nghiệp vi điện tử và bán dẫn dựa trên nền tảng silicon (Si) đã đạt đến các giới hạn vật lý cơ bản khi quá trình thu nhỏ kích thước linh kiện (scaling) theo định luật Moore tiến vào thang nanomet (dưới 22 nm và hướng tới 10 nm). Tại thang kích thước này, các hiệu ứng kênh ngắn (short-channel effects), sự gia tăng dòng rò do chui ngầm lượng tử, và khó khăn trong việc giải phóng nhiệt lượng đã tạo nên rào cản kỹ thuật lớn. Graphene—vật liệu hai chiều cấu tạo từ một lớp đơn nguyên tử carbon sắp xếp theo mạng tổ ong—sở hữu những tính chất vật lý đặc biệt như độ linh động hạt tải cao ($\mu_G \approx 200.000\text{ cm}^2/\text{V}\cdot\text{s}$), vận tốc Fermi lớn ($v_F \approx 10^6\text{ m/s}$, xấp xỉ $c/300$), độ bền cơ học cao và độ trong suốt quang học lớn. Tuy nhiên, việc ứng dụng graphene vào các linh kiện điện tử và quang điện tử đòi hỏi phải giải quyết các khoảng trống nghiên cứu lớn:

  1. Cơ chế biến điệu trạng thái điện tử và phản ứng quang học của graphene khi chịu tác động của các trường thế ngoài tuần hoàn (cấu trúc siêu mạng - GSLs).
  2. Bản chất tương tác liên kết điện tử và điện trở tiếp xúc nội tại tại mặt phân giới giữa graphene và các điện cực kim loại (M-G).
  3. Thiếu vắng các công cụ mô phỏng lượng tử tự nhất quán có khả năng mô tả chính xác linh kiện transistor hiệu ứng trường graphene (GFETs) có tính đến ảnh hưởng của tiếp xúc kim loại và cấu trúc cổng tự sắp xếp.

Mục tiêu nghiên cứu

Luận án xác định hai mục đích cơ bản:

  1. Nghiên cứu các tính chất cơ bản của điện tử bên trong các màng graphene dưới tác động của các điều kiện khác nhau từ bên ngoài để từ đó xem xét tiềm năng ứng dụng của loại vật liệu này trong lĩnh vực điện tử và quang điện tử.
  2. Thực hiện các nghiên cứu cho phép đóng góp tới việc khai thác và phát triển chương trình tính toán OPEDEVS (cụ thể là module GFETs).

Bài toán và câu hỏi nghiên cứu

Để thực hiện mục tiêu trên, luận án chia nội dung nghiên cứu thành hai bài toán lớn:

  • Bài toán 1: Nghiên cứu bức tranh cấu trúc vùng năng lượng của điện tử bên trong graphene trong các điều kiện nội tại và xử lý/chế tác khác nhau (khảo sát điện tử $2p_z$ trong graphene lý tưởng, dải nano graphene biên zigzag/armchair, và siêu mạng graphene chịu thế tuần hoàn một chiều GSLs).
  • Bài toán 2: Khảo sát các tính chất truyền dẫn của các trạng thái kích thích điện tử bên trong các màng graphene đóng vai trò kênh dẫn điện trong một số cấu trúc linh kiện điển hình (cơ chế tiêm điện tử qua mặt tiếp xúc kim loại - graphene và đặc trưng truyền dẫn dòng điện của linh kiện GFETs).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Mạng tinh thể graphene đơn lớp; dải nano graphene (GNRs) với hai dạng biên điển hình (zigzag và armchair); cấu trúc siêu mạng graphene chịu thế tuần hoàn 1D (A-GSLs và Z-GSLs); hệ tiếp xúc kim loại - graphene (M-G) với các kim loại mạng lập phương tâm mặt (FCC: Cu, Ag, Al, Ir, Pt, Au) và mạng lục giác xếp chặt (HCP: Co, Cd, Ru, Ti, Pd); linh kiện transistor trường GFETs có kênh dẫn dưới 100 nm sử dụng điện cực cổng dây nano GaN.
  • Phạm vi nghiên cứu: Tiếp cận bằng phương pháp lý thuyết giải tích và tính toán mô phỏng số lượng tử (gần đúng liên kết chặt Tight-Binding, phương trình Dirac, hình thức luận Kubo, phương pháp hàm Green không cân bằng NEGF kết hợp giải phương trình Poisson tự nhất quán). Luận án không tiến hành chế tạo thực nghiệm trực tiếp mà sử dụng các số liệu thực nghiệm công bố quốc tế để đối sánh mô hình.

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

Các hướng nghiên cứu trước trong tài liệu

  • Cấu trúc điện tử nền tảng của graphene: P. R. Wallace (1946) là người đầu tiên tính toán cấu trúc dải năng lượng của graphene từ gần đúng liên kết chặt để giải thích các thuộc tính của than chì. Trước năm 2004, các lý thuyết nhiệt động lực học (Peierls, Landau) cho rằng màng 2D không thể tồn tại độc lập do thăng giáng nhiệt phá hủy mạng tinh thể.
  • Thực nghiệm chế tạo và cô lập graphene: Ruoff và cộng sự (Đại học Northwestern, 1999) chế tạo đĩa mỏng than chì; Kim và cộng sự (Đại học Columbia) dùng đầu nhọn AFM bóc tách màng mỏng than chì khoảng 10 lớp. K. S. Novoselov, A. K. Geim và cộng sự (2004, 2005) đã bóc tách thành công đơn lớp graphene bằng băng dính cơ học, công bố trên Science, quan sát hiệu ứng trường và hiệu ứng Hall lượng tử dị thường ở nhiệt độ phòng. W. A. de Heer và cộng sự (Georgia Tech, 2004) phát triển kỹ thuật nung SiC ở $1300^\circ\text{C}$ để làm bay hơi Si và tạo màng graphene epitaxy.
  • Hiện tượng lượng tử và truyền dẫn Dirac: M. I. Katsnelson, K. S. Novoselov, A. K. Geim (2006) và sau đó là A. F. Young, P. Kim (2009) kiểm chứng nghịch lý chui ngầm Klein (Klein tunneling) của các fermion Dirac không khối lượng. J. Tworzydło và cộng sự (2006) đưa ra công thức Landauer xác định độ dẫn cực tiểu tổng quát của graphene $\sigma_{min} = 4e^2/(\pi h)$ thông qua các trạng thái sóng mờ (evanescent modes). R. R. Nair và cộng sự (2008) xác định độ hấp thụ ánh sáng của graphene liên hệ trực tiếp với hằng số cấu trúc tinh tế $\pi \alpha \approx 2,3%$.
  • Siêu mạng graphene (GSLs): C. X. Bai và X. D. Zhang (2007) khảo sát truyền qua rào thế và chui ngầm cộng hưởng; C. H. Park và cộng sự (2008) phát hiện tính truyền dẫn dị hướng mạnh của vận tốc nhóm và sự xuất hiện các điểm Dirac phụ (vệ tinh); J. H. Ho và cộng sự (2008), L. Brey và H. A. Fertig (2009) nghiên cứu thế tuần hoàn dạng cosin và Kronig-Penney.
  • Linh kiện GFETs và tiếp xúc kim loại: Nhóm IBM (2009, 2010) chế tạo GFETs đạt tần số cắt $f_T = 26\text{ GHz}$ và $100\text{ GHz}$ trên wafer $2 \times 2\text{ cm}^2$. Lei Liao và cộng sự (2010) chế tạo GFETs tự sắp xếp dùng dây nano GaN đạt độ hỗ dẫn $g_m > 2000\text{ S/m}$, dòng $I_{ON} \approx 5\text{ mA/mm}$ và $f_T \approx 300 - 1000\text{ GHz}$. J. Chauhan và J. Guo (2011) mô phỏng lượng tử NEGF cho GFETs kênh 20 nm. Khảo sát tiếp xúc M-G được thực hiện bởi P. A. Khomyakov (2009), C. Gong (2010) từ các tính toán nguyên lý đầu.
                  LỊCH SỬ PHÁT TRIỂN NGHIÊN CỨU GRAPHENE
1946 - 2004: Giai đoạn lý thuyết & tìm kiếm              2004: Khám phá đột phá
- Wallace (1946): Cấu trúc dải năng lượng              - Novoselov & Geim: Bóc tách đơn lớp
- Landau, Peierls: Giả thuyết bất ổn định 2D           - Quan sát hiệu ứng trường & Hall lượng tử
- Ruoff, Kim, de Heer: Nỗ lực tạo màng mỏng            - Khẳng định tính chất hạt Dirac 2D
2005 - 2015: Mở rộng nghiên cứu & Ứng dụng linh kiện
- Tworzydło: Độ dẫn cực tiểu 4e²/πh
- Katsnelson, Young, Kim: Chui ngầm Klein
- Park, Brey, Fertig: Siêu mạng graphene (GSLs)
- IBM, Lei Liao: Linh kiện GFETs cao tần GHz/THz
- Vị trí luận án: Làm rõ trạng thái định xứ GSLs,
  mô hình hóa tiếp xúc M-G và phát triển module GFETs (OPEDEVS)

Vị trí và khoảng trống luận án giải quyết

Kế thừa các kết quả nghiên cứu của nhóm nghiên cứu tại Đại học Bách Khoa Hà Nội (do TS. Đỗ Vân Nam dẫn dắt), luận án định vị vào việc làm rõ:

  1. Bản chất của các trạng thái định xứ dị thường trong siêu mạng GSLs xuất hiện bên ngoài giếng thế và nhúng trong phổ liên tục, cùng quy luật suy giảm độ dẫn quang trong dải năng lượng photon $(0, U_b)$.
  2. Xây dựng mô hình hiệu dụng liên kết chặt kết hợp dữ liệu DFT (VASP 4) để tính toán điện trở suất nội tại của lớp tiếp xúc giữa graphene và các cấu trúc kim loại FCC/HCP.
  3. Hoàn thiện module GFETs trong gói phần mềm OPEDEVS nhằm mô phỏng tự nhất quán Poisson-NEGF cho cấu trúc transistor cổng dây nano GaN có xét đến tự năng tiếp xúc nguồn/máng.

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

Cơ sở lý thuyết và khung phân tích

  • Mô hình gần đúng liên kết chặt (Tight-Binding Approximation - TBA): Mô tả trạng thái của các điện tử $2p_z$ tạo nên liên kết $\pi$ và $\pi^*$ trong mạng tinh thể graphene. Khoảng cách giữa hai nguyên tử carbon lân cận là $a_{CC} = 1,42\text{ \AA}$, độ dài vector cơ sở mạng thực $a = \sqrt{3} a_{CC} \approx 2,46\text{ \AA}$. Năng lượng nhảy bậc một giữa hai nguyên tử lân cận gần nhất là $t_{CC} \approx 2,7\text{ eV}$.
  • Phương trình Dirac cho hạt Fermion 2D không khối lượng: Khai triển Hamiltonian liên kết chặt lân cận hai điểm Dirac đối xứng $K$ và $K'$ trong vùng Brillouin thứ nhất dẫn đến hệ thức tán sắc tuyến tính dạng nón: $$E(k) = \pm \hbar v_F |\mathbf{k}|$$ Hamiltonian hiệu dụng viết dưới dạng ma trận Pauli: $$\hat{H}_K = v_F \boldsymbol{\sigma} \cdot \mathbf{p} = v_F (\sigma_x p_x + \sigma_y p_y)$$
  • Hình thức luận Kubo: Dùng để tính toán tenxơ độ dẫn quang học động học $\sigma_{\alpha\beta}(\omega)$ từ các phần tử ma trận xung lượng chuyển dời quang học giữa các dải năng lượng bị biến điệu trong siêu mạng.
  • Hình thức luận hàm Green không cân bằng (NEGF): Dùng để tính toán truyền dẫn lượng tử đạn đạo qua kênh dẫn và tiếp xúc, trong đó ma trận hàm Green được xác định bởi: $$G(E) = \left[ (E + i\eta)I - H - \Sigma_S - \Sigma_D \right]^{-1}$$ với $\Sigma_S, \Sigma_D$ là các ma trận tự năng mô tả tương tác của kênh dẫn với các điện cực nguồn (Source) và máng (Drain).

Phương pháp tính toán và nguồn dữ liệu

  • Tính toán cấu trúc điện tử và mật độ trạng thái (DOS): Chéo hóa ma trận Hamiltonian xây dựng từ hình thức luận lượng tử hóa lần 2. Sử dụng phần mềm nguyên lý đầu VASP 4 để tính $p_z$-DOS cho hệ tiếp xúc kim loại - graphene.
  • Mô phỏng tự nhất quán linh kiện: Thuật toán giải lặp tự nhất quán giữa phương trình tĩnh điện Poisson 2D/3D (tính thế tĩnh điện $V(\mathbf{r})$ từ phân bố điện tích hạt tải) và phương trình NEGF (tính nồng độ điện tử và lỗ trống từ hàm Green).
  • Công cụ tính toán: Tự xây dựng gói script giải ma trận dải nano GNRs và phát triển module GFETs tích hợp trong gói phần mềm mô phỏng linh kiện quang điện tử OPEDEVS (Opto-Electronic Devices Simulation).
  • Dữ liệu đối sánh: Số liệu đo đạc thực nghiệm đặc trưng $I_{DS}-V_{GS}$, $I_{DS}-V_{DS}$ và độ dẫn của nhóm Lei Liao (2010), cùng các công bố của nhóm IBM và Novoselov.

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa cơ sở vật lý của vật liệu carbon và graphene:

  • Phân tích cấu hình nguyên tử carbon ($1s^2 2s^2 2p^2$), sự lai hóa $sp^2$ tạo 3 liên kết $\sigma$ đồng phẳng tạo góc $120^\circ$ và orbital $2p_z$ vuông góc mặt phẳng tạo liên kết $\pi$ quyết định tính chất dẫn điện.
  • Thiết lập cấu trúc mạng tinh thể tổ ong gồm 2 mạng con tam giác lồng nhau A và B. Vector cơ sở mạng thực: $$\mathbf{a}1 = a{CC}\left(\frac{3}{2}, \frac{\sqrt{3}}{2}\right), \quad \mathbf{a}2 = a{CC}\left(\frac{3}{2}, -\frac{\sqrt{3}}{2}\right)$$ Diện tích ô cơ sở $S = \frac{3\sqrt{3}}{2} a_{CC}^2$. Vector cơ sở mạng đảo $\mathbf{b}1, \mathbf{b}2$, xác định vùng Brillouin thứ nhất dạng lục giác đều với diện tích $\Omega = \frac{8\pi^2}{3\sqrt{3}a{CC}^2}$ và các điểm đối xứng cao $\Gamma(0,0)$, $K\left(\frac{2\pi}{3a{CC}}, \frac{2\pi}{3\sqrt{3}a_{CC}}\right)$, $K'\left(\frac{2\pi}{3a_{CC}}, -\frac{2\pi}{3\sqrt{3}a_{CC}}\right)$, $M\left(\frac{2\pi}{3a_{CC}}, 0\right)$.
  • Nhóm đối xứng không gian tinh thể là nhóm $D_{6h} = D_3 \times {E, \sigma_h} \times {E, \sigma_i}$.
  • Giải bài toán trị riêng từ Hamiltonian liên kết chặt thu được hệ thức tán sắc tổng quát: $$E(\mathbf{k}) = \pm t_{CC} \sqrt{1 + 4\cos\left(\frac{3k_x a_{CC}}{2}\right)\cos\left(\frac{\sqrt{3}k_y a_{CC}}{2}\right) + 4\cos^2\left(\frac{\sqrt{3}k_y a_{CC}}{2}\right)}$$
  • Khảo sát cấu trúc vùng năng lượng của dải nano graphene (GNRs): dải biên zigzag luôn xuất hiện các trạng thái biên phẳng (edge states) tại mức năng lượng $E=0$ tập trung mật độ ở hai rìa; dải biên armchair thể hiện tính chất kim loại khi số chuỗi dimer $N = 3m - 1$ ($m$ nguyên) hoặc tính bán dẫn với độ rộng vùng cấm phụ thuộc nghịch đảo vào độ rộng dải. Tác giả đã lập trình đóng gói một công cụ phần mềm tính toán cấu trúc vùng cho đơn lớp và hai lớp GNRs.
  • Tổng quan về ứng dụng GFETs: phân tích nguyên nhân tỷ số đóng ngắt ($I_{ON}/I_{OFF}$) của GFETs thấp do không có vùng cấm (zero bandgap), dẫn đến sự chuyển đổi liên tục giữa kênh $p$ và kênh $n$ khi quét điện áp cổng $V_{GS}$ và tạo ra độ dẫn cực tiểu. Đặc tuyến $I_{DS}-V_{DS}$ thể hiện điểm uốn cong không bão hòa hoàn toàn do hiệu ứng chui ngầm Klein và trở vi phân âm. Đánh giá giải pháp dùng dây nano GaN của Lei Liao ($g_m > 2000\text{ S/m}$, $f_T \approx 1000\text{ GHz}$) và mô phỏng GFETs 20 nm của Jyotsna Chauhan ($C_{ins} = 486\text{ nF/cm}^2$).

Chương 2: Cấu trúc điện tử và tính chất quang của siêu mạng graphene

Chương 2 trình bày các tính toán cho cấu trúc siêu mạng graphene (GSLs) chịu tác động của thế năng tĩnh điện tuần hoàn 1D dạng bậc thang (Kronig-Penney) có chu kỳ $L$ và biên độ rào thế $U_b$:

  • Khảo sát hai cấu hình: Armchair-GSLs (A-GSLs, thế biến thiên dọc phương armchair) và Zigzag-GSLs (Z-GSLs, thế biến thiên dọc phương zigzag).
  • Ô đơn vị siêu mạng được mở rộng chứa $N$ dải nguyên tử ($N_1$ dải thuộc vùng giếng thế và $N - N_1$ dải thuộc vùng rào thế). Cấu trúc vùng năng lượng của GSLs được giải thích như sự chồng chập của các cấu trúc dải năng lượng dải nano tương ứng với 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 $U_b$.
  • Phát hiện tính chất định xứ kỳ lạ: Khảo sát hàm sóng điện tử $p_z$ cho thấy xuất hiện các trạng thái điện tử bị biến điệu mạnh, biên độ hàm sóng bị dập tắt gần như hoàn toàn trong vùng rào thế hoặc giếng thế. Hiện tượng định xứ này có hai đặc điểm dị thường:
    1. Xảy ra tại các mức năng lượng nằm bên ngoài miền thế năng ($E > U_b$), tức là điện tử bị định xứ ngay ở phía trên đỉnh rào thế.
    2. Các trạng thái định xứ không nằm tách biệt trong một cửa sổ năng lượng cấm mà bị 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 (extended states).
  • Tính chất quang học của GSLs: Tính toán độ dẫn quang động $\sigma(\omega)$ theo hình thức luận Kubo. Kết quả chỉ ra:
    1. Độ dẫn quang của GSLs phụ thuộc mạnh vào trạng thái phân cực của ánh sáng tới (dị hướng quang học) do thế ngoài 1D đã phá vỡ tính đối xứng quay của mạng tinh thể.
    2. Xuất hiện sự suy giảm rõ rệt của độ dẫn quang trong miền năng lượng photon thấp $(0, U_b)$. Nguyên nhân vật lý do 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ủ lưỡng cực chuyển mức quang học giữa các trạng thái dải hóa trị và dải dẫn.

Chương 3: Sự truyền dẫn điện tử qua bề mặt tiếp xúc kim loại-graphene

Chương 3 tập trung giải quyết bài toán tương tác và truyền dẫn qua mặt phân giới kim loại - graphene (M-G):

  • Phân loại các bề mặt tiếp xúc tinh thể:
    • Kim loại mạng lập phương tâm mặt (FCC) mặt (111): Cu (mô hình ô đơn vị graphene chứa 2 nguyên tử C); Ag, Al, Ir, Pt, Au (ô đơn vị chứa 8 nguyên tử C).
    • Kim loại mạng lục giác xếp chặt (HCP) mặt (0001): Co (ô đơn vị chứa 2 nguyên tử C); Cd, Ru, Ti (ô đơn vị chứa 8 nguyên tử C).
  • Sử dụng phần mềm VASP 4 tính toán cấu trúc dải điện tử và mật độ trạng thái $p_z$-DOS cho các tổ hợp G-Cu và G-Ti, làm rõ sự dịch chuyển mức Fermi và mức độ lai hóa giữa orbital $d$ của kim loại với orbital $p_z$ của graphene.
  • Xây dựng mô hình hiệu dụng liên kết chặt kết hợp hình thức luận NEGF để tính toán dòng truyền dẫn $I-V$ và xác suất truyền qua $T(E, k)$ cho các cấu trúc tiếp xúc hai đầu kim loại M-G-M (Cu-G-Cu, Au-G-Au, Pt-G-Pt, Pd-G-Pd, Ti-G-Ti).
          MÔ HÌNH TIẾP XÚC KIM LOẠI - GRAPHENE - KIM LOẠI (M-G-M)
   Tự năng tiếp xúc          Hamiltonian Dirac/TBA     Tự năng tiếp xúc
    Σ_S(E) = Re + i·Im            H_k (sp²-pz)          Σ_D(E) = Re + i·Im
                      Giải ma trận Hàm Green (NEGF):
                 G(E) = [(E + iη)I - H - Σ_S - Σ_D]⁻¹
                     Đặc trưng truyền dẫn & Dòng I-V:
               I = (2e/h) ∫ T(E) [f_S(E) - f_D(E)] dE

Bảng 1: Bảng tổng hợp các nhóm kim loại tiếp xúc với graphene khảo sát trong luận án

Nhóm mạng kim loại Mặt tinh thể tiếp xúc Kim loại khảo sát Số nguyên tử C trong ô cơ sở mô hình Đặc trưng tương tác điện tử
FCC (Lập phương tâm mặt) (1 1 1) Cu 2 nguyên tử Tương tác yếu, giữ cấu trúc nón Dirac, dịch Fermi
FCC (Lập phương tâm mặt) (1 1 1) Ag, Al, Ir, Pt, Au 8 nguyên tử Lai hóa orbital, dịch chuyển mức Fermi
HCP (Lục giác xếp chặt) (0 0 0 1) Co 2 nguyên tử Tương tác mạnh, mở rộng mức năng lượng
HCP (Lục giác xếp chặt) (0 0 0 1) Cd, Ru, Ti, Pd 8 nguyên tử Lai hóa mạnh $d-p_z$, làm biến dạng nón Dirac

Chương 4: Mô phỏng linh kiện GFETs

Chương 4 trình bày việc phát triển và ứng dụng công cụ mô phỏng linh kiện GFETs:

  • Cấu trúc mô phỏng: Transistor hiệu ứng trường graphene có độ dài kênh dẫn $L_c = 40\text{ nm}, 60\text{ nm}$ và $50 - 100\text{ nm}$, sử dụng điện cực cổng là dây nano GaN pha tạp $n$ cao ($n \approx 2 \times 10^{19}\text{ cm}^{-3}$), hằng số điện môi $\kappa \approx 10$, tạo tiếp xúc Schottky ngăn dòng rò với kênh dẫn graphene. Các điện cực nguồn và máng bằng kim loại Pt ($10\text{ nm}$) được chế tạo theo cơ chế tự sắp xếp.
  • Phát triển gói phần mềm: Thiết kế và tích hợp module GFETs vào package OPEDEVS. Nền tảng tính toán sử dụng ma trận thưa, thuật toán đệ quy giải hệ tự nhất quán 2 phương trình: phương trình Poisson xác định phân bố thế năng tĩnh điện $U(\mathbf{r})$ và phương trình NEGF xác định mật độ hạt tải và dòng truyền dẫn $I_{DS}$.
  • Đưa ảnh hưởng tiếp xúc kim loại vào mô hình: Tương tác M-G được đưa vào thông qua ma trận tự năng tiếp xúc $\Sigma_{G-M} = \text{Re}(\Sigma_{G-M}) + i,\text{Im}(\Sigma_{G-M})$, trong đó phần thực đại diện cho độ dịch thế năng hóa học tại tiếp xúc và phần ảo đại diện cho sự mở rộng mức do tán xạ tiêm hạt tải.
  • Kết quả tính toán:
    • Xác định biên dạng thế tĩnh điện và mật độ điện tử/lỗ trống dọc kênh dẫn tại các điện áp $V_{DS}, V_{GS}$ khác nhau.
    • Đường cong độ dẫn $G(V_{GS})$ và đặc tuyến $I_{DS}-V_{GS}$ giải thích rõ điểm Dirac thực nghiệm bị dịch chuyển do ảnh hưởng của $\text{Re}(\Sigma_{G-M})$.
    • Đặc tuyến $I_{DS}-V_{DS}$ tại các giá trị $V_{GS}$ khác nhau cho thấy sự tương đồng cao với kết quả đo đạc thực nghiệm của nhóm Lei Liao.
    • Làm rõ vai trò của phần ảo $\Gamma = \text{Im}(\Sigma_{G-M})$ làm mở rộng các trạng thái tán xạ và tăng cường độ dẫn cực tiểu trong kênh dẫn.

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

Đóng góp mới về mặt lý luận và khoa học

  1. Phát hiện quy luật định xứ dị thường trong siêu mạng graphene (GSLs): Chứng minh sự tồn tại của các trạng thái điện tử $p_z$ bị định xứ hoàn toàn trong một chu kỳ thế nhưng lại nằm ở mức năng lượng cao hơn bờ rào thế ($E > U_b$) và nhúng trong dải phổ liên tục của các trạng thái mở rộng.
  2. Lý giải cơ chế quang học vi mô của GSLs: Giải thích định lượng sự suy giảm độ dẫn quang trong miền năng lượng photon $(0, U_b)$ xuất phát từ sự suy giảm các yếu tố ma trận chuyển dời quang học do hàm sóng bị biến điệu, đồng thời xác lập tính dị hướng quang học phụ thuộc vào trạng thái phân cực photon của GSLs.
  3. Mô hình hóa thành công liên kết tiếp xúc kim loại - graphene: Đề xuất mô hình hiệu dụng mô tả liên kết điện tử giữa graphene và các bề mặt kim loại thuộc cấu trúc FCC và HCP, cho phép tính toán đặc trưng truyền dẫn $I-V$, xác suất truyền qua và đánh giá điện trở suất nội tại của lớp tiếp xúc M-G.

Đóng góp mới về mặt kỹ thuật và thực tiễn

  1. Xây dựng phần mềm mô phỏng cấu trúc dải nano graphene: Đóng gói chương trình tính toán cấu trúc vùng năng lượng cho dải nano đơn lớp và hai lớp với các cấu hình biên zigzag và armchair.
  2. Phát triển module GFETs trong gói phần mềm OPEDEVS: Mở rộng thành công package OPEDEVS (do TS. Đỗ Vân Nam sáng lập) bằng việc bổ sung module mô phỏng lượng tử tự nhất quán Poisson-NEGF dành riêng cho linh kiện GFETs, tích hợp tham số tự năng tiếp xúc kim loại $\Sigma_{G-M}$.

Bảng 2: Tổng hợp các phát hiện khoa học và công cụ phát triển trong luận án

Đối tượng nghiên cứu Phương pháp / Công cụ Phát hiện / Kết quả chính Công bố quốc tế liên quan
Siêu mạng GSLs (A-GSLs, Z-GSLs) TBA, Toán tử lượng tử 2, Hình thức luận Kubo - Trạng thái định xứ dị thường ($E > U_b$, nhúng trong phổ liên tục)
- Suy giảm độ dẫn quang trong miền $(0, U_b)$, phụ thuộc phân cực photon		 Appl. Phys. Lett.
J. Phys.: Condens. Matter
Tiếp xúc Kim loại - Graphene (M-G) Gần đúng TBA, DFT (VASP 4), NEGF - Mô hình hóa ô mạng tiếp xúc FCC (Cu, Au, Pt, Ag, Al, Ir) và HCP (Co, Ti, Cd, Ru, Pd)
- Xác định dòng $I-V$ và điện trở nội tại		 Appl. Phys. Lett. 101, 161605 (2012)
Linh kiện GFETs cổng GaN Nanowire Poisson - NEGF tự nhất quán (Module GFETs - OPEDEVS) - Mô phỏng biên dạng thế, mật độ hạt tải, đặc tuyến $I_{DS}-V_{GS}$, $I_{DS}-V_{DS}$ kênh $40-60\text{ nm}$
- Làm rõ ảnh hưởng của $\text{Re}(\Sigma_{G-M})$ và $\text{Im}(\Sigma_{G-M})$		 Đóng gói trong package OPEDEVS

Kiến nghị và đề xuất

  • Về mặt thiết kế linh kiện: Đề xuất lựa chọn các kim loại tiếp xúc có độ tích hợp tốt và kiểm soát thế hóa học tại tiếp xúc để giảm thiểu điện trở tiếp xúc ký sinh trong GFETs.
  • Về công nghệ chế tạo: Tận dụng các cấu trúc cổng tự sắp xếp và chất điện môi có hằng số điện môi cao ($\kappa \ge 10$) như dây nano GaN nhằm tối ưu hóa độ hỗ dẫn $g_m$ và nâng cao tần số làm việc của linh kiện RF.

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

Hạn chế của luận án

  • Mô phỏng GFETs ở Chương 4: Tác giả tự đánh giá các kết quả mô phỏng GFETs mới dừng lại ở các phân tích ban đầu (dáng điệu thế, nồng độ hạt tải, các đường đặc trưng volt-ampe cơ bản dưới sự thay đổi tham số tự năng tiếp xúc) mà chưa có điều kiện phân tích sâu sắc toàn diện mọi hiệu ứng vật lý để công bố trên các tạp chí quốc tế chuyên ngành trong thời hạn thực hiện luận án.
  • Giả định mô hình: Các tính toán truyền dẫn lượng tử trong luận án chủ yếu dựa trên gần đúng liên kết chặt một điện tử và truyền dẫn đạn đạo (ballistic transport), chưa tính toán đầy đủ ảnh hưởng của các cơ chế tán xạ phi đàn hồi phức tạp như tán xạ electron-phonon ở nhiệt độ cao, bẫy điện tích ngẫu nhiên trong đế điện môi hoặc các sai hỏng mạng tinh thể thực tế.

Hướng nghiên cứu tiếp

  • Hoàn thiện và tối ưu hóa module GFETs trong gói OPEDEVS để xử lý các bài toán động học hạt tải phức tạp hơn, bao gồm tán xạ phonon và hiệu ứng nhiệt độ.
  • Mở rộng nghiên cứu sang các cấu trúc tiếp xúc dị thể hai chiều khác (vd: graphene kết hợp với vật liệu chuyển tiếp dichalcogenide - TMDs) và các cấu trúc linh kiện spintronic/valleytronic trên nền graphene.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Vật lý kỹ thuật, Vật lý chất rắn, Khoa học và công nghệ nano, Vật liệu điện tử. Luận án cung cấp hệ thống giải tích toán học tường minh từ gần đúng liên kết chặt đến phương trình Dirac và hình thức luận Kubo.
  • Kỹ sư và chuyên gia nghiên cứu vi điện tử: Các nhóm thiết kế linh kiện nano và thiết bị cao tần (RF) tìm kiếm giải pháp vật lý để cải thiện tiếp xúc kim loại - bán dẫn và khai thác vật liệu 2D.

Nội dung đáng tham khảo

  • Phương pháp luận tính toán cấu trúc dải năng lượng GNRs và GSLs: Phần dẫn giải chi tiết cách thiết lập ma trận Hamiltonian và ma trận vận tốc trong không gian mạng đảo (trình bày tại Chương 1, Chương 2 và Phụ lục A, B).
  • Quy trình mô phỏng Poisson - NEGF: Thuật toán tự nhất quán kết hợp ma trận tự năng tiếp xúc kim loại $\Sigma_{G-M}$ trong phân tích linh kiện nano (trình bày tại Chương 4 và Phụ lục C).

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

1. Vì sao transistor hiệu ứng trường kênh dẫn graphene (GFETs) kích thước lớn lại có tỷ số đóng-ngắt dòng điện ($I_{ON}/I_{OFF}$) rất thấp?

Nguyên nhân vật lý xuất phát từ cấu trúc vùng năng lượng của graphene là chất bán kim không có vùng cấm (zero-bandgap), với hệ thức tán sắc tuyến tính chạm nhau tại các điểm Dirac. Khi thay đổi điện áp cổng $V_{GS}$, kênh dẫn chuyển đổi liên tục giữa hai loại hạt tải (điện tử và lỗ trống) mà không thể đóng hoàn toàn dòng điện, tạo ra một độ dẫn điện cực tiểu nội tại $\sigma_{min} = 4e^2/(\pi h)$. Do đó dòng $I_{OFF}$ luôn duy trì ở mức cao, khiến linh kiện không phù hợp cho mạch logic số CMOS nhưng thích hợp cho các ứng dụng khuếch đại cao tần (RF).

2. Hai đặc điểm dị thường của trạng thái định xứ điện tử trong siêu mạng graphene (GSLs) là gì?

Trong cấu trúc siêu mạng graphene chịu thế tuần hoàn 1D, các trạng thái định xứ của điện tử $p_z$ thể hiện hai đặc tính khác thường so với cơ học lượng tử truyền thống:

  1. Trạng thái định xứ xuất hiện ngay cả ở các mức năng lượng cao hơn biên độ rào thế ($E > U_b$), tức nằm phía trên mặt giếng thế thay vì bị giam giữ bên trong giếng.
  2. Các trạng thái này không tách biệt trong một vùng cấm năng lượng riêng biệt mà bị nhúng trực tiếp trong dải phổ năng lượng liên tục của các trạng thái truyền qua (trạng thái mở rộng).

3. Trong mô hình tiếp xúc kim loại - graphene (M-G), các kim loại FCC và HCP được phân loại như thế nào?

Luận án phân loại kim loại tiếp xúc thành hai họ cấu trúc tinh thể chính:

  • Họ mạng lập phương tâm mặt (FCC) tiếp xúc trên mặt (111): Cu (mô hình ô đơn vị graphene chứa 2 nguyên tử C); Ag, Al, Ir, Pt, Au (ô đơn vị chứa 8 nguyên tử C).
  • Họ mạng lục giác xếp chặt (HCP) tiếp xúc trên mặt (0001): Co (ô đơn vị chứa 2 nguyên tử C); Cd, Ru, Ti (ô đơn vị chứa 8 nguyên tử C).

4. Gói phần mềm OPEDEVS sử dụng phương pháp tính toán nền tảng nào và được phát triển thêm nội dung gì trong luận án?

OPEDEVS (Opto-Electronic Devices Simulation, do TS. Đỗ Vân Nam thiết kế) sử dụng phương pháp hàm Green không cân bằng (NEGF) kết hợp kỹ thuật tính toán ma trận thưa, đệ quy và tính toán song song. Trong luận án, tác giả đã phát triển và tích hợp thêm module GFETs, cho phép giải tự nhất quán hệ phương trình Poisson - NEGF để mô phỏng phân bố thế, mật độ hạt tải và đặc tuyến truyền dẫn $I-V$ của transistor graphene có tính đến tự năng tiếp xúc kim loại.

5. Cấu trúc GFETs của nhóm Lei Liao được mô phỏng trong Chương 4 có đặc điểm cấu tạo nào giúp đạt tần số cắt cao?

Linh kiện của Lei Liao sử dụng điện cực cổng bằng một dây nano GaN pha tạp loại $n$ cao ($n \approx 2 \times 10^{19}\text{ cm}^{-3}$) có hằng số điện môi cao ($\kappa \approx 10$), tiếp xúc với graphene tạo rào thế Schottky ngăn dòng rò cổng. Đồng thời, cấu trúc này kết hợp kỹ thuật tự sắp xếp (self-alignment) để tạo hai điện cực nguồn và máng Pt ($10\text{ nm}$) sát mép cổng dây nano, giúp triệt tiêu điện trở ký sinh lối vào và đạt độ hỗ dẫn $g_m > 2000\text{ S/m}$, đưa tần số ngưỡng nội tại lên tới cấp độ THz.

Kết luận

Luận án của NCS. Lê Hoàng Anh là công trình nghiên cứu hệ thống về vật lý lý thuyết và mô phỏng lượng tử các tính chất điện tử, quang học và truyền dẫn của vật liệu graphene cùng các cấu trúc linh kiện liên quan. Tác giả đã làm sáng tỏ hiện tượng định xứ dị thường và cơ chế suy giảm độ dẫn quang trong siêu mạng GSLs, đồng thời xây dựng thành công mô hình hiệu dụng tính toán tiếp xúc kim loại - graphene cho các họ mạng FCC và HCP. Bên cạnh các đóng góp lý thuyết được công bố trên các tạp chí quốc tế uy tín (Appl. Phys. Lett., J. Phys.: Condens. Matter), luận án đã phát triển thành công module GFETs trong gói phần mềm OPEDEVS, tạo tiền đề tính toán quan trọng cho các nghiên cứu chuyên sâu về linh kiện nano cao tần tiếp theo.