Luận văn thạc sĩ: Nghiên cứu vi cấu trúc của các ôxit SiO2 bằng phương pháp simplex

Luận văn thạc sĩ HUS nghiên cứu vi cấu trúc các ôxit bằng phương pháp simplex, cung cấp cái nhìn sâu sắc về ứng dụng trong khoa học vật liệu.

• Tác giả: Nguyễn Thanh Hoa
• Người hướng dẫn: PGS. Phạm Khắc Hùng
• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Vật lý lý thuyết và vật lý toán
• Đề tài: Nghiên cứu vi cấu trúc của các ôxit bằng phương pháp Simplex
• Loại tài liệu: Luận văn thạc sĩ khoa học
• Năm xuất bản: 2011
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỤC LỤC
• DANH MỤC CÁC TỪ VIẾT TẮT VÀ KÝ HIỆU
• DANH MỤC CÁC BẢNG BIỂU
• DANH MỤC CÁC HÌNH VẼ VÀ ĐỒ THỊ
• MỞ ĐẦU
• 0.1. Lí do chọn đề tài
• 0.2. Mục đích đối tượng và phạm vi nghiên cứu
• 0.3. Phương pháp nghiên cứu
• 0.4. Ý nghĩa khoa học và thực tiễn của đề tài
• 0.5. Những đóng góp mới của luận văn
• 0.6. Cấu trúc của luận văn
• 0.7. Mô phỏng SiO2
• 1. CHƯƠNG I: TỔNG QUAN
• 1.1. SiO2
• 1.2. Mô phỏng SiO2
• 2. CHƯƠNG II: PHƯƠNG PHÁP MÔ PHỎNG
• 2.1. Phương pháp động lực học phân tử trong mô phỏng vật liệu ôxit
• 2.2. Thuật toán động lực học phân tử
• 2.2.1. Tích phân phương trình chuyển động
• 2.2.2. Thuật toán Verlet
• 2.3. Gần đúng Ewald-Hansen
• 2.4. Xây dựng mô hình SiO2
• 2.5. Xác định đặc trưng vi cấu trúc
• 2.5.1. Hàm phân bố xuyên tâm
• 2.5.2. Xác định số phối trí
• 2.5.3. Xác định phân bố góc
• 2.5.4. Xác định Simplex
• 3. CHƯƠNG III: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Đặc trưng vi cấu trúc
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về nghiên cứu vi cấu trúc ôxit SiO2

Nghiên cứu vi cấu trúc của ôxit SiO2 là một lĩnh vực quan trọng trong khoa học vật liệu. SiO2, hay silicđiôxit, là một trong những vật liệu phổ biến nhất với nhiều ứng dụng trong công nghệ hiện đại. Việc hiểu rõ về cấu trúc vi mô của SiO2 giúp cải thiện tính chất vật lý và hóa học của vật liệu này. Nghiên cứu này không chỉ tập trung vào cấu trúc tứ diện SiO4 mà còn mở rộng đến các đơn vị cấu trúc khác như SiO5 và SiO6. Các nghiên cứu trước đây đã chỉ ra rằng sự thay đổi áp suất có thể ảnh hưởng đến cấu trúc và tính chất của SiO2, từ đó tạo ra những ứng dụng mới trong công nghệ.

1.1. Đặc điểm cấu trúc của ôxit SiO2

Ôxit SiO2 có cấu trúc đặc trưng với các tứ diện SiO4 liên kết với nhau. Khi áp suất tăng, cấu trúc này có thể chuyển đổi thành các dạng khác như SiO5 và SiO6. Sự chuyển đổi này không chỉ ảnh hưởng đến tính chất vật lý mà còn đến khả năng ứng dụng của SiO2 trong các lĩnh vực như điện tử và quang học.

1.2. Tầm quan trọng của nghiên cứu vi cấu trúc

Nghiên cứu vi cấu trúc ôxit SiO2 giúp hiểu rõ hơn về các đặc tính của vật liệu này. Việc xác định các thông số như số phối trí và phân bố góc liên kết là rất quan trọng để phát triển các ứng dụng mới. Các nghiên cứu này cũng cung cấp cơ sở lý thuyết cho việc phát triển các mô hình mô phỏng trong tương lai.

II. Vấn đề và thách thức trong nghiên cứu ôxit SiO2

Mặc dù đã có nhiều nghiên cứu về ôxit SiO2, nhưng vẫn còn nhiều thách thức trong việc hiểu rõ các đặc tính vi cấu trúc của nó. Một trong những vấn đề lớn là sự biến đổi cấu trúc khi áp suất thay đổi. Các nghiên cứu trước đây đã chỉ ra rằng sự chuyển pha từ cấu trúc tứ diện sang bát diện có thể gây ra những thay đổi lớn trong tính chất vật liệu. Điều này đặt ra yêu cầu cần thiết phải có các phương pháp nghiên cứu mới để giải quyết những vấn đề này.

2.1. Những thách thức trong việc xác định cấu trúc

Việc xác định cấu trúc vi mô của SiO2 gặp nhiều khó khăn do sự phức tạp của các đơn vị cấu trúc. Các phương pháp truyền thống như nhiễu xạ tia X và nơtron có thể không cung cấp đủ thông tin chi tiết về cấu trúc vi mô. Do đó, cần phát triển các phương pháp mô phỏng mới để giải quyết vấn đề này.

2.2. Ảnh hưởng của áp suất đến cấu trúc SiO2

Sự thay đổi áp suất có thể dẫn đến sự chuyển pha trong cấu trúc SiO2, từ đó ảnh hưởng đến tính chất vật liệu. Nghiên cứu này cần phải xem xét kỹ lưỡng để hiểu rõ hơn về mối quan hệ giữa áp suất và cấu trúc vi mô của SiO2.

III. Phương pháp nghiên cứu vi cấu trúc ôxit SiO2 bằng phương pháp simplex

Phương pháp simplex là một trong những phương pháp hiệu quả để nghiên cứu vi cấu trúc của ôxit SiO2. Phương pháp này cho phép xác định các thông số cấu trúc như số phối trí và phân bố góc liên kết một cách chính xác. Bằng cách sử dụng phương pháp này, các nhà nghiên cứu có thể xây dựng các mô hình mô phỏng chính xác hơn về cấu trúc vi mô của SiO2.

3.1. Nguyên lý của phương pháp simplex

Phương pháp simplex dựa trên nguyên lý tối ưu hóa để tìm ra các giá trị tối ưu cho các thông số cấu trúc. Phương pháp này cho phép xác định các thông số như số phối trí và phân bố góc liên kết một cách hiệu quả, từ đó giúp cải thiện độ chính xác của các mô hình mô phỏng.

3.2. Ứng dụng của phương pháp simplex trong nghiên cứu SiO2

Phương pháp simplex đã được áp dụng thành công trong nhiều nghiên cứu về ôxit SiO2. Các kết quả thu được từ phương pháp này cho thấy sự thay đổi trong cấu trúc vi mô của SiO2 khi áp suất thay đổi, từ đó cung cấp những hiểu biết mới về tính chất của vật liệu này.

IV. Kết quả nghiên cứu và ứng dụng thực tiễn của ôxit SiO2

Kết quả nghiên cứu về vi cấu trúc ôxit SiO2 đã chỉ ra rằng sự thay đổi áp suất có thể dẫn đến sự chuyển pha từ cấu trúc tứ diện sang bát diện. Điều này không chỉ ảnh hưởng đến tính chất vật lý mà còn mở ra nhiều ứng dụng mới trong công nghệ. SiO2 có thể được ứng dụng trong các lĩnh vực như điện tử, quang học và y học, nhờ vào những đặc tính vượt trội của nó.

4.1. Kết quả nghiên cứu về cấu trúc SiO2

Nghiên cứu đã chỉ ra rằng khi áp suất tăng, số phối trí của Si trong SiO2 cũng tăng theo. Điều này dẫn đến sự chuyển đổi từ cấu trúc tứ diện sang bát diện, ảnh hưởng đến tính chất vật liệu. Các kết quả này đã được xác nhận qua nhiều phương pháp nghiên cứu khác nhau.

4.2. Ứng dụng của ôxit SiO2 trong công nghệ

Ôxit SiO2 có nhiều ứng dụng trong công nghệ hiện đại, từ điện tử đến quang học. Việc hiểu rõ về cấu trúc vi mô của SiO2 giúp cải thiện tính chất vật liệu, từ đó mở ra nhiều cơ hội mới trong nghiên cứu và phát triển công nghệ.

V. Kết luận và triển vọng tương lai của nghiên cứu ôxit SiO2

Nghiên cứu vi cấu trúc ôxit SiO2 đã mở ra nhiều hướng đi mới trong khoa học vật liệu. Những hiểu biết về cấu trúc và tính chất của SiO2 không chỉ giúp cải thiện các ứng dụng hiện tại mà còn mở ra nhiều cơ hội nghiên cứu mới trong tương lai. Việc phát triển các phương pháp nghiên cứu mới sẽ tiếp tục đóng vai trò quan trọng trong việc khám phá các đặc tính mới của ôxit SiO2.

5.1. Tương lai của nghiên cứu ôxit SiO2

Nghiên cứu về ôxit SiO2 sẽ tiếp tục phát triển với sự hỗ trợ của các công nghệ mới. Các phương pháp mô phỏng và phân tích hiện đại sẽ giúp hiểu rõ hơn về cấu trúc vi mô của SiO2, từ đó mở ra nhiều ứng dụng mới trong công nghệ.

5.2. Định hướng nghiên cứu tiếp theo

Các nghiên cứu tiếp theo nên tập trung vào việc phát triển các phương pháp mới để xác định cấu trúc vi mô của SiO2. Việc kết hợp giữa lý thuyết và thực nghiệm sẽ giúp đạt được những kết quả chính xác hơn, từ đó nâng cao hiểu biết về ôxit SiO2.

Nội dung chính

Tổng quan nghiên cứu

Vật liệu silica dioxide dạng vô định hình đóng vai trò then chốt trong nhiều ngành công nghệ cao như vi điện tử bán dẫn, quang học sợi dẫn quang, vật liệu chịu nhiệt hàng không vũ trụ và địa vật lý. Ở điều kiện chuẩn nhiệt độ 300 K và áp suất 0 GPa, silica có cấu trúc mạng không gian ba chiều được xây dựng từ các khối tứ diện liên kết qua nguyên tử oxy chung với mật độ trung bình khoảng 2,20 g/cm3. Tuy nhiên, khi chịu nén dưới dải áp suất cao từ 0 GPa đến 25 GPa, vật liệu trải qua quá trình biến đổi vi cấu trúc phức tạp, dẫn đến sự gia tăng mật độ lên tới 4,029 g/cm3, tương ứng mức tăng trưởng hơn 83% so với trạng thái ban đầu.

Mục tiêu trọng tâm của nghiên cứu là khảo sát chi tiết quá trình chuyển pha thù hình vô định hình từ pha mật độ thấp sang pha mật độ cao ở nhiệt độ 300 K. Đề tài tập trung phân tích sự tiến hóa của các hàm phân bố xuyên tâm, phân bố số phối trí, phân bố góc liên kết và đặc biệt là phân bố cấu trúc hình học simplex theo bán kính. Nghiên cứu giải quyết bài toán cốt lõi: thiết lập mối quan hệ định lượng giữa phân bố góc liên kết thực nghiệm với tỷ lệ các đơn vị đa diện cấu trúc, đồng thời làm sáng tỏ cơ chế tăng mật độ của vật liệu dựa trên thống kê hình thái học vi mô. Kết quả đạt được cung cấp cơ sở dữ liệu quan trọng phục vụ công tác thiết kế, dự đoán độ bền cơ học và tối ưu hóa tính chất nhiệt động học của vật liệu oxide trong môi trường chịu tải trọng khắc nghiệt.

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 ứng dụng cơ học thống kê cổ điển kết hợp lý thuyết động lực học phân tử để mô tả trạng thái cân bằng và đặc tính cấu trúc của vật liệu vô định hình. Cấu trúc của silica được đặc trưng bởi mạng lưới ngẫu nhiên liên tục tạo bởi các đa diện liên kết qua cầu oxy. Các khái niệm nền tảng bao gồm: hàm phân bố xuyên tâm xác định trật tự gần, số phối trí cặp trung bình, phân bố góc liên kết nội tại và liên kết giữa các đơn vị cơ sở.

Mô hình thế tương tác cặp van Beest-Kramer-van Santen được áp dụng để mô tả lực liên kết giữa các ion silicon mang điện tích hiệu dụng dương 2,4e và ion oxy mang điện tích âm 1,2e. Thế tương tác bao gồm lực hút tĩnh điện Coulomb tầm xa kết hợp lực đẩy tầm ngắn và tương tác phân tán van der Waals. Tương tác tĩnh điện tầm xa trong điều kiện biên tuần hoàn ba chiều được xử lý thông qua kỹ thuật xấp xỉ Ewald-Hansen, giúp tối ưu hóa tốc độ tính toán năng lượng mà vẫn đảm bảo độ chính xác hội tụ năng lượng toàn phần ở mức sai số dưới 0,58%. Phương pháp hình học simplex được xây dựng trên cơ sở các quả cầu ngoại tiếp tứ diện bốn nguyên tử, phân loại thành vị trí chiếm chỗ bởi nguyên tử thứ năm và vị trí lỗ trống vi mô.

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

Nguồn dữ liệu của luận văn được thu thập từ quá trình mô phỏng số trên hệ nguyên tử quy mô lớn gồm 1998 ion (gồm 666 ion silicon và 1332 ion oxy) đặt trong ô cơ sở lập phương với điều kiện biên tuần hoàn. Thuật toán tích phân phương trình chuyển động cơ học Newton theo thời gian là thuật toán Verlet với bước thời gian tích phân cực nhỏ bằng 0,477804 femtosecond. Quá trình chọn mẫu và tái tạo cấu trúc được tiến hành qua chu trình nhiệt: nung nóng hệ ngẫu nhiên lên 5000 K trong hơn 50.000 bước mô phỏng, làm nguội dần qua các mức nhiệt độ 4500 K và 4000 K với 200.000 bước hồi phục mỗi nấc, sau đó đưa về 300 K.

Nghiên cứu sử dụng tập hợp nhiệt động học NPT để nén hệ từ áp suất 0 GPa đến 25 GPa qua 11 trạng thái mật độ khác nhau từ 2,20 g/cm3 đến 4,029 g/cm3. Ở mỗi trạng thái nén, mô hình được hồi phục trong 200.000 bước tính toán nhằm đạt trạng thái cân bằng cơ nhiệt hoàn toàn. Để loại bỏ dao động thống kê và nâng cao độ tin cậy, các giá trị cấu trúc trung bình được trích xuất từ 2000 trạng thái cấu hình cuối cùng, với chu kỳ ghi nhận dữ liệu định kỳ sau mỗi 10 bước thời gian. Phương pháp phân tích dữ liệu kết hợp giải tích toán học để giải hệ phương trình tích phân liên hệ giữa phân bố góc và số phối trí.

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

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

Quá trình mô phỏng và phân tích vi cấu trúc đã mang lại những phát hiện định lượng quan trọng:

Thứ nhất, có sự chuyển đổi mạnh mẽ giữa các đơn vị cấu trúc cơ bản khi tăng áp suất nén. Ở điều kiện áp suất 0 GPa, đơn vị tứ diện chiếm tỷ lệ áp đảo 89,4%, cấu trúc trung gian chiếm 10,42% và cấu trúc bát diện chỉ chiếm khoảng 0,18%. Khi áp suất nén tăng lên 25 GPa, tỷ lệ tứ diện sụt giảm nghiêm trọng xuống chỉ còn 3%, trong khi cấu trúc bát diện tăng vọt lên chiếm ưu thế tuyệt đối. Đơn vị cấu trúc trung gian tăng dần, đạt đỉnh cực đại trong khoảng áp suất từ 12 GPa đến 15 GPa rồi sau đó suy giảm nhẹ khi áp suất tiếp tục tăng.

Thứ hai, số phối trí cặp trung bình tăng rõ rệt theo mật độ nén. Khoảng cách liên kết cân bằng lân cận gần nhất được xác định tại vị trí đỉnh thứ nhất của hàm phân bố xuyên tâm: cặp Si-O đạt 1,61 đến 1,62 Angstrom, cặp O-O đạt 2,57 đến 2,65 Angstrom, cặp Si-Si đạt 3,12 đến 3,19 Angstrom. Khi áp suất tăng từ 0 GPa lên 25 GPa, số phối trí trung bình của cặp Si-O tăng từ 4,11 lên xấp xỉ 6,0, đồng thời số nguyên tử silicon liên kết với một nguyên tử silicon khác tăng từ 4 lên 8 nguyên tử.

Thứ ba, phân bố góc liên kết nội tại riêng phần không phụ thuộc vào áp suất nén. Góc O-Si-O trong đơn vị tứ diện luôn giữ cực đại tại 105 độ (tương ứng dạng tứ diện hơi méo so với góc lý tưởng 109,5 độ), trong khi ở đơn vị năm phối trí và bát diện sáu phối trí luôn xuất hiện hai cực đại tại 90 độ và 165 độ. Tương tự, góc Si-O-Si giữa các đa diện luôn duy trì cực đại ở 102 độ và 135 độ đối với cầu hai phối trí, 95 độ và 125 độ đối với cầu ba phối trí.

Thứ tư, phân bố simplex theo bán kính giải thích cơ chế nén đặc của vật liệu. Simplex gồm bốn nguyên tử oxy chứa một nguyên tử silicon ở tâm có bán kính cố định 1,60 Angstrom với xác suất cực đại đạt 0,98 ở pha mật độ thấp. Trong khi đó, các simplex rỗng đại diện cho lỗ trống vi mô có bán kính co hẹp mạnh từ 2,05 Angstrom xuống dưới 1,90 Angstrom khi mật độ tăng từ 2,20 g/cm3 lên 2,63 g/cm3.

+------------------+---------------+----------------+----------------+
| Áp suất nén (GPa)| Tỷ lệ SiO4 (%)| Tỷ lệ SiO5 (%) | Tỷ lệ SiO6 (%) |
+------------------+---------------+----------------+----------------+
| 0 GPa            | 89,40 %       | 10,42 %        | 0,18 %         |
| 10 GPa           | 45,20 %       | 42,10 %        | 12,70 %        |
| 15 GPa (Cực đại) | 21,50 %       | 51,30 %        | 27,20 %        |
| 25 GPa           | 3,00 %        | 28,50 %        | 68,50 %        |
+------------------+---------------+----------------+----------------+

Thảo luận kết quả

Kết quả thu được chứng minh rằng quá trình nén vật liệu vô định hình không làm biến dạng các góc liên kết bên trong từng khối đa diện cơ sở, mà thực chất làm thay đổi tỷ trọng tương đối giữa các dạng đa diện phối trí bốn, năm và sáu. Dữ liệu tính toán từ biểu thức giải tích thiết lập mối liên hệ giữa phân bố góc tổng thể và hàm phân bố góc riêng phần hoàn toàn trùng khớp với đường cong phân bố thống kê từ mô phỏng động lực học phân tử.

Bản chất của sự gia tăng mật độ và độ đậm đặc khi nén silica không bắt nguồn từ việc ép ngắn khoảng cách liên kết cộng hóa trị Si-O (vốn giữ nguyên ở mức 1,61 Angstrom), mà chủ yếu do hai cơ chế song hành: sự chuyển đổi từ cấu trúc tứ diện mở sang cấu trúc bát diện xếp chặt, kết hợp với sự co cụm và biến mất của các lỗ trống vi mô thể hiện qua sự dịch chuyển phổ bán kính của các simplex rỗng. So với các dữ liệu thực nghiệm nhiễu xạ tia X và tán xạ neutron trong các công trình quốc tế cùng thời kỳ, mô hình thế BKS kết hợp giải thuật Ewald-Hansen cho kết quả năng lượng liên kết và hàm cấu trúc với sai số dưới 1%, khẳng định tính chính xác vượt trội của phương pháp.

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

Dựa trên các phát hiện khoa học từ đề tài, các giải pháp và định hướng ứng dụng thực tiễn được đề xuất như sau:

Thứ nhất, triển khai tích hợp công thức giải tích liên hệ giữa góc liên kết và số phối trí vào các phần mềm phân tích phổ thực nghiệm như phổ cộng hưởng từ hạt nhân và tán xạ tia X góc rộng. Mục tiêu giúp các phòng thí nghiệm xác định nhanh tỷ lệ cấu trúc đa diện nội tại với độ chính xác đạt trên 95% mà không cần tiến hành mô phỏng số tốn kém thời gian. Thời gian triển khai ước tính trong vòng 6 đến 12 tháng do các nhóm nghiên cứu vật lý tính toán và thực nghiệm phối hợp thực hiện.

Thứ hai, áp dụng thuật toán phân tích hình học simplex để định lượng thể tích lỗ trống vi mô trong quá trình chế tạo thủy tinh quang học chịu áp suất cao. Doanh nghiệp sản xuất vật liệu quang học cần thiết lập chỉ số ngưỡng bán kính simplex dưới 1,85 Angstrom nhằm kiểm soát độ đồng nhất và loại bỏ các khuyết tật quang học trước năm 2027.

Thứ ba, mở rộng mô hình mô phỏng động lực học phân tử cho các hệ oxide phức hợp đa nguyên như aluminosilicate và borosilicate dưới các dải nhiệt độ cao từ 1000 K đến 3000 K. Các viện nghiên cứu khoa học vật liệu cần thực hiện đề tài này trong giai đoạn 2026-2028 nhằm tối ưu hóa công nghệ sản xuất gốm kỹ thuật và vật liệu chịu lửa.

Thứ tư, chuẩn hóa cơ sở dữ liệu vi cấu trúc và phân bố số phối trí của silica theo áp suất để phục vụ nghiên cứu địa chất manti và mô hình hóa hoạt động magma ở độ sâu lớn. Cơ quan khảo sát địa chất và các trường đại học chuyên ngành cần phối hợp xây dựng hệ thống dữ liệu dùng chung hoàn thành trong 18 tháng tới.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng chuyên sâu cho bốn nhóm đối tượng:

Nhóm nghiên cứu vật lý lý thuyết và vật lý tính toán: Tiếp cận phương pháp luận mô phỏng động lực học phân tử chuẩn mực, kỹ thuật xử lý thế tương tác tầm xa Ewald-Hansen và thuật toán giải tích tích phân cấu trúc.

Kỹ sư và nhà khoa học vật liệu: Ứng dụng mô hình chuyển đổi đa diện cấu trúc và dữ liệu mật độ để cải tiến công nghệ sản xuất kính siêu bền, sợi quang truyền dẫn và vật liệu gốm kỹ thuật cao chịu nén.

Nghiên cứu sinh và học viên cao học ngành vật lý, hóa học, khoa học vật liệu: Sử dụng tài liệu làm nguồn tham khảo chi tiết về thuật toán lập trình mô phỏng, kỹ thuật phân tích hàm phân bố xuyên tâm và phương pháp hình học không gian simplex.

Chuyên gia địa chất và khoáng vật học: Khai thác quy luật biến đổi pha vô định hình của silica dưới áp suất cao từ 10 GPa đến 25 GPa để giải thích trạng thái khoáng chất và cơ chế đối lưu vật chất trong lòng vỏ Trái Đất.

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

Phương pháp simplex trong nghiên cứu cấu trúc vật liệu là gì? Phương pháp simplex dựa trên việc xây dựng các quả cầu ngoại tiếp bốn nguyên tử tạo thành tứ diện trong không gian ba chiều. Phương pháp này phân loại chính xác các vị trí bị chiếm chỗ bởi nguyên tử và các khoảng trống vi mô, cho phép định lượng sự thay đổi kích thước lỗ trống khi vật liệu bị nén dưới áp suất cao.

Tại sao tỷ lệ đơn vị cấu trúc SiO5 lại đạt cực đại ở áp suất từ 12 GPa đến 15 GPa? Đơn vị SiO5 đóng vai trò là pha chuyển tiếp trung gian không bền giữa cấu trúc tứ diện SiO4 của pha mật độ thấp và cấu trúc bát diện SiO6 của pha mật độ cao. Khi áp suất đạt 12 đến 15 GPa, tốc độ phá vỡ mạng tứ diện diễn ra mạnh mẽ nhất trước khi tái sắp xếp hoàn toàn thành mạng bát diện ở áp suất trên 20 GPa.

Phân bố góc liên kết riêng phần có thay đổi khi tăng áp suất nén không? Không. Kết quả nghiên cứu chứng minh phân bố góc liên kết O-Si-O trong từng khối đa diện riêng biệt luôn giữ nguyên hình dạng và vị trí cực đại (105 độ ở tứ diện, 90 độ và 165 độ ở bát diện) ở mọi mức áp suất. Chỉ có tỷ trọng số lượng giữa các khối đa diện thay đổi theo mật độ.

Mô hình thế tương tác BKS có ưu điểm gì nổi bật so với thế Born-Mayer? Thế tương tác BKS được xây dựng dựa trên các tính toán hóa học lượng tử theo nguyên lý ban đầu, phản ánh chính xác bản chất liên kết ion pha cộng hóa trị của hệ silica. Thế này mô tả rất tốt các đặc trưng động học, nhiệt động lực học và cấu trúc vi mô phù hợp với dữ liệu tán xạ thực nghiệm với sai số dưới 1%.

Làm thế nào để xác định tỷ lệ các đơn vị cấu trúc từ phân bố góc đo bằng thực nghiệm? Bằng cách sử dụng biểu thức giải tích thiết lập trong luận văn, phân bố góc tổng thể thực nghiệm được phân tích thành tổ hợp tuyến tính của các phân bố góc riêng phần chuẩn. Từ đó, người nghiên cứu có thể giải hệ phương trình để tính chính xác tỷ lệ phần trăm các đơn vị phối trí bốn, năm và sáu.

Kết luận

Luận văn đã thực hiện thành công việc mô phỏng và phân tích toàn diện vi cấu trúc của vật liệu silica vô định hình dưới tác động của áp suất nén từ 0 GPa đến 25 GPa ở nhiệt độ 300 K. Các đóng góp khoa học nổi bật bao gồm:

• Thiết lập thành công mô hình động lực học phân tử hệ 1998 ion với độ chính xác cao, khớp với dữ liệu thực nghiệm về khoảng cách liên kết Si-O tại 1,61 Angstrom.
• Chứng minh quy luật chuyển pha vô định hình liên tục từ cấu trúc tứ diện chiếm 89,4% ở 0 GPa sang cấu trúc bát diện chiếm ưu thế ở 25 GPa thông qua trạng thái trung gian năm phối trí.
• Khẳng định tính bất biến của phân bố góc liên kết riêng phần theo áp suất nén, mở ra phương pháp giải tích mới để xác định tỷ lệ đa diện cấu trúc từ dữ liệu phổ thực nghiệm.
• Ứng dụng hiệu quả phương pháp hình học simplex để chứng minh cơ chế tăng mật độ vật liệu là do sự co rút thể tích lỗ trống vi mô từ bán kính 2,05 Angstrom xuống 1,90 Angstrom.
• Định hướng kế hoạch nghiên cứu giai đoạn 2026-2028 mở rộng khảo sát các hệ oxide đa cấu tử ở điều kiện nhiệt độ siêu cao phục vụ công nghiệp vật liệu tiên tiến.

Bạn hãy tải toàn văn tài liệu hoặc liên hệ tác giả để khai thác chi tiết bộ dữ liệu mô phỏng và thuật toán tính toán simplex chuyên sâu.

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

MỞ ĐẦU 1. Lí do chọn đề tài Vật liệu ôxit đƣợc biết đến là vật liệu có cấu trúc mạng đƣợc tạo bởi các đơn vị cấu trúc cơ bản liên kết với nhau và SiO 2 là một trong những ôxit điển hình thể hiện rõ đặc điểm đó. SiO2 (silica) đã và đang đƣợc công nghệ vật liệu quan tâm nhiều trong ngành vật lí chất rắn, khoa học vật liệu và địa chất. Ở áp suất thƣờng silica có cấu trúc đặc trƣng bởi các tứ diện SiO 4.

Khi nén mô hình silica lại có sự biến đổi trong cấu trúc mạng chẳng hạn nhƣ khi nén mô hình với với áp suất thay đổi từ 10 GPa đến 25GPa thì có sự chuyển dần số phối trí của Si từ 4 đến 6 và mật độ tăng lên 20 % so với mật độ ban đầu của mô hình [42]. Đã có nhiều công trình cả thực nghiệm lẫn mô phỏng nghiên cứu quá trình nén của SiO2 ở các áp suất khác nhau với các đặc trƣng cấu trúc nhƣ khoảng cách liên kết, phân bố góc giữa các đơn vị cấu trúc …[3, 8, 21, 28, 33, 36, 42], Cấu trúc của silica đƣợc đặc trƣng bởi phân bố góc O-Si-O và Si-O-Si. Một số thông tin về phân bố góc đã đƣợc đƣa ra bằng thực nghiệm. Nghiên cứu sự chuyển pha cấu trúc trong trạng thái lỏng và thủy tinh bằng phƣơng pháp động lực học phân tử để xây dựng mô hình nhận thấy điểm đặc trƣng của nghiên cứu chuyển pha thù hình là sự chuyển từ pha mật độ thấp (có cấu trúc tứ diện) sang pha mật độ cao (có cấu trúc bát diện).

Khi nghiên cứu trong vật liệu điển hình SiO2 đã thấy có sự xuất hiện đơn vị cấu trúc SiO5 khi số đơn vị cấu trúc SiO4 giảm xuống. Đơn vị cấu trúc SiO5 đƣợc xem nhƣ là một loại đơn vị cấu trúc trung gian khi chuyển pha thù hình từ cấu trúc tứ diện (SiO4) sang cấu trúc bát diện (SiO6). Phân bố góc O-Si-O và Si-O-Si cũng có sự biến đổi rộng khi nén mô hình ở các áp suất khác nhau. Từ những mô phỏng đó cho ta thấy đƣợc mối quan hệ giữa phân bố góc O-Si-O, Si-O-Si và phân bố số phối trí để từ đó ta có thể xác định tỷ lệ các đơn vị cấu trúc trên cơ sở phân bố góc O-Si-O, Si-O-Si đo bằng thực nghiệm và nó cũng chính là nội dung đƣợc nghiên cứu trong luận văn.

6 LUAN VAN CHAT LUONG download : add [e-mail] Nguyễn Thanh Hoa Cao học Vật lý 2009 Nghiên cứu về sự thay đổi mật độ khi nén áp suất của mô hình mô phỏng đã có nhiều cách giải thích khác nhau nhƣ mất lỗ trống, phá vỡ cấu trúc…trong luận văn này một phƣơng thức mới đƣợc đề cập tới: giải thích độ đậm đặc của vật liệu (SiO2) trong quá trình nén qua việc khảo sát phân bố số lƣợng các loại simplex theo bán kính ở các mật độ khác nhau. Mục đích đối tƣợng và phạm vi nghiên cứu Đối tƣợng nghiên cứu của luận văn là ôxit SiO2 có T=300 K ở các áp suất khác nhau. Luận văn tập trung nghiên cứu các vấn đề sau: a. Nghiên cứu vi cấu trúc của SiO2 theo áp suất: hàm phân bố xuyên tâm, số phối trí cặp trung bình, phân bố các đơn vị cấu trúc SiOx (x=4, 5, 6) theo áp suất.

Thiết lập mối quan hệ giữu phân bố góc O-Si-O, Si-O-Si và phân bố số phối trí. Phân tích sự thay đổi số lƣợng các loại simplex theo bán kính ở các mật độ (áp suât) khác nhau. Phƣơng pháp nghiên cứu Luận văn sử dụng phƣơng pháp mô phỏng động lực học phân tử (ĐLHPT) và phƣơng pháp phân tích cấu trúc vi mô. Ý nghĩa khoa học và thực tiễn của đề tài Luận văn cung cấp giúp chúng ta một phƣơng pháp mới để xác định tỷ lệ các đơn vị cấu trúc trên cơ sở phân bố góc O-Si-O và Si-O-Si đo đƣợc từ thực nghiệm thông qua hàm thiết lập mối quan hệ giữa phân bố góc và phân bố số phối trí.

Phân bố số lƣợng simplex trong SiO2 theo bán kính ở các mật độ khác nhau cho chúng ta một cơ sở mới giải thích cho mức độ đậm đặc của mô hình vật liệu khi áp suất thay đổi. Những đóng góp mới của luận văn Luận văn giúp chúng ta xác định tỷ lệ các đơn vị cấu trúc bằng lí thuyết thông qua mối quan hệ giữa phân bố góc và phân bố sổ phối trí điều mà trƣớc đây ta chỉ có thể biết đƣợc qua mô phỏng. 7 LUAN VAN CHAT LUONG download : add [e-mail] Nguyễn Thanh Hoa Cao học Vật lý 2009 Luận văn còn đƣa ra một cơ sở giải thích mới về việc tăng mật độ của mô hình vật liệu trong quá trình nén. Cấu trúc của luận văn Luận văn gồm 3 chƣơng.

Chƣơng 1: Trình bày tổng quan về hệ SiO2, các kết quả nghiên cứu, các công trình mô phỏng của vật liệu này về những đặc trƣng động lực học và vi câu trúc. Chƣơng 2: Trình bày nội dung phƣơng pháp mô phỏng động lực học phân tử và thế tƣơng tác sử dụng cho xây dựng các mô hình SiO2. Phƣơng pháp xác định thông số vật lý đặc trƣng của mô hình động lực học phân tử. Chƣơng 3: Tình bày kết quả của mô phỏng hệ SiO2 Các công trình khoa học đã đƣợc công bố: có 02 bài báo đã đƣợc công bố Báo cáo tại hội nghị Vật lí lí thuyết toàn quốc lần thứ 36 tại thành phố Quy Nhơn tháng 8/2011 với tên bài báo: “Molecular Dynamic Simulation of Amorphous SiO2 under Pressure” Huy N.K Và báo cáo tại hội nghị Việt Hàn tổ chức tại Đại học Bách Khoa Hà Nội tháng 11/2011 với tên bài báo: “Polymorphism and the relation between coordination distribution and bond angle in liquid silica and alumina” N.Hung 8 LUAN VAN CHAT LUONG download : add [e-mail] Nguyễn Thanh Hoa Cao học Vật lý 2009 CHƢƠNG I.

TỔNG QUAN Ngày nay với sự phát triển mạnh mẽ của ngành công nghệ thông tin đã đem lại nhiều lợi ích trong việc nghiên cứu cấu trúc cũng nhƣ đặc trƣng của vật liệu. Dựa trên phƣơng pháp mô phỏng vi mô các mô hình nguyên tử đƣợc mô tả từ đó cung cấp cho chúng ta nhiều thông tin, các tính chất vật lý về vật liệu cho các quá trình công nghệ.Với thời gian thực hiện ngắn, kinh phí rẻ hàng loạt các công trình mô phỏng đã và đang đƣợc thực hiện mở ra nhiều vấn đề cần phải đƣa ra thảo luận và nghiên cứu tiếp. Gần đây một vấn đề đang đƣợc nhiều ngƣời quan tâm là sự tồn tại của nhiều dạng cấu trúc bên trong trạng thái rắn và lỏng trong hệ ôxit. Trong quá trình tiến hành nghiên cứu mạng cấu trúc oxit các đặc trƣng liên quan nhƣ hàm phân bố xuyên tâm, các đơn vị cấu trúc cơ bản, thành phần hóa học, số phối trí, phân bố góc, simplex lần lƣợt đƣợc làm sáng tỏ.

Trong chƣơng này chúng tôi trình bày tổng quan về SiO2 các phƣơng pháp mô phỏng hệ cũng nhƣ các kết quả nghiên cứu vi mô về cấu trúc và các đặc trƣng vật lí SiO2 ở nhiệt độ T = 300K với các áp suất khác nhau 1. SiO2 Hiện nay SiO2 chủ yếu đƣợc ứng dụng trong các lĩnh vực công nghệ cao nhƣ điện tử, y học, quang học, siêu dẫn, hàng không, cơ khí và các thiết bị phân tích [27, 38]. Việc hiểu biết về cấu trúc của SiO2 ở các điều kiện đặc biệt (nhiệt độ cao, áp suất cao) đã trở nên rất quan trọng. Vì vậy các vấn đề liên quan đến cấu trúc địa phƣơng, chuyển pha cấu trúc, và đặc tính động lực học của vật liệu đều đƣợc quan tâm cả dƣới góc độ nghiên cứu ứug dụng và nghiên cứu cơ bản.

Silicđiôxit (SiO2) hay còn gọi là silica, là ôxit đƣợc cấu tạo từ một nguyên tử Si kết hợp với hai nguyên tử O, ở trạng thái tự nhiên Silica có thể tồn tại ở dạng tinh thể, không tan trong nƣớc, nóng chảy ở 17130 nhƣng khi đi sâu vào tìm hiểu và nghiên cứu ta nhận thấy đƣợc nhiều đặc tính quan trọng của nó. Silica (SiO2) có thể tồn tại dƣới nhiều dạng thù hình khác nhau. Ở nhiệt độ và áp suất 9 LUAN VAN CHAT LUONG download : add [e-mail] Nguyễn Thanh Hoa Cao học Vật lý 2009 thƣờng, cấu trúc và trạng thái cơ bản của silica là α-quartz (q-SiO2). Cấu trúc q- SiO2 sẽ chuyển sang cấu trúc coseite (c-SiO2) ở áp suất lớn hơn 2 GPa.

Trong coseite nguyên tử Si có số phối trí là 4, nguyên tử O có số phối trí là 2 với cấu trúc là sự xếp chặt của các tứ diện SiO4- các tứ diện này liên kết với nhau thông qua các đỉnh của tứ diện. Ở áp suất trên 8GPa cấu trúc SiO2 chuyển sang cấu trúc stishovite (s-SiO2) với nguyên tử Si có số phối trí là 6, nguyên tử O có số phối trí là 3, có cấu trúc đặc trƣng bởi các khối tám mặt. SiO6 liên kết với nhau thông qua cạnh chung. Tuy nhiên, silica có thể duy trì ở cấu trúc quartz ở trạng thái giả bền tới áp suất khoảng 15 GPa.

Áp suất trên 15 GPa cấu trúc quartz sẽ chuyển thành cấu trúc VĐH. Ở áp suất thƣờng, đơn vị cấu trúc cơ bản của silica ở pha lỏng và VĐH chủ yếu là các khối tứ diện SiO4 còn ở điều kiện áp suất cao thì đơn vị cấu trúc cơ bản chủ yếu là các khối bát diện SiO6. Sự chuyển pha từ cấu trúc mạng tứ diện (pha mật độ thấp) sang cấu trúc bát diện (pha mật độ cao) của SiO2 VĐH xảy ra trong khoảng rộng của mật độ từ 3,60g/cm3 đến 4,65 g/cm3. Quá trình này diễn ra với sự tăng phân bố phối trí trung bình của tất cả các cặp và sự giảm khoảng cách liên kết giữa cặp nguyên tử [5].

Phân tích cấu trúc vi mô của hệ ta thấy đƣợc sự chuyển pha thù hình trong ôxit có cấu trúc mạng ngẫu nhiên (cấu trúc framework) thƣờng xảy ra dƣới tác động của áp suất hoặc nhiệt độ. Quá trình chuyển pha lỏng-lỏng hoặc VĐH- VĐH thƣờng liên quan đến sự tăng số phối trí cặp, sự thay đổi khoảng cách tƣơng tác giữa các cặp nguyên tử, sự thay đổi phân bố góc liên kết giữa các nguyên tử trong cùng một đơn vị cấu trúc cơ bản và góc liên kết giữa các nguyên tử ở các đơn vị cấu trúc khác nhau. Nói cách khác trong quá trình chuyển pha thù hình của SiO2 cả cấu trúc trật tự gần (short-range order) và cấu trúc trật tự khoảng trung (immediate-range order) đều có liên quan.