Luận văn: Nghiên cứu mô hình hai trạng thái trong hệ GeO2 bằng mô phỏng - ĐH Sư phạm Thái Nguyên

Luận văn nghiên cứu mô hình hai trạng thái trong hệ GeO2 lỏng bằng phương pháp mô phỏng động lực học phân tử, phân tích vi cấu trúc và tính chất động học.

• Tác giả: Nammavong Thidakham
• Người hướng dẫn: TS. Luyện Thị San
• Trường học: Trường Đại học Sư phạm - Đại học Thái Nguyên
• Chuyên ngành: Vật lý chất rắn
• Đề tài: Nghiên cứu mô hình hai trạng thái trong hệ GeO2 bằng phương pháp mô phỏng
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2023
• Địa điểm: Thái Nguyên

Nội dung chính

Tổng quan nghiên cứu

Nghiên cứu cấu trúc vi mô của vật liệu ôxít đóng vai trò then chốt trong vật lý chất rắn và khoa học vật liệu hiện đại, đặc biệt đối với germani đioxit (GeO2) – một hệ vật liệu tương đồng quan trọng với silic đioxit nhưng thể hiện các biến đổi pha ở điều kiện áp suất ôn hòa hơn. Nghiên cứu thực hiện mô phỏng hệ 5499 nguyên tử trong dải áp suất nén cực hạn từ 0,1 GPa đến 119,6 GPa tại nhiệt độ 3500 K nhằm giải quyết bài toán còn nhiều tranh cãi về bản chất chuyển pha lỏng - lỏng và cơ chế vi mô của hiện tượng đa thù hình. Mục tiêu cốt lõi của công trình là làm sáng tỏ mô hình hai trạng thái cấu trúc, xác định sự biến đổi số phối trí và phân tích hiện tượng động học không đồng nhất diễn ra trong mạng lưới liên kết. Bằng cách thiết lập mối tương quan định lượng giữa các đơn vị cấu trúc cơ bản và tỷ phần liên kết ôxy cầu, công trình chứng minh sự cùng tồn tại của hai pha gồm pha mật độ thấp 3,62 g/cm³ và pha mật độ cao 6,89 g/cm³. Ý nghĩa khoa học của nghiên cứu thể hiện qua việc làm rõ bước chuyển đổi số phối trí trung bình Ge-O từ 3,7 lên 6,7, cung cấp bức tranh vật lý hoàn chỉnh về sự cạnh tranh không gian giữa các cấu hình tứ diện và bát diện dưới áp suất siêu cao, đồng thời tạo tiền đề cho việc tối ưu hóa quy trình tổng hợp vật liệu quang học chịu áp lực.

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ủa mô hình hai trạng thái cấu trúc vốn bắt nguồn từ các khám phá thực nghiệm về nước đá và silic vô định hình. Theo mô hình này, trạng thái lỏng của hệ ôxít tứ diện không phải là một pha đồng nhất duy nhất mà là sự đan xen giữa pha mật độ thấp, đặc trưng bởi các đơn vị tứ diện GeO4 nối qua liên kết hai phối trí OGe2 với khối lượng riêng xấp xỉ 3,62 g/cm³, và pha mật độ cao gồm các đơn vị GeO5, GeO6 nối qua liên kết ba phối trí OGe3 với khối lượng riêng đạt 6,89 g/cm³. Bên cạnh đó, nghiên cứu tích hợp lý thuyết tái cấu trúc phối hợp của Adam-Gibbs và lý thuyết thể tích tự do nhằm lý giải mối liên hệ giữa sự sụt giảm entropy cấu hình với tính linh động khuếch tán. Ba khái niệm trung tâm bao gồm: đơn vị cấu trúc cơ bản GeOx, liên kết cầu ôxy OGey và tham số trật tự chuyển pha η(P). Hàm phân bố góc liên kết O-Ge-O trong GeO4 đạt cực đại tại 105 độ, trong khi GeO5 và GeO6 xuất hiện hai đỉnh đặc trưng tại 90 độ và 165 độ, phản ánh sự biến đổi sâu sắc của trật tự tầm gần và trật tự tầm trung.

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

Nghiên cứu sử dụng phương pháp mô phỏng động lực học phân tử cổ điển với thế tương tác cặp Oeffner-Elliott, một dạng thế năng lượng kết hợp tương tác tĩnh điện Coulomb và hàm thế đẩy Born-Mayer-Huggins có độ tin cậy vượt trội khi so sánh với thực nghiệm. Cỡ mẫu nghiên cứu gồm 5499 nguyên tử, được cấu thành từ 1833 nguyên tử Ge và 3666 nguyên tử O nhằm đảm bảo đúng tỉ lệ hóa học hợp thức 1:2. Phương pháp chọn mẫu ban đầu được tiến hành bằng cách phân bố ngẫu nhiên các hạt vào khối lập phương sao cho khoảng cách giữa hai hạt bất kỳ luôn lớn hơn 1,0 Å, đạt mật độ chuẩn 3,5 g/cm³ tại 0 GPa. Lý do lựa chọn mô phỏng động lực học phân tử xuất phát từ việc các phương pháp thực nghiệm nhiễu xạ tia X hay tán xạ nơtron rất khó thực hiện và tốn kém ở nhiệt độ nóng chảy 3500 K và áp suất lên đến 120 GPa. Quy trình mô phỏng trải qua các mốc thời gian chặt chẽ: mẫu được nung đến 5000 K trong 500.000 bước tích phân với bước thời gian 0,48 fs, làm lạnh về 3500 K với tốc độ 1 K/ps, sau đó nén đẳng nhiệt đến 13 áp suất đích trong 1.000.000 bước và hồi phục cân bằng thêm 100.000 bước. Dữ liệu được trích xuất qua hàm phân bố xuyên tâm g(r), phân bố góc liên kết và kỹ thuật phân tích cụm link-cluster.

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

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

Thứ nhất, số phối trí trung bình của cặp Ge-O gia tăng đơn điệu và liên tục từ 3,7 tại áp suất 0,1 GPa lên mức 4,9 tại 14,5 GPa, rồi đạt 6,7 tại áp suất 119,6 GPa, bác bỏ hoàn toàn giả thuyết chuyển pha đột ngột loại một. Thứ hai, sự phân bố các đơn vị cấu trúc biến thiên mạnh mẽ: tỷ phần đơn vị tứ diện GeO4 giảm dốc từ 95% ở 0,1 GPa xuống 44% tại 14,5 GPa, chỉ còn 16% ở 18,6 GPa và triệt tiêu hoàn toàn khi áp suất vượt 79,1 GPa; ngược lại, đơn vị GeO6 tăng từ 0% lên 87% tại 119,6 GPa. Thứ ba, đơn vị phối trí năm GeO5 đóng vai trò trạng thái trung gian chủ đạo, tăng vọt từ 5% ở 0,1 GPa lên cực đại 43% tại 14,5 GPa trước khi giảm dần về 8% ở 119,6 GPa. Thứ tư, khoảng cách liên kết trung bình Ge-O thể hiện hành vi phi tuyến dị thường: dãn nở từ 1,77 Å ở 0,1 GPa lên giá trị cực đại 1,89 Å tại 18,6 GPa (tăng khoảng 6,8%), sau đó co ngắn lại mức 1,82 Å khi áp suất đạt 119,6 GPa do lực nén cơ học chiếm ưu thế.

Thảo luận kết quả

Sự tiến triển bất thường của độ dài liên kết Ge-O phản ánh cuộc cạnh tranh vật lý trực tiếp giữa hai cơ chế vi mô: sự gia tăng số phối trí trung bình làm dãn rộng vỏ cầu phối trí thứ nhất và hiệu ứng nén thủy tĩnh cưỡng bức ép các nguyên tử lại gần nhau. Kết quả tính toán hoàn toàn tương thích với các dữ liệu phổ hấp thụ và tán xạ tia X thực nghiệm của Spiekermann và Kono, khẳng định chuyển động trơn tru của quá trình tái sắp xếp cấu trúc.

Để minh họa trực quan, mối tương quan giữa khối lượng riêng và áp suất có thể được trình bày rõ nét qua biểu đồ đường cong hàm số hai trạng thái kết hợp bảng dữ liệu phân bố phối trí, trong đó đường cong mô phỏng bám sát hoàn hảo công thức lý thuyết kết hợp giữa tỷ phần liên kết OGe2 (giảm từ 85% xuống 16%) và OGe3 (tăng từ 8% lên 68%). Các lát cắt đồ họa không gian ba chiều cho thấy các liên kết OGe2 tụ tập tạo thành các đám vi mô mật độ thấp, trong khi OGe3 hình thành mạng lưới mật độ cao bao quanh. Động học không đồng nhất được xác nhận khi phân tích 300 hạt linh động nhất (nhóm A1) và 300 hạt kém linh động nhất (nhóm A2); các hạt chuyển động nhanh có xu hướng kết tụ thành 208 đến 282 cụm, phân tách không gian rõ rệt so với 302 đến 331 cụm của nhóm hạt ngẫu nhiên A0 tại áp suất 0 GPa.

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

Thứ nhất, hoàn thiện và chuẩn hóa bộ tham số thế tương tác Oeffner-Elliott cải biên trong các phần mềm tính toán vật liệu chuyên dụng; hướng đến mục tiêu giảm sai số dự báo cấu trúc tinh thể và vô định hình xuống dưới 2,0% trong chu kỳ 6 tháng, do các viện nghiên cứu vật lý tính toán và trường đại học chuyên ngành chủ trì thực hiện.

Thứ hai, tối ưu hóa công nghệ chế tạo thủy tinh quang học chịu áp lực cao bằng cách kiểm soát vùng chuyển pha trung gian từ 12 GPa đến 20 GPa; đặt chỉ tiêu gia tăng độ bền nén cơ học và khả năng chống nứt vỡ lên 25% trong vòng 12 tháng, do các kỹ sư trưởng tại trung tâm nghiên cứu và phát triển sản xuất sợi quang và kính kỹ thuật triển khai.

Thứ ba, xây dựng thuật toán phân tích Link-Cluster tự động tích hợp học máy để sàng lọc và định vị các vùng động học không đồng nhất; hướng tới mục tiêu nâng tốc độ xử lý quỹ đạo nguyên tử lên 40% với độ nhạy phân tách cụm đạt trên 95% trong thời gian 9 tháng, giao cho nhóm nghiên cứu công nghệ thông tin khoa học vật liệu.

Thứ tư, mở rộng mô hình nghiên cứu sang hệ ôxít hỗn hợp đa thành phần GeO2-SiO2 dưới dải nhiệt độ 2000 K đến 4000 K và áp suất 0 GPa đến 50 GPa; nhắm đến mục tiêu tối ưu hóa suy hao truyền dẫn cáp quang thế hệ mới giảm thêm 15% trong lộ trình 18 tháng, do các tập đoàn công nghệ viễn thông và phòng thí nghiệm vật liệu trọng điểm phối hợp tiến hành.

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

Thứ nhất, nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn, Hóa lý thuyết và Khoa học vật liệu. Trường hợp áp dụng: Kế thừa quy trình mô phỏng động lực học phân tử cho 14 mẫu áp suất từ 0,1 GPa đến 119,6 GPa và phương pháp phân tích hàm phân bố góc liên kết O-Ge-O nhằm hoàn thành các đề tài chuyển pha vô định hình.

Thứ hai, kỹ sư phát triển vật liệu tại các doanh nghiệp sản xuất linh kiện bán dẫn và thiết bị quang học cao cấp. Trường hợp áp dụng: Vận dụng số liệu thực nghiệm về tỷ phần đơn vị GeO4, GeO5 và độ dài liên kết để hiệu chỉnh quá trình ủ nhiệt và nén ép thủy tinh germani có chỉ số khúc xạ tối ưu trong dải áp suất 12 đến 20 GPa.

Thứ ba, nhà khoa học địa chất và vật lý địa cầu nghiên cứu cấu trúc sâu của vỏ Trái Đất. Trường hợp áp dụng: Sử dụng pha mật độ cao 6,89 g/cm³ và bước chuyển phối trí bát diện của GeO2 như một mô hình chất lỏng tương đương để mô phỏng dòng magma silicat ở manti sâu.

Thứ tư, chuyên gia phân tích dữ liệu mô phỏng và kỹ sư phần mềm khoa học tính toán. Trường hợp áp dụng: Ứng dụng thuật toán phân loại 5 nhóm tập hợp hạt vi mô từ A0 đến A4 và độ dịch chuyển bình phương trung bình vào việc phát triển các công cụ mô phỏng động học nguyên tử trên quy mô 5499 hạt.

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

Mô hình hai trạng thái trong hệ GeO2 lỏng giải thích sự biến đổi mật độ như thế nào? Mô hình hai trạng thái xem GeO2 lỏng là hỗn hợp của pha mật độ thấp 3,62 g/cm³ và pha mật độ cao 6,89 g/cm³. Khối lượng riêng tổng thể ở mỗi áp suất từ 0,1 GPa đến 119,6 GPa được xác định chính xác thông qua tỷ lệ liên kết ôxy cầu OGe2 và OGe3, cho thấy sự chuyển đổi trơn tru giữa hai pha mà không có sự đột biến thể tích gián đoạn.

Vì sao khoảng cách liên kết Ge-O lại tăng từ 1,77 Å lên 1,89 Å trước khi co rút ở áp suất cao? Hiện tượng này bắt nguồn từ sự cạnh tranh vi mô: từ 0,1 GPa đến 18,6 GPa, số phối trí Ge-O tăng mạnh từ 3,7 lên 5,3 làm vỏ phối trí nở rộng ra 1,89 Å. Khi áp suất vượt quá 20 GPa, số phối trí tăng chậm lại khiến tác động nén cơ học chiếm ưu thế tuyệt đối, ép khoảng cách liên kết co ngắn dần về 1,82 Å tại 119,6 GPa.

Hiện tượng động học không đồng nhất được lượng hóa bằng những chỉ số cụ thể nào? Hiện tượng này được định lượng qua độ dịch chuyển bình phương trung bình của 300 hạt linh động nhất thuộc nhóm A1 so với 300 hạt chậm nhất thuộc nhóm A2. Tại 0 GPa, các hạt nhóm A1 liên kết thành các cụm lớn có từ 11 đến 22 hạt, lớn hơn đáng kể so với cụm từ 5 đến 11 hạt của nhóm ngẫu nhiên A0.

Vai trò của đơn vị cấu trúc GeO5 trong việc làm sáng tỏ cơ chế chuyển pha là gì? Đơn vị GeO5 đóng vai trò cầu nối cấu trúc trung gian quyết định. Tỷ phần GeO5 tăng vọt từ 5% lên đỉnh 43% tại 14,5 GPa đã giải thích bản chất chuyển pha liên tục của hệ lỏng, thay thế cho quan điểm chuyển pha loại một đột ngột từng được giả định trong khoảng áp suất hẹp 7 đến 9 GPa.

Kết quả mô phỏng trong luận văn có độ tương đồng ra sao so với các dữ liệu thực nghiệm quốc tế? Các thông số mô phỏng từ 0,1 đến 119,6 GPa có độ tương thích cao với dữ liệu nhiễu xạ tia X của Spiekermann và Kono. Tỷ phần liên kết OGe2 giảm từ 85% xuống 16% và OGe3 tăng từ 8% lên 68% hoàn toàn khớp với quy luật biến thiên mật độ thực nghiệm, khẳng định độ tin cậy của mô hình tính toán.

Kết luận

• Xây dựng thành công 14 mô hình GeO2 lỏng ở nhiệt độ 3500 K trong dải áp suất rộng từ 0,1 GPa đến 119,6 GPa với cỡ mẫu chuẩn xác 5499 nguyên tử bằng phương pháp động lực học phân tử.
• Chứng minh tính xác thực của mô hình hai trạng thái cấu trúc với sự cùng tồn tại của pha mật độ thấp 3,62 g/cm³ và pha mật độ cao 6,89 g/cm³ thông qua mạng lưới liên kết OGe2 và OGe3.
• Khẳng định bản chất chuyển pha lỏng - lỏng diễn ra liên tục nhờ sự hiện diện trung gian của đơn vị phối trí năm GeO5 đạt cực đại 43% tại 14,5 GPa, loại bỏ giả thuyết chuyển pha đột ngột loại một.
• Phát hiện quy luật phi tuyến dị thường của độ dài liên kết Ge-O đạt giá trị cực đại 1,89 Å tại 18,6 GPa do sự cạnh tranh giữa gia tăng số phối trí từ 3,7 lên 6,7 và lực nén thủy tĩnh.
• Định lượng rõ nét hiện tượng động học không đồng nhất thông qua sự kết đám không gian của 300 hạt linh động cao trong toàn bộ tiến trình mô phỏng.

Đóng góp cốt lõi của công trình là thiết lập hệ phương trình toán học liên hệ trực tiếp giữa mật độ thực tế và tỷ phần vi cấu trúc ôxy cầu, mở ra hướng tiếp cận định lượng mới cho vật lý chất rắn. Kế hoạch tiếp theo trong 12 tháng tới tập trung mở rộng mô hình cho hệ pha tạp GeO2-SiO2 và kiểm chứng thực nghiệm phổ tán xạ tia X nhiệt độ cao. Các viện nghiên cứu và doanh nghiệp công nghệ quang điện tử có thể ứng dụng trực tiếp cơ sở dữ liệu vi cấu trúc này để nâng cao hiệu năng vật liệu thế hệ mới.

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

Chương 1 VẬT LIỆU Ô XÍT GeO2 Nội dung trong chương này, luận văn trình bày tổng quan các kết quả nghiên cứu về cấu trúc cũng như động học của hệ ôxít GeO2 bằng nhiều phương pháp thực nghiệm, lý thuyết và mô phỏng. Tổng quan về một số hệ ôxít Các hệ ôxít MxOy (M là Si, Mg, Al, Pb,.; x, y là các số nguyên) là nhóm các vật liệu có nhiều ứng dụng quan trọng trong khoa học công nghệ cũng như trong đời sống. Vì vậy việc cung cấp thêm các thông tin và đặc tính vi cấu trúc, các đặc tính về động học của hệ ôxít MxOy ở các nhiệt độ khác nhau và áp suất khác nhau là quan trọng và cần thiết. Cho đến nay, theo hiểu biết của chúng tôi, đã có nhiều nghiên cứu thực nghiệm (nhiễu xạ nơtrôn, nhiễu xạ tia-X, phổ Raman, cộng hưởng từ hạt nhân, phổ hấp thụ tia-X) và nhiều nghiên cứu mô phỏng (Thống kê hồi phục, động lực học phân tử, Monte Carlo,.) đã được sử dụng nhằm phát hiện các tính chất cấu trúc cũng như các đặc tính động học của hệ ôxít này.

Mặc dù các kết quả và số liệu nhận được từ nghiên cứu thực nghiệm và mô phỏng cho các hệ ôxít MxOy đã chỉ ra nhiều tính chất thú vị về vi cấu trúc và cơ chế chuyển động. Tuy nhiên, vẫn còn nhiều khía cạnh về cấu trúc và đặc tính động trong xảy ra trong các hệ ôxít MxOy chưa được hiểu biết đầy đủ và cần phải có thêm các nghiên cứu tiếp theo. Về mặt cấu trúc, các nghiên cứu thực nghiệm và mô phỏng thường sử dụng HPBXT, hàm phân bố góc, số phối trí, khoảng cách liên kết để cho thấy sự thay đổi cấu trúc của vật liệu khi áp suất hoặc nhiệt độ biến đổi. Từ đó, họ đã phát hiện các hiện tượng đa thù hình trong các hệ ôxít MxOy.

Về mặt động học, một số hiện tượng thường được tập trung nghiên cứu như cơ chế khuếch tán dị thường, động học không đồng nhất, tính chất dẫn điện và hệ số nhớt. Các cơ chế khuếch tán dị thường được quan sát trong một số chất lỏng cấu trúc mạng tứ diện như MgO, H2O, Al2O3, SiO2 và một số hơp chất của các ôxít này. Khuếch tán dị 4 thường là hiện tượng hệ số khuếch tán tăng khi áp suất tăng (trái với qui luật thông thường). Ngoài khuếch tán dị thường, các chất lỏng có cấu trúc mạng tứ diện còn biểu hiện các tính dị thường trong các đặc trưng nhiệt động học.

Trong vật liệu ôxít tồn tại các vùng hạt chuyển động nhanh hơn hoặc các vùng hạt chuyển động chậm hơn so với những vùng hạt khác trong hệ được gọi là động học không đồng nhất. Các vùng hạt chuyển động nhanh và các vùng hạt chuyển động chậm có thể biến đổi trong không gian theo thời gian. Sự phân bố tỉ phần các đơn vị Qn trong vật liệu thủy tinh (MgO)x-(SiO2)100-x Với các hệ ôxít đa thành phần, ví dụ như các hệ silicate (MgO-SiO2, PbO- SiO2, Al2O3-SiO2, Na2O-SiO2, .), sự biến đổi áp suất và nhiệt độ của hệ cũng dẫn tới nhiều hiện tượng thú vị. Ví dụ như, mức độ polymer hóa của SiO2 ảnh hưởng đến các đặc trưng của vật liệu.

Mức độ polymer hóa của các hệ SiO2 phân tích thông qua số lượng của các đại lượng Qn, trong đó Q là ĐVCT SiO4, n là số hạt O liên kết cầu trong ĐVCT SiO4. Ôxy cầu kí hiệu là BO, là những O liên kết đồng thời với 2 hạt Si. Sự thay đổi hàm lượng Si, áp suất và nhiệt độ làm thay đổi sự phân bố của các ĐVCT Qn.