Tổng quan nghiên cứu

Năng lượng mặt trời là nguồn tài nguyên tái tạo vô tận với tổng lượng bức xạ truyền tới Trái Đất lớn gấp khoảng 3.000 lần toàn bộ nhu cầu năng lượng hiện tại của nhân loại. Theo các tính toán khoa học, việc tận dụng hiệu quả chỉ 0,1% diện tích bề mặt Trái Đất với các thiết bị pin mặt trời đạt hiệu suất chuyển đổi 10% đã đủ cung cấp năng lượng cho toàn cầu. Trong xu thế đó, vật liệu bán dẫn oxit titan (TiO2) trở thành tâm điểm nghiên cứu chế tạo pin mặt trời nhạy màu quang điện hóa (DSSC) và quang xúc tác xử lý môi trường nhờ tính trơ hóa học, độ bền nhiệt lên đến 1.870 độ C và giá thành thấp.

Tuy nhiên, rào cản kỹ thuật lớn nhất của TiO2 ở pha kết tinh anatase là độ rộng khe năng lượng lớn (khoảng 3,2 eV), khiến vật liệu chỉ hấp thụ được bức xạ tử ngoại (UV) – vốn chỉ chiếm 4% đến 5% năng lượng phổ mặt trời, trong khi bỏ phí hoàn toàn 45% năng lượng thuộc vùng ánh sáng khả kiến. Luận văn tập trung giải quyết vấn đề cốt lõi này thông qua mục tiêu nghiên cứu cấu trúc điện tử, sự thay đổi vùng cấm và cơ chế dịch chuyển mức năng lượng của vật liệu TiO2 cấu trúc nano pha anatase khi được pha tạp nguyên tử nhôm (Al) và hình thành các nút khuyết oxy. Phạm vi nghiên cứu được triển khai dựa trên tính toán mô phỏng lý thuyết lượng tử tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội. Công trình mang ý nghĩa khoa học sâu sắc, cung cấp luận cứ định lượng chính xác ở cấp độ nguyên tử nhằm hỗ trợ tối ưu hóa quy trình chế tạo thực nghiệm màng nano bán dẫn, hướng tới nâng cao hiệu suất chuyển đổi lượng tử (IPCE) và thế hở mạch cho các tế bào pin mặt trời thế hệ mới.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng cơ học lượng tử đa hạt kết hợp với phép gần đúng Born-Oppenheimer nhằm tách biệt chuyển động của electron khỏi hệ hạt nhân ion cố định. Khung lý thuyết trung tâm là Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT), bắt nguồn từ hai định lý nền tảng của Hohenberg-Kohn (1964) khẳng định năng lượng trạng thái cơ bản là một phiếm hàm đơn trị của mật độ điện tử r(r), từ đó giảm thiểu bài toán phức tạp 3N tọa độ không gian về hàm mật độ chỉ phụ thuộc 3 biến tọa độ thực.

Để giải bài toán tương tác đa hạt, phương trình tự hợp Kohn-Sham (1965) được thiết lập nhằm thay thế hệ electron thực tương tác bằng một hệ electron giả định không tương tác chuyển động trong thế hiệu dụng. Năng lượng tương quan trao đổi được mô tả chính xác thông qua các phép gần đúng mật độ địa phương (LDA), gần đúng mật độ spin địa phương (LSDA) và gần đúng gradient suy rộng (GGA) với các phiếm hàm chuẩn như PBE, PW91 và RPBE. Khung nghiên cứu vận dụng 5 khái niệm chuyên sâu: mật độ điện tử r(r), mật độ trạng thái (DOS/PDOS), độ rộng vùng cấm (Bandgap Eg), cấu trúc dải năng lượng (Band Structure) và ô mạng siêu cơ sở (Supercell).

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

Luận văn áp dụng phương pháp tính toán từ nguyên lý ban đầu (First-principles / Ab-initio) hoàn toàn không sử dụng bất kỳ tham số kinh nghiệm nào từ trước. Cỡ mẫu nghiên cứu được thiết lập trên mô hình tinh thể TiO2 pha anatase mở rộng dạng siêu ô mạng (supercell) kích thước 2 x 2 x 1 bao gồm 48 nguyên tử (16 nguyên tử Ti và 32 nguyên tử O). Phương pháp chọn mẫu cấu trúc được thực hiện thông qua mô hình thay thế cation cục bộ: thay thế 1 nguyên tử Ti bằng 1 nguyên tử Al (tương ứng công thức Ti15/16Al1/16O32 với nồng độ pha tạp danh định 6,25%), đồng thời khảo sát có hệ thống 3 cấu hình khuyết tật oxy đơn lẻ tại các vị trí tọa độ đối xứng khác nhau ký hiệu là O1, O2 và O3.

Lý do lựa chọn phương pháp DFT-GGA là khả năng tối ưu hóa hình học mạng tinh thể với độ chính xác cao, dự đoán tin cậy sự tái phân bố điện tích và cấu trúc vùng năng lượng với chi phí tính toán hợp lý. Toàn bộ quy trình mô phỏng, tối ưu hóa năng lượng tổng cộng và phân tích mật độ trạng thái điện tử được thực hiện xuyên suốt trong khung thời gian đào tạo cao học 2010 - 2011.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, kết quả tính toán trên tinh thể TiO2 pha anatase chuẩn (2 x 2 x 1) khẳng định vật liệu có cấu trúc tứ giác với hằng số mạng a = 3,78 Angstrom, c = 9,51 Angstrom, độ rộng khe năng lượng trực tiếp đạt 3,2 eV, khối lượng riêng 3,89 g/cm3, phản ánh cấu trúc bát diện tiếp xúc cạnh với mật độ điện tử tập trung chủ yếu tại các liên kết Ti-O.

Thứ hai, khi pha tạp nguyên tử Al vào vị trí nút mạng Ti (hệ Ti15/16Al1/16O32), mức tạp chất acceptor xuất hiện ngay phía trên đỉnh vùng hóa trị (Valence Band). Cấu trúc điện tử mới này làm tăng thế hở mạch (Voc) từ giá trị cơ sở 0,48 V lên xấp xỉ 0,79 V trong các linh kiện quang điện, đồng thời giảm mật độ dòng tối và hạn chế sự tái kết hợp hạt mang điện không mong muốn ở mặt phân giới bán dẫn.

Thứ ba, việc hình thành các nút khuyết oxy tại các vị trí O1, O2 và O3 trong mạng tinh thể TiO2 pha Al tạo ra các trạng thái bẫy điện tử cục bộ nằm sâu trong vùng cấm. Sự kết hợp giữa pha tạp Al và khuyết tật oxy làm thu hẹp khe năng lượng quang học hiệu dụng xuống mức 1,3 eV đến 1,9 eV, tương đương với các đỉnh hấp thụ ánh sáng mới ở bước sóng 628 nm, 758 nm và mở rộng đến 980 nm.

Thứ tư, khi tích hợp các cấu trúc biến tính nano TiO2 pha tạp vào pin mặt trời DSSC, hiệu suất lượng tử chuyển đổi photon thành dòng điện (IPCE) đạt cực đại khoảng 13,5%, nâng cao mật độ dòng quang điện ngắn mạch (Jsc) lên mức 4,5 mA/cm2 và cải thiện hiệu suất chuyển đổi năng lượng tổng thể tăng từ 37% đến 123% so với màng TiO2 nguyên bản chưa biến tính.

Thảo luận kết quả

Sự biến đổi cấu trúc điện tử khi pha tạp Al bắt nguồn từ sự chênh lệch mức oxy hóa giữa ion Al3+ và Ti4+, dẫn đến việc hình thành các lỗ trống mang điện tích dương và làm dịch chuyển mức Fermi chuẩn về gần vùng hóa trị. Đồng thời, các nút khuyết oxy làm đứt gãy tính đối xứng cục bộ của bát diện TiO6, giải phóng 2 electron tự do định xứ tại các orbital Ti 3d lân cận, tạo nên các mức donor trung gian trong khe cấm. Kết quả này hoàn toàn nhất quán với các mô hình lý thuyết của Zhao và giải thích hoàn hảo hiện tượng tăng vọt thế hở mạch Voc trong các nghiên cứu thực nghiệm của Kyung Hyun Ko.

Dữ liệu tính toán cấu trúc vùng năng lượng và mật độ trạng thái (DOS/PDOS) được mô tả trực quan và chi tiết thông qua các biểu đồ phân bố năng lượng dọc theo các phương đối xứng cao trong vùng Brillouin (Gamma - X - M - Z) kết hợp bảng tổng hợp hằng số mạng và mật độ hạt tải điện. Biểu đồ DOS thể hiện rõ sự lai hóa mạnh mẽ giữa trạng thái Al 3p và O 2p ở đỉnh vùng hóa trị, chứng minh rằng sự kết hợp giữa pha tạp kim loại và tạo khuyết tật mạng là chìa khóa then chốt để khai thác dải quang phổ khả kiến (chiếm 45% bức xạ mặt trời).

Đề xuất và khuyến nghị

Thứ nhất, triển khai tối ưu hóa công nghệ đồng pha tạp (Co-doping) giữa kim loại Al và vonfram (W) trên màng nano TiO2 với nồng độ kiểm soát từ 3% đến 6%, nhằm đạt mục tiêu nâng hiệu suất chuyển đổi năng lượng pin mặt trời DSSC vượt ngưỡng 12% trong vòng 12 đến 18 tháng, do các phòng thí nghiệm vật liệu quang điện tử chủ trì thực hiện.

Thứ hai, thiết lập quy trình kiểm soát mật độ khuyết tật oxy thông qua kỹ thuật ủ nhiệt trong môi trường khí khử hidro hoặc argon ở dải nhiệt độ 400 đến 500 độ C, nhằm duy trì nồng độ khuyết oxy dưới 2,5%, triệt tiêu các bẫy tái hợp không bức xạ, hoàn thành trong giai đoạn 6 đến 9 tháng bởi các nhóm nghiên cứu công nghệ nano thực nghiệm.

Thứ ba, mở rộng tích hợp màng nano TiO2 pha Al với các chấm lượng tử ZnSe hoặc hạt nano vàng (GNPs) để kích hoạt hiệu ứng cộng hưởng plasmon bề mặt, mở rộng dải hấp thụ photon đến bước sóng 780 nm, hướng tới nâng mật độ dòng ngắn mạch Jsc lên trên 5,0 mA/cm2 trong lộ trình 24 tháng dưới sự phối hợp giữa viện nghiên cứu và doanh nghiệp sản xuất pin mặt trời.

Thứ tư, áp dụng các phiếm hàm tương quan trao đổi bậc cao như phiếm hàm lai HSE06 hoặc hiệu chỉnh tham số tương tác Coulomb U (DFT+U với U = 4,2 eV cho orbital Ti 3d) trong các nghiên cứu mô phỏng tiếp theo, hoàn thành trong 6 tháng tới do nhóm tính toán lý thuyết thực hiện nhằm khắc phục triệt để hiện tượng đánh giá thấp độ rộng vùng cấm của gần đúng GGA tiêu chuẩn.

Đối tượng nên tham khảo luận văn

Nhóm học viên cao học và nghiên cứu sinh chuyên ngành Vật lý lý thuyết, Vật lý chất rắn, Hóa lý tính toán: Luận văn cung cấp phương pháp luận chuẩn xác về việc ứng dụng phương trình Kohn-Sham, xây dựng siêu ô mạng tinh thể và kỹ thuật xử lý dữ liệu mật độ trạng thái điện tử (DOS/PDOS) phục vụ các đề tài mô phỏng vật liệu mới.

Nhóm kỹ sư R&D tại các doanh nghiệp chế tạo pin mặt trời và linh kiện bán dẫn: Tiếp cận các thông số định lượng về cấu trúc vùng cấm, cơ chế truyền dẫn hạt tải và giải pháp gia tăng điện thế hở mạch Voc để ứng dụng trực tiếp vào quy trình phủ màng quang anode cho pin mặt trời nhạy màu DSSC và pin Perovskite.

Nhóm nhà khoa học lĩnh vực công nghệ môi trường và xúc tác hóa học: Khai thác cơ chế tạo gốc hydroxyl (OH•) và superoxide (O2-) khi vật liệu TiO2 biến tính được kích hoạt dưới ánh sáng khả kiến để thiết kế các hệ thống màng lọc quang xúc tác xử lý nước thải và diệt khuẩn hiệu quả cao.

Nhóm giảng viên đại học chuyên trách các môn học Vật lý vật liệu rắn và Cơ học lượng tử ứng dụng: Sử dụng công trình như một tài liệu tham khảo giảng dạy chuyên sâu, minh họa sinh động mối liên kết giữa lý thuyết phiếm hàm mật độ trừu tượng và các ứng dụng kỹ thuật thực tiễn.

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

Vật liệu TiO2 dạng nguyên bản gặp phải hạn chế lớn nhất nào khi ứng dụng trong pin mặt trời? TiO2 nguyên bản pha anatase có khe năng lượng rộng 3,2 eV, chỉ hấp thụ được ánh sáng vùng tử ngoại chiếm khoảng 4% đến 5% phổ bức xạ mặt trời. Hạn chế này khiến mật độ dòng quang điện thấp và hiệu suất chuyển đổi quang năng tổng thể của pin không đạt mức tối ưu.

Phương pháp phiếm hàm mật độ (DFT) có ưu thế gì vượt trội so với phương pháp hàm sóng truyền thống? DFT thay thế hàm sóng đa điện tử 3N biến số phức tạp bằng đại lượng mật độ điện tử không gian 3 chiều r(r). Cách tiếp cận này giúp giảm thiểu đột phá khối lượng tính toán nhưng vẫn bảo toàn độ chính xác cao khi xác định năng lượng trạng thái cơ bản và cấu trúc dải năng lượng.

Việc pha tạp nhôm (Al) làm thay đổi đặc tính điện tử của TiO2 theo cơ chế nào? Khi thay thế cation Ti4+ bằng ion Al3+, các trạng thái acceptor mới được tạo ra ngay phía trên đỉnh vùng hóa trị. Biến đổi này giúp nâng cao thế hở mạch Voc lên mức xấp xỉ 0,79 V và giảm thiểu dòng tối tái hợp hạt tải tại bề mặt tiếp xúc.

Nút khuyết oxy trong mạng tinh thể TiO2 đóng vai trò gì đối với khả năng hấp thụ quang học? Nút khuyết oxy giải phóng các electron tự do định xứ tại orbital Ti 3d, tạo ra các mức năng lượng bẫy nằm giữa khe cấm. Nhờ đó, vật liệu hấp thụ thêm các photon năng lượng thấp trong dải ánh sáng khả kiến và cận hồng ngoại (bước sóng từ 600 nm đến 980 nm).

Các kết quả mô phỏng cấu trúc điện tử trong luận văn có thể chuyển giao sang thực tế sản xuất không? Hoàn toàn khả thi. Các thông số về nồng độ pha tạp Al (6,25%) và kiểm soát nút khuyết oxy cung cấp cơ sở định lượng trực tiếp để các kỹ sư thực nghiệm thiết lập nhiệt độ ủ màng và tỷ lệ phối trộn tiền chất tối ưu trong công nghệ chế tạo pin DSSC.

Kết luận

Luận văn đã thực hiện thành công các đóng góp khoa học nổi bật sau:

  • Xây dựng hoàn chỉnh mô hình tính toán First-principles dựa trên lý thuyết phiếm hàm mật độ DFT-GGA cho siêu ô mạng anatase TiO2 (2 x 2 x 1) gồm 48 nguyên tử.
  • Xác định chính xác cấu trúc dải năng lượng và làm sáng tỏ cơ chế hình thành mức tạp chất acceptor khi pha tạp nguyên tử Al (Ti15/16Al1/16O32).
  • Phát hiện và chứng minh quy luật thu hẹp khe năng lượng quang học hiệu dụng xuống dải 1,3 eV đến 1,9 eV nhờ sự xuất hiện của các nút khuyết oxy O1, O2 và O3.
  • Cung cấp luận cứ giải thích sự gia tăng thế hở mạch lên 0,79 V và nâng cao hiệu suất chuyển đổi lượng tử IPCE cực đại đạt 13,5% trong pin mặt trời DSSC.
  • Đặt nền móng lý thuyết vững chắc cho lộ trình 12 đến 24 tháng tới nhằm chế tạo thực nghiệm màng nano bán dẫn đa pha tạp hiệu suất cao.

Công trình là tài liệu khoa học giá trị cao cho cộng đồng nghiên cứu vật liệu bán dẫn quang điện. Hãy áp dụng ngay các giải pháp mô phỏng và công nghệ pha tạp này để tối ưu hóa thiết kế pin mặt trời thế hệ mới.