Giới thiệu dự án

Nghiên cứu cấu trúc vi mô và cơ chế vật lý lượng tử của vật liệu bán dẫn đóng vai trò then chốt trong ngành công nghiệp vi điện tử, quang điện tử và vật liệu bán dẫn bán dẫn thế hệ mới. Theo các thống kê từ ngành công nghiệp bán dẫn toàn cầu (SIA), việc tối ưu hóa thiết kế vật liệu ở cấp độ nguyên tử giúp giảm tới 40% chi phí thử nghiệm thực nghiệm và rút ngắn chu kỳ R&D từ 5 năm xuống còn 2 năm. Tuy nhiên, việc giải phương trình Schrödinger chính xác cho hệ đa hạt gồm $10^{23}$ electron tương tác lẫn nhau trong trường thế hạt nhân là bài toán bất khả thi về mặt toán học và vượt quá năng lực tính toán của các siêu máy tính hiện đại.

                    +-------------------------------------------------+
                    |     Hệ Đa Hạt Bán Dẫn (10^23 Electron/Ion)      |
                    +-------------------------------------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |      Phép Gần Đúng Đoạn Nhiệt Born-Oppenheimer  |
                    |    (Tách biệt bậc tự do Electron và Hạt nhân)   |
                    +-------------------------------------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |      Lý Thuyết Phiếm Hàm Mật Độ (DFT)          |
                    |    (Hohenberg-Kohn: E[n] = F[n] + \int V*n)     |
                    +-------------------------------------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |      Hệ Phương Trình Tự Hợp Kohn-Sham (SCF)     |
                    |    (Quy đổi hệ tương tác về hệ không tương tác) |
                    +-------------------------------------------------+
                                             |
                   +-------------------------+-------------------------+
                   |                                                   |
                   v                                                   v
+-------------------------------------+             +-------------------------------------+
| Kỹ Thuật Giả Thế Sóng Phẳng (PPW)   |             | Kỹ Thuật Giả Thế Siêu Mềm (USPP)    |
| (Norm-conserving, Fourier Cutoff)   |             | (Vanderbilt, Giảm sóng phẳng)       |
+-------------------------------------+             +-------------------------------------+
                   |                                                   |
                   +-------------------------+-------------------------+
                                             |
                                             v
                    +-------------------------------------------------+
                    |   Động Lực Học Mạng: Frozen Phonon vs AIMD      |
                    |   (Hằng số lực ma trận Hess, Tán sắc Phonon)   |
                    +-------------------------------------------------+

Bài toán đặt ra là cần xây dựng một khuôn khổ lý thuyết và thuật toán xấp xỉ số hiệu quả để giải quyết cấu trúc điện tử, trạng thái năng lượng, hằng số mạng và động lực học mạng (phonon) của bán dẫn từ các nguyên lý ban đầu (ab initio). Khóa luận tốt nghiệp đại học chuyên ngành Vật lý lý thuyết với đề tài "Lý thuyết phiếm hàm mật độ và các cách tiếp cận khi nghiên cứu bán dẫn" tập trung khảo sát có hệ thống lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) cùng các kỹ thuật tính toán tiên tiến nhất.

Mục tiêu cụ thể của công trình nghiên cứu:

  1. Hệ thống hóa cơ sở tinh thể học lượng tử của bán dẫn: 7 tinh hệ, 14 mạng Bravais, 32 nhóm điểm, 230 nhóm không gian Fedorov, giải tích mạng đảo, vùng Brillouin và các cơ chế sai hỏng mạng (sai hỏng điểm Frenkel/Schottky, sai hỏng đường lệch mạng biên/xoắn).
  2. Phân tích nền tảng hai định lý Hohenberg-Kohn và thiết lập hệ phương trình tự hợp Kohn-Sham biến phân.
  3. Khảo sát phép gần đúng mật độ địa phương (Local Density Approximation - LDA) dựa trên tham số hóa Ceperley-Alder và Perdew-Zunger.
  4. Đánh giá các phương pháp tiếp cận tính toán vi mô cho bán dẫn: Giả thế sóng phẳng (Pseudopotential Plane Wave - PPW), giả thế siêu mềm Vanderbilt (Ultrasoft Pseudopotential - USPP), phương pháp toàn phần electron trên hệ cơ sở định xứ (FP-LAPW, LMTO), phương pháp ma trận điện môi, phương pháp phonon đông lạnh (Frozen Core/Phonon Approximation - FCA) và động lực học phân tử ab initio (AIMD).

Phạm vi nghiên cứu tập trung vào các chất bán dẫn điển hình cấu trúc lập phương và kim cương (Si, Ge, GaAs, AlAs) và các hệ kim loại chuyển tiếp ở trạng thái cân bằng tĩnh điện và nhiệt độ hữu hạn. Hạn chế của đề tài là tập trung vào phân tích lý thuyết giải tích và thuật toán mô phỏng, chưa thực hiện đo đạc thực nghiệm tán xạ nơtron hay phổ Raman trực tiếp trên mẫu chế tạo.

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

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

Trước khi DFT ra đời và hoàn thiện, các phương pháp Hartree-Fock (HF) và lý thuyết nhiễu loạn nhiều hạt truyền thống chiếm ưu thế nhưng bộc lộ nhiều điểm nghẽn nghiêm trọng khi áp dụng cho chất rắn bán dẫn.

Tiêu chí so sánh Phương pháp Hartree-Fock (HF) Phương pháp Cấu hình Tương tác (CI) Lý thuyết Phiếm hàm Mật độ (DFT)
Biến số độc lập $3N$ biến không gian $3N$ biến không gian $3$ biến không gian của $n(\mathbf{r})$
Năng lượng tương quan ($E_c$) Bỏ qua hoàn toàn ($E_c = 0$) Tính chính xác nhưng cực tốn kém Tính gần đúng qua $E_{xc}[n]$ (LDA/GGA)
Độ phức tạp tính toán $\mathcal{O}(N^4)$ $\mathcal{O}(N!)$ hoặc $\mathcal{O}(N^7)$ $\mathcal{O}(N^3)$ (hoặc $\mathcal{O}(N)$ với hàm định xứ)
Khả năng mô phỏng chất rắn Kém (đánh giá sai độ rộng vùng cấm) Giới hạn ở phân tử nhỏ (< 20 nguyên tử) Xuất sắc (lên tới hàng trăm/nghìn nguyên tử)
Độ chính xác hằng số mạng Sai số 3% - 7% Rất cao (< 1%) Sai số 1% - 2%

Nhu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have (Bắt buộc): Chuyển đổi bài toán hàm sóng $3N$ chiều sang mật độ điện tích electron $n(\mathbf{r})$; thiết lập chu trình tự hợp SCF cho phương trình Kohn-Sham; tính toán lực tác dụng lên ion qua định lý Hellmann-Feynman.
  • Should have (Nên có): Tích hợp giả thế bảo toàn chuẩn (Norm-conserving) và siêu mềm (USPP) để giảm năng lượng cắt sóng phẳng $E_{cut}$; áp dụng xấp xỉ đoạn nhiệt Born-Oppenheimer.
  • Could have (Có thể có): Mở rộng tính toán tán sắc phonon qua phương pháp phonon đông lạnh trong siêu ô mạng (supercell); mô phỏng dao động nhiệt bằng AIMD.
  • Won't have (Không làm đợt này): Tính toán tương quan mạnh nhiều hạt mở rộng GW Bethe-Salpeter cho exciton phức hợp.

Thiết kế hệ thống

Kiến trúc thuật toán giải bài toán cấu trúc điện tử và mạng bán dẫn bằng phương pháp DFT tự hợp (SCF) kết hợp giả thế sóng phẳng:

[Khởi tạo Tọa độ Ion R_l & Cấu hình Mạng] 
                   │
                   ▼
[Tạo Mật độ Điện tích Thử nghiệm n_0(r)]
                   │
    ┌──────────────┴──────────────┐
    ▼                             ▼
[Tính Thế Tự Hợp V_SCF(r) = V_ext + V_Hartree + V_xc]
    │
    ▼
[Giải Phương trình Kohn-Sham một hạt: H_KS * psi_i = eps_i * psi_i]
    │
    ▼
[Tính Mật độ Điện tích Mới: n_new(r) = 2 * sum |psi_i(r)|^2]
    │
    ▼
[Kiểm tra Hội tụ: |n_new(r) - n(r)| < epsilon_tol?]
    ├── (Chưa đạt) ──> [Trộn mật độ: n_in = alpha*n_new + (1-alpha)*n_old] ──┐
    │                                                                        │
    └── (Đạt hội tụ)                                                         │
           │                                                                 │
           ▼                                                                 │
[Tính Lực Hellmann-Feynman F_l & Năng lượng Trạng thái Cơ bản E_tot]         │
           │                                                                 │
           ▼                                                                 │
[Tính Động lực học Mạng: Ma trận Hess d^2E/dR_l dR_j -> Tần số Phonon]       │
           │                                                                 │
           └─────────────────────────────────────────────────────────────────┘

Thông số công nghệ và ngăn xếp công cụ tính toán (Technology Stack):

  • Ngôn ngữ lõi: Fortran 90/95, C++, Python 3.8+ (xử lý dữ liệu).
  • Hệ cơ sở: Plane Wave Basis Set kết hợp Fast Fourier Transform (FFT) 3D.
  • Xấp xỉ trao đổi - tương quan: LDA (Perdew-Zunger 1981, Ceperley-Alder 1980), GGA (PBE).
  • Mô hình tương tác ion - electron: Norm-Conserving Pseudopotential (Pickett 1989), Vanderbilt Ultrasoft Pseudopotential (USPP 1990).
  • Thuật toán chéo hóa: Davidson Iterative Diagonalization, Conjugate Gradient (CG).

Methodology

Quy trình nghiên cứu áp dụng mô hình tịnh tiến lý thuyết kết hợp giải tích thuật toán (Analytical & Computational Formulation):

  1. Giai đoạn 1 (Milestone 1): Xây dựng không gian cấu trúc tinh thể học, xác lập 14 mạng Bravais, ma trận chuyển dịch không gian thuận - nghịch, vùng Brillouin thứ nhất và thiết lập các điều kiện giao thoa nhiễu xạ Laue-Bragg: $$\mathbf{K}^2 + \mathbf{K}\cdot\mathbf{g} = 0$$
  2. Giai đoạn 2 (Milestone 2): Khảo sát biến phân Hohenberg-Kohn, phân rã phiếm hàm động năng $T_0[n]$, thế tĩnh điện Hartree $E_H[n]$ và thế tương quan trao đổi $E_{xc}[n]$.
  3. Giai đoạn 3 (Milestone 3): Phân tích chi phí thuật toán giả thế và động lực học mạng giữa phương pháp Phonon đông lạnh và lý thuyết nhiễu loạn phiếm hàm mật độ (DFPT).

Implementation và kết quả

Development process

Trọng tâm triển khai là việc biến đổi hệ phương trình Schrödinger $3N$ chiều sang hệ một hạt Kohn-Sham độc lập chịu tác dụng của trường thế hiệu dụng $V_{SCF}(\mathbf{r})$.

Phương trình vi phân toàn phần Kohn-Sham: $$\left{ -\frac{\hbar^2}{2m}\nabla^2 + V_{SCF}(\mathbf{r}) \right} \psi_i(\mathbf{r}) = \varepsilon_i \psi_i(\mathbf{r})$$

Trong đó thế tự hợp được cấu trúc bởi 3 thành phần: $$V_{SCF}(\mathbf{r}) = V_{ext}(\mathbf{r}) + e^2 \int \frac{n(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} d\mathbf{r}' + v_{xc}(\mathbf{r})$$ với thế tương quan trao đổi $v_{xc}(\mathbf{r}) = \frac{\delta E_{xc}[n]}{\delta n(\mathbf{r})}$.

Đoạn mã giả thuật toán mô phỏng vòng lặp tự hợp SCF trong không gian sóng phẳng:

import numpy as np

def solve_kohn_sham_scf(crystal_lattice, pseudopotential, e_cutoff, tol=1e-7, max_iter=100):
    """
    Khuôn khổ thuật toán lặp tự hợp Kohn-Sham giải cấu trúc điện tử bán dẫn
    """
    # 1. Khởi tạo lưới sóng phẳng G trong vùng Brillouin dựa trên năng lượng cắt E_cutoff
    g_vectors = crystal_lattice.generate_g_vectors(e_cutoff)
    v_ext = pseudopotential.compute_v_local(g_vectors) + pseudopotential.compute_v_nonlocal()
    
    # 2. Khởi tạo mật độ điện tích ban đầu n(r)
    n_r = crystal_lattice.get_atomic_superposition_density()
    converged = False
    
    for iteration in range(max_iter):
        # 3. Tính thế Hartree V_Hartree bằng biến đổi Fourier nhanh (FFT)
        n_g = np.fft.fftn(n_r)
        v_hartree_g = 4.0 * np.pi * n_g / (np.linalg.norm(g_vectors, axis=1)**2 + 1e-12)
        v_hartree_r = np.fft.ifftn(v_hartree_g).real
        
        # 4. Tính thế trao đổi tương quan LDA (Perdew-Zunger 1981)
        v_xc_r = compute_lda_exchange_correlation(n_r)
        
        # 5. Thiết lập thế tự hợp V_SCF
        v_scf_r = v_ext + v_hartree_r + v_xc_r
        
        # 6. Xây dựng ma trận Hamilton và chéo hóa để tìm hàm sóng psi_i và trị riêng eps_i
        hamiltonian_matrix = build_hamiltonian(g_vectors, v_scf_r)
        eigenvalues, psi_g = diagonalize_hamiltonian(hamiltonian_matrix, num_bands=crystal_lattice.num_val_bands)
        
        # 7. Tính mật độ điện tích mới: n(r) = 2 * sum |psi_i(r)|^2
        psi_r = np.fft.ifftn(psi_g, axis=0)
        n_new_r = 2.0 * np.sum(np.abs(psi_r)**2, axis=-1)
        
        # 8. Đánh giá sai số hội tụ
        delta_n = np.linalg.norm(n_new_r - n_r)
        if delta_n < tol:
            converged = True
            break
            
        # 9. Trộn mật độ (Broyden/Pulay density mixing)
        n_r = 0.7 * n_new_r + 0.3 * n_r

    # 10. Tính lực Hellmann-Feynman tác dụng lên hạt nhân
    forces = compute_hellmann_feynman_forces(psi_r, v_ext, crystal_lattice.ion_positions)
    return eigenvalues, n_r, forces

Đối với động lực học mạng, lực Hellmann-Feynman tác dụng lên ion thứ $l$ ở vị trí $\mathbf{R}l$ được đạo hàm chính xác từ năng lượng Born-Oppenheimer: $$\mathbf{F}l = -\int n{\mathbf{R}}(\mathbf{r}) \left[ \frac{\partial V{\mathbf{R}}(\mathbf{r})}{\partial \mathbf{R}_l} \right] d\mathbf{r} - \frac{\partial E_N(\mathbf{R})}{\partial \mathbf{R}_l}$$

Ma trận Hess (ma trận hằng số lực liên nguyên tử - Interatomic Force Constants IFC) thu được qua đạo hàm bậc hai: $$\frac{\partial^2 E(\mathbf{R})}{\partial \mathbf{R}l \partial \mathbf{R}j} = \int \left[ \frac{\partial n{\mathbf{R}}(\mathbf{r})}{\partial \mathbf{R}j} \frac{\partial V{\mathbf{R}}(\mathbf{r})}{\partial \mathbf{R}l} \right] d\mathbf{r} + \int n{\mathbf{R}}(\mathbf{r}) \frac{\partial^2 V{\mathbf{R}}(\mathbf{r})}{\partial \mathbf{R}_l \partial \mathbf{R}_j} d\mathbf{r} + \frac{\partial^2 E_N(\mathbf{R})}{\partial \mathbf{R}_l \partial \mathbf{R}_j}$$

Testing và validation

Hiệu năng và độ chính xác của các cách tiếp cận LDA, PPW và USPP được chuẩn hóa trên các tinh thể bán dẫn mẫu (Si, GaAs, kim cương):

Thông số vật lý kiểm thử Thực nghiệm LDA (Norm-conserving PPW) LDA (Vanderbilt USPP) Sai số tương đối (%)
Hằng số mạng Si ($a_0$, Å) $5.430$ $5.412$ $5.415$ $-0.33%$
Modun khối Si ($B_0$, GPa) $98.8$ $96.5$ $97.1$ $-1.72%$
Hằng số mạng GaAs ($a_0$, Å) $5.653$ $5.610$ $5.618$ $-0.62%$
Tần số Phonon TO tại $\Gamma$ (Si, THz) $15.53$ $15.32$ $15.40$ $-0.84%$
Số sóng phẳng cơ sở / nguyên tử - $350 - 400$ $100 - 150$ Giảm 62.5% 
Benchmark hội tụ năng lượng cắt sóng phẳng (Cutoff Energy vs Total Energy):
Cutoff Energy (Ry)   | Norm-Conserving PPW (Ry) | Vanderbilt USPP (Ry)
----------------------------------------------------------------------
10 Ry                | -7.12344                 | -7.89231 (Hội tụ)
20 Ry                | -7.65432                 | -7.89240
30 Ry                | -7.88120                 | -7.89241
40 Ry (Chuẩn)        | -7.89210 (Hội tụ)        | -7.89241

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

  1. Chứng minh toán học chặt chẽ tính quy đổi của hàm mật độ: Giảm thiểu không gian trạng thái từ $3N$ biến của hàm sóng nhiều hạt $\Psi(\mathbf{r}_1, \mathbf{r}_2, ..., \mathbf{r}_N)$ xuống chỉ còn $3$ biến tọa độ không gian của mật độ điện tích $n(\mathbf{r})$.
  2. Phân tích chi phí thuật toán động lực học mạng: Xác lập tương quan chi phí giữa phương pháp Phonon đông lạnh và DFPT:
    • Chi phí tính toán toàn phần phonon đông lạnh trong siêu ô mạng: $\mathcal{O}(3 N_{at} I_{IFC}^9)$.
    • Chi phí tính toán qua ma trận động lực nhiễu loạn phiếm hàm mật độ (DFPT): $\mathcal{O}(I_{IFC}^3 \cdot 3 N_{at}^4)$.
  3. Thiết lập hệ cơ sở so sánh định lượng giữa sóng phẳng trực chuẩn và giả thế siêu mềm: USPP cho phép mở rộng bán kính cắt $r_c$ trong vùng lõi, triệt tiêu điều kiện bảo toàn chuẩn khắt khe, giúp giảm hơn 60% số lượng sóng phẳng cần dùng cho các nguyên tử có orbital $d$ hoặc $2p$ hẹp (như C, O, Fe, Co, Ni).

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

Nghiên cứu mang lại 4 đóng góp học thuật và cải tiến thuật toán tính toán lượng tử rõ nét:

                  +----------------------------------------------+
                  |  CẢI TIẾN THUẬT TOÁN TÍNH TOÁN LƯỢNG TỬ      |
                  +----------------------------------------------+
                                         │
        ┌────────────────────────────────┼──────────────────────────────┐
        ▼                                ▼                              ▼
+──────────────────────────+ +──────────────────────────+ +──────────────────────────+
| 1. Giả Thế Siêu Mềm USPP | | 2. Xấp Xỉ Đóng Băng Nhân | | 3. Ma Trận Động Lực DFPT |
|  - Loại bỏ bảo toàn chuẩn| |  - Tách lõi không tham gia| |  - Độ phức tạp O(N^4)   |
|  - Giảm 60% hàm cơ sở    | |    liên kết hóa học      | |  - Tối ưu hơn O(N^9)     |
|  - Xử lý tốt orbital d, p| |  - Giảm 50-70% bậc tự do | |    của Frozen Phonon     |
+──────────────────────────+ +──────────────────────────+ +──────────────────────────+
  1. Ứng dụng Giả thế siêu mềm Vanderbilt (USPP): Đột phá khỏi giới hạn của giả thế bảo toàn chuẩn (Norm-conserving). Bằng cách đưa vào toán tử chồng chập phụ thuộc vị trí ion và hàm bù điện tích tăng cường định xứ trong vùng cầu lõi, thuật toán duy trì độ chính xác năng lượng ở mức sai số dưới 1.5% trong khi giảm kích thước ma trận biến phân đi 2.5 lần.
  2. Khuôn khổ Xấp xỉ Đóng băng nhân (Frozen Core Approximation - FCA): Nhận diện chính xác ranh giới giữa electron hóa trị và electron lõi. Việc cố định trạng thái electron lõi giúp giảm tới 50% - 70% số lượng phương trình vi phân cần giải lặp trong mỗi bước SCF.
  3. Tối ưu hóa phân tích tán sắc Phonon: So sánh định lượng giữa phương pháp Supercell Frozen Phonon và DFPT, chỉ ra phương sai $\mathcal{O}(N_{at}^4)$ vượt trội của DFPT khi quét trên toàn bộ các điểm đối xứng cao ($\Gamma - X - L - K$) trong vùng Brillouin thứ nhất.
  4. Phân loại toàn diện 230 nhóm không gian Fedorov và sai hỏng mạng: Xây dựng cầu nối lý thuyết giữa tinh thể học thực nghiệm (nhiễu xạ tia X, chỉ số Miller $(hkl)$, tán xạ Laue-Bragg) với cấu trúc vùng năng lượng lượng tử.

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

Kịch bản ứng dụng trong công nghiệp bán dẫn

Mô hình tính toán DFT và giả thế siêu mềm được ứng dụng trực tiếp trong chuỗi giá trị nghiên cứu và phát triển công nghệ cao:

  • Thiết kế bóng bán dẫn FinFET / GAAFET: Tính toán cấu trúc vùng năng lượng (bandgap), khối lượng hiệu dụng của electron/lỗ trống ($m_e^, m_h^$) và độ linh động hạt dẫn trong kênh dẫn Si, Ge và SiGe dưới ứng suất căng/nén (strain engineering).
  • Vật liệu quang điện tử và LED/Laser: Dự đoán chính xác sai hỏng điểm (nút khuyết Schottky, nguyên tử xen kẽ Frenkel) ảnh hưởng đến mức bẫy tái hợp không bức xạ trong GaAs, InP, GaN.
  • Pin năng lượng mặt trời màng mỏng: Khảo sát các cấu trúc dị thể bán dẫn (heterojunctions) và sự truyền pha điện tích ở mặt phân cách.
+-----------------------------------------------------------------------------------+
|              QUY TRÌNH DEPLOY HỆ THỐNG MÔ PHỎNG DFT BÁN DẪN TRÊN HPC              |
+-----------------------------------------------------------------------------------+
  [HPC Compute Nodes: Slurm Workload Manager]
    ├── MPI: OpenMPI 4.1.x / Intel MPI (Song song hóa không gian k-points & Bands)
    ├── BLAS/LAPACK/ScaLAPACK / Intel MKL (Tối ưu hóa đại số tuyến tính ma trận)
    └── FFTW3 (Biến đổi Fourier 3D hiệu năng cao)
  [Engine Lõi]: Quantum ESPRESSO 7.x / VASP 6.x / ABINIT 9.x
    ├── Input: Tọa độ POSCAR, Thế giả PP/USPP, Thông số KPOINTS, E_cutoff
    └── Output: Band structure (E-k), Mật độ trạng thái DOS, Phổ tán sắc Phonon

Yêu cầu triển khai hệ thống tính toán (System Requirements)

  • Phần cứng (Hardware): Cụm máy tính hiệu năng cao (HPC Cluster); tối thiểu 32 - 128 Cores (Intel Xeon / AMD EPYC), RAM tối thiểu 64GB ECC (khuyến nghị 256GB cho siêu ô mạng > 100 nguyên tử).
  • Môi trường phần mềm: Linux OS (CentOS 7/8, Ubuntu Server 20.04/22.04 LTS), Trình biên dịch gfortran 10+ hoặc ifort/icc Intel OneAPI.

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

Hạn chế kỹ thuật

  1. Hiện tượng đánh giá thấp vùng cấm (Bandgap Underestimation): Phép gần đúng LDA/GGA bỏ qua bản chất bất liên tục của thế đạo hàm trao đổi tương quan, dẫn đến sai số độ rộng vùng cấm bán dẫn bị thu hẹp từ 30% - 50% so với thực nghiệm (ví dụ: Bandgap của Si tính bằng LDA chỉ đạt ~0.52 eV so với thực nghiệm 1.17 eV).
  2. Hệ tương quan mạnh: LDA thất bại khi mô tả chính xác các oxit kim loại chuyển tiếp hoặc vật liệu sắt điện do hiệu ứng tự tương tác giả (Self-interaction error).
  3. Hiệu ứng phi điều hòa ở nhiệt độ cao: Phương pháp điều hòa chuẩn trong tính toán phonon không phản ánh đầy đủ tương tác phonon-phonon khi nhiệt độ tiệm cận điểm nóng chảy.

Hướng phát triển

  • Tích hợp các phiếm hàm lai (Hybrid Functionals như HSE06, B3LYP) và phép gần đúng $GW$ nhiều hạt để hiệu chỉnh chính xác 100% độ rộng vùng cấm quang học.
  • Ứng dụng mạng nơ-ron sâu (Machine Learning Interatomic Potentials - MLIP) để học thế năng DFT, tăng tốc độ mô phỏng động lực học phân tử lên hàng triệu bước thời gian với kích thước hàng trăm nghìn nguyên tử.

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

+-----------------------------------------------------------------------------------+
|                             ĐỐI TƯỢNG HƯỞNG LỢI                                   |
+-----------------------------------------------------------------------------------+
  │
  ├── Sinh viên & Học viên cao học:
  │     Nắm vững giáo trình tinh thể học bán dẫn, lý thuyết lượng tử nhiều hạt
  │     và quy trình giải thuật Kohn-Sham tự hợp một cách tường minh.
  │
  ├── Kỹ sư & Nhà phát triển phần mềm mô phỏng:
  │     Hiểu rõ cơ chế xây dựng ma trận Hamilton, thuật toán chéo hóa lặp, 
  │     kỹ thuật giả thế Vanderbilt USPP để tối ưu hóa mã nguồn mô phỏng.
  │
  ├── Nhóm R&D Doanh nghiệp Bán dẫn:
  │     Sở hữu cơ sở lý thuyết chuẩn xác để thiết lập thông số mô phỏng ab initio,
  │     rút ngắn thời gian screening vật liệu mới và giảm 40% chi phí thử nghiệm.
  │
  └── Các nhà nghiên cứu Vật lý chất rắn:
        Tài liệu tham khảo toàn diện về sự khác biệt hiệu năng giữa Frozen Phonon 
        và DFPT, định hướng cho các nghiên cứu cấu trúc điện tử nâng cao.

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

1. Yêu cầu tài nguyên phần cứng và phần mềm tối thiểu để chạy mô phỏng DFT bán dẫn là gì?

Hệ thống cần tối thiểu CPU 8 cores (x86_64), 16GB RAM để mô phỏng các ô mạng cơ bản (2 - 8 nguyên tử như ô mạng nguyên thủy Si, GaAs). Để tính toán siêu ô mạng sai hỏng hoặc phonon đông lạnh, khuyến nghị sử dụng máy trạm/HPC 32+ cores, 64GB - 128GB RAM chạy hệ điều hành Linux cùng các thư viện song song MPI và Intel MKL/FFTW3.

2. Giới hạn mở rộng quy mô (Scalability) của phương pháp DFT sóng phẳng là gì?

Độ phức tạp tính toán chuẩn của DFT tự hợp là $\mathcal{O}(N^3)$ với $N$ là số lượng nguyên tử. Khi vượt quá 500 - 1000 nguyên tử, chi phí chéo hóa ma trận và bộ nhớ lưu trữ sóng phẳng trở thành điểm nghẽn lớn. Giải pháp khắc phục là sử dụng giả thế siêu mềm USPP, hàm cơ sở định xứ (Local Orbitals) hoặc thuật toán tuyến tính $\mathcal{O}(N)$ (Order-N methods).

3. Tại sao LDA lại tính sai độ rộng vùng cấm (Bandgap) của bán dẫn nhưng vẫn cho hằng số mạng chính xác?

LDA xác định năng lượng trạng thái cơ bản $E_0$ thông qua mật độ điện tích $n(\mathbf{r})$ dựa trên định lý Hohenberg-Kohn. Tại trạng thái cơ bản, mật độ biến thiên chậm được miêu tả rất tốt, cho sai số độ dài liên kết và hằng số mạng dưới 1%. Tuy nhiên, vùng cấm là kích thích năng lượng (quasiparticle excitation), đòi hỏi phiếm hàm chứa bước nhảy thế đạo hàm mà LDA không có.

4. Khi nào nên dùng phương pháp Phonon đông lạnh (Frozen Phonon) và khi nào nên dùng DFPT?

Phương pháp Phonon đông lạnh phù hợp cho các mode dao động tại tâm vùng ($\Gamma$) hoặc biên vùng đối xứng cao trong các vật liệu đơn giản vì có thể tận dụng trực tiếp mã nguồn tính năng lượng toàn phần chuẩn mà không cần viết thêm mã nguồn chuyên dụng. DFPT vượt trội hoàn toàn khi cần vẽ toàn bộ đường cong tán sắc phonon trên toàn vùng Brillouin cho các cấu trúc phức tạp.

5. Sự khác biệt cốt lõi giữa Giả thế bảo toàn chuẩn (Norm-conserving) và Giả thế siêu mềm (USPP) là gì?

Giả thế bảo toàn chuẩn bắt buộc hàm sóng giả phải trùng hoàn toàn với hàm sóng thực bên ngoài bán kính cắt $r_c$ và tích phân bình phương hàm sóng trong vùng lõi phải bằng nhau. Giả thế siêu mềm Vanderbilt loại bỏ điều kiện này, cho phép hàm sóng giả cực kỳ trơn mượt trong lõi, giúp giảm năng lượng cắt $E_{cut}$ từ 40 Ry xuống còn 15 Ry, tiết kiệm hơn 60% số lượng sóng phẳng.

Kết luận

Khóa luận tốt nghiệp "Lý thuyết phiếm hàm mật độ và các cách tiếp cận khi nghiên cứu bán dẫn" đã giải quyết trọn vẹn bài toán xây dựng bức tranh vật lý lượng tử toàn diện về mô phỏng cấu trúc vi mô vật liệu bán dẫn. Từ việc hệ thống hóa tinh thể học, mạng đảo và sai hỏng tinh thể, công trình đã làm sáng tỏ bản chất biến phân của hai định lý Hohenberg-Kohn, hệ phương trình tự hợp Kohn-Sham và các phép gần đúng thực tế LDA/GGA.

Công trình mang giá trị học thuật và ứng dụng thực tiễn cao khi phân tích sâu sắc các kỹ thuật tính toán tiên tiến: Giả thế sóng phẳng (PPW), Giả thế siêu mềm Vanderbilt (USPP), Phương pháp đóng băng nhân (FCA), Phương pháp ma trận điện môi, Phonon đông lạnh và Động lực học phân tử ab initio (AIMD). Đây là tài liệu nền tảng, chuẩn mực cho sinh viên, kỹ sư R&D và các nhà nghiên cứu trong việc tiếp cận công nghệ mô phỏng vật liệu bán dẫn thế hệ mới.