Nghiên cứu mô hình hóa và mô phỏng cấu trúc nano xốp ZnO

Trong phần này, tài liệu cung cấp cái nhìn tổng quan về cấu trúc nano và vật liệu nano, đặc biệt là ZnO. 157 trang · Trường Đại học Bách Khoa Hà Nội · 2023.

• Tác giả: Nguyễn Thị Thảo
• Người hướng dẫn: PGS. Vũ Ngọc Tước
• Trường học: Trường Đại học Bách Khoa Hà Nội
• Chuyên ngành: Vật lý kỹ thuật
• Đề tài: Mô Hình Hóa Và Mô Phỏng Cấu Trúc Nano Xốp ZnO
• Loại tài liệu: luận án
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• 1. CHƯƠNG 1: TỔNG QUAN VỀ VẬT LIỆU VÀ CÁC CẤU TRÚC NANO BÁN DẪN THẤP CHIỀU
• 1.1. Tổng quan về vật liệu và cấu trúc nano bán dẫn thấp chiều
• 1.2. Định nghĩa vật liệu nano
• 1.3. Phân loại vật liệu nano
• 1.4. Chế tạo vật liệu nano
• 1.5. Tổng quan về các vật liệu nghiên cứu
• 1.6. Sơ lược về vật liệu bán dẫn
• 1.7. Phân loại vật liệu bán dẫn theo cấu trúc nguyên tử
• 1.8. Vật liệu ôxit kẽm (ZnO)
• 1.9. Đặc điểm cấu trúc và các thuộc tính: Ứng dụng tiềm năng
• 1.10. Cấu trúc nano xốp
• 1.11. Phân loại nano xốp
• 1.12. Zeolite - Nano xốp vô cơ
• 1.13. Khung kim loại hữu cơ - MOF
• 1.14. Siêu vật liệu
• 1.15. Các tinh thể Fullerite
• 1.16. Các khoáng sét nanoclay
• 2. CHƯƠNG 2: CƠ SỞ LÝ THUYẾT PHIẾM HÀM MẬT ĐỘ VÀ GẦN ĐÚNG LIÊN KẾT CHẶT DỰA TRÊN DFT
• 2.1. Phương trình Schrödinger độc lập thời gian
• 2.2. Gần đúng Born-Oppenheimer
• 2.3. Lý thuyết phiếm hàm mật độ DFT. Các Định lý Höhenberg-Kohn
• 2.4. Phương trình Kohn-Sham
• 2.5. Thế hiệu dụng Kohn-Sham
• 2.6. Phiếm hàm trao đổi tương quan
• 2.7. Gần đúng mật độ cục bộ
• 2.8. Gần đúng Gradient tổng quát
• 2.9. Phương pháp trường tự hợp
• 2.10. Phương pháp phiếm hàm mật độ kết hợp gần đúng liên kết chặt tự hợp điện tích SCC-DFTB
• 2.11. Mô hình gần đúng liên kết chặt
• 2.12. Phương pháp SCC-DFTB
• 2.13. Sắp xếp lại công thức năng lượng tổng Kohn-Sham
• 2.14. Các gần đúng SCC-DFTB
• 3. CHƯƠNG 3: NGHIÊN CỨU CÁC CẤU TRÚC NANO XỐP MẬT ĐỘ THẤP BẰNG PHƯƠNG PHÁP TIẾP CẬN TỪ DƯỚI LÊN
• 3.1. Phương pháp dự đoán cấu trúc bằng cách tiếp cận từ dưới lên
• 3.2. Tính toán chi tiết
• 3.3. Năng lượng liên kết, độ bền vững của cấu trúc
• 3.4. Cấu trúc vùng năng lượng của điện tử
• 3.5. Các thông số cấu trúc
• 4. CHƯƠNG 4: NGHIÊN CỨU CÁC CẤU TRÚC NANO XỐP KÊNH RỖNG DẠNG LỤC GIÁC BẰNG PHƯƠNG PHÁP TIẾP CẬN TỪ TRÊN XUỐNG
• 4.1. Phương pháp thiết kế cấu trúc bằng phương pháp từ trên xuống
• 4.2. Chi tiết tính toán
• 4.3. Năng lượng dao động tự do
• 4.4. Tính mật độ trạng thái phonon bằng lý thuyết DFT
• 4.5. Các cấu trúc nano xốp kênh rỗng dạng lục giác
• 4.6. Mô tả về cách phân loại cấu trúc
• 4.7. Năng lượng liên kết, độ bền vững của pha và phương trình trạng thái
• 4.8. Cấu trúc vùng điện tử
• 4.9. Kết quả và thảo luận
• 5. CHƯƠNG 5: NGHIÊN CỨU CÁC CẤU TRÚC NANO XỐP KÊNH RỖNG DẠNG TAM GIÁC VÀ THOI BẰNG PHƯƠNG PHÁP TIẾP CẬN TỪ TRÊN XUỐNG
• 5.1. Thiết kế các cấu trúc nano xốp kênh rỗng dạng thoi và tam giác
• 5.2. Chi tiết tính toán
• 5.3. Đánh giá độ bền vững của cấu trúc
• 5.4. Mô phỏng ảnh nhiễu xạ tia X
• 5.5. Tính chất cơ học của các cấu trúc
• 5.6. Tính chất điện tử
• KẾT LUẬN VÀ KIẾN NGHỊ
• TÀI LIỆU THAM KHẢO
• DANH MỤC CÁC CÔNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN

Tóm tắt

I. Tổng quan về vật liệu và cấu trúc nano

Trong phần này, tài liệu cung cấp cái nhìn tổng quan về cấu trúc nano và vật liệu nano, đặc biệt là ZnO. Mô hình hóa và mô phỏng là những công cụ quan trọng trong nghiên cứu vật liệu nano, cho phép dự đoán các tính chất vật lý và hóa học của chúng. Cấu trúc nano xốp ZnO được phân tích từ góc độ cấu trúc tinh thể và các thuộc tính của nó. Các loại vật liệu nano được phân loại, từ zeolite đến MOF, cho thấy sự đa dạng trong ứng dụng và tính chất. Đặc biệt, ZnO nano nổi bật với các ứng dụng trong lĩnh vực quang học và điện tử, nhờ vào tính chất quang học độc đáo của nó.

1.1. Định nghĩa và phân loại vật liệu nano

Vật liệu nano được định nghĩa là các vật liệu có kích thước trong khoảng từ 1 đến 100 nm. Phân loại vật liệu nano bao gồm vật liệu vô cơ, hữu cơ và các cấu trúc lai. Cấu trúc nano xốp là một trong những loại vật liệu nano quan trọng, với khả năng hấp thụ và lưu trữ năng lượng cao. ZnO là một ví dụ điển hình cho vật liệu nano với tính chất quang học và điện tử vượt trội, được ứng dụng trong các thiết bị quang điện và cảm biến. Việc hiểu rõ về cấu trúc nano và các loại vật liệu này là cần thiết để phát triển các ứng dụng mới trong công nghệ nano.

II. Cơ sở lý thuyết phiếm hàm mật độ

Phần này trình bày các nguyên lý cơ bản của lý thuyết phiếm hàm mật độ (DFT), một công cụ quan trọng trong mô phỏng và mô hình hóa vật liệu. DFT cho phép tính toán các tính chất của hệ thống ở cấp độ nguyên tử, từ đó dự đoán các đặc tính của cấu trúc nano. Các định lý của Hohenberg-Kohn và phương trình Kohn-Sham là nền tảng cho việc áp dụng DFT trong nghiên cứu vật liệu. Việc sử dụng DFT trong nghiên cứu ZnO nano giúp hiểu rõ hơn về các tính chất điện tử và quang học của nó, từ đó mở ra hướng đi mới cho các ứng dụng trong công nghệ nano.

2.1. Các định lý Hohenberg Kohn

Định lý Hohenberg-Kohn khẳng định rằng mọi tính chất của một hệ nhiều electron có thể được xác định hoàn toàn từ mật độ electron. Điều này có nghĩa là việc tính toán mật độ electron là chìa khóa để hiểu các tính chất vật lý của vật liệu. Trong nghiên cứu cấu trúc nano xốp ZnO, việc áp dụng DFT giúp dự đoán chính xác các đặc tính như năng lượng liên kết và độ bền của cấu trúc. Sự phát triển của các phương pháp gần đúng trong DFT, như GGA và LDA, đã cải thiện độ chính xác trong mô phỏng các tính chất của vật liệu nano.

III. Nghiên cứu các cấu trúc nano xốp mật độ thấp

Nghiên cứu này tập trung vào việc thiết kế và mô phỏng các cấu trúc nano xốp mật độ thấp bằng phương pháp tiếp cận từ dưới lên. Các tính toán chi tiết về năng lượng liên kết và độ bền của cấu trúc được thực hiện để đánh giá tính khả thi của các cấu trúc này. Kết quả cho thấy rằng các cấu trúc nano xốp có thể đạt được độ bền cao và tính chất quang học tốt, mở ra khả năng ứng dụng trong các lĩnh vực như cảm biến và lưu trữ năng lượng.

3.1. Phương pháp dự đoán cấu trúc

Phương pháp dự đoán cấu trúc từ dưới lên sử dụng các thuật toán tối ưu hóa để tìm kiếm cấu trúc có năng lượng thấp nhất. Các thông số cấu trúc như năng lượng liên kết và mật độ trạng thái được tính toán để đánh giá tính ổn định của các cấu trúc nano xốp. Kết quả cho thấy rằng các cấu trúc này không chỉ ổn định mà còn có tính chất quang học và điện tử vượt trội, cho thấy tiềm năng ứng dụng trong công nghệ nano.

IV. Nghiên cứu các cấu trúc nano xốp kênh rỗng

Phần này nghiên cứu các cấu trúc nano xốp kênh rỗng dạng lục giác và tam giác, sử dụng phương pháp tiếp cận từ trên xuống. Các tính toán chi tiết về năng lượng dao động tự do và mật độ trạng thái phonon được thực hiện để đánh giá tính chất của các cấu trúc này. Kết quả cho thấy rằng các cấu trúc này có độ bền cao và tính chất quang học tốt, mở ra khả năng ứng dụng trong các lĩnh vực như cảm biến và lưu trữ năng lượng.

4.1. Thiết kế cấu trúc kênh rỗng

Thiết kế các cấu trúc nano xốp kênh rỗng yêu cầu sự kết hợp giữa lý thuyết và thực nghiệm. Các thông số như kích thước hốc rỗng và độ dày vách được tối ưu hóa để đạt được tính chất mong muốn. Kết quả cho thấy rằng các cấu trúc này không chỉ ổn định mà còn có tính chất quang học và điện tử vượt trội, cho thấy tiềm năng ứng dụng trong công nghệ nano.

V. Kết luận và kiến nghị

Tài liệu kết luận rằng việc mô hình hóa và mô phỏng các cấu trúc nano xốp ZnO là một lĩnh vực nghiên cứu đầy tiềm năng. Các kết quả đạt được không chỉ cung cấp cái nhìn sâu sắc về tính chất của vật liệu mà còn mở ra hướng đi mới cho các ứng dụng trong công nghệ nano. Đề xuất nghiên cứu tiếp theo bao gồm việc mở rộng nghiên cứu sang các loại vật liệu khác và phát triển các phương pháp mô phỏng mới để cải thiện độ chính xác và hiệu quả.

5.1. Đề xuất nghiên cứu tiếp theo

Nghiên cứu tiếp theo nên tập trung vào việc mở rộng các mô hình mô phỏng để bao gồm các loại vật liệu khác ngoài ZnO. Việc phát triển các phương pháp mô phỏng mới có thể giúp cải thiện độ chính xác và hiệu quả trong việc dự đoán các tính chất của vật liệu nano. Hơn nữa, việc kết hợp giữa lý thuyết và thực nghiệm sẽ giúp xác thực các kết quả mô phỏng và mở rộng ứng dụng của vật liệu nano trong thực tế.

Nội dung chính

Tổng quan luận án

Luận án tiến sĩ với đề tài "Nghiên cứu mô hình hóa và mô phỏng cấu trúc nano xốp ZnO" do nghiên cứu sinh Nguyễn Thị Thảo thực hiện dưới sự hướng dẫn của PGS. Vũ Ngọc Tước tại Viện Vật lý kỹ thuật, Trường Đại học Bách Khoa Hà Nội phối hợp cùng Trường Đại học Hồng Đức. Công trình thuộc chuyên ngành Vật lý kỹ thuật, tập trung vào lĩnh vực Khoa học tính toán vật liệu (Computational Materials Science - CMS).

Tính cấp thiết và khoảng trống nghiên cứu

Trong khoa học và công nghệ nano, khi kích thước vật liệu giảm xuống thang nanomet, các hiệu ứng lượng tử và hiệu ứng bề mặt xuất hiện làm thay đổi đặc tính quang, điện, từ và cơ học mà không cần thay đổi thành phần hóa học. Tuy nhiên, các phương pháp chế tạo thực nghiệm truyền thống theo quy trình "thử - loại" đối mặt với chi phí cao và giới hạn công nghệ khi tiếp cận kích thước nguyên tử. Khoa học tính toán vật liệu đóng vai trò như một "thí nghiệm ảo", cho phép thiết kế và dự đoán các cấu trúc vật liệu mới từ các nguyên lý ban đầu trước khi tiến hành thực nghiệm.

Kẽm oxit (ZnO) là vật liệu bán dẫn nhóm II-VI có độ rộng vùng cấm thẳng 3,37 eV và năng lượng liên kết exciton lớn (60 meV) ở nhiệt độ phòng. Dù các cấu trúc nano 0D, 1D, 2D của ZnO đã được nghiên cứu rộng rãi, các cấu trúc tích hợp tuần hoàn ba chiều như vật liệu nano xốp (nanoporous) và siêu vật liệu (metamaterials) vẫn tồn tại nhiều khoảng trống lý thuyết:

• Thiếu các mô hình dự đoán cấu trúc tinh thể nano xốp mật độ thấp bền vững được xây dựng từ các cụm nguyên tử ma thuật cơ bản.
• Chưa có nghiên cứu hệ thống về các họ cấu trúc nano xốp kênh rỗng với các dạng topo hình học khác nhau, cũng như quy luật biến đổi phi tuyến của tính chất cơ học (module khối) và tính chất điện tử (vùng cấm) theo độ xốp và bề dày vách.

Mục tiêu nghiên cứu

Mục tiêu tổng quát của luận án là ứng dụng các phương pháp mô hình hóa và mô phỏng lượng tử để thiết kế, dự đoán lý thuyết các pha đa hình nano xốp tinh thể mới từ bán dẫn ZnO, đồng thời khảo sát các đặc trưng vật lý hướng tới ứng dụng trong nano quang-spin điện tử và nano sinh học. Các mục tiêu cụ thể gồm:

1. Thiết kế và dự đoán lý thuyết các pha đa hình nano xốp tinh thể dạng lồng rỗng mật độ thấp từ các cụm nguyên tử ma thuật bằng phương pháp tiếp cận từ dưới lên (bottom-up).
2. Thiết kế và khảo sát các họ cấu trúc nano xốp dạng kênh rỗng với các dạng hình học khác nhau (lục giác, thoi, tam giác) bằng phương pháp tiếp cận từ trên xuống (top-down).
3. Đánh giá độ bền nhiệt động, độ ổn định động học, cấu trúc vùng năng lượng, mật độ trạng thái điện tử (DOS/PDOS), phổ phonon, ảnh nhiễu xạ tia X (XRD) và tính chất cơ học (module khối) của các cấu trúc được thiết kế.
4. Rút ra các quy luật tương quan giữa cấu trúc hình học topo, kích thước hốc xốp, bề dày vách với module khối và độ rộng vùng cấm năng lượng.

Đối tượng và phạm vi nghiên cứu

• Đối tượng nghiên cứu: Các cấu trúc nano bán dẫn ZnO được tích hợp thành các mạng tuần hoàn tinh thể 3D, bao gồm các pha đa hình nano xốp dạng lồng rỗng và các họ nano xốp dạng kênh rỗng.
• Phạm vi nghiên cứu: Giới hạn trong các tính toán mô phỏng lý thuyết lượng tử (DFT và SCC-DFTB) đối với các cấu trúc tuần hoàn lý tưởng, không xét đến các sai hỏng bất trật tự vô định hình quy mô lớn vượt quá năng lực tính toán.

---

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước

Luận án điểm lại các mốc phát triển chính trong lĩnh vực tính toán vật liệu và cấu trúc nano:

• Dự đoán cấu trúc tinh thể (Crystal Structure Prediction): M. Neumann (2007) trên tạp chí Nature đã chứng minh khả năng tiên đoán chính xác cấu trúc tinh thể phân tử hữu cơ từ các nguyên lý ban đầu mà không cần dữ liệu thực nghiệm.
• Khái niệm siêu vật liệu (Metamaterials): Thuật ngữ được Rodger M. Walser (Đại học Texas tại Austin) đặt ra năm 1999 để chỉ các vật liệu nhân tạo tuần hoàn có tính chất vĩ mô phụ thuộc vào cấu trúc hình học hơn là thành phần. Các hướng nghiên cứu tiếp nối gồm chiết suất âm (V. Veselago; Pendry và Smith; A. Grbic và G.V. Eleftheriades tại Đại học Toronto) và siêu vật liệu cơ học tổ ong mật độ thấp chịu lực cao của Oliver Kraft (KIT, CHLB Đức).
• Công nghệ chế tạo và lắp ráp nano: K.A. Geim (2004) với phương pháp bóc tách vi cơ chế tạo graphene; Y. Tian et al. (2016) phát triển kỹ thuật tự lắp ráp mạng hạt nano sử dụng khung đa diện DNA.
• Nghiên cứu các pha nano xốp ZnO: Các công trình của Carrasco et al. (Nanoscale, 2014), Sponza et al. (Physical Review B, 2015) và Zagorac et al. (2015) đã đề xuất 7 pha đa hình nano xốp ZnO dạng zeolite gồm: SOD, FAU, LTA, EMT, CAN, GME và $\text{Zn}{16}\text{O}{16}\text{-cub}$.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù 7 pha đa hình cấu trúc zeolite của ZnO đã được các nhóm tác giả quốc tế công bố, không gian cấu hình của các pha lồng rỗng mật độ thấp từ các cụm nguyên tử ma thuật nhỏ hơn vẫn chưa được khai phá đầy đủ. Đồng thời, đối với cấu trúc nano xốp kênh rỗng 3D, chưa có công trình nào phân loại toàn diện 4 họ topo hình học (lục giác, thoi, hai dạng tam giác) và thiết lập mô hình định lượng giải thích sự phụ thuộc phi tuyến của module khối cũng như độ rộng vùng cấm vào độ dày vách và kích thước kênh. Luận án được định vị nhằm giải quyết trọn vẹn các bài toán này.

---

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

Khung lý thuyết và mô hình toán học

Nghiên cứu sử dụng các phương pháp tính toán lượng tử từ các nguyên lý ban đầu (ab-initio) và phương pháp bán kinh nghiệm dựa trên DFT:

• Phương trình Schrödinger độc lập thời gian và Gần đúng Born-Oppenheimer: Tách biệt chuyển động của hạt nhân nặng và hệ điện tử nhẹ, cho phép giải bài toán trạng thái dừng của hệ nhiều hạt trên tọa độ hạt nhân cố định.
• Lý thuyết phiếm hàm mật độ (DFT): Dựa trên hai định lý Hohenberg-Kohn chứng minh mật độ điện tử $n(\mathbf{r})$ là biến số cơ bản quyết định mọi tính chất trạng thái cơ bản, và hệ phương trình tự hợp Kohn-Sham biến đổi hệ electron tương tác phức tạp thành hệ đơn hạt chuyển động trong thế hiệu dụng $V_{eff}[n(\mathbf{r})]$.
• Phiếm hàm trao đổi - tương quan (XC): Sử dụng các gần đúng mật độ cục bộ (LDA), gần đúng gradient tổng quát (GGA với các tham số PBE, PBESol) và phiếm hàm lai HSE.
• Phương pháp phiếm hàm mật độ kết hợp gần đúng liên kết chặt tự hợp điện tích (SCC-DFTB): Khai triển năng lượng toàn phần Kohn-Sham đến bậc hai của thăng giáng mật độ điện tích quanh mật độ nguyên tử chồng chập, sử dụng bộ tích phân Slater-Koster không chứa tham số thực nghiệm. Phương pháp này cho phép tính toán các hệ tuần hoàn chứa hàng trăm đến hàng ngàn nguyên tử với tốc độ vượt trội nhưng vẫn duy trì độ chính xác của DFT.

Phương pháp tiếp cận thiết kế cấu trúc

1. Tiếp cận từ dưới lên (Bottom-up): Chọn các cụm nguyên tử ma thuật đứng độc lập $\text{Zn}_k\text{O}_k$ đã tối ưu hóa làm thành tố cơ bản, sau đó liên kết tuần hoàn trong không gian ba chiều để tạo thành các mạng tinh thể lồng rỗng.
2. Tiếp cận từ trên xuống (Top-down): Xuất phát từ mạng tinh thể ZnO khối dạng wurtzite lý tưởng, tiến hành tạo các kênh rỗng tuần hoàn chạy dọc theo trục tinh thể với các hình học xác định, sau đó hồi phục hoàn toàn cấu trúc hình học và ô mạng đơn vị.

Công cụ tính toán và kỹ thuật phân tích

• Phần mềm thực thi: DFTB+, VASP, SIESTA; phần mềm biểu diễn cấu trúc tinh thể và mật độ điện tử VESTA.
• Kỹ thuật phân tích:
  • Khớp đường cong năng lượng - thể tích $E(V)$ bằng phương trình trạng thái Birch-Murnaghan bậc ba để xác định năng lượng liên kết riêng, thể tích cân bằng và module khối $B_0$.
  • Tính toán năng lượng dao động tự do Helmholtz $F(T)$ và mật độ trạng thái phonon bằng DFT để kiểm tra tính ổn định động học.
  • Phân tích cấu trúc dải năng lượng (Band structure), mật độ trạng thái tổng cộng (DOS), mật độ trạng thái thành phần (PDOS), orbital phân tử HOMO/LUMO và hiệu mật độ điện tích.
  • Mô phỏng phổ nhiễu xạ tia X (XRD) để cung cấp dấu vết nhận diện cho thực nghiệm.

---

Nội dung chính theo từng chương

Chương 1: Tổng quan về vật liệu và các cấu trúc nano bán dẫn thấp chiều

Chương này hệ thống hóa các khái niệm nền tảng về vật liệu nano, hiện tượng giam hãm lượng tử và các công nghệ chế tạo:

• Trình bày công thức lượng tử hóa năng lượng theo các chiều giam hãm không gian (0D, 1D, 2D) và sự biến đổi của mật độ trạng thái năng lượng tương ứng.
• Phân tích ưu nhược điểm của các kỹ thuật chế tạo từ trên xuống (quang khắc, ăn mòn hóa học, bóc tách vi cơ MC) và từ dưới lên (lắng đọng pha hơi hóa học CVD, epitaxy chùm phân tử MBE, tự lắp ráp DNA).
• Tổng quan các tính chất vật lý của bán dẫn ZnO: cấu trúc tinh thể wurtzite ($P6_3mc$, $a = 0,325\text{ nm}$, $c/a \approx 1,600$) và zincblende; tính chất áp điện mạnh do phân cực liên kết ion $\text{Zn}^{2+}-\text{O}^{2-}$; độ rộng vùng cấm 3,37 eV và năng lượng liên kết exciton 60 meV.
• Phân loại vật liệu nano xốp theo tiêu chuẩn IUPAC (vi xốp < 2 nm, xốp trung bình 2–50 nm, xốp vĩ mô > 50 nm); tổng quan về Zeolite, khung kim loại hữu cơ (MOF), siêu vật liệu điện từ và siêu vật liệu cơ học.

Chương 2: Cơ sở lý thuyết phiếm hàm mật độ và gần đúng liên kết chặt dựa trên DFT

Nội dung chương trình bày chi tiết các công cụ toán học và gần đúng lượng tử:

• Thiết lập phương trình Schrödinger nhiều hạt, áp dụng gần đúng Born-Oppenheimer.
• Trình bày dạng tường minh của hai định lý Hohenberg-Kohn và phương trình Kohn-Sham, thế hiệu dụng trao đổi tương quan $V_{xc}$.
• Chi tiết hóa sơ đồ lặp tự hợp (SCF) giải phương trình ma trận Kohn-Sham.
• Xây dựng mô hình SCC-DFTB: phân tích công thức năng lượng tổng, thế đẩy ion-ion hai tâm, ma trận Hamilton Slater-Koster và điều kiện tự hợp điện tích bậc hai $\Delta q_\alpha$.

Chương 3: Nghiên cứu các cấu trúc nano xốp mật độ thấp bằng phương pháp tiếp cận từ dưới lên

Chương này thực hiện thiết kế và mô phỏng 10 pha đa hình nano xốp dạng lồng rỗng từ các cụm nguyên tử ma thuật $\text{Zn}k\text{O}k$ đứng độc lập ($\text{Zn}9\text{O}9$, $\text{Zn}{12}\text{O}{12}\text{-a}$, $\text{Zn}{12}\text{O}{12}\text{-b}$, $\text{Zn}{16}\text{O}{16}\text{-cub}$):

• Tái lập và kiểm chứng 7 pha cấu trúc tham chiếu: SOD, FAU, LTA, EMT, CAN, GME và $\text{Zn}{16}\text{O}{16}\text{-cub}$.
• Dự đoán 3 pha đa hình tinh thể mới: AST, SAT và SBT.
• Kết quả tính toán năng lượng liên kết riêng theo thể tích tỷ đối cho thấy 3 pha mới có độ ổn định nhiệt động cao, năng lượng liên kết tương đương với các pha tham chiếu quốc tế.
• Phân tích cấu trúc dải năng lượng cho thấy tất cả các pha nano xốp đều duy trì vùng cấm thẳng hoặc gần thẳng, với độ rộng vùng cấm mở rộng đáng kể so với ZnO wurtzite khối do hiệu ứng giam hãm lượng tử trên thành lồng rỗng.

Chương 4: Nghiên cứu các cấu trúc nano xốp kênh rỗng dạng lục giác bằng phương pháp tiếp cận từ trên xuống

Chương 4 tập trung nghiên cứu họ cấu trúc nano xốp có kênh rỗng hình lục giác dạng tổ ong chạy dọc trục $c$ của mạng wurtzite:

• Thiết kế 8 chuỗi cấu trúc với độ dày vách biến thiên: Half Wall (HW), Single Wall (SW), 1.5 Wall, Double Wall (DW), 2.5 Wall, Triple Wall (TW), 3.5 Wall, Quadruple Wall (QW) cùng các kích thước hốc $n = 2, 3, 4, 5...$
• Đánh giá độ bền động học: Tính toán phổ phonon và năng lượng tự do Helmholtz $F(T)$ khẳng định cấu trúc không xuất hiện mode dao động mềm (tần số ảo), duy trì độ ổn định nhiệt động từ 0 K đến nhiệt độ cao.
• Kiểm chứng phương pháp: So sánh trực tiếp kết quả giữa DFTB+ và DFT trên các cấu trúc cơ sở cho thấy sự trùng khớp cao về hằng số mạng và năng lượng.
• Phát hiện sự phụ thuộc phi tuyến của module khối $B_0$ vào mật độ tỷ đối và bán kính hốc: Module khối giảm nhanh khi độ xốp tăng, sau đó đạt trạng thái bão hòa. Độ rộng vùng cấm tăng tỷ lệ nghịch với bán kính hốc và độ dày vách xốp.

Chương 5: Nghiên cứu các cấu trúc nano xốp kênh rỗng dạng tam giác và thoi bằng phương pháp tiếp cận từ trên xuống

Chương 5 mở rộng phạm vi nghiên cứu sang ba họ cấu trúc rỗng mới:

• Họ A: Kênh rỗng hình thoi.
• Họ B: Kênh rỗng hình tam giác hướng lên.
• Họ C: Kênh rỗng hình tam giác hướng xuống.
• Tổng số cấu trúc khảo sát trên toàn bộ 4 họ kênh rỗng đạt 119 cấu trúc.
• Tính toán cơ học và liên kết: So sánh module khối giữa 4 họ chỉ ra quy luật phân hạng độ cứng theo topo hình học: ở cùng một giá trị độ xốp, kênh lục giác có độ cứng lớn nhất, tiếp theo lần lượt là tam giác hướng lên (Họ B), hình thoi (Họ A), và thấp nhất là tam giác hướng xuống (Họ C).
• Xác định cấu hình tối ưu: Bề dày vách 3 lớp ô lục giác (Triple Wall - TW) là ngưỡng tối ưu giúp vật liệu tiệm cận đường chuẩn độ cứng trên tỷ trọng.
• Mô phỏng phổ XRD: Cung cấp phổ nhiễu xạ đặc trưng của từng họ cấu trúc, làm rõ sự dịch chuyển và tách đỉnh nhiễu xạ so với pha wurtzite nguyên bản.

Đặc trưng / Pha tinh thể	ZnO Wurtzite (WZ)	ZnO Rocksalt (B1)	ZnO Zincblende (ZB)	Cấu trúc Nano xốp (Dự đoán)
Nhóm không gian	$P6_3mc$	$Fm\bar{3}m$	$F\bar{4}3m$	Đa dạng (lục giác, thoi, tam giác, zeolite)
Hằng số mạng	$a = 0,325\text{ nm}, c/a \approx 1,600$	Biến đổi theo áp suất	$a_0 = 0,464\text{ nm}$	Phụ thuộc kích thước hốc và bề dày vách
Độ dài liên kết Zn-O	$0,234\text{ nm}$	-	$0,281\text{ nm}$	Biến thiên theo độ cong và phối vị cục bộ
Độ rộng vùng cấm ($E_g$)	$3,37\text{ eV}$ (thẳng, 300 K)	-	$3,21\text{ eV}$ (thẳng, 300 K)	Mở rộng $> 3,37\text{ eV}$ (tùy kích thước hốc)
Module khối ($B_0$) DFT	$113,7\text{ GPa}$ (0 K, GGA)	$175\text{ GPa}$ (0 K, GGA)	$109,8\text{ GPa}$ (0 K, GGA)	Biến thiên phi tuyến theo độ xốp
Module khối Thực nghiệm	$139,6\text{ GPa}$ (300 K)	$228\text{ GPa}$ (300 K)	-	-
Module khối DFTB+	$142\text{ GPa}$ (0 K)	-	-	Thấp hơn WZ khối, bão hòa theo độ xốp

---

Kết quả và những đóng góp mới

Luận án đạt được các kết quả và đóng góp mới sau:

1. Đề xuất 3 pha đa hình nano xốp lồng rỗng mới: Bằng phương pháp tiếp cận từ dưới lên xuất phát từ các cụm nguyên tử ma thuật $\text{Zn}_k\text{O}_k$, luận án đã dự đoán thành công 3 pha tinh thể mới gồm AST, SAT và SBT. Các tính toán nhiệt động khẳng định 3 pha này có độ bền vững sánh ngang với 7 pha cấu trúc zeolite tham chiếu do các tác giả quốc tế đề xuất trước đó.
2. Xây dựng cơ sở dữ liệu 4 họ cấu trúc kênh rỗng (119 cấu trúc): Bằng phương pháp tiếp cận từ trên xuống, luận án đã thiết kế và tính toán chi tiết 119 cấu trúc thuộc 4 họ topo kênh rỗng (lục giác, thoi, tam giác hướng lên, tam giác hướng xuống), chứng minh các cấu trúc này hoàn toàn ổn định cơ học và nhiệt động mà không bị sụp đổ mạng.
3. Phát hiện quy luật topo tổng quan về cơ học: Luận án chỉ ra mối quan hệ phi tuyến giữa module khối và độ xốp. Ở cùng một tỷ đối thể tích rỗng, độ cứng của vật liệu phụ thuộc chặt chẽ vào dạng topo hình học theo thứ tự: Kênh lục giác > Kênh tam giác B > Kênh thoi A > Kênh tam giác C.
4. Xác lập tiêu chí thiết kế vách tối ưu: Xác định cấu trúc có bề dày vách bằng 3 lớp ô lục giác (Triple Wall) là cấu hình tối ưu để đạt tỷ lệ độ cứng trên trọng lượng cao nhất, tiệm cận đường chuẩn cơ học của siêu vật liệu.
5. Cung cấp dữ liệu phổ XRD và cấu trúc dải năng lượng: Thiết lập hệ thống dữ liệu phổ nhiễu xạ tia X mô phỏng và quy luật mở rộng vùng cấm theo kích thước hình học, phục vụ việc định danh và "may đo" tính chất vật liệu trong thực nghiệm chế tạo linh kiện quang - điện tử và cảm biến.

---

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

• Các mô hình tính toán được thực hiện trên các ô mạng tinh thể tuần hoàn lý tưởng, chưa xét đến ảnh hưởng của sai hỏng mạng thực tế như nút khuyết nguyên tử ngẫu nhiên, sai lệch biên hạt hoặc pha tạp vô định hình quy mô lớn.
• Tính toán cấu trúc điện tử và phonon chủ yếu thực hiện ở trạng thái 0 K hoặc trong gần đúng điều hòa nhiệt độ, chưa mô phỏng động học phân tử đầy đủ ở điều kiện nhiệt độ và áp suất môi trường thực tế phức tạp.
• Giới hạn về tài nguyên tính toán khiến các cấu trúc kênh rỗng chỉ khảo sát đến kích thước vài nanomet, chưa thể mở rộng sang thang vi mô.

Hướng nghiên cứu tiếp

• Mở rộng mô hình nghiên cứu ảnh hưởng của pha tạp kim loại chuyển tiếp (Mn, Fe, Co) vào các vách nano xốp ZnO nhằm khảo sát tính chất sắt từ ở nhiệt độ phòng phục vụ lĩnh vực spin tử (spintronics).
• Nghiên cứu cơ chế hấp phụ phân tử khí ($\text{H}_2$, $\text{CO}_2$) và tương tác sinh học (protein, DNA) trên bề mặt vách kênh nano xốp ZnO phục vụ cảm biến y sinh và lưu trữ năng lượng.
• Phối hợp với các nhóm nghiên cứu thực nghiệm sử dụng phương pháp ăn mòn chọn lọc hoặc tự lắp ráp hóa học để hiện thực hóa các pha cấu trúc đã dự đoán.

---

Giá trị tham khảo

Luận án mang lại giá trị tham khảo khoa học cho nhiều đối tượng:

• Nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn, Khoa học vật liệu và Vật lý tính toán: Luận án là tài liệu tham khảo chi tiết về quy trình kết hợp giữa DFT và gần đúng SCC-DFTB, phương pháp thiết kế cấu trúc từ dưới lên và từ trên xuống, cũng như kỹ thuật khớp phương trình trạng thái Birch-Murnaghan.
• Giảng viên và nhà nghiên cứu lý thuyết: Cung cấp hệ thống dữ liệu chuẩn về thông số mạng, năng lượng liên kết, phổ phonon và module khối của các pha đa hình ZnO.
• Các kỹ sư và nhà thực nghiệm công nghệ nano: Bảng số liệu phổ XRD mô phỏng và các quy luật biến đổi vùng cấm là công cụ hữu ích để đối chiếu, phân tích cấu trúc sản phẩm màng xốp hoặc vật liệu nano xốp tổng hợp trong phòng thí nghiệm.

---

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

1. Ba pha nano xốp lồng rỗng mới do luận án đề xuất là những pha nào và được xây dựng từ đâu?

Ba pha đa hình mới gồm AST, SAT và SBT. Chúng được thiết kế bằng phương pháp tiếp cận từ dưới lên (bottom-up), sử dụng các cụm nguyên tử ma thuật bền vững $\text{Zn}9\text{O}9$ và $\text{Zn}{12}\text{O}{12}\text{-a}$ làm đơn vị cấu trúc cơ bản để liên kết tuần hoàn trong không gian 3 chiều.

2. Module khối của 4 họ nano xốp kênh rỗng phụ thuộc vào dạng topo hình học như thế nào?

Ở cùng một giá trị độ xốp (tỷ đối thể tích rỗng), module khối (đại diện cho độ cứng chống nén) biến đổi có quy luật theo hình học kênh: kênh lục giác (tổ ong) có độ cứng cao nhất, tiếp theo là kênh tam giác hướng lên (họ B), kênh hình thoi (họ A), và thấp nhất là kênh tam giác hướng xuống (họ C).

3. Vì sao luận án kết hợp phương pháp SCC-DFTB bên cạnh DFT truyền thống?

DFT nguyên lý ban đầu có độ chính xác cao nhưng đòi hỏi chi phí tính toán rất lớn, bị giới hạn đối với các ô mạng chứa hàng trăm hoặc hàng ngàn nguyên tử. SCC-DFTB sắp xếp lại năng lượng Kohn-Sham và sử dụng tích phân Slater-Koster tự hợp điện tích, giúp tăng tốc độ tính toán lên nhiều bậc, cho phép khảo sát toàn diện 119 cấu trúc nano xốp kích thước lớn mà vẫn đảm bảo độ tin cậy lượng tử.

4. Hiệu ứng kích thước lượng tử ảnh hưởng thế nào đến độ rộng vùng cấm của nano xốp ZnO?

Khi kích thước hốc rỗng tăng lên hoặc bề dày vách xốp giảm đi (tương đương với mức độ giam hãm không gian tăng), độ rộng vùng cấm năng lượng của các cấu trúc nano xốp ZnO có xu hướng mở rộng đáng kể so với giá trị 3,37 eV của ZnO wurtzite khối.

5. Bề dày vách xốp nào được xác định là tối ưu cho tính chất cơ học của nano xốp ZnO?

Kết quả tính toán chỉ ra bề dày vách tương ứng với 3 lớp ô lục giác (Triple Wall - TW) là giá trị tối ưu, giúp cấu trúc đạt được tỷ lệ module khối trên mật độ khối lượng cao nhất, tiệm cận đường chuẩn độ cứng của siêu vật liệu cơ học.

---

Kết luận

Luận án "Nghiên cứu mô hình hóa và mô phỏng cấu trúc nano xốp ZnO" của tác giả Nguyễn Thị Thảo đã ứng dụng các phương pháp tính toán lượng tử tiên tiến DFT và SCC-DFTB để dự đoán thành công 3 pha lồng rỗng mới (AST, SAT, SBT) và 119 cấu trúc thuộc 4 họ nano xốp kênh rỗng. Công trình đã làm sáng tỏ quy luật phụ thuộc phi tuyến của module khối và độ rộng vùng cấm vào độ xốp, bề dày vách và dạng hình học topo. Hệ thống kết quả tính toán nhiệt động, phổ phonon và phổ XRD mô phỏng cung cấp cơ sở dữ liệu lý thuyết đáng tin cậy, định hướng cho các nghiên cứu chế tạo thực nghiệm vật liệu nano xốp ZnO trong lĩnh vực quang điện tử và siêu vật liệu.

Trích đoạn nội dung tài liệu

CHƯƠNG 1: TỔNG QUAN VỀ VẬT LIỆU VÀ CÁC CẤU TRÚC NANO BÁN DẪN THẤP CHIỀU Trong vòng chỉ một hai thập kỷ gần đây, các tiến bộ trong lĩnh vực KH&CN nano đã thúc đẩy mạnh mẽ các hiểu biết nền tảng của khoa học vật liệu (KHVL) cũng như tăng cường cầu nối giữa lý thuyết cơ bản và các quan sát thực nghiệm. Cùng với lý thuyết và thực nghiệm, KH tính toán (computational science) đang trở thành ‘trụ cột thứ ba’ (pillar) của khoa học, cho phép các nhà nghiên cứu xây dựng và thử nghiệm mô hình của các hiện tượng phức tạp, đặc biệt là trong lĩnh vực tính toán KHVL (Computational Materials Science - CMS). Gần đây, CMS đã nổi lên như một nhánh liên ngành của KH&CN mà sự quan trọng và tính xác đáng của nó đến từ sự khao khát để hiểu biết vi mô về các hiện tượng và vật liệu ở kích cỡ nano. Năng lực dự đoán của các kỹ thuật CMS trên máy tính không những giảm bớt đáng kể các nhiệm vụ của các thí nghiệm đắt tiền mà còn mở ra một chân trời mới cho các khả năng “dự đoán các vật liệu mới”.

Sẽ không còn xa cái ngày mà các vật liệu với các đặc tính “có chủ định’’ được thiết kế trên máy tính trước khi tiến hành trên các thí nghiệm thực [106]. Vật liệu nano là một trong những hướng nghiên cứu được quan tâm đặc biệt của hầu hết các quốc gia và các trường đại học hàng đầu trên thế giới, trong đó có Việt Nam. Chỉ trong hơn một thập kỷ gần đây, các cấu trúc nano bán dẫn, chẳng hạn hạt nano, ống nano, dây nano, dải nano, nano xốp, nano-composite… đã và đang thu hút được rất nhiều quan tâm do tính xác đáng và tầm quan trọng của nó trong cả nghiên cứu cơ bản lẫn các tiềm năng ứng dụng CN. Chúng không chỉ là các hệ lý tưởng để nghiên cứu các hiệu ứng giam hãm lượng tử dẫn đến các đặc tính cơ, điện tử và quang phụ thuộc vào kích thước mà còn bộc lộ một loạt các chức năng thiết bị sẵn sàng trở thành các thành tố cơ bản cho nano điện tử và nano quang-spin tử (photonics, spintronics) tương lai.

Chúng vừa có thể làm chức năng thiết bị vừa có thể làm chức năng dây kết nối, có thể được kết hợp thành transistor hiệu ứng trường (FET), các đảo mạch, bộ phận cảm quang, cảm biến nano, diode phát quang và laser cũng như trong các cấu trúc mới có tiềm năng ứng dụng trong năng lượng bền vững như lưu trữ hydro và tế bào quang điện [3]. Trong luận án này chúng tôi tập trung nghiên cứu, mô phỏng tính toán các cấu trúc và vật liệu nano bán dẫn dạng tích hợp trong các cấu trúc ba chiều - các pha tinh thể nano xốp, nano composite. Các vật liệu quan tâm của luận án là các vật liệu thuộc nhóm bán dẫn 1 II-VI như ZnO, ZnS, ZnSe, CdSe, CdTe …. Trong chương này, chúng tôi sẽ tóm lược các kiến thức tổng quan về các vật liệu và cấu trúc nano nghiên cứu của luận án.Tổng quan về vật liệu và cấu trúc nano bán dẫn thấp chiều 1.

Định nghĩa vật liệu nano Vật liệu nano là vật liệu trong đó ít nhất một chiều của các thành tố cơ bản tạo nên nó có kích thước ở cỡ nanomet, tức là từ 1-1000nm nhưng thường là chỉ trong khoảng từ 1- 100nm (1nm=10-9m). Đây là đối tượng nghiên cứu của KH nano và CN nano, nó liên kết hai lĩnh vực với nhau. KH&CN nano [123] dựa trên ba cơ sở là: (1) Sự chuyển tiếp từ đặc tính cổ điển sang đặc tính lượng tử.