Tổng quan về luận án

Nhu cầu chuyển dịch năng lượng và khai thác nguồn nhiệt thải công nghiệp đặt ra yêu cầu cấp thiết về việc phát triển các vật liệu nhiệt điện (thermoelectric - TE) hiệu suất cao, thân thiện với môi trường và bền vững ở nhiệt độ trung - cao. Luận án tiến sĩ chuyên ngành Khoa học Vật liệu (Mã số: 62440122) của nghiên cứu sinh Phạm Thanh Tuấn Anh, dưới sự hướng dẫn của PGS. TS. Trần Cao Vinh và GS. TS. Phan Bách Thắng tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia TP. Hồ Chí Minh (VNU-HCM), mang tiêu đề "Studying thermoelectric properties of Ga and In co-doped ZnO thin films". Công trình đại diện cho bước đột phá tiên phong trong việc giải quyết mâu thuẫn nhiệt điện cốt lõi giữa hệ số Seebeck ($S$), độ dẫn điện ($\sigma$) và độ dẫn nhiệt ($\kappa_{\text{tot}}$) trên hệ vật liệu nền oxit kẽm kẽm đồng pha tạp Indium và Gallium (IGZO).

Khoảng trống nghiên cứu (research gap) lớn nhất trong lĩnh vực nhiệt điện oxit là độ dẫn nhiệt mạng tinh thể ($\kappa_{\text{lat}}$) của ZnO dạng khối truyền thống rất cao ($\sim 50\text{ W/m}\cdot\text{K}$ ở nhiệt độ phòng do liên kết ion-cộng hóa trị mạnh và dao động mạng điều hòa) và hệ số công suất ($PF = S^2\sigma$) bị giới hạn bởi nồng độ hạt tải ($n_e$) và độ suy biến vùng đơn ($N_v = 1$). Mặc dù các kỹ thuật đơn pha tạp (Al, Ga hoặc In) đã cải thiện độ dẫn điện, nhưng giới hạn hòa tan rắn (solubility limit) thấp thường dẫn đến hiện tượng phân tách pha thứ cấp không kiểm soát, làm suy giảm độ linh động hạt tải ($\mu_H$). Trong khi đó, các nghiên cứu màng mỏng trước đây chủ yếu tập trung vào tính chất điện quang mà bỏ qua việc định lượng chính xác độ dẫn nhiệt màng mỏng siêu nhanh hoặc thiếu sự tối ưu đồng thời cả ba đại lượng nhiệt điện.

Luận án thiết lập hệ thống 4 câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  1. Hypothesis 1 ($H_1$): Cơ chế thiêu kết một giai đoạn (one-step sintering) ở $1400^\circ\text{C}$ đối với vật liệu khối IGZO sẽ kiểm soát mật độ phân tách pha spinel $\text{Ga}_2\text{ZnO}4$, tạo ra các tâm tán xạ phonon chọn lọc giúp giảm $\kappa{\text{tot}}$ mà vẫn duy trì hệ số phẩm chất nhiệt điện $B$ cao.
  2. Hypothesis 2 ($H_2$): Chiến lược đồng pha tạp cation bất đối xứng kích thước ($\text{Ga}^{3+} \sim 0.062\text{ nm} < \text{Zn}^{2+} \sim 0.074\text{ nm} < \text{In}^{3+} \sim 0.081\text{ nm}$) kết hợp xử lý nhiệt sau lắng đọng sẽ kỹ thuật hóa khuyết tật mạng (lattice-defect engineering), giảm mạnh $\kappa_{\text{lat}}$ màng mỏng xuống dưới $1.0\text{ W/m}\cdot\text{K}$.
  3. Hypothesis 3 ($H_3$): Tinh chỉnh tỷ lệ In/Ga sẽ làm phẳng độ cong đáy vùng dẫn (band flattening), gia tăng khối lượng hiệu dụng trạng thái ($m_d^*$) và điều biến mức Fermi ($E_F$), từ đó phá vỡ quan hệ tỷ lệ nghịch Pisarenko truyền thống giữa $S$ và $n_e$.
  4. Hypothesis 4 ($H_4$): Màng mỏng IGZO cấu trúc nano lắng đọng bằng phương pháp phún xạ magnetron DC có thể đạt hệ số phẩm chất không thứ nguyên $ZT \ge 0.20$ ở dải nhiệt độ $573\text{ K}$, tiệm cận hiệu suất của các vật liệu khối nano tiên tiến nhất.

Khung lý thuyết của nghiên cứu dựa trên mô hình đơn vùng parabolic (Single Parabolic Band - SPB), lý thuyết tán xạ phonon âm học, định luật Wiedemann–Franz và tương quan năng lượng phân tán $E(k)$. Phạm vi thực nghiệm bao gồm tổng hợp hệ mẫu khối 25 g ($\text{In}{0.05}\text{Ga}{0.05}\text{Zn}{0.90}\text{O}$) qua 3 quy trình thiêu kết khác nhau, chế tạo bia phún xạ 3 inch (khối lượng 140 g) và lắng đọng màng mỏng $\text{In}x\text{Ga}{0.05-x}\text{Zn}{0.95}\text{O}$ ($x = 0 \div 1.5\text{ at}%$) với dải khảo sát nhiệt độ rộng từ $300\text{ K}$ đến $673\text{ K}$.

Tác động định lượng mang tính đột phá của luận án được khẳng định qua việc màng mỏng $\text{In}{0.5}\text{Ga}{4.5}\text{Zn}_{95}\text{O}$ đạt hệ số công suất kỷ lục $PF = 745.2\ \mu\text{W/m}\cdot\text{K}^2$ tại $573\text{ K}$, độ dẫn nhiệt giảm sâu xuống $0.95\text{ W/m}\cdot\text{K}$, dẫn đến giá trị $ZT = 0.20$. Đây là minh chứng thực nghiệm xuất sắc cho việc chế tạo màng nhiệt điện oxit hiệu năng cao phục vụ vi thiết bị tự cấp nguồn (self-powered IoT) và cảm biến đeo không độc hại.

Literature Review và Positioning

Bản đồ nghiên cứu nhiệt điện oxit kẽm thế giới phân hóa thành hai luồng học thuật chính: tối ưu hóa cấu trúc nano để tán xạ phonon và kỹ thuật hóa vùng năng lượng (band engineering) để tăng hệ số công suất.

Luồng thứ nhất tập trung vào việc giảm độ dẫn nhiệt mạng ($\kappa_{\text{lat}}$) thông qua cấu trúc nano đa cấp và tạo khuyết tật điểm. Han và cộng sự (2012) khi nghiên cứu vật liệu khối $\text{Zn}{0.98}\text{Al}{0.02}\text{O}$ đã chứng minh việc chuyển đổi hình thái hạt từ thanh (rods), phiến (platelets) sang hạt nano (nanoparticles) giúp giảm $\kappa_{\text{tot}}$ từ $36\text{ W/m}\cdot\text{K}$ xuống $8.5\text{ W/m}\cdot\text{K}$ ở $373\text{ K}$ và đạt $3.21\text{ W/m}\cdot\text{K}$ ở $1223\text{ K}$. Tương tự, Takemoto và đồng nghiệp (2006) chỉ ra rằng việc đồng pha tạp Ga và In tạo ra các lỗi xếp chồng (stacking faults) trong mạng lưới ZnO, giảm $\kappa_{\text{tot}}$ khối xuống $1.2\text{ W/m}\cdot\text{K}$. Tuy nhiên, hạn chế cố hữu của luồng nghiên cứu này là sự gia tăng ranh giới hạt và pha thứ cấp thường gây tán xạ mạnh điện tử, làm suy giảm nghiêm trọng độ linh động Hall ($\mu_H$).

Luồng thứ hai tập trung vào việc gia tăng nồng độ hạt tải ($n_e$) và khối lượng hiệu dụng ($m_d^$). Jood và cộng sự (2011) đạt được mức cải thiện độ dẫn điện gấp 1000 lần trên mẫu khối ZnO pha tạp 0.25 at% In, nâng $PF$ lên $\sim 90\ \mu\text{W/m}\cdot\text{K}^2$ ở nhiệt độ phòng. Jung và các cộng sự (2017) đã tăng cường $m_d^$ trong vật liệu khối $\text{Zn}{0.98}\text{Ga}{0.02}\text{O}$ nhờ hiện tượng nén hóa học (chemical compression) do độ hòa tan của Ga mở rộng, đạt $PF = 1250\ \mu\text{W/m}\cdot\text{K}^2$ tại $1273\text{ K}$.

Công trình quốc tế Hệ vật liệu & Cấu trúc Độ dẫn nhiệt $\kappa_{\text{tot}}\ (\text{W/m}\cdot\text{K})$ Hệ số công suất $PF\ (\mu\text{W/m}\cdot\text{K}^2)$ Hệ số phẩm chất $ZT$ Hạn chế chính
Teehan et al. (2017) Màng mỏng $(\text{ZnO})\text{Al}_{0.02}$ Không công bố $2210$ (tại $975\text{ K}$) Không xác định Không kiểm soát và định lượng được $\kappa_{\text{tot}}$
Seo et al. (2018) Siêu mạng $\text{InGaO}_3(\text{ZnO})_m$ / Sapphire $7.00$ (tại $300\text{ K}$) $80$ (tại $375\text{ K}$) Rất thấp ($< 0.05$) Độ kết tinh kém do sai mạng nền; $PF$ thấp
Zheng et al. (2019) Màng mỏng In-doped AZO Không công bố $222$ (tại $300\text{ K}$) Không xác định Thiếu dữ liệu phân tích $ZT$ và phonon
Liu et al. (2020) Màng nano xốp AZO/CNT $1.16$ (tại $300\text{ K}$) $120$ (tại $300\text{ K}$) $0.10$ Quy trình phức tạp, khó mở rộng quy mô
Luận án này (2024) Màng mỏng $\text{In}{0.5}\text{Ga}{4.5}\text{Zn}_{95}\text{O}$ (IGZO) $0.95$ (tại $573\text{ K}$) $745.2$ (tại $573\text{ K}$) $0.20$ (tại $573\text{ K}$) Tối ưu hóa đồng thời $S$, $\sigma$, $\kappa$ qua biến điệu cấu trúc vùng

Tranh luận học thuật cốt lõi tồn tại giữa quan điểm cho rằng "pha thứ cấp và ranh giới hạt nano luôn gây hại cho vận chuyển điện tử" và quan điểm đối lập khẳng định "phân tách pha có kiểm soát là chìa khóa duy nhất để triệt tiêu độ dẫn nhiệt mạng". Luận án đã định vị chính xác vị thế khoa học của mình bằng cách chứng minh: ở quy mô vật liệu khối, thiêu kết một giai đoạn trực tiếp ở $1400^\circ\text{C}$ tạo pha spinel phân tán nano giúp tối ưu hóa hệ số phẩm chất $B$; đồng thời khi chuyển giao công nghệ sang dạng màng mỏng phún xạ magnetron, tính đồng nhất pha rắn được khôi phục hoàn toàn, cho phép điều biến cấu trúc vùng năng lượng (band structure engineering) mà không làm suy giảm độ dẫn điện.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và thách thức các giới hạn của lý thuyết vận chuyển điện tử và phonon trong chất bán dẫn suy biến:

  1. Mở rộng lý thuyết Pisarenko và mô hình SPB (Pisarenko, 1953; Snyder et al., 2008): Mô hình đơn vùng parabolic cổ điển thiết lập mối quan hệ tuyến tính nghịch đảo giữa $S$ và nồng độ hạt tải $n_e$: $$S = \frac{8\pi^2 k_B^2 T}{3 e h^2} m_d^* \left(\frac{\pi}{3 n_e}\right)^{2/3}$$ Trong hệ ZnO có độ suy biến thung lũng $N_v = 1$, việc tăng $n_e$ luôn làm suy giảm $S$. Luận án chứng minh rằng thông qua việc đồng pha tạp In và Ga, sự đóng góp của các orbital định vị bất đối xứng tạo nên hiện tượng làm phẳng đáy vùng dẫn (band flattening), trực tiếp nâng cao khối lượng hiệu dụng vùng trung bình $m_b^$, từ đó gia tăng $m_d^ = N_v^{2/3} m_b^*$ độc lập với $n_e$.
           MÔ HÌNH VÙNG NĂNG LƯỢNG TRONG MÀNG MỎNG IGZO
           
       Undoped ZnO (Vùng dốc)            IGZO 0.5 (Làm phẳng đáy vùng dẫn)
              E                                    E
              |       /                            |    \________/
              |      /                             |    |        |
              |     /  m_b* nhỏ                    |    | m_b* lớn (Band flattening)
              |    /   S thấp                      |    | S tăng cao
              |___/________ k                      |____|________ k
                 k=0                                  k=0
  1. Phát triển mô hình độ linh động trọng số (Weighted Mobility $\mu_w$) và Hệ số phẩm chất $B$ (Goldsmid & Sharp, 1999; Snyder et al., 2020): Luận án áp dụng hệ phương thức trích xuất $\mu_w$ trực tiếp từ dữ liệu thực nghiệm $S$ và $\sigma$ mà không phụ thuộc vào phép đo Hall phức tạp: $$\mu_w = \frac{3 h^3 \sigma}{8 \pi e (2 m_0 k_B T)^{3/2}} \left[ \frac{\exp\left(\frac{|S|}{k_B/e} - 2\right)}{1 + \exp\left[-5\left(\frac{|S|}{k_B/e} - 1\right)\right]} + \frac{\frac{3}{\pi^2}\frac{|S|}{k_B/e}}{1 + \exp\left[5\left(\frac{|S|}{k_B/e} - 1\right)\right]} \right]$$ Từ đó, hệ số phẩm chất nhiệt điện dưới cơ chế tán xạ phonon âm học $B = \left(\frac{k_B}{e}\right)^2 \frac{2 e (2 m_0 k_B T)^{3/2}}{3 \pi \hbar^3} \frac{\mu_w}{\kappa_{\text{lat}}}$ được định lượng chính xác, chứng minh quy trình thiêu kết một bước đạt giá trị $\mu_w/\kappa_{\text{lat}}$ cao nhất.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết nền tảng:

  • Lý thuyết khuyết tật điểm bán dẫn (Kröger–Vink notation): Phân tích tương tác giữa các tâm donor thay thế ($\text{Ga}{\text{Zn}}^{\bullet}$, $\text{In}{\text{Zn}}^{\bullet}$), khuyết tật xen kẽ ($\text{Zn}i^{\bullet\bullet}$, $\text{Ga}i^{\bullet\bullet\bullet}$) và vacancy oxy ($V{\text{O}}^{\bullet\bullet}$, $V{\text{O}}^{\times}$).
  • Định luật Wiedemann–Franz: Tách biệt chính xác thành phần điện tử $\kappa_e = L\sigma T$ (với số Lorentz $L = 2.45 \times 10^{-8}\ \text{W}\cdot\Omega/\text{K}^2$) và thành phần mạng $\kappa_{\text{lat}} = \kappa_{\text{tot}} - \kappa_e$.
  • Lý thuyết biến điệu độ cong vùng năng lượng (Energy Dispersion Curvature): $$m_b^* = \hbar^2 \left( \frac{\partial^2 E(k)}{\partial k^2} \right)^{-1}$$ Việc pha tạp In tạo ra bán kính bất đối xứng lớn ($\text{In-O} \approx 0.230\text{ nm}$ so với $\text{Zn-O} \approx 0.193\text{ nm}$), gây biến dạng mạng cục bộ (local lattice strain), làm giảm đạo hàm bậc hai $\partial^2 E / \partial k^2$, dẫn tới tăng vọt $m_b^*$.

Điều kiện biên lý thuyết (boundary conditions) được xác định rõ: Khung phân tích giả định cơ chế tán xạ hạt tải chiếm ưu thế ở nhiệt độ cao là tán xạ phonon âm học và tán xạ tạp chất ion hóa, trong điều kiện bán dẫn suy biến nồng độ hạt tải đạt ngưỡng $10^{19} \div 10^{21}\text{ cm}^{-3}$.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng nghiêm ngặt (positivism), sử dụng phương pháp tiếp cận từ trên xuống (top-down dimensional reduction) kết hợp đa phương thức (multi-method experimental design). Thiết kế nghiên cứu chia làm hai cấp độ liên kết chặt chẽ:

  • Cấp độ 1 (Vật liệu khối): Nghiên cứu cơ chế hòa tan rắn, tách pha spinel $\text{Ga}_2\text{ZnO}_4$ và thiết lập điều kiện thiêu kết tối ưu cho bia phún xạ.
  • Cấp độ 2 (Màng mỏng 2D): Kế thừa thành phần tối ưu từ khối để phún xạ màng mỏng, khảo sát ảnh hưởng của hiệu ứng giảm chiều không gian và kỹ thuật hóa cấu trúc vùng.
       QUY TRÌNH NGHIÊN CỨU VÀ TỔNG HỢP VẬT LIỆU TOÀN DIỆN
       
 [Bột oxit: ZnO (3N), Ga2O3 (4N), In2O3 (4N)]
                      |
                      v (Nghiền bi hành tinh 5h + Sấy chân không 120°C)
       +--------------+--------------+
       |                             |
       v                             v
[Quy trình 1: Ép 14 MPa]    [Quy trình 2 & 3: Nung sơ bộ 1000°C]
       |                             |
       v                             v
[Thiêu kết trực tiếp 1400°C]   [Nghiền/Ép lại -> Thiêu kết 1400°C]
       |                             |
  (IGZO-1 Bulk: Tối ưu)       (IGZO-2, IGZO-3 Bulks)
       |
       +---> [Chế tạo bia phún xạ 3-inch (D > 90%)]
                      |
                      v (Phún xạ DC Magnetron - Leybold Univex-450)
            [Màng mỏng IGZO x = 0 - 1.5 at%]
                      |
                      v (Ủ nhiệt 400°C - 500°C trong chân không/khí quyển)
            [Đo đạc TDTR, ZEM-3, LSR-3, HMS-5500, XPS, XRD, TEM, PL]

Quy trình nghiên cứu rigorous

  1. Tổng hợp vật liệu khối: Bột ZnO (độ tinh khiết 99.9%, Merck), $\text{Ga}_2\text{O}_3$ (99.99%, Sigma-Aldrich) và $\text{In}_2\text{O}_3$ (99.99%, Sigma-Aldrich) được cân bằng cân phân tích điện tử Sartorius BP211D (độ chính xác 0.01 mg). Hỗn hợp được nghiền bi hành tinh (Ceramic Instruments srl) trong 5 giờ với bi gốm $\text{Al}_2\text{O}_3$ (đường kính 10 mm và 18.5 mm), sấy khô ở $120^\circ\text{C}$ trong 24 giờ. Ba quy trình thiêu kết được áp dụng:
    • Quy trình 1 (IGZO-1): Ép thủy lực 14 MPa tạo viên $30 \times 30\text{ mm}^2$, thiêu kết trực tiếp ở $1400^\circ\text{C}$ trong 3 giờ (tốc độ gia nhiệt $2.5^\circ\text{C/phút}$).
    • Quy trình 2 (IGZO-2): Ép 14 MPa, nung sơ bộ (calcination) $1000^\circ\text{C}$ trong 3 giờ, nghiền mịn bằng cối mã não, ép lại và thiêu kết ở $1400^\circ\text{C}$.
    • Quy trình 3 (IGZO-3): Nung sơ bộ bột ở $1000^\circ\text{C}$, ép viên và thiêu kết ở $1400^\circ\text{C}$.
  2. Chế tạo màng mỏng: Chế tạo bia phún xạ đường kính 75 mm, độ dày 5 mm (đạt khối lượng riêng tương đối $> 90%$ so với mật độ lý thuyết $5.606\text{ g/cm}^3$). Màng mỏng được lắng đọng bằng hệ phún xạ magnetron DC thương mại Leybold Univex-450 trong môi trường khí Ar độ tinh khiết cao ($5N$), áp suất cơ bản $< 10^{-5}\text{ mbar}$. Quá trình ủ nhiệt sau lắng đọng được thực hiện ở $400^\circ\text{C}$ và $500^\circ\text{C}$.
  3. Triangulation & Độ tin cậy thực nghiệm: Kiểm chứng chéo cấu trúc tinh thể qua XRD (nhiễu xạ tia X góc rộng và góc quét hẹp), TEM/SAED và HR-TEM. Phân tích trạng thái liên kết hóa học qua phổ quang điện tử tia X (XPS) độ phân giải cao cho các mức năng lượng lõi Zn 2p, Ga 2p, In 3d, O 1s và đường Auger Zn $L_3M_{4,5}M_{4,5}$. Độ phát quang quang học được đo bằng hệ PL ở hai bước sóng kích thích 266 nm và 355 nm.

Data và phân tích

  • Hệ số Seebeck và Điện trở suất: Đo đồng thời bằng hệ thống chuyên dụng thương mại Linseis LSR-3 (tại Viện INOMAR) và ULVAC ZEM-3 (tại Viện KERI, Hàn Quốc) từ $300\text{ K}$ đến $673\text{ K}$.
  • Thông số hạt tải Hall: Hệ thống Ecopia HMS-3000 và HMS-5500 tích hợp bộ gá mẫu gia nhiệt HTSST3 (đo từ nhiệt độ phòng đến $300^\circ\text{C}$) theo cấu hình Van der Pauw trên mẫu vuông có tiếp xúc In 4 góc.
  • Đo độ dẫn nhiệt màng mỏng: Đây là thách thức thực nghiệm lớn nhất do màng mỏng có độ trong suốt cao và độ khuếch tán nhiệt siêu nhanh, vượt quá giới hạn đo của phương pháp phân tích chớp laser (Laser Flash Analysis - LFA) thông thường (như Netzsch LFA-457 hay Linseis LFA500/1000). Luận án đã áp dụng kỹ thuật phản xạ nhiệt miền thời gian (Time-Domain Thermoreflectance - TDTR) và hệ Linseis TF-LFA chuyên dụng, sử dụng hai chùm xung laser pump-probe tách biệt với độ phân giải thời gian pico-giây, cho phép xác định chính xác độ khuếch tán nhiệt $\alpha$ và độ dẫn nhiệt $\kappa_{\text{tot}} = \alpha \cdot C_p \cdot \rho$ (với nhiệt dung riêng $C_p$ đo từ DSC Netzsch 204 F1 và khối lượng riêng $\rho$ từ cân tỉ trọng Sartorius BSA224S-CW).

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Khám phá vai trò của phân tách pha spinel $\text{Ga}_2\text{ZnO}_4$ trong mẫu khối IGZO-1: Quy trình 1 (thiêu kết trực tiếp ở $1400^\circ\text{C}$) tạo ra sự phân tách pha spinel nano đồng đều, đóng vai trò như các rào cản phân tán phonon hiệu quả. Kết quả thực nghiệm cho thấy mẫu IGZO-1 đạt độ dẫn nhiệt thấp nhất và hệ số công suất cao nhất ở nhiệt độ cao, đồng thời duy trì độ linh động trọng số $\mu_w$ và hệ số phẩm chất $B$ vượt trội so với các quy trình hai bước (IGZO-2 và IGZO-3).

  2. Kỹ thuật hóa khuyết tật mạng triệt tiêu $\kappa_{\text{tot}}$ màng mỏng xuống $0.95\text{ W/m}\cdot\text{K}$: Trích dẫn trực tiếp từ văn bản luận án: "The combination of the dual-doping effect and post-thermal treatment is found to reduce significantly the thermal conductivity of the IGZO thin films to 0.95 W/mK, leading to the highest figure of merit ZT = 0.20." Sự kết hợp giữa biến dạng mạng do sai khác bán kính ion ($\text{In}^{3+}/\text{Ga}^{3+}$) và ủ nhiệt ở $500^\circ\text{C}$ đã tối ưu hóa mật độ vacancy oxy ($V_{\text{O}}$) và triệt tiêu $\kappa_{\text{lat}}$ tới mức cực hạn của giới hạn vô định hình (amorphous limit).

  3. Hiện tượng làm phẳng đáy vùng dẫn (Conduction Band Flattening) và biến điệu mức Fermi: Trích dẫn trực tiếp từ văn bản luận án: "Tuning In and Ga contents is found to coupling modification of Fermi level and band flattening, leading to the ZT of 0.2 and PF of 745.2 uW/mK2 at 573 K for the IGZO films deposited from the In0.5Ga4.5Zn95O target." Phổ XPS vùng hóa trị và đồ thị Tauc chỉ ra rằng mẫu $\text{In}{0.5}\text{Ga}{4.5}\text{Zn}_{95}\text{O}$ tạo ra sự dịch chuyển mức Fermi $E_F$ vào sâu trong vùng dẫn đồng thời làm tăng $m_d^*$, giúp $PF$ tăng vọt lên $745.2\ \mu\text{W/m}\cdot\text{K}^2$ ở $573\text{ K}$, cao gấp hơn 3.5 lần so với các màng mỏng ZnO đơn pha tạp thông thường ($\sim 200\ \mu\text{W/m}\cdot\text{K}^2$).

  4. Phân tích trạng thái hóa học và cơ chế phát quang PL: Phổ PL ở bước sóng kích thích 266 nm và 355 nm kết hợp phân rã đường phổ XPS O 1s chứng minh tỷ lệ diện tích đỉnh liên kết $\text{O}{\text{vacancy}}$/$\text{O}{\text{lattice}}$ đạt giá trị tối ưu ở mẫu ủ $500^\circ\text{C}$, giải phóng nồng độ điện tử tự do trong khoảng lý tưởng $10^{20}\text{ cm}^{-3}$ mà không tạo ra tán xạ tạp chất trung hòa quá mức.

+---------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP CÁC PHÁT HIỆN THỰC NGHIỆM ĐỘT PHÁ              |
+===========================================================================+
| 1. Khối IGZO-1: Thiêu kết trực tiếp 1400°C tạo pha spinel Ga2ZnO4 nano    |
|    -> Giảm mạnh độ dẫn nhiệt mạng, tối đa hóa hệ số phẩm chất B.         |
+---------------------------------------------------------------------------+
| 2. Màng mỏng IGZO 0.5@500: Đạt độ dẫn nhiệt cực tiểu κ_tot = 0.95 W/m·K   |
|    -> Phá vỡ giới hạn truyền nhiệt tinh thể wurtzite thông thường.        |
+---------------------------------------------------------------------------+
| 3. Hiện tượng Band Flattening: Nâng cao m_d* và điều chỉnh vị trí E_F     |
|    -> Đạt PF kỷ lục = 745.2 µW/m·K² tại 573 K.                            |
+---------------------------------------------------------------------------+
| 4. Hệ số phẩm chất ZT = 0.20 tại 573 K                                    |
|    -> Tiệm cận đẳng cấp vật liệu khối nano tiên tiến nhất thế giới.       |
+---------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm rõ ràng nhất trên hệ màng mỏng oxit chứng minh rằng việc biến điệu độ cong vùng dẫn ($m_b^*$) thông qua pha tạp cation bất đối xứng bán kính ion có thể phá vỡ thế lưỡng nan đánh đổi (trade-off) giữa $S$ và $\sigma$.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn mực trong việc đo độ dẫn nhiệt màng mỏng oxit bán dẫn trong suốt bằng kỹ thuật TDTR và TF-LFA có độ nhạy pico-giây, giải quyết triệt để sự thiếu hụt dữ liệu $\kappa_{\text{tot}}$ trong các y văn quốc tế suốt nhiều thập kỷ.
  • Về mặt ứng dụng thực tiễn: Mở ra lộ trình công nghệ khả thi để chế tạo các module nhiệt điện màng mỏng (Thin-Film Thermoelectric Generators - TFTEGs) hoạt động ở công suất từ nW đến mW, trực tiếp tích hợp vào chip vi điện tử, cảm biến IoT và thiết bị đeo y tế thu hồi nhiệt cơ thể người.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 4 hạn chế mang tính ranh giới kỹ thuật:

  1. Giới hạn pha tạp Indium: Khi hàm lượng In vượt quá 0.5 at% (mẫu $\text{In}{1.0}$ và $\text{In}{1.5}$), ứng suất dư ($\epsilon$) và tán xạ ranh giới hạt tăng đột biến, làm suy giảm độ linh động Hall $\mu_H$ và kéo theo sự sụt giảm của hệ số công suất $PF$.
  2. Nhiệt độ đo giới hạn: Các phép đo màng mỏng hiện mới khảo sát tối đa đến $573\text{ K} - 673\text{ K}$ do giới hạn chịu nhiệt của đế thủy tinh/silicon và hệ gá tiếp xúc điện cực In, chưa đánh giá hết tiềm năng của vật liệu ở vùng nhiệt độ rất cao ($> 900\text{ K}$).
  3. Độ ổn định oxy hóa dài hạn: Sự suy thoái tiềm tàng của các khuyết tật vacancy oxy trong môi trường không khí ở nhiệt độ cao sau hàng ngàn giờ vận hành liên tục chưa được khảo sát bằng các thử nghiệm lão hóa gia tốc (accelerated aging tests).
  4. Hiệu ứng tiếp xúc tiếp giáp (Interface resistance): Luận án tập trung vào tính chất nội tại của màng mỏng mà chưa đánh giá điện trở tiếp xúc giữa màng IGZO và các lớp điện cực kim loại trong module hoàn chỉnh.

Chương trình nghiên cứu 5-10 năm tới mở ra các hướng đi cụ thể:

  • Hướng 1: Phát triển kỹ thuật tính toán phiếm hàm mật độ (DFT) kết hợp học máy (Machine Learning) để mô phỏng chính xác cấu trúc vi mô vùng dẫn và mặt tán xạ phonon dưới tác động của các cặp đồng pha tạp nhóm III-A khác (như Al-In, Ga-Sc).
  • Hướng 2: Chế tạo thiết bị nhiệt điện màng mỏng đa tầng (multi-junction thin-film TE module) kết hợp lớp p-type tương thích (như $\text{CuCrO}_2$ hoặc $\text{SnO}$) trên đế dẻo (flexible substrates) chịu nhiệt cao như polyimide hoặc mica.
  • Hướng 3: Tích hợp hệ màng mỏng IGZO vào các cấu trúc lai ghép nhiệt điện - quang điện (hybrid Thermoelectric-Photovoltaic devices) nhằm thu hồi nhiệt tổn hao từ pin mặt trời.

Tác động và ảnh hưởng

Luận án của NCS. Phạm Thanh Tuấn Anh mang lại giá trị học thuật và ứng dụng sâu rộng:

  • Ảnh hưởng học thuật: Kết quả nghiên cứu đã được công bố tại các diễn đàn khoa học uy tín hàng đầu như Hội nghị Nhiệt điện Quốc tế (ICT), Hội nghị Nhiệt điện Châu Á (ACT), Hội nghị Nhiệt điện Đông Nam Á (SACT) và các tạp chí chuyên ngành vật liệu ISI/Scopus. Công trình cung cấp hệ thống dữ liệu toàn diện về thông số tán xạ, $m_d^*$, $\mu_w$, và độ dẫn nhiệt TDTR, ước tính sẽ thu hút số lượng trích dẫn học thuật cao trong cộng đồng nhiệt điện oxit thế giới.
  • Chuyển đổi công nghiệp & Năng lượng xanh: Đóng góp trực tiếp vào chiến lược phát triển năng lượng quốc gia Việt Nam đến năm 2030, tầm nhìn đến năm 2045 (Nghị quyết số 55-NQ/TW của Bộ Chính trị ban hành ngày 11/02/2020) về việc nâng cao hiệu quả sử dụng năng lượng và phát triển công nghệ tái tạo.
  • Lợi ích môi trường và xã hội: Thay thế hoàn toàn các vật liệu nhiệt điện truyền thống chứa nguyên tố độc hại, quý hiếm hoặc đắt đỏ (như $\text{Bi}_2\text{Te}_3$, $\text{PbTe}$, $\text{CoSb}_3$) bằng hệ oxit kẽm thân thiện sinh học, trữ lượng dồi dào trong vỏ trái đất, giá thành thấp và an toàn tuyệt đối với sức khỏe con người.

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

  1. Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận phương pháp luận mẫu mực về phân tích nhiệt điện hiện đại, công thức tính toán độ linh động trọng số $\mu_w$, trích xuất $m_d^*$ theo SPB và kỹ thuật xử lý phổ XPS/PL/TDTR phân giải cao.
  2. Chuyên gia & Giáo sư Vật liệu: Thụ hưởng khung lý thuyết hoàn chỉnh về cơ chế làm phẳng vùng dẫn và tương tác khuyết tật mạng trong hệ oxit kẽm đồng pha tạp.
  3. Kỹ sư R&D Công nghệ Bán dẫn & Vi cơ điện tử (MEMS): Sở hữu thông số công nghệ thực chứng (áp suất phún xạ, tỷ lệ pha tạp, chế độ ủ nhiệt) để chế tạo màng nhiệt điện tích hợp trong vi mạch tự cấp nguồn, hệ thống IoT công nghiệp và cảm biến y tế đeo tay.
  4. Các nhà hoạch định chính sách môi trường & Quỹ khoa học: Có cơ sở khoa học vững chắc từ nguồn tài trợ của Bộ Khoa học & Công nghệ, Quỹ NAFOSTED, ĐHQG-HCM, Quỹ Đổi mới sáng tạo Vingroup (VINIF) và Học bổng Vallet để tiếp tục nhân rộng các dự án năng lượng xanh bền vững.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc chứng minh thực nghiệm thành công cơ chế đồng biến điệu vị trí mức Fermi và làm phẳng đáy vùng dẫn (Fermi-level & Band Flattening Co-engineering) trong màng mỏng bán dẫn oxit kẽm suy biến. Luận án đã mở rộng mô hình đơn vùng parabolic (Single Parabolic Band - SPB) và tương quan Pisarenko bằng cách chỉ ra rằng việc đồng pha tạp cặp cation bất đối xứng kích thước ion ($\text{In}^{3+}/\text{Ga}^{3+}$) cho phép gia tăng khối lượng hiệu dụng trạng thái $m_d^$ từ việc tăng khối lượng dải trung bình $m_b^$, giúp hệ số Seebeck $S$ duy trì ở mức cao ngay cả khi nồng độ hạt tải $n_e$ tăng mạnh.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với nghiên cứu của Teehan et al. (2017) và Zheng et al. (2019) – những công trình chỉ đo được hệ số công suất $PF$ mà hoàn toàn bế tắc trong việc xác định độ dẫn nhiệt $\kappa_{\text{tot}}$, và nghiên cứu của Seo et al. (2018) trên cấu trúc siêu mạng phức tạp nhưng hiệu suất $ZT$ rất thấp – luận án đã thiết lập quy trình đo đạc toàn diện và chính xác bậc nhất bằng phương pháp phản xạ nhiệt miền thời gian (TDTR) và hệ TF-LFA. Điều này cho phép phân tách rạch ròi thành phần $\kappa_e$ và $\kappa_{\text{lat}}$, lượng hóa chính xác $ZT$ màng mỏng mà không dựa trên bất kỳ phép ngoại suy lý thuyết thiếu tin cậy nào.

3. Phát hiện thực nghiệm nào bất ngờ nhất và lời giải thích vật lý?

Phát hiện bất ngờ nhất là ở quy trình thiêu kết khối IGZO: mẫu thiêu kết trực tiếp một bước ở $1400^\circ\text{C}$ (IGZO-1) lại cho hiệu suất nhiệt điện và hệ số phẩm chất $B$ cao vượt trội so với các quy trình hai bước phức tạp có nung sơ bộ (IGZO-2 và IGZO-3). Về mặt vật lý, quá trình thiêu kết trực tiếp ở nhiệt độ cao kích thích sự phân tách pha spinel $\text{Ga}_2\text{ZnO}4$ phân tán kích thước nanomet đồng đều, vừa tạo ranh giới tán xạ mạnh phonon dải tần trung - cao để giảm $\kappa{\text{lat}}$, vừa hạn chế bẫy hạt tải ở các ranh giới hạt kích thước lớn, bảo toàn độ linh động trọng số $\mu_w$.

4. Luận án có cung cấp đầy đủ giao thức tái lập thực nghiệm (Replication Protocol) không?

Có, toàn bộ giao thức thực nghiệm được mô tả chi tiết với độ chính xác tuyệt đối: từ tỷ lệ phối liệu oxit theo bảng tính khối lượng mol chi tiết (Bảng 2.1 và 2.2), thông số nghiền bi (5 giờ, tỷ lệ bi gốm $\text{Al}_2\text{O}_3$), áp lực ép thủy lực ($14\text{ MPa}$), tốc độ gia nhiệt ($2.5^\circ\text{C/phút}$) trong lò Nabertherm 1550, cho đến thông số lắng đọng phún xạ magnetron DC trên hệ Leybold Univex-450 và nhiệt độ ủ chân không ($400^\circ\text{C}, 500^\circ\text{C}$).

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?

Chương trình nghiên cứu 10 năm định hướng: (1) Ứng dụng mô phỏng DFT bậc cao để tính toán cấu trúc vùng năng lượng các hệ đồng pha tạp mới; (2) Chế tạo module nhiệt điện màng mỏng hoàn chỉnh (TFTEG) công suất mW trên đế dẻo; (3) Tích hợp thiết bị thu hồi nhiệt thải vào pin mặt trời và cảm biến y sinh IoT tự cấp nguồn thương mại hóa.

Kết luận

  1. Luận án đã làm sáng tỏ toàn diện cơ chế vận chuyển điện tử và phonon trên hệ vật liệu ZnO đồng pha tạp In và Ga ở cả hai dạng cấu trúc khối (bulk) và màng mỏng (thin film).
  2. Phát hiện và chứng minh cơ chế phân tách pha spinel $\text{Ga}_2\text{ZnO}_4$ qua quy trình thiêu kết một giai đoạn ở $1400^\circ\text{C}$, tối ưu hóa hệ số phẩm chất nhiệt điện $B$ và độ dẫn nhiệt khối.
  3. Kỹ thuật hóa khuyết tật mạng thành công trên màng mỏng IGZO, triệt tiêu độ dẫn nhiệt màng mỏng xuống mức kỷ lục $0.95\text{ W/m}\cdot\text{K}$.
  4. Hiện thực hóa chiến lược làm phẳng đáy vùng dẫn (band flattening) và điều biến mức Fermi, đạt hệ số công suất $PF = 745.2\ \mu\text{W/m}\cdot\text{K}^2$ tại $573\text{ K}$.
  5. Thiết lập giá trị hệ số phẩm chất $ZT = 0.20$ cho màng mỏng $\text{In}{0.5}\text{Ga}{4.5}\text{Zn}_{95}\text{O}$, chính thức đưa màng mỏng nhiệt điện oxit kẽm bước vào nhóm vật liệu nhiệt điện nano hiệu năng cao.
  6. Mở ra 3 luồng nghiên cứu học thuật mới: Kỹ thuật biến điệu cấu trúc vùng bằng pha tạp đa nguyên tố bất đối xứng; Chuẩn hóa phương pháp đo TDTR cho màng mỏng nhiệt điện trong suốt; và Phát triển các vi thiết bị thu hồi năng lượng nhiệt tự cấp nguồn thân thiện với môi trường.