Bối cảnh và vấn đề nghiên cứu
Nhu cầu sử dụng năng lượng trên toàn cầu ngày càng gia tăng trong khi các nguồn nhiên liệu hóa thạch truyền thống đang dần cạn kiệt, đồng thời gây ra tình trạng ô nhiễm môi trường và những tác động tiêu cực đến sức khỏe con người. Việc chuyển dịch sang khai thác các nguồn năng lượng tái tạo sạch và bền vững, đặc biệt là năng lượng mặt trời thông qua hiệu ứng quang điện, trở thành giải pháp kỹ thuật tất yếu.
Trong số các công nghệ quang điện, pin mặt trời dị thể vô định hình/tinh thể Silic (a-Si:H/c-Si) thu hút sự quan tâm lớn nhờ sở hữu nhiều ưu điểm vượt trội: hiệu suất chuyển đổi cao, nhiệt độ chế tạo thấp (dưới hoặc bằng 200 °C) và hệ số suy giảm do nhiệt độ thấp. Dù cấu trúc tiếp xúc mặt sau dạng răng lược (IBC) có thể đạt hiệu suất lên đến 26,7%, quy trình chế tạo của cấu trúc này rất phức tạp, đòi hỏi chi phí sản xuất cao và khó triển khai ở quy mô thương mại.
Đối với dòng pin dị thể a-Si:H/c-Si hai mặt tiêu chuẩn, thiết kế emitter đặt ở phía trước (Front Emitter - FE) bộc lộ hạn chế cố hữu: khó cân bằng giữa tính chất điện và tính chất quang của lớp a-Si:H(p). Cụ thể, độ dẫn của lớp a-Si:H(p) vốn thấp hơn lớp a-Si:H(n); việc tăng nồng độ pha tạp để cải thiện độ dẫn sẽ gia tăng tốc độ tái hợp cặp điện tử - lỗ trống làm suy giảm thế hở mạch ($V_{oc}$), trong khi tăng độ dày lớp màng lại dẫn đến hiện tượng hấp thụ ký sinh làm giảm mật độ dòng ngắn mạch ($J_{sc}$).
Cấu trúc emitter đặt ở phía sau (Rear Emitter - RE) khắc phục được những hạn chế trên: lớp a-Si:H(n) ở mặt trước có hệ số hấp thụ thấp hơn lớp a-Si:H(p) nên cho phép nhiều ánh sáng đi vào lớp hấp thụ c-Si(n) hơn, giảm tổn thất quang học; hạt tải điện tử được thu thập theo phương thẳng đứng về lưới điện cực kim loại phía trước giúp giảm yêu cầu dẫn điện đối với lớp oxit dẫn điện trong suốt (TCO) và cho phép linh hoạt hơn trong thiết kế cấu trúc bề mặt. Mặc dù vậy, đa phần các công trình công bố về cấu trúc RE trước đây tập trung vào thực nghiệm tốn kém thời gian và chi phí, thiếu vắng các nghiên cứu mô phỏng khảo sát đồng thời và toàn diện ảnh hưởng của lớp TCO, lớp bán dẫn mặt trước a-Si:H(n+) và lớp hấp thụ c-Si(n).
Mục tiêu của đề tài là khảo sát, phân tích ảnh hưởng của các lớp vật liệu đến các thông số quang điện của pin mặt trời dị thể cấu trúc RE bằng phương pháp mô phỏng số, từ đó xác lập bộ thông số tối ưu để nâng cao hiệu suất chuyển đổi năng lượng. Khóa luận triển khai 4 nhiệm vụ nghiên cứu cụ thể:
- Khảo sát hiệu suất chuyển đổi năng lượng của pin mặt trời dị thể có cấu trúc emitter phía sau khi thay đổi công thoát điện tử của lớp TCO.
- Khảo sát ảnh hưởng khi thay đổi nồng độ pha tạp và độ dày của lớp bán dẫn mặt trước a-Si:H(n+).
- Khảo sát ảnh hưởng khi thay đổi mật độ oxy ($D_{od}$), điện trở suất và độ dày của lớp hấp thụ c-Si(n).
- Phân tích, đánh giá các kết quả đạt được từ đó rút ra điều kiện tối ưu để chế tạo pin mặt trời có hiệu suất cao nhất.
Đối tượng nghiên cứu là pin mặt trời dị thể vô định hình/tinh thể Silic cấu trúc emitter phía sau. Phạm vi nghiên cứu giới hạn trong tính toán mô phỏng số một chiều (1D) trên phần mềm AFORS-HET dưới điều kiện thử nghiệm tiêu chuẩn (STC: phổ bức xạ mặt trời AM1.5G, cường độ chiếu xạ 1000 W/m² và nhiệt độ hoạt động 25 °C).
Cơ sở lý thuyết và phương pháp
Khóa luận vận dụng lý thuyết vật lý chất bán dẫn và chuyển tiếp dị thể bán dẫn (heterojunction). Nguyên lý hoạt động dựa trên sự phân tách cặp điện tử - lỗ trống quang sinh nhờ điện trường nội tiếp xúc và sự khác biệt về độ rộng vùng cấm giữa Silic vô định hình ($E_g \approx 1,72\text{ eV}$) và Silic tinh thể ($E_g \approx 1,124\text{ eV}$). Sự chênh lệch này tạo ra độ lệch vùng dẫn ($\Delta E_c$) và độ lệch vùng hóa trị ($\Delta E_v$) tại các mặt tiếp xúc, hình thành rào thế ngăn cản hạt tải thiểu số khuếch tán sai hướng và cho phép hạt tải đa số di chuyển về điện cực tương ứng qua cơ chế phát xạ nhiệt hoặc xuyên hầm. Các cơ chế tái hợp được xem xét bao gồm tái hợp Shockley-Read-Hall (SRH) qua các mức bẫy khuyết tật trong vùng cấm và tái hợp Auger.
Phương pháp nghiên cứu chính là mô phỏng số 1D sử dụng phần mềm học thuật chuyên dụng AFORS-HET (Automat FOR Simulation of HETerostructures do Viện Helmholtz-Zentrum Berlin phát triển). Phần mềm chia quá trình tính toán thành hai phân hệ:
- Mô hình quang: Áp dụng khi năng lượng photon lớn hơn hoặc bằng năng lượng vùng cấm ($E_{ph} \ge E_g$). Tốc độ sinh cặp điện tử - lỗ trống $G(x,t)$ được tính toán dựa trên định luật hấp thụ Lambert-Beer đối với bề mặt gồ ghề và mô hình phản xạ liên tục/không liên tục đối với bề mặt phẳng.
- Mô hình điện: Áp dụng khi $E_{ph} \le E_g$. Phần mềm giải hệ phương trình vi phân Poisson kết hợp với các phương trình dòng liên tục cho điện tử và lỗ trống dưới các điều kiện biên tiếp xúc kim loại - bán dẫn (tiếp xúc Schottky hoặc Ohmic):
$$\frac{d^2\psi(x)}{dx^2} = \frac{q}{\epsilon}(n' - p + N_A^- - N_D^+)$$
Nguồn dữ liệu mô phỏng được thiết lập từ các thông số vật liệu thực tế đo đạc từ các nghiên cứu trước đây và cấu hình vật liệu chuẩn trong phần mềm.
| Thông số vật liệu | Lớp c-Si(n) | Lớp c-Si(n) bề mặt | Lớp a-Si:H(p) | Lớp a-Si:H(n+) |
|---|---|---|---|---|
| Độ dày lớp | Biến thiên | $3 \times 10^{-7}\text{ cm}$ | $7 \times 10^{-7}\text{ cm}$ | Biến thiên |
| Hằng số điện môi ($\epsilon$) | 11,9 | 11,9 | 11,9 | 11,9 |
| Ái lực điện tử $\chi$ (eV) | 4,05 | 4,05 | 3,9 | 3,9 |
| Độ rộng vùng cấm $E_g$ (eV) | 1,124 | 1,124 | 1,72 | 1,72 |
| Mật độ trạng thái vùng dẫn $N_c$ ($\text{cm}^{-3}$) | $2,843 \times 10^{19}$ | $2,846 \times 10^{19}$ | $1,0 \times 10^{20}$ | $1,0 \times 10^{21}$ |
| Mật độ trạng thái vùng hóa trị $N_v$ ($\text{cm}^{-3}$) | $2,682 \times 10^{19}$ | $2,685 \times 10^{19}$ | $1,0 \times 10^{20}$ | $1,0 \times 10^{21}$ |
| Độ linh động điện tử $\mu_n$ ($\text{cm}^2/\text{V}\cdot\text{s}$) | 582,7 | 1111 | 5 | 5 |
| Độ linh động lỗ trống $\mu_p$ ($\text{cm}^2/\text{V}\cdot\text{s}$) | 273,7 | 421,6 | 1 | 1 |
| Nồng độ pha tạp lỗ trống $N_A$ ($\text{cm}^{-3}$) | 0 | 0 | $2,08 \times 10^{20}$ | 0 |
| Nồng độ pha tạp điện tử $N_D$ ($\text{cm}^{-3}$) | $1,8 \times 10^7$ | $1,5 \times 10^{16}$ | 0 | Biến thiên |
| Hệ số tái hợp Auger điện tử ($\text{cm}^6/\text{s}$) | 0 | $2,2 \times 10^{-31}$ | 0 | 0 |
| Hệ số tái hợp Auger lỗ trống ($\text{cm}^6/\text{s}$) | 0 | $9,9 \times 10^{-32}$ | 0 | 0 |
Nội dung chính theo từng chương
Chương 1. Tổng quan về pin mặt trời
Chương này hệ thống hóa tiến trình phát triển của pin mặt trời từ năm 1839, phân loại chi tiết 4 thế hệ công nghệ quang điện:
- Thế hệ thứ nhất: Pin Silic dạng khối (đơn tinh thể đạt hiệu suất thương mại 20,4%, đa tinh thể đạt 12 - 14%), hạn chế lớn là giá thành nguyên liệu cao và quy trình xử lý nhiệt độ cao (800 - 900 °C).
- Thế hệ thứ hai: Pin màng mỏng (a-Si/μc-Si, CdTe, CIGS) giảm chi phí nhiệt độ chế tạo (~600 °C), nhưng hiệu suất thấp hơn (6 - 9%) và dễ suy giảm phẩm chất theo thời gian.
- Thế hệ thứ ba và thứ tư: Pin nano tinh thể, pin chất màu, polymer và pin lai vô cơ - hữu cơ, có ưu điểm giá thành rẻ nhưng độ bền kém và hiệu suất không ổn định.
Tác giả phân tích các phương trình xác lập đặc trưng dòng - thế ($J\text{ - }V$), định nghĩa các thông số cốt lõi: mật độ dòng ngắn mạch $J_{sc}$, thế hở mạch $V_{oc}$, hệ số lấp đầy $FF$ và hiệu suất chuyển đổi quang điện $\eta$. Nội dung trọng tâm của chương tập trung làm rõ cấu trúc pin dị thể HIT (Heterojunction with Intrinsic Thin layer) do Sanyo phát triển, trong đó lớp màng siêu mỏng a-Si:H(i) thụ động hóa các liên kết chưa bão hòa tại bề mặt đế c-Si, làm giảm các tâm tái hợp.
Phần so sánh giữa cấu trúc Emitter mặt trước (FE) và Emitter mặt sau (RE) chỉ ra: ở cấu trúc FE, lỗ trống di chuyển theo phương ngang qua lớp TCO có điện trở lớn; ở cấu trúc RE, phần lớn điện tử di chuyển theo phương thẳng đứng về lưới điện cực kim loại phía trước qua lớp c-Si(n) và lớp TCO ghép song song, giảm điện trở nối tiếp và tổn thất quang học. Chương 1 điểm qua các công trình tiêu biểu: Taguchi (2014, FE đạt $\eta = 24,7%$), Bivour (2012, RE đạt $\eta = 22,8%$, $FF = 81,5%$), Renaud Varache (2014), Sen Li (2018, RE đạt $\eta = 21,71%$), Sangho Kim (2018, $\eta = 22,34%$), Duy Phong Pham (2019, $\eta = 21,8%$), Janne-Petteri Niemelä (2019, $\eta > 21%$) và Park (2019, mô phỏng RE đạt $\eta = 24,01\text{ - }25,35%$).
Chương 2. Phần mềm mô phỏng AFORS - HET
Chương 2 mô tả chi tiết công cụ mô phỏng AFORS-HET phiên bản 1D. Cấu trúc mô phỏng được thiết lập gồm các lớp tuần tự: $$\text{Điện cực Ag trước} \ / \ \text{TCO trước (80 nm)} \ / \ \text{a-Si:H}(n^+) \ / \ \text{c-Si}(n)\text{ bề mặt} \ / \ \text{c-Si}(n)\text{ hấp thụ} \ / \ \text{a-Si:H}(p)\text{ (7 nm)} \ / \ \text{TCO sau} \ / \ \text{Điện cực Ag sau}$$
Tác giả trình bày phương pháp giải các phương trình vi phân bán dẫn một chiều dưới điều kiện biên tiếp xúc Schottky và Ohmic, tích hợp mô hình thống kê tái hợp Shockley-Read-Hall (SRH) và tái hợp Auger. Khung mô phỏng sử dụng phổ chuẩn toàn cầu AM1.5G, nhiệt độ 300 K. Các ưu điểm của AFORS-HET được tổng hợp: cho phép thiết lập linh hoạt thông số mật độ bẫy khuyết tật trong khối và tại bề mặt, mô hình hóa rào thế tiếp xúc tại các mặt phân cách dị thể, xuất các đường đặc trưng $J\text{ - }V$, hiệu suất lượng tử ($QE$), phân bố nồng độ hạt mang điện và giản đồ dải năng lượng.
Chương 3. Kết quả và thảo luận
Chương này trình bày các khảo sát mô phỏng và phân tích cơ chế vật lý:
- Ảnh hưởng của công thoát điện tử lớp TCO ($\Phi_{TCO}$): Khảo sát trong dải $\Phi_{TCO}$ từ 3,7 eV đến 4,3 eV (dải thực tế của màng mỏng ZnO:Al chế tạo bằng phương pháp phún xạ). Khi $\Phi_{TCO}$ tăng từ 3,7 eV đến 3,9 eV, $V_{oc}$ và $J_{sc}$ giảm nhẹ trong khi $FF$ giữ nguyên. Khi tiếp tục tăng $\Phi_{TCO}$ lên 4,3 eV, tất cả các thông số đều giảm mạnh: $V_{oc}$ giảm từ 740,6 mV xuống 739,5 mV; $J_{sc}$ giảm từ 39,50 mA/cm² xuống 37,85 mA/cm²; $FF$ giảm từ 85,16% xuống 85,14%; hiệu suất $\eta$ giảm từ 24,92% xuống 23,83%.
Giản đồ cấu trúc vùng năng lượng cho thấy khi $\Phi_{TCO}$ tăng, dải năng lượng bị uốn cong lên trên; thế tiếp xúc $V_D$ tại mặt phân cách Schottky TCO/a-Si:H(n+) tăng thêm 0,3 V, làm tăng rào thế $\Delta E_c$, cản trở điện tử di chuyển về phía TCO. Mật độ điện tử tại mặt tiếp xúc giảm từ $3,5 \times 10^{19}\text{ cm}^{-3}$ (ở 3,7 eV) xuống còn $2 \times 10^{14}\text{ cm}^{-3}$ (ở 4,3 eV). Đồng thời, rào thế $\Delta E_v$ bị hạ thấp khiến các lỗ trống tích tụ tại mặt phân cách a-Si:H(n+)/c-Si(n) đạt mật độ $5,7 \times 10^{10}\text{ cm}^{-3}$, làm tăng mạnh dòng bão hòa ngược $J_0$ do tái hợp, trực tiếp kéo tụt $V_{oc}$ và $FF$. Phổ hiệu suất lượng tử ($QE$) phản ánh rõ: ở bước sóng ngắn $\lambda = 300\text{ nm}$, $QE$ giảm sâu từ 0,66 xuống 0,23 do sự sụt giảm hấp thu hạt tải ở mặt trước; trong khi ở bước sóng dài $\lambda = 800\text{ nm}$, $QE$ không đổi (đạt 0,80) do ánh sáng hấp thụ sâu bên trong khối c-Si.
-
Ảnh hưởng của lớp a-Si:H(n+): Khảo sát nồng độ pha tạp $N_d$ từ $5 \times 10^{14}\text{ cm}^{-3}$ đến $5 \times 10^{18}\text{ cm}^{-3}$ (quanh giá trị mật độ sai hỏng mạng $5 \times 10^{16}\text{ cm}^{-3}$) và độ dày lớp $T_{a\text{-Si:H}(n^+)}$ (từ 5 nm đến 10 nm). Khi $N_d$ thấp, vùng nghèo điện tích không hình thành đủ điện trường, làm gia tăng tái hợp tại bề mặt. Lớp a-Si:H(n+) quá mỏng dẫn đến nguy cơ lớp TCO tiếp xúc trực tiếp với c-Si(n) gây lỗi đoản mạch cục bộ, trong khi lớp quá dày gây tổn thất hấp thụ quang học.
-
Ảnh hưởng của lớp hấp thụ c-Si(n): Khảo sát mật độ bẫy oxy $D_{od}$ ($1 \times 10^{10}\text{ cm}^{-3}$ và $1 \times 10^{11}\text{ cm}^{-3}$), điện trở suất $\rho$ (0,5 $\Omega\cdot\text{cm}$ và 5 $\Omega\cdot\text{cm}$) và độ dày lớp hấp thụ $T_{c\text{-Si}(n)}$ (150 μm và 250 μm). Khi mật độ oxy tăng, tốc độ tái hợp SRH trong khối c-Si tăng rõ rệt, làm suy giảm thời gian sống của hạt tải $\tau_{bulk}$ và độ dài khuếch tán hiệu dụng của lỗ trống $L_{p,eff}$, dẫn đến sự suy giảm của đường đặc trưng $J\text{ - }V$ và phổ $QE$. Việc tối ưu hóa độ dày ở mức thích hợp vừa đảm bảo hấp thụ trọn vẹn photon bước sóng dài, vừa hạn chế chiều dài quãng đường lỗ trống phải di chuyển xuyên qua đế về emitter mặt sau.
Kết quả và đóng góp
Khóa luận đã xác lập bộ số liệu mô phỏng định lượng chi tiết về sự phụ thuộc của các thông số pin dị thể RE vào các lớp vật liệu:
| Giá trị $\Phi_{TCO}$ (eV) | $V_{oc}$ (mV) | $J_{sc}$ ($\text{mA/cm}^2$) | $FF$ (%) | Hiệu suất $\eta$ (%) |
|---|---|---|---|---|
| 3,7 | 740,6 | 39,50 | 85,16 | 24,92 |
| 3,9 | 740,4 | 39,45 | 85,16 | 24,88 |
| 4,1 | 740,1 | 38,90 | 85,15 | 24,51 |
| 4,3 | 739,5 | 37,85 | 85,14 | 23,83 |
Sau khi kết hợp tối ưu hóa đồng thời các thông số: công thoát điện tử TCO ($\Phi_{TCO} = 3,7\text{ eV}$), nồng độ pha tạp và độ dày lớp bán dẫn mặt trước a-Si:H(n+), độ dày, điện trở suất và kiểm soát mật độ sai hỏng oxy lớp c-Si(n), hiệu suất chuyển đổi năng lượng cực đại của pin đạt tới 27,33% (được ghi nhận trong biên bản đánh giá của giảng viên hướng dẫn và đường đặc trưng $J\text{ - }V$ ở điều kiện tối ưu).
Đóng góp khoa học của đề tài:
- Khắc phục khoảng trống nghiên cứu mô phỏng đa thông số: Đây là công trình tính toán mô phỏng khảo sát đồng thời cả ba lớp cấu trúc then chốt (TCO, a-Si:H(n+) và c-Si(n)) của pin dị thể emitter phía sau.
- Cung cấp cơ chế giải thích biến thiên quang điện dựa trên giản đồ vùng năng lượng: Làm rõ bản chất vật lý của sự suy giảm $V_{oc}$ và $J_{sc}$ theo $\Phi_{TCO}$ thông qua sự thay đổi thế tiếp xúc Schottky $V_D$, mật độ hạt tải tích tụ và phổ hiệu suất lượng tử $QE$ ở dải bước sóng ngắn (300 - 700 nm).
- Cung cấp bộ tham số định hướng thực nghiệm: Kết quả tính toán chỉ ra điều kiện tối ưu để phòng thí nghiệm FM&D (Đại học Duy Tân) và các nhóm nghiên cứu thực nghiệm áp dụng trực tiếp vào quy trình phún xạ TCO và lắng đọng hóa học pha hơi tăng cường plasma (RF-PECVD), giúp tiết kiệm thời gian và chi phí vật liệu.
Hạn chế và hướng nghiên cứu tiếp
Khóa luận tồn tại các hạn chế kỹ thuật gắn liền với phạm vi đề tài:
- Toàn bộ kết quả đều thu được từ phương pháp mô phỏng số một chiều (1D) trên phần mềm AFORS-HET, chưa thực hiện kiểm chứng thực nghiệm bằng việc chế tạo linh kiện vật lý thực tế.
- Mô hình 1D chưa phản ánh đầy đủ các hiệu ứng không gian 2D/3D trong thực tế sản xuất, như hình thái bề mặt tạo vân kim tự tháp (texture) bất đối xứng, điện trở tiếp xúc giữa màng kim loại Ag và lớp TCO, cũng như điện trở nối tiếp của lưới kim loại ngón tay (finger grid).
- Theo nhận xét của giảng viên hướng dẫn, khâu xử lý và phân tích số liệu ở một số đồ thị chuyên sâu chưa đạt mức độ sâu sắc toàn diện.
Hướng nghiên cứu phát triển tiếp theo được xác định:
- Chế tạo thực nghiệm linh kiện pin mặt trời dị thể cấu trúc RE tại phòng thí nghiệm FM&D dựa trên bộ thông số tối ưu thu được từ mô phỏng.
- Mở rộng khảo sát các vật liệu oxit dẫn điện trong suốt thay thế không chứa Indium quý hiếm, cụ thể là màng ZnO:Al (AZO) lắng đọng bằng phương pháp lắng đọng lớp nguyên tử (ALD).
- Nghiên cứu thay thế lớp bán dẫn mặt trước bằng các màng mỏng hợp chất vi tinh thể như $\mu\text{c-SiO:H}(n)$ và $\mu\text{c-Si:H}(n)$ nhằm cải thiện hơn nữa độ rộng vùng cấm quang học, giảm thiểu tối đa hiện tượng hấp thụ ký sinh ở dải bước sóng ngắn.
Giá trị tham khảo
Bản khóa luận là tài liệu tham khảo hữu ích cho sinh viên, học viên cao học và cán bộ nghiên cứu thuộc các chuyên ngành Công nghệ Vật liệu, Vật lý kỹ thuật, Năng lượng tái tạo và Kỹ thuật vi điện tử quang điện.
Những nội dung có giá trị ứng dụng và tham khảo cao gồm:
- Bảng thông số đầu vào chi tiết của phần mềm AFORS-HET (Bảng 2 trong văn bản), bao gồm các thông số điện tử, thông số quang học, mật độ trạng thái dải dẫn/dải hóa trị và hệ số tái hợp Auger cho cả vật liệu a-Si:H và c-Si.
- Quy trình chuẩn mô phỏng số pin mặt trời dị thể: từ khâu thiết lập mô hình quang/mô hình điện, cách khai báo tiếp xúc Schottky tại biên TCO/bán dẫn, đến phương pháp trích xuất và đối chiếu giản đồ vùng năng lượng, phân bố mật độ hạt tải, phổ $QE$ và đường đặc trưng $J\text{ - }V$.
Câu hỏi thường gặp
1. Tại sao cấu trúc Emitter phía sau (RE) lại giảm được tổn thất quang học so với cấu trúc Emitter phía trước (FE)? Trong cấu trúc FE thông thường, lớp emitter a-Si:H(p) nằm ở mặt đón sáng; vì a-Si:H(p) có độ dẫn kém nên nếu tăng độ dày hoặc nồng độ pha tạp sẽ gây tái hợp mạnh và hấp thụ ký sinh ánh sáng. Ngược lại, ở cấu trúc RE, lớp emitter a-Si:H(p) được đưa về mặt sau; lớp mặt trước là a-Si:H(n) có hệ số hấp thụ quang học thấp hơn, cho phép tỷ lệ ánh sáng truyền qua đi vào lớp hấp thụ c-Si(n) cao hơn đáng kể.
2. Công thoát điện tử của lớp TCO ($\Phi_{TCO}$) tác động như thế nào đến rào thế năng lượng của tiếp xúc mặt trước? Khi $\Phi_{TCO}$ tăng từ 3,7 eV lên 4,3 eV, giản đồ vùng năng lượng bị uốn lên trên, thế tiếp xúc Schottky $V_D$ tại mặt phân cách TCO/a-Si:H(n+) tăng thêm 0,3 V. Điều này làm tăng rào thế $\Delta E_c$ cản trở điện tử thoát ra cực trước (mật độ điện tử giảm từ $3,5 \times 10^{19}\text{ cm}^{-3}$ xuống $2 \times 10^{14}\text{ cm}^{-3}$), đồng thời hạ thấp rào thế $\Delta E_v$ khiến lỗ trống tràn về mặt phân cách a-Si:H(n+)/c-Si(n) (đạt $5,7 \times 10^{10}\text{ cm}^{-3}$), thúc đẩy tái hợp cặp điện tử - lỗ trống và làm suy giảm $V_{oc}$.
3. Vì sao phổ hiệu suất lượng tử ($QE$) chỉ nhạy cảm với sự thay đổi của $\Phi_{TCO}$ ở bước sóng ngắn (300 - 700 nm)? Photon có bước sóng ngắn (năng lượng cao) bị hấp thụ chủ yếu ở các lớp bán dẫn mỏng sát bề mặt mặt trước (lớp TCO và a-Si:H(n+)), do đó sự thay đổi rào thế và mật độ hạt tải tại tiếp xúc mặt trước do $\Phi_{TCO}$ gây ra sẽ trực tiếp làm thay đổi $QE$ ở vùng này. Ngược lại, photon bước sóng dài ($\lambda \ge 700\text{ nm}$) có độ xuyên sâu lớn, được hấp thụ sâu bên trong khối c-Si(n), nên không chịu ảnh hưởng bởi biến thiên điện trường tại mặt tiếp xúc TCO mặt trước.
4. Khóa luận sử dụng các phương trình vật lý cơ bản nào trong mô hình tính toán của AFORS-HET? Phần mềm giải phương trình Poisson một chiều để xác định phân bố điện thế $\psi(x)$ theo mật độ điện tích không gian (điện tử $n'$, lỗ trống $p$, ion $N_A^-$, $N_D^+$); giải các phương trình dòng liên tục kết hợp mô hình tái hợp Shockley-Read-Hall (SRH) và tái hợp Auger; đồng thời áp dụng định luật hấp thụ Lambert-Beer để tính toán tốc độ sinh cặp hạt tải $G(x,t)$ dưới tác xạ của thông lượng photon $\Phi(\lambda)$ từ phổ AM1.5G.
5. Hiệu suất chuyển đổi quang điện cao nhất đạt được trong nghiên cứu này là bao nhiêu và vượt trội ra sao so với các nghiên cứu trước? Bằng việc tối ưu hóa phối hợp giữa $\Phi_{TCO} = 3,7\text{ eV}$, lớp a-Si:H(n+) và lớp c-Si(n), hiệu suất chuyển đổi của pin đạt mức cực đại 27,33%. Kết quả mô phỏng này cao hơn các công bố thực nghiệm trước đó của Bivour (2012, 22,8%), Watahiki (2015, 23,43%), Sen Li (2018, 21,71%), Duy Phong Pham (2019, 21,8%) và nghiên cứu mô phỏng của Park (2019, 24,01% - 25,35%).
Kết luận
Khóa luận tốt nghiệp của sinh viên Hồ Công Oanh Kiều đã thực hiện thành công việc mô phỏng và phân tích toàn diện ảnh hưởng của các lớp vật liệu TCO, a-Si:H(n+) và c-Si(n) lên pin mặt trời dị thể vô định hình/tinh thể Silic cấu trúc emitter phía sau. Nghiên cứu đã làm sáng tỏ cơ chế vật lý chi phối sự suy giảm quang điện thông qua giản đồ vùng năng lượng, thế tiếp xúc Schottky và phổ hiệu suất lượng tử ở dải bước sóng ngắn. Với việc xác lập bộ thông số tối ưu mang lại hiệu suất chuyển đổi tiềm năng lên tới 27,33%, công trình là tài liệu học thuật giá trị, cung cấp cơ sở dữ liệu quan trọng phục vụ công tác chế tạo thực nghiệm pin quang điện dị thể hiệu suất cao.