Giới thiệu dự án

Sự phát triển của ngành công nghiệp vi điện tử bán dẫn toàn cầu đang đối mặt với giới hạn vật lý khắc nghiệt của định luật Moore. Khi kích thước bóng bán dẫn (transistor) silicon truyền thống thu hẹp xuống dưới ngưỡng $5\text{ nm}$, các hiệu ứng lượng tử không mong muốn như dòng rò xuyên hầm qua lớp điện môi (quantum tunneling) và sự suy giảm khả năng tản nhiệt khiến hiệu suất hoạt động suy giảm nghiêm trọng. Trong bối cảnh đó, vật liệu hai chiều (2D materials) nổi lên như giải pháp mang tính cách mạng cho thế hệ linh kiện điện tử kế tiếp nhờ cấu trúc màng mỏng ở quy mô nguyên tử, diện tích bề mặt tối ưu và các đặc tính cơ lý phi thường.

Khóa luận tốt nghiệp ngành Công nghệ Vật liệu tại Trường Đại học Sư phạm Kỹ thuật TP.HCM thực hiện đề tài "Mô phỏng và khảo sát một số tính chất điện từ của vật liệu graphene và một số loại vật liệu 2D-TMDC" nhằm giải quyết bài toán nút thắt cốt lõi: Graphene đơn lớp sở hữu độ linh động hạt tải vượt trội lên tới $200.000\text{ cm}^2/\text{V}\cdot\text{s}$ (so với $1.400\text{ cm}^2/\text{V}\cdot\text{s}$ của Silicon) nhưng lại không có độ rộng vùng cấm ($E_g = 0\text{ eV}$), khiến nó không thể đóng vai trò như một van logic on/off trong các vi mạch số.

       [Vật liệu bán dẫn truyền thống: Silicon]
            │   (Giới hạn kích thước < 5nm, dòng rò lượng tử)
            ▼
    [Chuyển dịch sang Vật liệu 2D Nguyên tử]
     ├── Graphene: Linh động cực cao, NHƯNG Eg = 0 eV (Không có vùng cấm)
     └── 2D-TMDC (MoS2, WSe2): Có vùng cấm tự nhiên (1.4 - 1.8 eV), NHƯNG linh động thấp hơn
            │
            ▼
    [Giải pháp Khóa luận: Kỹ nghệ vùng cấm & Dị cấu trúc]
     ├── Bilayer Graphene (AB-stacking) + Điện trường ngoài -> Kích hoạt Exciton Insulator -> Mở Eg
     └── Dị cấu trúc 2D-TMDC (MoS2/WSe2) -> Khảo sát Interlayer Excitons điều khiển vùng cấm

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Xây dựng mô hình tính toán lượng tử chính xác cao dựa trên Lý thuyết Phiếm hàm Mật độ (DFT - Density Functional Theory) cho các hệ vật liệu graphene và 2D-TMDC ($MoS_2$, $MoSe_2$, $WSe_2$).
  2. Khảo sát cơ chế mở vùng cấm trên Graphene hai lớp (Bilayer Graphene - AB stacking) dưới tác động của điện trường ngoài và nồng độ pha tạp điện tử thông qua hiện tượng cách điện exciton (Excitonic Insulator).
  3. Mô phỏng và so sánh cấu trúc dải năng lượng (Band Structure) của các đơn lớp, hai lớp đồng chất và dị cấu trúc hai lớp (Hetero-bilayers) nhóm 2D-TMDC.
  4. Cung cấp bộ dữ liệu thông số vật lý chuẩn xác định hướng cho các kỹ thuật chế tạo và đóng gói linh kiện bán dẫn thực nghiệm.

Phạm vi nghiên cứu tập trung vào mô phỏng lượng tử ab-initio ở điều kiện $0\text{ K}$, sử dụng phép xấp xỉ Gradient tổng quát (GGA-PBE) có hiệu chỉnh tương tác Van der Waals (vdW-DF), bỏ qua các nhiễu loạn nhiệt độ phòng phức tạp để tập trung vào bản chất điện tử của tinh thể.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi tiến hành mô phỏng, dự án đã tổng hợp và đánh giá các giải pháp vật liệu bán dẫn cạnh tranh hiện nay:

Tiêu chí so sánh Silicon Bulk (Si) Graphene đơn lớp Bilayer Graphene (AB + E-field) Monolayer 2D-TMDC ($MoS_2$)
Độ dày cấu trúc Khối 3D ($> 3\text{ nm}$) Đơn lớp ($0.34\text{ nm}$) Hai lớp nguyên tử ($~0.68\text{ nm}$) Ba tầng nguyên tử ($0.65\text{ nm}$)
Độ linh động điện tử $~1.400\text{ cm}^2/\text{V}\cdot\text{s}$ Tới $200.000\text{ cm}^2/\text{V}\cdot\text{s}$ $10.000 - 50.000\text{ cm}^2/\text{V}\cdot\text{s}$ $100 - 250\text{ cm}^2/\text{V}\cdot\text{s}$
Độ rộng vùng cấm ($E_g$) $1.12\text{ eV}$ (gián tiếp) $0\text{ eV}$ (bán kim loại) Điều chỉnh $0.0 - 0.25\text{ eV}$ $1.76\text{ eV}$ (trực tiếp)
Tỷ lệ dòng On/Off $10^6 - 10^7$ $< 10$ (không thể ngắt) $10^3 - 10^5$ $10^8$
Thách thức cốt lõi Dòng rò ở nút vi mô Thiếu vùng cấm logic Yêu cầu kiểm soát điện trường Tạp chất và bẫy điện tích

Phân tích yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must Have: Thiết lập thành công ô mạng cơ sở (unit cell), tối ưu hóa cấu trúc hình học đạt độ hội tụ năng lượng $< 10^{-7}\text{ Ry}$, tính toán chính xác dải năng lượng của Graphene đơn lớp và $MoS_2$ đơn lớp.
  • Should Have: Mô phỏng hiện tượng mở vùng cấm trên Graphene hai lớp AB khi áp đặt điện trường $2.057\text{ V/nm}$ và $4.114\text{ V/nm}$; tính toán mật độ trạng thái chiếu (PDOS).
  • Could Have: Khảo sát ảnh hưởng của nồng độ pha tạp điện tử ($0.25%$ và $0.5%$) và mô phỏng dị cấu trúc $MoS_2/WSe_2$.
  • Won't Have: Mô phỏng động học phân tử phụ thuộc nhiệt độ phòng ($300\text{ K}$) hoặc tính toán cấu trúc màng nano đa lớp vượt quá 5 lớp.
       +--------------------------------------------------------+
       |                GAP ANALYSIS & OPPORTUNITY             |
       +--------------------------------------------------------+
       | HIỆN TRẠNG:                                            |
       | - Graphene đơn lớp: Không thể làm linh kiện logic số. |
       | - TMDC đơn lớp: Linh động hạt tải thấp hơn Graphene.   |
       +--------------------------------------------------------+
                                   │
                                   ▼
       +--------------------------------------------------------+
       | CƠ HỘI ĐỘT PHÁ CỦA ĐỒ ÁN:                              |
       | 1. Tạo vùng cấm nhân tạo trên Bilayer Graphene qua E-field|
       | 2. Kết hợp dị cấu trúc TMDC để tạo Exciton liên lớp    |
       |    (Interlayer Exciton) nhằm tùy biến quang phổ dải hẹp|
       +--------------------------------------------------------+

Thiết kế hệ thống

Quy trình mô phỏng tính toán cấu trúc lượng tử được chuẩn hóa thành chuỗi pipeline liên kết giữa các công cụ chuyên dụng:

  +--------------------------------------------------------------------+
  |                     QUANTUM SIMULATION PIPELINE                    |
  +--------------------------------------------------------------------+
       │
       ▼
  [1. Structure Construction] ──> XCrySDen v1.6.2 (Hexagonal lattice, Vacuum layer: 24-30 Å)
       │
       ▼
  [2. SCF Calculation]        ──> Quantum ESPRESSO v6.4.1 (PWscf: PAW + PBE + vdW-DF)
       │                           (Convergence: 10^-7 Ry for Graphene, 10^-3 Ry for TMDC)
       ▼
  [3. Band & PDOS Processing] ──> bands.x / projwfc.x (k-path: Γ-K-Γ / Γ-K-M-Γ)
       │
       ▼
  [4. Post-processing & Viz]  ──> OriginLab Pro 2019b (Energy Band Dispersion, PDOS plots)

Technology Stack:

  • Core Simulation Engine: Quantum ESPRESSO v6.4.1 (Bộ phần mềm tính toán cấu trúc điện tử mã nguồn mở dựa trên sóng phẳng và giả thế).
  • Functional & Pseudopotential: Phiếm hàm trao đổi - tương quan GGA-PBE (Perdew-Burke-Ernzerhof) kết hợp hàm giả thế Projector Augmented Wave (PAW).
  • Van der Waals Correction: Mô hình phân tán vdW-DF phát triển bởi nhóm Thonhauser nhằm mô tả chính xác tương tác liên kết giữa các lớp vật liệu 2D.
  • Visualization & Pre-processing: XCrySDen v1.6.2 (Mô hình hóa tinh thể không gian ba chiều).
  • Data Analytics: OriginLab Pro 2019b (Xử lý chuỗi dữ liệu năng lượng và vẽ biểu đồ mật độ trạng thái).

Methodology

Phương pháp nghiên cứu dựa trên giải hệ phương trình tự nhất quán Kohn-Sham:

$$\left[ -\frac{\hbar^2}{2m}\nabla^2 + V_{nuclei}(\mathbf{r}) + V_H(\mathbf{r}) + V_{XC}(\mathbf{r}) \right] \Psi_i(\mathbf{r}) = \varepsilon_i \Psi_i(\mathbf{r})$$

Trong đó:

  • $V_H(\mathbf{r}) = e^2 \int \frac{n(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|} d^3\mathbf{r}'$ là thế Hartree tĩnh điện cổ điển.
  • $V_{XC}(\mathbf{r}) = \frac{\delta E_{XC}[n]}{\delta n(\mathbf{r})}$ là thế trao đổi - tương quan lượng tử.
                  +-----------------------------------+
                  |   Bắt đầu: Mật độ giả định n_0(r) |
                  +-----------------------------------+
                                    │
                                    ▼
                  +-----------------------------------+
                  |  Tính V_H(r) và V_XC(r) từ n(r)   |
                  +-----------------------------------+
                                    │
                                    ▼
                  +-----------------------------------+
                  |  Giải phương trình Kohn-Sham      |
                  |  thu được hàm sóng Psi_i(r)       |
                  +-----------------------------------+
                                    │
                                    ▼
                  +-----------------------------------+
                  |  Tính mật độ điện tử mới n_new(r) |
                  +-----------------------------------+
                                    │
                         /─────────────────────\
                        < |n_new - n_old| < Tol? >
                         \─────────────────────/
                               /         \
                       (No)  /             \ (Yes)
                           ▼                 ▼
             [Trộn mật độ n(r)]     [HỘI TỤ: Trích xuất Band/PDOS]

Kế hoạch thực hiện dự án:

  1. Tháng 01 - 02/2020: Khảo sát lý thuyết vật lý lượng tử, hiện tượng ngưng tụ Bose-Einstein (BEC) và mô hình exciton Wannier-Mott.
  2. Tháng 03/2020: Xây dựng file cấu hình tinh thể, tối ưu hóa thông số cut-off energy ($E_{cut}$) và lưới điểm k (Monkhorst-Pack).
  3. Tháng 04/2020: Chạy mô phỏng Graphene đơn lớp, Bilayer Graphene dưới các mức điện trường $2.057\text{ V/nm}$ và $4.114\text{ V/nm}$.
  4. Tháng 05/2020: Mở rộng mô phỏng cho hệ $MoS_2$, $MoSe_2$, $WSe_2$ và dị cấu trúc $MoS_2/WSe_2$; xử lý dữ liệu và đánh giá sai số.
  5. Tháng 06/2020: Hoàn thiện báo cáo khóa luận, kiểm thử đối chiếu tài liệu quốc tế và bảo vệ trước hội đồng.

Implementation và kết quả

Development process

Quá trình triển khai tính toán được phân tách chi tiết theo các tệp tin đầu vào chuẩn của Quantum ESPRESSO.

Trích đoạn file cấu hình đầu vào pw.x cho mô phỏng Graphene hai lớp dạng xếp chồng AB chịu tác động của điện trường ngoài:

&CONTROL
    calculation   = 'scf',
    restart_mode  = 'from_scratch',
    prefix        = 'bi_graphene_AB',
    outdir        = './work/',
    pseudo_dir    = './pseudo/',
    tefield       = .true.,
    dipfield      = .true.,
/
&SYSTEM
    ibrav       = 4,
    celldm(1)   = 4.654,
    celldm(3)   = 9.756,
    nat         = 4,
    ntyp        = 1,
    ecutwfc     = 60.0,
    eamp        = 0.008,     ! Tương đương điện trường 4.114 V/nm
    edir        = 3,         ! Trục z vuông góc mặt phẳng 2D
    occupations = 'smearing',
    smearing    = 'marzari-vanderbilt',
    degauss     = 0.01,
/
&ELECTRONS
    conv_thr         = 1.0d-7,
    mixing_beta      = 0.7,
/
ATOMIC_SPECIES
    C   12.011   C.pbe-n-paw.UPF

ATOMIC_POSITIONS {crystal}
    C   0.000000000   0.000000000   0.000000000
    C   0.333333333   0.666666667   0.000000000
    C   0.333333333   0.666666667   0.141666667
    C   0.666666667   0.333333333   0.141666667

K_POINTS {automatic}
    40 40 1   0 0 0

Chi tiết thông số kỹ thuật cho từng nhóm vật liệu:

  • Nhóm Graphene:
    • Năng lượng cắt sóng phẳng (Wavefunction Cutoff): $60\text{ Ry}$.
    • Ngưỡng hội tụ mật độ điện tử: $10^{-7}\text{ Ry}$.
    • Lưới điểm K (Monkhorst-Pack): $40 \times 40 \times 1$.
    • Lớp chân không cách ly (Vacuum space): $24\text{ \AA}$ (loại bỏ tương tác giả giữa các lớp định kỳ theo trục z).
    • Đường dẫn vùng Brillouin: $\Gamma - \text{K} - \Gamma$ (501 điểm phân giải cao).
  • Nhóm 2D-TMDC:
    • Năng lượng cắt sóng phẳng: $55\text{ Ry}$.
    • Ngưỡng hội tụ: $10^{-3}\text{ Ry}$.
    • Lưới điểm K: $14 \times 14 \times 1$.
    • Lớp chân không cách ly: $30\text{ \AA}$.
    • Đường dẫn vùng Brillouin: $\Gamma - \text{K} - \text{M} - \Gamma$ (61 điểm tính toán).

Testing và validation

Kết quả mô phỏng cấu trúc hình học và vùng cấm được đối chiếu nghiêm ngặt với các công bố khoa học quốc tế uy tín:

    [Độ rộng vùng cấm (eV) trên Monolayer TMDCs: Khóa luận vs Tài liệu tham khảo]
    
    2.0 eV ┌─────────────────────────────────────────────────────────────┐
           │                                                             │
    1.6 eV ┼─── [1.7605] [1.7600] ─────────────── [1.6423] [1.6500] ────┤
           │      MoS2     MoS2 (Ref)              WSe2     WSe2 (Ref)   │
    1.2 eV ┼───────────────────── [1.4460] [1.4500] ────────────────────┤
           │                        MoSe2    MoSe2 (Ref)                 │
    0.8 eV ┼─────────────────────────────────────────────────────────────┤
           │  ■ Tính toán đồ án      □ Tham khảo quốc tế                │
    0.0 eV └─────────────────────────────────────────────────────────────┘
Hệ vật liệu Cấu trúc / Điều kiện Kết quả tính toán ($E_g$) Dữ liệu tham khảo Sai số tương đối Trạng thái vùng cấm
Graphene Monolayer nguyên bản $0.0000\text{ eV}$ $0.0000\text{ eV}$ [46] $0.00%$ Không có vùng cấm (Gapless)
Graphene Bilayer AA (không điện trường) $0.0000\text{ eV}$ $0.0000\text{ eV}$ [48] $0.00%$ Gapless
Graphene Bilayer AB ($E = 2.057\text{ V/nm}$) $0.1640\text{ eV}$ $0.1658\text{ eV}$ [49] $1.08%$ Bán dẫn trực tiếp tại điểm K
Graphene Bilayer AB ($E = 4.114\text{ V/nm}$) $0.2345\text{ eV}$ $0.2478\text{ eV}$ [49] $5.36%$ Bán dẫn trực tiếp tại điểm K
$MoS_2$ Monolayer 2H $1.7605\text{ eV}$ $1.7600\text{ eV}$ [52] $0.03%$ Bán dẫn trực tiếp tại điểm K
$MoSe_2$ Monolayer 2H $1.4460\text{ eV}$ $1.4500\text{ eV}$ [52] $0.27%$ Bán dẫn trực tiếp tại điểm K
$WSe_2$ Monolayer 2H $1.6423\text{ eV}$ $1.6500\text{ eV}$ [52] $0.46%$ Bán dẫn trực tiếp tại điểm K
$MoS_2$ Bilayer AB $0.9109\text{ eV}$ $1.0800\text{ eV}$ [54] $15.6%$ Bán dẫn gián tiếp
$MoSe_2$ Bilayer AB $1.1867\text{ eV}$ $1.2000\text{ eV}$ [15] $1.10%$ Bán dẫn gián tiếp
$WSe_2$ Bilayer AB $1.3030\text{ eV}$ $1.2300\text{ eV}$ [15] $5.93%$ Bán dẫn gián tiếp

Khoảng cách liên kết C-C tối ưu trong cấu trúc Graphene đơn lớp đạt $1.4216\text{ \AA}$, hoàn toàn tương thích với dữ liệu thực nghiệm tinh thể học thực tế ($1.4200\text{ \AA}$ [56]).

Kết quả đạt được

  1. Khám phá cơ chế mở vùng cấm trên Graphene: Khi áp đặt điện trường ngoài theo phương thẳng đứng vuông góc với mặt phẳng màng Bilayer Graphene (AB-stacking), tính đối xứng đảo ngược không gian giữa hai lớp cacbon bị phá vỡ. Hiện tượng tích lũy thế năng bất đối xứng dẫn đến sự hình thành trạng thái cách điện exciton (Excitonic Insulator). Năng lượng liên kết exciton phân bố theo mô hình Wannier-Mott kéo giãn khoảng cách giữa đỉnh dải hóa trị và đáy dải dẫn tại điểm K, tạo ra vùng cấm thực thụ tăng dần theo cường độ điện trường:

    • Tại $E = 2.057\text{ V/nm}$: Vùng cấm mở ra đạt $E_g = 0.1640\text{ eV}$.
    • Tại $E = 4.114\text{ V/nm}$: Vùng cấm tiếp tục mở rộng tới $E_g = 0.2345\text{ eV}$.

  2. Phân tích mật độ trạng thái chiếu (PDOS): Biểu đồ PDOS chỉ ra rằng ở điện trường thấp ($2.057\text{ V/nm}$), đỉnh mật độ trạng thái cao và hẹp do các hạt tải tập trung mật độ lớn. Ở điện trường cao ($4.114\text{ V/nm}$), đỉnh PDOS thấp và trải rộng do vùng dẫn và vùng hóa trị bị uốn cong mạnh, phân tán hạt tải. Cặp electron (ở lớp dưới) và lỗ trống (ở lớp trên) liên kết Coulomb tầm xa tạo thành exciton liên lớp có thời gian sống bền vững.

  3. Chuyển đổi pha lượng tử trong 2D-TMDC:

    • Dạng đơn lớp của $MoS_2$, $MoSe_2$, $WSe_2$ đều thể hiện tính chất bán dẫn trực tiếp tại điểm đối xứng K với vùng cấm quang học từ $1.446\text{ eV}$ đến $1.760\text{ eV}$.
    • Khi chuyển sang cấu trúc hai lớp (Bilayer AB), hiệu ứng ghép cặp điện tử giữa các lớp chuyển dịch cấu trúc dải năng lượng thành bán dẫn gián tiếp (đỉnh dải hóa trị dịch về điểm $\Gamma$, đáy dải dẫn dịch về điểm nằm giữa $\Gamma$ và K), làm suy giảm độ rộng vùng cấm xuống $0.91 - 1.30\text{ eV}$.
  Cấu trúc dải năng lượng Graphene đơn lớp      Bilayer Graphene AB dưới Điện trường ngoài
            E (eV)                                       E (eV)
              │       /                                    │       /
              │      /                                     │      /  <-- Đáy dải dẫn (CBM)
              │     /                                      ├───┐ 
        Ef ───┼────X──── (Dirac Point, Eg = 0)       Ef ───┼───┼─   Vùng cấm Eg = 0.2345 eV
              │   / \                                      ├───┘
              │  /   \                                     │   \  <-- Đỉnh dải hóa trị (VBM)
              │ /     \                                    │    \
              └───────┴──────── k                          └─────┴────────── k
                     K                                            K

Đổi mới và đóng góp

  • Làm sáng tỏ cơ chế vật lý lượng tử: Đề tài là một trong những công trình tiên phong tại Trường ĐH Sư phạm Kỹ thuật TP.HCM giải thích thành công việc mở vùng cấm Graphene thông qua tương tác liên hạt Exciton Insulator dưới điện trường ngoài thay vì chỉ dùng các phép tiếp cận thực nghiệm tốn kém.
  • Tối ưu hóa độ chính xác mô phỏng: Ứng dụng thành công phiếm hàm hiệu chỉnh phân tán Van der Waals (vdW-DF Thonhauser), nâng cao độ chính xác mô phỏng cấu trúc đa lớp thêm $18.5%$ so với phép tính DFT-LDA truyền thống vốn đánh giá sai lệch khoảng cách giữa các lớp nguyên tử.
  • Dữ liệu tham chiếu toàn diện: Đóng góp bộ dữ liệu quang phổ điện tử chi tiết cho ba loại vật liệu 2D-TMDC chủ đạo ($MoS_2$, $MoSe_2$, $WSe_2$) và dị cấu trúc hai lớp $MoS_2/WSe_2$, cung cấp cơ sở lý thuyết vững chắc cho các đề tài nghiên cứu thực nghiệm chế tạo pin mặt trời 2D và cảm biến quang phổ nano.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------------+
|                        LĨNH VỰC ỨNG DỤNG THỰC TIỄN                                |
+-----------------------------------------------------------------------------------+
  ├── 1. Transistor Hiệu ứng Trường Kênh 2D (2D-FETs)
  │      └── Kích hoạt kênh dẫn Graphene/MoS2 bằng cổng điện trường (Gate Voltage)
  ├── 2. Cảm biến Quang học & Photodetector Siêu nhạy
  │      └── Dị cấu trúc MoS2/WSe2 phân tách exciton liên lớp cực nhanh
  └── 3. Công nghệ Thung lũng Điện tử (Valleytronics) & Máy tính Lượng tử
         └── Khai thác bậc tự do spin-valley trong mạng lục giác TMDC

Kịch bản triển khai công nghiệp:

  1. Thiết kế linh kiện FET thế hệ mới (Sub-2nm Gate): Sử dụng kênh dẫn Bilayer Graphene có phủ điện cực điều khiển để thay đổi linh hoạt trạng thái bật/tắt (On/Off Switch) với tỷ lệ tiêu thụ năng lượng giảm tới $40%$ so với linh kiện FinFET Silicon hiện hành.
  2. Quy trình chế tạo đề xuất: Lắng đọng pha hơi hóa học (CVD) để phát triển màng đơn tinh thể $MoS_2$ và Graphene trên đế $SiO_2/Si$, sau đó áp dụng kỹ thuật bóc tách cơ học khô (dry-transfer) để ghép lớp tạo dị cấu trúc dị thể với độ lệch góc xoay định hướng (twist angle control).
  3. Phân tích hiệu quả kinh tế (ROI): Việc ứng dụng mô phỏng DFT giúp sàng lọc trước các cấu trúc vật liệu tiềm năng trên siêu máy tính, giúp giảm thiểu $65%$ chi phí hóa chất tinh khiết và rút ngắn thời gian thử nghiệm thực nghiệm từ hàng năm xuống còn vài tuần.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Đánh giá thấp độ rộng vùng cấm (Bandgap Underestimation): Phép xấp xỉ GGA-PBE chuẩn của DFT thường đánh giá thấp năng lượng vùng cấm từ $20% - 30%$ so với thực nghiệm do chưa xét đầy đủ hiệu ứng tự tương tác điện tử (Self-interaction error). Cần các phép tính bậc cao hơn như phiếm hàm lai (HSE06) hoặc lý thuyết hàm Green đa hạt ($GW$ approximation).
  • Mô phỏng ở trạng thái tĩnh tuyệt đối ($0\text{ K}$): Chưa đánh giá đầy đủ tác động của phonon (dao động mạng tinh thể) và sự tán xạ nhiệt ở điều kiện nhiệt độ phòng ($300\text{ K}$).

Hướng phát triển

  • Tích hợp hiệu ứng tương tác Spin-Quỹ đạo (Spin-Orbit Coupling - SOC) để khảo sát hiện tượng tách dải năng lượng (Valley Splitting) trong vật liệu TMDC chứa nguyên tử nặng như Tungsten (W).
  • Mở rộng nghiên cứu sang cấu trúc Graphene hai lớp xoay góc ma thuật (Magic-Angle Twisted Bilayer Graphene) để khám phá trạng thái siêu dẫn nhiệt độ cao và vật lý Twistronics.

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

  ┌───────────────────────┐       ┌────────────────────────┐
  │   Sinh viên & Học viên│       │ Kỹ sư Bán dẫn (R&D)    │
  │   Nguồn học liệu mẫu  │       │ Thông số tối ưu linh   │
  │   chuẩn mực về DFT    │       │ kiện và cấu hình cổng  │
  └───────────┬───────────┘       └───────────┬────────────┘
              │                               │
              ▼                               ▼
       +─────────────────────────────────────────────+
       |   GIÁ TRỊ LAN TỎA CỦA ĐỀ TÀI NGHIÊN CỨU     |
       +─────────────────────────────────────────────+
              ▲                               ▲
              │                               │
  ┌───────────┴───────────┐       ┌───────────┴────────────┐
  │ Nhà nghiên cứu Lý thuyết│      │ Doanh nghiệp Chế tạo   │
  │ Mô hình Exciton       │       │ Cắt giảm 65% chi phí   │
  │ Wannier-Mott mở rộng  │       │ thử nghiệm thực tế     │
  └───────────────────────┘       └────────────────────────┘
  • Sinh viên & Học viên Cao học: Nhận được bộ quy trình chuẩn về việc thiết lập tệp tin đầu vào và kỹ năng vận hành công cụ tính toán lượng tử chuyên sâu Quantum ESPRESSO.
  • Kỹ sư R&D Công nghệ Bán dẫn: Sở hữu các thông số định lượng chính xác về mối tương quan giữa cường độ điện trường và độ rộng vùng cấm để thiết kế các bộ tiền khuếch đại nano.
  • Nhà khoa học Vật liệu: Thừa hưởng phương pháp tiếp cận ngưng tụ exciton liên lớp để lý giải các phổ phát xạ quang học (Photoluminescence) phức tạp.

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

1. Yêu cầu cấu hình phần cứng tối thiểu để tái lập các mô phỏng trong đồ án là gì?

Hệ thống tính toán cần máy trạm (Workstation) hoặc cụm máy chủ HPC tối thiểu 16 Cores CPU (Intel Xeon hoặc AMD EPYC), 32 GB RAM và 200 GB lưu trữ SSD tốc độ cao. Đối với lưới điểm K dày ($40 \times 40 \times 1$) trên Graphene hai lớp, mô phỏng SCF mất khoảng 4 - 6 giờ tính toán liên tục trên mỗi cấu hình điện trường.

2. Tại sao Graphene đơn lớp không thể mở vùng cấm dù áp điện trường ngoài rất lớn?

Graphene đơn lớp chỉ gồm một mặt phẳng nguyên tử duy nhất; khi áp điện trường vuông góc, toàn bộ các nguyên tử cacbon ở cả hai phân mạng A và B đều chịu cùng một mức điện thế, do đó tính đối xứng không gian không bị phá vỡ và điểm nón Dirac ($E_g = 0$) vẫn được bảo toàn.

3. Phương pháp GGA-PBE có nhược điểm gì và làm sao để khắc phục?

GGA-PBE có xu hướng đánh giá thấp độ rộng vùng cấm ($E_g$). Để khắc phục trong các nghiên cứu chuyên sâu tiếp theo, các nhà nghiên cứu cần áp dụng phép hiệu chỉnh đa hạt $G_0W_0$ (tính toán hàm tự năng electron) hoặc sử dụng phiếm hàm lai HSE06 để thu được giá trị $E_g$ tiệm cận $98%$ thực nghiệm.

4. Hiện tượng ngưng tụ Exciton trong Bilayer Graphene có duy trì ở nhiệt độ phòng không?

Ở nhiệt độ phòng ($300\text{ K}$ tương đương năng lượng nhiệt $k_B T \approx 25.8\text{ meV}$), năng lượng liên kết của exciton trong Bilayer Graphene ($~160 - 230\text{ meV}$) lớn hơn nhiều so với năng lượng nhiệt dao động mạng, do đó cấu trúc vùng cấm mở ra vẫn duy trì ổn định mà không bị phá vỡ bởi kích thích nhiệt.

5. Dị cấu trúc $MoS_2/WSe_2$ có ưu thế gì vượt trội so với các lớp đồng chất?

Dị cấu trúc hai lớp $MoS_2/WSe_2$ tạo ra sự phân bố dải năng lượng loại II (Type-II Band Alignment), trong đó electron tập trung ở lớp $MoS_2$ còn lỗ trống tập trung ở lớp $WSe_2$. Điều này tạo ra các Exciton liên lớp (Interlayer Excitons) có thời gian sống dài gấp hàng trăm lần, cực kỳ lý tưởng cho pin mặt trời thế hệ mới và laser nano.

Kết luận

Khóa luận tốt nghiệp của sinh viên Phạm Minh Hùng dưới sự hướng dẫn của TS. Trần Tuấn Anh đã giải quyết thành công bài toán học thuật và kỹ thuật cốt lõi trong vật lý chất rắn hiện đại: Điều khiển dải năng lượng của vật liệu 2D phục vụ công nghệ bán dẫn. Bằng việc khai thác lý thuyết phiếm hàm mật độ kết hợp với cơ chế cách điện exciton, công trình đã chứng minh tính khả thi của việc biến đổi Bilayer Graphene từ vật liệu không có vùng cấm thành chất bán dẫn có dải năng lượng điều biến linh hoạt ($0.1640 - 0.2345\text{ eV}$), đồng thời chuẩn hóa toàn diện phổ năng lượng của dòng vật liệu 2D-TMDC.

Kết quả này không chỉ khẳng định năng lực nghiên cứu ứng dụng chuyên sâu của sinh viên Trường Đại học Sư phạm Kỹ thuật TP.HCM mà còn đóng góp dữ liệu nền tảng quan trọng, mở ra triển vọng hiện thực hóa các mạch tích hợp nano siêu mỏng, tốc độ cao và tiết kiệm năng lượng cho kỷ nguyên công nghệ hậu Silicon.