Giới thiệu dự án
Ngành công nghiệp vi mạch bán dẫn toàn cầu đang chứng kiến sự chuyển dịch mang tính bước ngoặt khi quy trình thu nhỏ kích thước bóng bán dẫn tiệm cận các giới hạn vật lý lượng tử của vật liệu Silicon truyền thống. Theo các báo cáo từ Hiệp hội Công nghiệp Bán dẫn (SIA), nhu cầu tối ưu hóa hiệu suất, giảm thiểu suy hao năng lượng và tích hợp mật độ linh kiện cao đòi hỏi các giải pháp đột phá về cả khoa học vật liệu lẫn công nghệ thiết kế hỗ trợ máy tính (TCAD). Các chất bán dẫn oxide kim loại dạng màng mỏng (TCO) và bán dẫn hai chiều đang trở thành tâm điểm nghiên cứu nhằm thay thế dần các kênh dẫn Silicon và Germanium thế hệ cũ.
Vấn đề cốt lõi đặt ra là phần lớn các oxide kim loại tự nhiên đều thể hiện tính dẫn loại n ($n\text{-type}$), trong khi các vật liệu bán dẫn oxide loại p ($p\text{-type}$) có độ linh động hạt tải cao và độ rộng vùng cấm lớn rất hiếm và khó tổng hợp do sự định vị sâu của orbital $O\ 2p$ tại đỉnh vùng hóa trị (Valence Band Maximum - VBM). Thiếc monoxide ($\text{SnO}$) nổi lên như một ứng viên tiềm năng bậc nhất với độ rộng vùng cấm trực tiếp $E_g \approx 3.6\text{ eV}$ và độ linh động lỗ trống $\mu_p \approx 1.851 - 5.27\text{ cm}^2/\text{Vs}$. Tuy nhiên, chi phí hóa chất tiền chất tinh khiết cao, quy trình lắng đọng màng mỏng phức tạp và rủi ro hỏng hóc trong phòng thí nghiệm sạch (Cleanroom) tạo nên rào cản lớn cho việc thử nghiệm trực tiếp. Do đó, việc ứng dụng mô phỏng số đa vật lý nhằm dự báo chính xác các đặc tính điện trước khi chế tạo thực nghiệm là một yêu cầu cấp thiết.
+-----------------------------------------------------------------------------------+
| QUY TRÌNH TỔNG QUAN DỰ ÁN |
+-----------------------------------------------------------------------------------+
| [Thiết kế 3D CAD & Vật liệu] --> [Mô phỏng Đa vật lý FEM (COMSOL)] |
| | |
| v |
| [Thực nghiệm Phún xạ DC Magnetron] --> [Đo đạc Đặc tuyến I-V & Đối chuẩn Sai số] |
+-----------------------------------------------------------------------------------+
Dự án khóa luận tập trung giải quyết bài toán trên thông qua 5 mục tiêu cụ thể:
- Xây dựng mô hình hình học 3D hoàn chỉnh của diode chuyển tiếp dị thể $p\text{-SnO} / n\text{-a-Si:H}$ trên nền tảng phần mềm mô phỏng đa vật lý COMSOL Multiphysics.
- Ứng dụng phương pháp phần tử hữu hạn (Finite Element Method - FEM) giải hệ phương trình bảo toàn dòng điện và phân bố thế điện tĩnh trong miền không gian ba chiều dưới chế độ quét thế phân cực thuận (Forward Bias Sweep từ $0\text{V}$ đến $2\text{V}$, bước nhảy $0.05\text{V}$).
- Đánh giá định lượng sự suy biến của mật độ dòng điện $J$ và điện trường $E$ khi biến thiên bề dày lớp $p\text{-SnO}$ ở 3 mức: $100\text{ nm}$, $150\text{ nm}$ và $200\text{ nm}$.
- Chế tạo mẫu thực nghiệm diode màng mỏng cấu trúc $\text{Cu} / p\text{-SnO} / n\text{-a-Si:H} / \text{Cu}$ bằng kỹ thuật phún xạ Magnetron DC tại Phòng Thí nghiệm Bộ môn Vật lý Chất rắn (Đại học Khoa học Tự nhiên, ĐHQG-HCM).
- Phân tích đối chuẩn đặc tuyến $I-V$ giữa mô phỏng và thực nghiệm, bóc tách nguyên nhân gây ra sự chênh lệch điện áp ngưỡng ($V_{th}$).
Kết quả kỳ vọng của công trình là thiết lập thành công mô hình tham số hóa kiểm chứng được phân bố điện trường cục bộ tại bề mặt tiếp giáp dị thể, giải thích cơ chế chuyển vận hạt tải và rút ngắn chu kỳ thử nghiệm vật liệu bán dẫn mới. Phạm vi nghiên cứu giới hạn ở chế độ hoạt động tĩnh (Stationary), điều kiện nhiệt độ phòng chuẩn ($T_{ref} = 25^\circ\text{C} - 32^\circ\text{C}$) và cấu trúc tiếp giáp diện tích lớn ($4.0\text{ cm} \times 2.5\text{ cm}$).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trong lĩnh vực chế tạo linh kiện điện tử màng mỏng, việc kết hợp các chất bán dẫn vô cơ truyền thống và oxide kim loại vô định hình luôn tồn tại sự đánh đổi giữa hiệu năng dẫn điện, độ phức tạp công nghệ và chi phí gia công.
| Tiêu chí so sánh |
Diode Silicon PN truyền thống |
Diode Bán dẫn Hữu cơ (OLED/OPV) |
Diode Dị thể $p\text{-SnO} / n\text{-a-Si:H}$ (Dự án) |
| Độ rộng vùng cấm ($E_g$) |
$1.1\text{ eV}$ (Hẹp, hấp thụ quang thấp) |
$1.8 - 2.5\text{ eV}$ (Trung bình) |
$3.6\text{ eV}$ (Rộng, truyền qua quang cao) |
| Độ linh động hạt tải |
Rất cao ($\mu_e \approx 1400, \mu_h \approx 500$) |
Rất thấp ($\mu < 0.1\text{ cm}^2/\text{Vs}$) |
Cao đối với Oxide p ($\mu_p \approx 1.851 - 5.27\text{ cm}^2/\text{Vs}$) |
| Nhiệt độ chế tạo |
Rất cao ($>900^\circ\text{C}$ khuếch tán) |
Thấp ($<150^\circ\text{C}$) |
Thấp đến trung bình ($100 - 200^\circ\text{C}$) |
| Khả năng tích hợp |
Đế cứng Si Wafer |
Đế dẻo linh hoạt |
Đế thủy tinh/đế dẻo đa dụng |
| Chi phí chế tạo mẫu |
Cực cao (Yêu cầu Fab tiêu chuẩn) |
Trung bình |
Tối ưu nhờ mô phỏng FEM định hướng |
Yêu cầu kỹ thuật của hệ thống mô phỏng và thực nghiệm được phân loại theo ma trận MoSCoW:
- Must have (Bắt buộc): Mô hình hóa cấu trúc 3D chính xác kích thước màng mỏng; thiết lập bộ tham số vật liệu thực ($\sigma$, $\varepsilon_r$, $\rho_0$); giải phương trình dòng điện theo định luật Ohm mở rộng; đo đạc thực nghiệm đường cong $I-V$.
- Should have (Nên có): Chia lưới Sweep Mesh tùy biến cho tỷ lệ kích thước biên lớn ($Aspect\ Ratio > 10^5$); mô phỏng quét tham số bề dày ($Auxiliary\ Sweep$).
- Could have (Có thể có): Tích hợp hiệu ứng điện trở tiếp xúc nhiệt độ tuyến tính $\rho(T) = \rho_0(1 + \alpha(T - T_{ref}))$.
- Won't have (Chưa thực hiện): Mô phỏng đầy đủ động học tái hợp Shockley-Read-Hall (SRH) và tiếp giáp Schottky lượng tử.
Thách thức kỹ thuật lớn nhất nằm ở tỷ lệ hình học cực hạn: diện tích bề mặt tiếp xúc là $4.0\text{ cm} \times 2.5\text{ cm}$ ($10^{-2}\text{ m}$) trong khi độ dày các lớp chỉ từ $150\text{ nm}$ đến $500\text{ nm}$ ($10^{-7}\text{ m}$). Sự chênh lệch bậc độ lớn lên tới $10^5$ lần khiến việc chia lưới tự động thông thường làm tràn bộ nhớ RAM hoặc tạo ra các phần tử lưới có tỷ số co giãn méo mó (skewness), dẫn đến hiện tượng không hội tụ ma trận độ cứng FEM.
Thiết kế hệ thống
Kiến trúc mô phỏng được phân rã thành các miền không gian hình khối 3D ghép lớp. Do điện cực Anode phía trên chỉ có kích thước $1.0\text{ cm} \times 1.0\text{ cm}$ đặt tại tâm, miền hình học được phân tách thành 5 khối con đối xứng nhằm đồng nhất mặt chiếu từ trên xuống, tối ưu hóa quá trình quét lưới theo phương thẳng đứng.
+---------------------------------------+
| Anode Cu (Domain 13): 350nm (1x1cm) |
+------+---------------------------------------+------+
| p-SnO (Domains 3,6,9,12,16): 150nm |
+-----------------------------------------------------+
| n-a-Si:H (Domains 2,5,8,11,15): 200nm |
+-----------------------------------------------------+
| Cathode Cu (Domains 1,4,7,10,14): 500nm |
+-----------------------------------------------------+
Cấu hình ngăn xếp công nghệ và thông số vật liệu thiết lập trong COMSOL Multiphysics:
COMSOL Multiphysics (v5.6a / v6.0)
├── Physics Module: AC/DC -> Electric Currents (ec)
├── Study Type: Stationary Solver with Auxiliary Continuation Sweep
└── Material Database Configuration:
├── Copper (Cu - Anode/Cathode):
│ ├── Relative Permittivity (εr): 1.0
│ └── Reference Resistivity (ρ0): 1.72e-8 Ω·m
├── p-type Tin Monoxide (SnO):
│ ├── Relative Permittivity (εr): 8.86
│ ├── Electrical Conductivity (σ): 0.10 S/cm (10 S/m)
│ └── Reference Resistivity (ρ0): 9.81 Ω·m
└── n-type Hydrogenated Amorphous Silicon (a-Si:H):
├── Relative Permittivity (εr): 11.7
├── Electrical Conductivity (σ): 7.84 S/cm (784 S/m)
└── Reference Resistivity (ρ0): 1.272e-3 Ω·m
Hệ phương trình vi phân đạo hàm riêng (PDE) chủ đạo chi phối tương tác vật lý:
$$\nabla \cdot \mathbf{J} = Q_{j,v}$$
$$\mathbf{J} = \sigma \mathbf{E} + \mathbf{J}_e$$
$$\mathbf{E} = -\nabla V$$
Trong đó $\mathbf{J}$ là mật độ dòng điện tổng ($\text{A/m}^2$), $\mathbf{E}$ là cường độ điện trường ($\text{V/m}$), $V$ là thế điện tĩnh ($\text{V}$), $\sigma$ là độ dẫn điện ($\text{S/m}$), và nguồn dòng cục bộ bên ngoài được xác định $Q_{j,v} = 0$, $\mathbf{J}_e = 0$.
Methodology
Nghiên cứu áp dụng phương pháp luận phối hợp song song giữa tính toán mô phỏng phần tử hữu hạn và chế tạo thực nghiệm vật lý chất rắn:
[BƯỚC 1: TIỀN XỬ LÝ MÔ PHỎNG]
├── Dựng Geometry 3D 5 miền
├── Gán điều kiện biên V_Anode = VDD, V_Cathode = GND
└── Tạo lưới quét Swept Mesh độ mịn Extra Fine
│
v
[BƯỚC 2: GIẢI SỐ TỰ ĐỘNG]
├── Bộ giải giải tích Stationary
└── Quét tham số Auxiliary Sweep (0V -> 2V, step 0.05V)
│
v
[BƯỚC 3: GIA CÔNG THỰC NGHIỆM]
├── Tẩy rửa siêu âm đế thủy tinh -> Sấy chân không
├── Phún xạ DC Magnetron màng SnO (Áp suất 10^-3 mTorr)
└── Lắng đọng màng n-a-Si:H và bốc bay điện cực Cu
│
v
[BƯỚC 4: HẬU XỬ LÝ & ĐỐI CHUẨN]
├── Tích phân miền không gian tính I_total = ∬ J dS
├── Đo đặc tuyến thực nghiệm bằng hệ đo Probe Station
└── Bóc tách sai số và đánh giá cơ chế rào thế tiếp xúc
Implementation và kết quả
Development process
Quá trình xây dựng mô hình trong COMSOL Multiphysics được hiện thực hóa thông qua việc thiết lập trực tiếp các toán tử tích phân miền không gian và cấu hình biên vật lý:
% Định nghĩa các toán tử tích phân hình học trong COMSOL Model Java/API
model.component("comp1").cparam.set("VDD", "2[V]");
model.component("comp1").cparam.set("d_p", "150[nm]");
% Thiết lập tích phân biên phục vụ trích xuất dữ liệu dòng điện tổng
model.component("comp1").cxtrn.create("P_volume", "Integration");
model.component("comp1").cxtrn("P_volume").set("entitydim", 3);
model.component("comp1").cxtrn("P_volume").selection.set([3 6 9 12 16]);
model.component("comp1").cxtrn("N_volume", "Integration");
model.component("comp1").cxtrn("N_volume").set("entitydim", 3);
model.component("comp1").cxtrn("N_volume").selection.set([2 5 8 11 15]);
% Cấu hình điều kiện biên dẫn điện (Electric Currents Interface)
model.component("comp1").physics("ec").feature("pot1").set("V0", "VDD");
model.component("comp1").physics("ec").feature("gnd1").selection.set([1 4 7 10 14]);
Kỹ thuật chia lưới thích ứng (Custom Swept Meshing) được áp dụng: Mặt phẳng 2D trên cùng của lớp Anode được tạo lưới tam giác tự do (Free Triangular) với chất lượng siêu mịn (Extra Fine), sau đó thực hiện phép quét (Swept) tịnh tiến xuyên qua độ dày của từng lớp màng mỏng ($150\text{ nm} - 500\text{ nm}$). Kỹ thuật này giới hạn số lượng phần tử ở mức $1.2 \times 10^5$ bậc tự do (DOFs), loại bỏ hoàn toàn lỗi tràn bộ nhớ và đảm bảo sai số phần dư chuẩn hóa nhỏ hơn $10^{-6}$.
+---------------------------------------------------------------+
| CẤU TRÚC LƯỚI QUÉT SWEPT MESH CHI TIẾT |
+---------------------------------------------------------------+
| [Free Triangular Mesh trên bề mặt Anode/Cathode] |
| │ |
| ▼ (Quét tịnh tiến đa lớp theo trục Z) |
| [Lớp 1: Cu Anode (350nm) -> 5 lớp phần tử phân bố] |
| [Lớp 2: p-SnO (150nm) -> 8 lớp phần tử phân bố cao độ] |
| [Lớp 3: n-a-Si:H (200nm) -> 8 lớp phần tử phân bố cao độ] |
| [Lớp 4: Cu Cathode (500nm) -> 5 lớp phần tử phân bố] |
+---------------------------------------------------------------+
Testing và validation
Quá trình quét tham số phân cực thuận được thực hiện trên toàn dải điện áp $0\text{V} \le V_{Anode} \le 2.0\text{V}$ với bước nhảy $\Delta V = 0.05\text{V}$ (tổng cộng 41 bước tính toán trên mỗi kịch bản độ dày).
Dòng điện I (A)
^
| * Mô phỏng (100nm)
| * + Mô phỏng (150nm)
| * + # Mô phỏng (200nm)
| * + # - Thực nghiệm (150nm)
| * + #
| * + #
| * + #
| * + #
| * + # ------------------ Thực nghiệm (Vth = 0.92V)
+--------------------+-------------------> Điện áp cấp V (Volt)
0 0.92 1.52 2.0
Các kết quả đo đạc phân bố trường cục bộ cho thấy:
- Mật độ dòng điện ($J$): Tập trung cao độ tại khu vực diện tích đối xứng dưới chân đế cực Anode ($1.0\text{ cm} \times 1.0\text{ cm}$) và suy giảm phi tuyến dạng hàm mũ về phía các mép biên tự do của màng bán dẫn.
- Phân bố điện trường ($E$): Gradient điện trường đạt cực đại tại mặt tiếp xúc dị thể $p\text{-SnO} / n\text{-a-Si:H}$ với độ lớn xấp xỉ $1.3 \times 10^7\text{ V/m}$ ở điện áp phân cực $2\text{V}$.
Kết quả đạt được
Bảng đối chuẩn định lượng giữa mô phỏng số và thực nghiệm chế tạo:
| Đại lượng đo lường |
Mô phỏng COMSOL ($d_p = 150\text{ nm}$) |
Thực nghiệm chế tạo ($d_p = 150\text{ nm}$) |
Sai lệch tuyệt đối / Ghi chú |
| Điện áp ngưỡng mở ($V_{th}$) |
$1.52\text{ V}$ |
$0.92\text{ V}$ |
$\Delta V = 0.60\text{ V}$ ($39.4%$) |
| Dòng thuận tại $V = 2.0\text{ V}$ |
$12.4\text{ mA}$ |
$8.7\text{ mA}$ |
Do điện trở tiếp xúc ký sinh |
| Dòng rò nghịch ($V = -2.0\text{ V}$) |
$\approx 10^{-7}\text{ A}$ |
$4.2 \times 10^{-6}\text{ A}$ |
Ảnh hưởng bởi dòng rò bề mặt |
| Tỷ số chỉnh lưu ($I_{on}/I_{off}$) |
$> 10^4$ |
$\approx 2.1 \times 10^3$ |
Đặc tính diode rõ nét |
| Độ dẫn điện màng SnO ($\sigma$) |
$0.10\text{ S/cm}$ |
$0.086\text{ S/cm}$ (Đo Hall) |
Phù hợp với màng ủ nhiệt |
+------------------------------------------------------------------------------------+
| ẢNH HƯỞNG CỦA BỀ DÀY LỚP P-SnO LÊN DÒNG THUẬN (V = 2.0V) |
+------------------------------------------------------------------------------------+
| Độ dày d_p = 100 nm | Mật độ dòng J = 14.2 mA/cm² (Điện trở khối thấp nhất) |
| Độ dày d_p = 150 nm | Mật độ dòng J = 12.4 mA/cm² (Điểm tối ưu thực nghiệm) |
| Độ dày d_p = 200 nm | Mật độ dòng J = 10.1 mA/cm² (Tổn hao điện trở khối tăng) |
+------------------------------------------------------------------------------------+
Sự chênh lệch giữa điện áp ngưỡng mô phỏng ($1.52\text{ V}$) và thực nghiệm ($0.92\text{ V}$) bắt nguồn từ việc mô hình dẫn điện tuyến tính trong module AC/DC xem các lớp tiếp xúc kim loại - bán dẫn ($\text{Cu}/p\text{-SnO}$ và $n\text{-a-Si:H}/\text{Cu}$) là tiếp xúc Ohmic lý tưởng không có điện trở tiếp xúc bề mặt. Trong thực tế chế tạo màng mỏng, các trạng thái bẫy bề mặt (interface trap states), sai hỏng mạng tinh thể (nút khuyết thiếc $V_{Sn}$ và oxy xen kẽ $O_i$) tạo ra rào thế Schottky phi lý tưởng và điện trở nối tiếp phân tán, dẫn đến hiện tượng sụt áp tiếp xúc thực tế.
Đổi mới và đóng góp
Nghiên cứu mang lại những đóng góp học thuật và kỹ thuật quan trọng cho lĩnh vực linh kiện bán dẫn màng mỏng:
Orbital O 2p (Sâu, định vị cao) <--- Lai hóa ---> Orbital Sn 5s (Hình cầu, mở rộng)
│
▼
[VBM mở rộng, giảm khối lượng hiệu dụng lỗ trống m_h*]
│
▼
[Độ linh động lỗ trống μ_p tăng vượt bậc: ~5.27 cm²/Vs]
- Làm sáng tỏ cơ chế dẫn điện của vật liệu mới $p\text{-SnO}$: Giải thích bản chất tính dẫn loại p tự nhiên xuất phát từ năng lượng hình thành thấp của các nút khuyết thiếc ($V_{Sn}$). Sự lai hóa giữa orbital bất đối xứng $\text{Sn}\ 5s$ và $\text{O}\ 2p$ giải phóng đỉnh vùng hóa trị (VBM), làm giảm khối lượng hiệu dụng của lỗ trống ($m_h^*$), tạo ra độ linh động $\mu_p \approx 1.851 - 5.27\text{ cm}^2/\text{Vs}$, vượt trội so với các oxide kim loại loại p truyền thống như $\text{NiO}$ ($\mu_p < 0.1\text{ cm}^2/\text{Vs}$) và $\text{Cu}_2\text{O}$ ($\mu_p \approx 1.0\text{ cm}^2/\text{Vs}$).
- Phương pháp chia lưới swept mesh giải quyết tỷ lệ hình học $10^5$: Đưa ra kỹ thuật phân tách 5 khối miền không gian để xử lý bài toán hội tụ trường điện tĩnh trong cấu trúc màng mỏng có bề dày cỡ nanomet trên bề mặt centimet mà không làm bùng nổ tài nguyên tính toán.
- Quy trình kết nối TCAD - Thực nghiệm tinh gọn: Xây dựng cầu nối trực tiếp giữa tính toán số đa vật lý và quy trình phún xạ màng mỏng Magnetron DC, giúp định hình thông số chế tạo tối ưu (bề dày $p\text{-SnO} = 150\text{ nm}$ cho dòng dẫn ổn định nhất).
Ứng dụng thực tế và triển khai
Cấu trúc diode dị thể $p\text{-SnO} / n\text{-a-Si:H}$ mở ra tiềm năng thương mại hóa và ứng dụng kỹ thuật trong nhiều lĩnh vực:
+----------------------------------------------------------------------------------+
| HỆ SINH THÁI ỨNG DỤNG CỦA MÀNG MỎNG p-SnO |
+----------------------------------------------------------------------------------+
| [Thiết bị Hiển thị Trong suốt] ---> TFT suốt quang học, màn hình Micro-LED |
| [Cảm biến Quang & Khí sinh học] ---> Cảm biến phát hiện UV, màng hấp phụ NO2/CO |
| [Năng lượng Tái tạo] ---> Diode tách quang pin mặt trời thế hệ mới |
| [Mạch Tích hợp 3D Phía sau] ---> Lớp Diode giải mã chọn lọc bộ nhớ ReRAM |
+----------------------------------------------------------------------------------+
Lộ trình triển khai ứng dụng công nghiệp và hiệu quả kinh tế:
Giai đoạn 1: Chuẩn hóa Model TCAD (Tháng 1 - 3)
├── Tích hợp module Drift-Diffusion & điều kiện bẫy rào Schottky
└── Sai số mô phỏng - thực nghiệm < 10%
│
v
Giai đoạn 2: Tối ưu Quy trình Phún xạ RF/DC (Tháng 4 - 8)
├── Kiểm soát nồng độ dòng khí Ar/O2 ở mức 5 - 10 SCCM
└── Gia nhiệt đế ổn định 150°C - 200°C tăng độ tinh thể màng
│
v
Giai đoạn 3: Chế tạo Module Mạch Màng mỏng Thử nghiệm (Tháng 9 - 12)
├── Tích hợp diode trên đế thủy tinh kích thước 100mm x 100mm
└── Thử nghiệm độ bền hoạt động liên tục 1000 giờ
- Phân tích chi phí - lợi ích (Cost-Benefit Analysis): Việc sử dụng COMSOL Multiphysics để quét tối ưu hóa bề dày và diện tích điện cực trước khi chế tạo giúp giảm thiểu $65%$ số lượng ca thực nghiệm phún xạ trong phòng sạch, tiết kiệm trung bình $3,000 - $5,000 USD chi phí bia vật liệu tinh khiết và khí phân tích cho mỗi chu kỳ R&D linh kiện mới.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Mô hình vật lý dòng điện giản lược: Nghiên cứu hiện tại sử dụng module Electric Currents (
ec) dựa trên định luật Ohm mở rộng và điện trở suất tuyến tính hóa thay vì giải hệ phương trình tự nhất quán Poisson kết hợp trôi - khuếch tán (Poisson - Drift-Diffusion) của module Semiconductor chuyên sâu.
- Bỏ qua hiệu ứng tiếp xúc mặt ngoài: Chưa tích hợp rào thế tiếp xúc dị thể (heterojunction band offsets), thế chắn Schottky tại bề mặt tiếp xúc kim loại dẫn $\text{Cu}$ và các mức bẫy tái hợp bề mặt.
- Giới hạn kích thước chế tạo: Diện tích linh kiện thực nghiệm còn lớn ($4.0\text{ cm} \times 2.5\text{ cm}$) do giới hạn của hệ thống mặt nạ che (shadow mask) cơ khí sẵn có, chưa đạt tới kích thước micro/nano của vi mạch tích hợp.
Hướng phát triển
-
Nâng cấp mô hình bán dẫn toàn diện: Chuyển dịch toàn bộ mô hình sang module Semiconductor (semi) của COMSOL, giải đồng thời phương trình liên tục cho điện tử và lỗ trống:
$$\nabla \cdot \mathbf{J}n = q R{SRH}, \quad \nabla \cdot \mathbf{J}p = -q R{SRH}$$
-
Tối ưu hóa công nghệ lắng đọng: Áp dụng hệ phún xạ Magnetron RF kết hợp kiểm soát lưu lượng khí Oxy tinh vi bằng bộ điều khiển dòng khối (Mass Flow Controller - MFC) ở dải $1 - 2\text{ SCCM}$ nhằm triệt tiêu các donor bù trừ tự sinh, nâng cao độ tinh khiết pha của màng $p\text{-SnO}$.
-
Mở rộng cấu trúc linh kiện: Phát triển mô hình bóng bán dẫn màng mỏng kênh dẫn kép (Complementary TFT) và cảm biến khí oxide kim loại hoạt động ở nhiệt độ phòng.
Đối tượng hưởng lợi
+-----------------------------------------------------------------------------------+
| GIÁ TRỊ MANG LẠI CHO CÁC BÊN |
+-----------------------------------------------------------------------------------+
| Sinh viên & Học viên | Nắm vững phương pháp mô phỏng TCAD và công nghệ màng mỏng |
| Kỹ sư Vi mạch Back-end| Hiểu sâu cơ chế trường tĩnh tầng vật lý và tiếp xúc kim loại |
| Viện & Nhóm nghiên cứu| Rút ngắn quy trình khảo sát vật liệu bán dẫn oxide 2D mới |
| Doanh nghiệp Bán dẫn | Cắt giảm 65% chi phí chế tạo thử nghiệm R&D phòng sạch |
+-----------------------------------------------------------------------------------+
- Sinh viên & Học viên ngành Bán dẫn: Tiếp cận tài liệu chuẩn hóa về quy trình xây dựng mô hình FEM 3D màng mỏng, nắm bắt kỹ thuật xử lý lưới swept mesh trên các miền kích thước chênh lệch cực lớn.
- Kỹ sư Thiết kế Vi mạch & Thiết bị (TCAD Engineers): Sở hữu bộ tham số vật liệu chi tiết của $\text{SnO}$ và $a\text{-Si:H}$ được đối chuẩn trực tiếp với thực nghiệm để tích hợp vào các thư viện mô phỏng linh kiện.
- Nhà nghiên cứu Khoa học Vật liệu: Phương pháp luận rõ ràng trong việc khảo sát đặc tính chuyển tiếp dị thể $p-n$ của các vật liệu bán dẫn oxide mới nổi.
- Doanh nghiệp Sản xuất Linh kiện Điện tử: Tiết kiệm tài nguyên, thời gian và chi phí vận hành máy lắng đọng chân không thông qua việc ứng dụng công cụ mô phỏng số tiền khả thi.
Câu hỏi thường gặp
1. Yêu cầu cấu hình phần cứng máy tính để chạy mô hình FEM màng mỏng trong COMSOL là gì?
Mô hình yêu cầu tối thiểu CPU 4 nhân (khuyến nghị Intel Core i7/AMD Ryzen 7 thế lên), RAM tối thiểu 16 GB (khuyến nghị 32 GB để xử lý các phép quét tham số mật độ cao) và ổ cứng SSD NVMe để lưu trữ các tập tin ma trận độ cứng trung gian có kích thước từ $2 - 5\text{ GB}$.
2. Tại sao lưới tự động (Physics-controlled Mesh) không thể áp dụng cho mô hình này?
Do tỷ lệ kích thước biên hình học ($Aspect\ Ratio$) giữa chiều dài ($4\text{ cm}$) và độ dày màng ($150\text{ nm}$) vượt quá $10^5$ lần. Lưới tự động sẽ tạo ra các phần tử tứ diện vô cùng dày đặc làm tràn bộ nhớ RAM (vượt quá $128\text{ GB}$) hoặc tạo ra các phần tử quá méo mó gây lỗi không hội tụ toán tử giải ma trận vi phân.
3. Làm thế nào để khắc phục hiện tượng sụt áp và sai lệch $V_{th}$ giữa mô phỏng và thực tế?
Cần bổ sung tính chất rào thế kim loại - bán dẫn (Schottky Barrier Height) tại các mặt biên cực Anode/Cathode và khai báo lớp điện trở màng tiếp xúc bề mặt (Contact Resistance Interface) trong phần mềm COMSOL để bù trừ cho điện trở tiếp xúc thực tế của hệ đo.
4. Tại sao SnO lại có tính dẫn điện loại p tự nhiên thay vì loại n như đa số các oxide kim loại khác?
Tính dẫn loại p của $\text{SnO}$ được quy định bởi năng lượng hình thành cực thấp của các nút khuyết thiếc ($V_{Sn}$) trong mạng tinh thể bốn phương. Đồng thời, sự lai hóa mạnh mẽ giữa orbital $\text{Sn}\ 5s$ và $\text{O}\ 2p$ giúp giải tỏa sự định vị cục bộ ở đỉnh vùng hóa trị, cho phép các lỗ trống tự do di chuyển dễ dàng.
5. Chi phí triển khai mô phỏng COMSOL so với chế tạo thực nghiệm màng mỏng chênh lệch như thế nào?
Chi phí trang bị bản quyền phần mềm và máy trạm tính toán chỉ cần đầu tư một lần, trong khi một ca chế tạo màng mỏng trong phòng sạch tiêu tốn từ $500 - $1,500 USD cho khí phân tích, bia kim loại độ tinh khiết cao ($99.99%$), hóa chất tẩy rửa quang khắc và điện năng chân không siêu cao. Ứng dụng mô phỏng giúp giảm thiểu từ $60 - 70%$ số lần thử nghiệm thất bại.
Kết luận
Đề tài khóa luận tốt nghiệp "Ứng dụng phần mềm COMSOL Multiphysics vào nghiên cứu và mô phỏng Diode" đã hoàn thành xuất sắc việc xây dựng mô hình số 3D khảo sát thành công tiếp giáp dị thể giữa vật liệu mới $p\text{-SnO}$ và $n\text{-a-Si:H}$. Công trình không chỉ làm chủ kỹ thuật mô phỏng phần tử hữu hạn (FEM) xử lý tỷ lệ hình học cực hạn mà còn kiểm chứng thành công các đặc tính chỉnh lưu $I-V$ thông qua quy trình chế tạo thực nghiệm màng mỏng phún xạ Magnetron DC.
Sự kết hợp giữa công nghệ TCAD tiên tiến và thực nghiệm vật lý chất rắn đóng vai trò nền tảng, mở đường cho việc phát triển các thế hệ linh kiện bán dẫn oxide trong suốt, pin mặt trời hiệu suất cao và cảm biến thông minh. Các kỹ sư, nhà nghiên cứu và sinh viên có thể kế thừa trực tiếp mô hình tham số này để mở rộng phát triển cho các kiến trúc linh kiện phức tạp hơn như Transistor màng mỏng (TFT) hay cấu trúc nhớ bán dẫn phi bay hơi trong tương lai.
