Luận văn: Cấu trúc và cơ chế khuếch tán trong SiO2 lỏng bằng phương pháp mô phỏng

Nghiên cứu cấu trúc và cơ chế khuếch tán trong SiO2 lỏng bằng phương pháp mô phỏng. Phân tích chuyên sâu các đặc tính vật lý và hóa học của SiO2 lỏng.

• Tác giả: Phạm Xuân Trường
• Người hướng dẫn: TS. Phạm Hữu Kiên
• Trường học: Đại học Sư phạm Thái Nguyên
• Chuyên ngành: Vật lí chất rắn
• Đề tài: Nghiên cứu cấu trúc và cơ chế khuếch tán trong SiO2 lỏng 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: 2017
• Địa điểm: Thái Nguyên

Tóm tắt

I. Tổng quan về nghiên cứu cấu trúc và khuếch tán SiO2 lỏng

Silicon dioxide (SiO2) lỏng là vật liệu quan trọng trong nhiều ngành công nghiệp. Từ sản xuất thủy tinh đến luyện kim và địa chất học, hiểu biết về cấu trúc vi mô của SiO2 lỏng có ý nghĩa khoa học lớn. Ở trạng thái lỏng, SiO2 tồn tại dưới dạng mạng ngẫu nhiên ba chiều. Các nguyên tử silicon liên kết với nguyên tử oxy tạo thành các đơn vị cấu trúc tứ diện SiO4. Các tứ diện này nối với nhau qua nguyên tử cầu oxy ở đỉnh. Phương pháp mô phỏng động lực học phân tử là công cụ mạnh mẽ để nghiên cứu hệ vật liệu này. Phương pháp cho phép theo dõi chuyển động từng nguyên tử theo thời gian. Từ đó xác định được các đặc trưng vi cấu trúc như hàm phân bố xuyên tâm, số phối trí. Đồng thời tính toán hệ số khuếch tán và cơ chế vận chuyển vật chất trong hệ. Nghiên cứu SiO2 lỏng bằng mô phỏng giúp giải thích nhiều hiện tượng vật lý phức tạp. Kết quả mô phỏng cung cấp thông tin chi tiết mà thực nghiệm khó đạt được.

1.1. Đặc điểm cấu trúc vi mô của SiO2 lỏng

Cấu trúc vi mô của SiO2 lỏng được đặc trưng bởi mạng liên kết ngẫu nhiên. Nguyên tử silicon đóng vai trò trung tâm, liên kết với bốn nguyên tử oxy tạo thành tứ diện SiO4. Các tứ diện này nối với nhau qua nguyên tử oxy cầu. Liên kết Si-O có tính cộng hóa trị mạnh với độ dài khoảng 1,6 angstrom. Góc liên kết Si-O-Si trung bình vào khoảng 150 độ. Số phối trí trung bình của nguyên tử silicon thường là 4 ở điều kiện áp suất thường. Khi áp suất tăng, cấu trúc tứ diện có thể chuyển sang bát diện SiO6.

1.2. Tầm quan trọng của nghiên cứu SiO2 lỏng trong khoa học vật liệu

SiO2 lỏng là thành phần cơ bản của nhiều loại thủy tinh silicate. Nghiên cứu trạng thái lỏng giúp hiểu quá trình tạo thủy tinh và kết tinh. Trong luyện kim, xỉ silicate lỏng ảnh hưởng đến chất lượng kim loại. Địa chất học sử dụng dữ liệu magma silicate để mô hình hóa hoạt động núi lửa. Hiểu biết cơ chế khuếch tán trong SiO2 lỏng giúp dự đoán tính chất vận chuyển. Kiến thức này áp dụng trong thiết kế vật liệu chịu lửa và gốm sứ kỹ thuật. Nghiên cứu còn đóng góp vào lý thuyết về trạng thái lỏng của oxide.

II. Phân tích các vấn đề trong nghiên cứu SiO2 lỏng bằng mô phỏng

Nghiên cứu SiO2 lỏng gặp nhiều thách thức về mặt kỹ thuật và lý thuyết. Thực nghiệm trực tiếp quan sát cấu trúc vi mô ở nhiệt độ cao rất khó khăn. SiO2 lỏng có nhiệt độ nóng chảy cao khoảng 1986 K. Điều này gây khó khăn cho phép đo cấu trúc bằng nhiễu xạ tia X hay neutron. Hệ số tự khuếch tán trong SiO2 lỏng gần điểm nóng chảy rất nhỏ. Giá trị vào khoảng 10⁻⁶ cm²/s, đòi hỏi thời gian mô phỏng rất dài. Số bước mô phỏng phải đạt hàng trăm ngàn bước mới có kết quả tin cậy. Nếu hệ số khuếch tán nhỏ hơn 10⁻⁷ cm²/s thì gần như không thể xác định. Tính đa thù hình của SiO2 lỏng cũng là vấn đề phức tạp. Hệ có thể chuyển pha cấu trúc từ tứ diện sang bát diện dưới áp suất cao. Việc mô tả chính xác các đơn vị cấu trúc SiOx cần phương pháp tính toán phù hợp. Hiệu ứng nhiệt độ và áp suất tác động phức tạp lên vi cấu trúc hệ.

2.1. Thách thức trong xác định vi cấu trúc SiO2 lỏng

Xác định vi cấu trúc SiO2 lỏng đòi hỏi độ chính xác cao về vị trí nguyên tử. Hàm phân bố xuyên tâm g(r) là công cụ chính để phân tích cấu trúc. Tuy nhiên, việc tách riêng các cặp nguyên tử Si-Si, Si-O, O-O cần dữ liệu mô phỏng đủ lớn. Số phối trí trung bình có thể thay đổi theo nhiệt độ và áp suất. Phân bố góc liên kết Si-O-Si có dạng rộng, phản ánh tính ngẫu nhiên của mạng. Các đơn vị cấu trúc SiO4, SiO5, SiO6 cùng tồn tại tạo thành hệ phức tạp. Cần phương pháp thống kê mạnh để xử lý dữ liệu mô phỏng.

2.2. Vấn đề mô tả cơ chế khuếch tán trong hệ oxide lỏng

Cơ chế khuếch tán trong SiO2 lỏng khác biệt so với chất lỏng đơn giản. Chuyển động khuếch tán có thể xảy ra theo cơ chế nhảy đơn hoặc tập thể. Nguyên tử có thể đổi chỗ trực tiếp với lân cận hoặc đi đường vòng. Bình phương trung bình độ dịch chuyển gồm thành phần dao động và khuếch tán. Công thức ⟨r²⟩ = a(t) + 6Dt mô tả sự phụ thuộc vào thời gian. Hệ số D được tính từ độ dốc của đường tiệm cận. Với hệ nhớt cao như SiO2, thời gian đạt tiệm cận rất dài. Điều này tạo thách thức lớn cho tính toán mô phỏng.

III. Phương pháp mô phỏng động lực học phân tử nghiên cứu SiO2 lỏng

Phương pháp động lực học phân tử (MD) là công cụ hiệu quả để nghiên cứu SiO2 lỏng. MD dựa trên việc giải phương trình chuyển động Newton cho từng nguyên tử. Thế tương tác giữa các nguyên tử được mô tả bằng thế BKS (van Beest, Kramer, van Santen). Thế này bao gồm tương tác Coulomb và tương tác Lennard-Jones. Gần đúng Ewald-Hansen được sử dụng để xử lý tương tác điện tích dài tầm. Phương pháp này chia tương tác thành thành phần tầm ngắn và tầm dài. Thành phần tầm ngắn tính trực tiếp trong không gian thực. Thành phần tầm dài tính trong không gian Fourier. Hệ mô phỏng thường chứa vài trăm đến vài nghìn nguyên tử. Thời gian mô phỏng từ vài picosecond đến nanosecond. Bước thời gian điển hình là 1-2 femtosecond. Từ轨迹 nguyên tử, có thể tính hàm phân bố xuyên tâm g(r). Số phối trí được xác định từ tích phân của g(r). Phân bố góc liên kết được tính từ vị trí ba nguyên tử liên tiếp. Hệ số khuếch tán tính từ độ dốc bình phương dịch chuyển.

3.1. Thế tương tác BKS và gần đúng Ewald Hansen

Thế BKS là thế tương tác phổ biến cho hệ SiO2. Thế gồm thành phần Coulomb giữa các ion và thành phần Lennard-Jones. Các hệ số được hiệu chỉnh để tái tạo tính chất thực nghiệm của thạch anh. Gần đúng Ewald-Hansen xử lý hiệu quả tương tác điện tích dài tầm. Phương pháp chia tổng tương tác thành hai phần riêng biệt. Phần tầm ngắn tính trực tiếp khi khoảng cách nhỏ hơn bán kính cắt. Phần tầm dài chuyển đổi bằng biến đổi Fourier nhanh FFT. Cách tiếp cận này đảm bảo hội tụ nhanh và tiết kiệm thời gian tính toán.

3.2. Xác định đặc trưng vi cấu trúc từ mô phỏng

Hàm phân bố xuyên tâm g(r) là đại lượng cơ bản mô tả cấu trúc. g(r) cho biết xác suất tìm thấy nguyên tử ở khoảng cách r so với nguyên tử gốc. Số phối trí trung bình tính bằng cách tích phân g(r) đến cực tiểu đầu tiên. Phân bố góc liên kết Si-O-Si được xác định từ ba nguyên tử liên tiếp. Góc liên kết phản ánh mức độ mở của mạng oxide. Phương pháp Voronoi cũng được sử dụng để phân tích không gian. Từ轨迹 nguyên tử, tính toán bình phương dịch chuyển theo thời gian. Độ dốc đường tiệm cận cho hệ số khuếch tán D.

IV. Kết quả và ứng dụng của nghiên cứu SiO2 lỏng bằng mô phỏng

Kết quả mô phỏng cho thấy cấu trúc SiO2 lỏng thay đổi rõ rệt theo áp suất. Ở áp suất 0,1 GPa và nhiệt độ 3200 K, tỷ lệ đơn vị SiO4 chiếm 92%. Đơn vị SiO5 khoảng 8% và SiO6 gần như không tồn tại. Khi áp suất tăng lên 25,2 GPa, tỷ lệ SiO4 giảm xuống chỉ còn 5%. Tỷ lệ SiO5 đạt khoảng 36% và SiO6 chiếm 58%. Đặc biệt trong khoảng 12-15 GPa, tỷ lệ SiO5 đạt cực đại. Điều này chứng tỏ sự chuyển pha cấu trúc từ tứ diện sang bát diện. Cơ chế khuếch tán trong SiO2 lỏng là khuếch tán tập thể. Nguyên tử di chuyển bằng cách đổi chỗ gián tiếp với nguyên tử lân cận. Hệ số khuếch tán phụ thuộc mạnh vào nhiệt độ theo luật Arrhenius. Kết quả áp dụng trong nghiên cứu magma núi lửa và thiết kế vật liệu thủy tinh. Hiểu biết về chuyển pha áp suất cao có ý nghĩa trong địa chất học hành tinh.

4.1. Chuyển pha cấu trúc từ tứ diện sang bát diện dưới áp suất

Khi áp suất tăng từ 0 đến 25 GPa, SiO2 lỏng trải qua chuyển pha cấu trúc. Đơn vị cấu trúc SiO4 chiếm ưu thế ở áp suất thấp. Khi nén áp suất, nguyên tử oxy bị đẩy vào gần silicon hơn. Điều này tạo điều kiện hình thành liên kết Si-O bổ sung. Đơn vị SiO5 xuất hiện như trung gian chuyển pha. Ở áp suất 12-15 GPa, tỷ lệ SiO5 đạt cực đại khoảng 40%. Trên 15 GPa, SiO5 chuyển dần thành SiO6. Ở 25 GPa, cấu trúc bát diện SiO6 chiếm ưu thế với 58% tỷ lệ. Chuyển pha này liên quan đến hiện tượng lỏng-lỏng trong oxide.

4.2. Ứng dụng kết quả nghiên cứu trong khoa học và công nghệ

Kết quả nghiên cứu SiO2 lỏng có nhiều ứng dụng quan trọng. Trong địa chất học, dữ liệu giúp mô hình hóa quá trình magma silicate. Hiểu biết cấu trúc áp suất cao áp dụng cho nghiên cứu lớp vỏ hành tinh. Trong công nghiệp thủy tinh, kết quả tối ưu hóa quy trình tạo phôi. Kiến thức về cơ chế khuếch tán giúp kiểm soát tốc độ kết tinh. Vật liệu chịu lửa silicate được thiết kế dựa trên tính chất nhiệt độ cao. Kết quả còn đóng góp vào lý thuyết trạng thái lỏng oxide cơ bản. Nghiên cứu mở đường cho phát triển vật liệu nano silicate mới.

Nội dung chính

Tổng quan nghiên cứu

Hợp chất điôxít silic chiếm hơn 60% khối lượng vỏ Trái Đất và giữ vai trò cốt lõi trong địa chất học, vật lý chất rắn cũng như công nghệ vật liệu bán dẫn và silicat. Tuy nhiên, việc khảo sát trực tiếp cấu trúc vi mô và động học khuếch tán của vật liệu này ở trạng thái lỏng dưới điều kiện khắc nghiệt gặp rất nhiều rào cản kỹ thuật. Tại nhiệt độ cao khoảng 3200 K và dải áp suất từ 0 đến 25 GPa tương đương điều kiện sâu trong lòng Trái Đất, các phép đo thực nghiệm nhiễu xạ tia X hay tán xạ nơtron rất khó thực hiện với độ phân giải nguyên tử đầy đủ.

Nghiên cứu tập trung giải quyết bài toán mô hình hóa vi cấu trúc, làm sáng tỏ tính đa thù hình và cơ chế khuếch tán phi đồng nhất của hệ ôxít lỏng. Mục tiêu cụ thể là xây dựng mô hình mô phỏng gồm 1998 nguyên tử (bao gồm 666 nguyên tử Silic và 1332 nguyên tử Ôxi) ở nhiệt độ 3200 K trong dải áp suất nén từ 0.10 GPa đến 25.20 GPa, đồng thời mở rộng khảo sát cơ chế khuếch tán động học tại các mốc nhiệt độ 2600 K, 3000 K, 3200 K và 3500 K.

Ý nghĩa học thuật của đề tài thể hiện qua việc cung cấp các số liệu định lượng chuẩn xác về sự chuyển đổi mạng liên kết. Số phối trí trung bình của cặp Si-O tăng từ 4.07 lên 5.50 khi nén áp suất, đi kèm sự thay đổi mật độ mẫu từ 2.634 g/cm3 lên 3.991 g/cm3. Những kết quả này đóng góp cơ sở dữ liệu quan trọng cho việc dự báo các đặc tính nhiệt động, dẫn truyền địa chấn và hỗ trợ định hướng tổng hợp vật liệu vô định hình chịu áp suất cao.

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 lý thuyết mạng vô định hình liên tục và mô hình đa diện phối trí, trong đó cấu trúc của mạng chất lỏng được hình thành từ các khối đa diện cơ bản có tâm là Silic và đỉnh là Ôxi. Sự biến đổi trật tự gần và trật tự trung bình được giải thích thông qua sự sắp xếp của các đơn vị tứ diện, tam giác kép và bát diện.

Lý thuyết động học không đồng nhất không gian do Adam và Gibbs khởi xướng được ứng dụng để mô tả động lực học chất lỏng siêu lạnh và chất lỏng nhiệt độ cao. Theo lý thuyết này, quá trình dịch chuyển cấu trúc không diễn ra đồng nhất trên toàn hệ mà phân bố theo các vùng độc lập, tạo thành các đám nguyên tử có độ linh động khác nhau.

Định luật khuếch tán Fick thứ nhất và thứ hai kết hợp với quy luật nhiệt động Arrhenius được sử dụng làm cơ sở xác định hệ số tự khuếch tán thông qua độ dịch chuyển bình phương trung bình của các nguyên tử theo thời gian. Hàm thế tương tác cặp van Beest, Kramer and van Santen được áp dụng với các thông số điện tích hiệu dụng chính xác, trong đó điện tích của Ôxi là -1.2e và Silic là +2.4e, đảm bảo tính ổn định của liên kết cộng hóa trị lai ion trong mạng cấu trúc.

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

Hệ mô phỏng sử dụng cỡ mẫu chuẩn hóa gồm 1998 nguyên tử (666 nguyên tử Silic và 1332 nguyên tử Ôxi) được đặt trong ô cơ sở lập phương với điều kiện biên tuần hoàn ba chiều nhằm triệt tiêu hoàn toàn các sai số do hiệu ứng bề mặt.

Quy trình chọn mẫu và tạo cấu hình ban đầu được thực hiện bằng cách gieo ngẫu nhiên các nguyên tử, sau đó nung nóng hệ lên 5000 K và giữ trong hơn 50.000 bước thời gian tích phân để xóa bỏ trạng thái nhớ ban đầu. Mẫu tiếp tục được làm nguội đẳng nhiệt qua các nấc 4500 K và 4000 K (mỗi nấc hồi phục trong 200.000 bước) trước khi đưa về trạng thái cân bằng tại 3200 K với thời gian hồi phục 1.000.000 bước tích phân. Bước thời gian mô phỏng được cố định ở mức 0.47 femtosecond.

Phương pháp động lực học phân tử cổ điển tích hợp thuật toán Verlet được lựa chọn nhờ khả năng theo dõi chuyển động chi tiết của từng hạt theo thời gian thực. Đối với tương tác tĩnh điện Coulomb tầm xa, nghiên cứu sử dụng kỹ thuật gần đúng Ewald-Hansen kết hợp ngắt tương tác tầm ngắn để tối ưu hóa thời gian tính toán mà vẫn bảo toàn độ chính xác tuyệt đối. Quá trình tạo mẫu ở các mức áp suất khác nhau sử dụng mô hình đẳng nhiệt - đẳng áp kết hợp kỹ thuật điều chỉnh tọa độ với bước tăng áp suất rất nhỏ, sau đó chuyển sang mô hình bảo toàn năng lượng để trích xuất hàm phân bố xuyên tâm và phân bố góc liên kết.

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 trên dải áp suất từ 0.10 GPa đến 25.20 GPa đã làm sáng tỏ sự tiến hóa vi cấu trúc rõ rệt của hệ vật liệu:

1. Sự chuyển dịch tỷ lệ các đơn vị cấu trúc: Ở áp suất thường 0.10 GPa, đơn vị tứ diện chiếm tỷ lệ áp đảo lên tới 92%, đơn vị 5 phối trí chiếm khoảng 8% và đơn vị bát diện 6 phối trí hầu như không xuất hiện. Khi áp suất nén tăng lên 25.20 GPa, tỷ lệ tứ diện giảm mạnh 87% xuống chỉ còn 5%, trong khi đơn vị bát diện tăng vọt lên 58% và đơn vị 5 phối trí chiếm 36%.

2. Sự xuất hiện cực đại của cấu trúc trung gian: Trong khoảng áp suất chuyển tiếp từ 12 GPa đến 15 GPa, đơn vị cấu trúc 5 phối trí đạt nồng độ cực đại trước khi suy giảm dần dưới áp lực nén cao hơn để nhường chỗ cho cấu trúc bát diện 6 phối trí.

3. Biến thiên số phối trí và khoảng cách liên kết: Số phối trí cặp trung bình của Si-O tăng 35.1% từ giá trị 4.07 lên 5.50. Độ dài liên kết Si-O tăng nhẹ từ 1.60 Angstrom lên 1.64 Angstrom. Đồng thời, khoảng cách liên kết O-O co ngắn từ 2.60 Angstrom xuống 2.44 Angstrom, khiến số phối trí trung bình O-O tăng vọt 86.0% từ 8.17 lên 15.20.

4. Cơ chế khuếch tán theo đám phản ứng: Các nguyên tử phân tách thành hai miền rõ rệt gồm đám không phản ứng và đám phản ứng tuần hoàn. Sự khuếch tán diễn ra chủ yếu thông qua các phản ứng bẻ gãy và tái tạo liên kết có hiệu quả giữa các đơn vị 4 phối trí và 5 phối trí lân cận.

Thảo luận kết quả

Khi áp suất nén tăng từ 0.10 GPa lên 25.20 GPa, mật độ của hệ tăng từ 2.634 g/cm3 lên 3.991 g/cm3. Việc thu hẹp thể tích tự do buộc các nguyên tử Ôxi phải sắp xếp lại gần nhau hơn, làm biến dạng các góc liên kết O-Si-O và Si-O-Si. Quá trình chuyển pha đa thù hình từ pha tứ diện mật độ thấp sang pha bát diện mật độ cao không diễn ra đột ngột mà là một quá trình liên tục thông qua trạng thái trung gian là đơn vị 5 phối trí.

Dữ liệu mô phỏng có thể được minh họa trực quan thông qua các đường cong hàm phân bố xuyên tâm thành phần với các đỉnh thứ nhất sắc nét và bảng đối sánh số phối trí theo từng nấc áp suất. Sự dịch chuyển của đỉnh hàm phân bố cặp Si-O từ 1.60 Angstrom sang 1.64 Angstrom phản ánh sự gia tăng độ dài liên kết do lực đẩy tĩnh điện giữa các nguyên tử Ôxi tăng lên khi số phối trí tăng từ 4 lên 6.

So sánh với các nghiên cứu thực nghiệm nhiễu xạ tia X công bố trước đây ở áp suất 0 GPa, kết quả tính toán độ dài liên kết Si-O đạt 1.60 Angstrom, O-O đạt 2.60 Angstrom và Si-Si đạt 3.10 Angstrom hoàn toàn trùng khớp với sai số nhỏ hơn 0.01 Angstrom. Điều này khẳng định thế tương tác BKS và phương pháp tính toán gần đúng Ewald-Hansen có độ tin cậy vượt trội trong việc mô tả vật liệu silicat lỏng.

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

1. Mở rộng quy mô tính toán hệ nguyên tử: Tiến hành nâng cấp kích thước mô hình từ 1998 nguyên tử lên hơn 50.000 nguyên tử nhằm khảo sát chi tiết hơn các hiệu ứng ranh giới và hiện tượng không đồng nhất động học quy mô lớn. Giải pháp này cần được các nhóm nghiên cứu vật lý tính toán triển khai trong thời gian từ 6 đến 12 tháng tới.

2. Hoàn thiện hàm thế tương tác đa thành phần: Nghiên cứu tích hợp các hiệu ứng phân cực nguyên tử vào thế tương tác nhằm giảm sai số năng lượng kích hoạt khuếch tán xuống dưới mức 5%, thực hiện định kỳ hàng quý bởi các chuyên gia mô hình hóa vật liệu.

3. Chuẩn hóa cơ sở dữ liệu vi cấu trúc silicat: Xây dựng hệ thống bảng tra cứu số phối trí, phân bố góc liên kết và hệ số khuếch tán trên dải nhiệt độ từ 2600 K đến 3500 K trong vòng 18 tháng, do các viện nghiên cứu vật liệu phối hợp với các trường đại học thực hiện.

4. Ứng dụng số liệu vào mô hình địa vật lý: Ứng dụng trực tiếp các hàm mật độ phụ thuộc áp suất từ 2.634 đến 3.991 g/cm3 vào phần mềm mô phỏng truyền sóng địa chấn tại đới chuyển tiếp manti Trái Đất, hướng tới mục tiêu nâng cao độ chính xác dự báo địa chấn trên 90% trong giai đoạn 2 năm.

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

1. Các nhà nghiên cứu và chuyên gia vật lý chất rắn: Nắm vững quy trình thiết lập mô phỏng động lực học phân tử với thế BKS, kỹ thuật tính gần đúng tương tác Coulomb và các phương pháp giải tích hàm phân bố xuyên tâm để ứng dụng cho các vật liệu oxit vô định hình khác.

2. Các nhà địa vật lý và khoáng vật học: Sử dụng dữ liệu định lượng về vi cấu trúc và tính đa thù hình của silicat nóng chảy tại nhiệt độ 3200 K để giải mã trạng thái macma và động lực học sâu trong lòng Trái Đất.

3. Giảng viên, nghiên cứu sinh và học viên cao học: Tiếp cận tài liệu tham khảo chuẩn mực về phương pháp phân tích số phối trí riêng biệt cho từng loại đa diện, phương pháp kiểm soát nhiệt độ và áp suất trong các mô hình tính toán hiện đại.

4. Kỹ sư công nghệ sản xuất thủy tinh và vật liệu chịu lửa: Vận dụng cơ chế khuếch tán theo đám và phản ứng bẻ gãy liên kết để tối ưu hóa quy trình nấu luyện, làm nguội và kiểm soát độ nhớt của silicat lỏng trong công nghiệp.

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

1. Tại sao nghiên cứu lại chọn mốc nhiệt độ 3200 K để khảo sát cấu trúc mạng lỏng? Nhiệt độ 3200 K tương ứng với điều kiện nhiệt động thực tế tại đới manti trong lòng Trái Đất, nơi hợp chất silicat tồn tại dưới dạng lỏng nhớt. Việc khảo sát ở 3200 K giúp tái hiện chính xác các chuyển biến cấu trúc từ mật độ 2.634 g/cm3 lên 3.991 g/cm3, cung cấp bằng chứng vật lý quan trọng cho các mô hình địa chất học.

2. Quá trình chuyển pha đa thù hình dưới áp suất cao diễn ra theo cơ chế nào? Dưới áp suất nén từ 0.10 GPa đến 25.20 GPa, hệ chuyển đổi liên tục từ pha tứ diện sang pha bát diện thông qua cấu trúc trung gian 5 phối trí. Tỷ lệ tứ diện giảm từ 92% xuống 5%, trong khi bát diện tăng từ 0% lên 58%, và cấu trúc trung gian đạt cực đại ở dải 12 đến 15 GPa.

3. Thế tương tác BKS mang lại ưu điểm gì vượt trội trong bài toán mô phỏng này? Thế BKS mô tả chính xác tương tác cặp Coulomb tầm xa và thế đẩy Born-Mayer tầm ngắn với điện tích hiệu dụng Silic (+2.4e) và Ôxi (-1.2e). Mô hình này tái tạo chính xác độ dài liên kết Si-O đạt 1.60 Angstrom, hoàn toàn phù hợp với dữ liệu thực nghiệm nhiễu xạ tia X.

4. Cơ chế khuếch tán trong mạng chất lỏng nhớt tuân theo quy luật gì? Quá trình khuếch tán không diễn ra đơn lẻ mà vận hành theo cơ chế chuyển động đám không đồng nhất. Sự dịch chuyển của các hạt được kích hoạt chủ yếu thông qua các phản ứng trao đổi liên kết có hiệu quả giữa các cấu trúc 4 phối trí và 5 phối trí trên dải nhiệt độ 2600 K đến 3500 K.

5. Cỡ mẫu 1998 nguyên tử có đảm bảo độ tin cậy và tính đại diện cho vật liệu thực? Hệ 1998 hạt gồm 666 Silic và 1332 Ôxi kết hợp điều kiện biên tuần hoàn đã triệt tiêu hoàn toàn hiệu ứng biên bề mặt. Các kiểm chứng đối sánh cho thấy các đặc trưng hàm phân bố xuyên tâm và số phối trí của hệ hoàn toàn tương đồng với các mô hình lớn hàng trăm nghìn hạt.

Kết luận

• Luận văn đã xây dựng thành công mô hình động lực học phân tử cho hệ 1998 nguyên tử ở nhiệt độ 3200 K trong dải áp suất nén từ 0.10 GPa đến 25.20 GPa.
• Xác định định lượng sự gia tăng số phối trí cặp trung bình Si-O từ 4.07 lên 5.50 và mật độ mẫu từ 2.634 g/cm3 lên 3.991 g/cm3.
• Làm sáng tỏ tính đa thù hình với quá trình chuyển pha liên tục từ tứ diện (92%) sang bát diện (58%) thông qua đơn vị trung gian 5 phối trí đạt đỉnh tại 12 đến 15 GPa.
• Khám phá cơ chế khuếch tán phi đồng nhất theo đám và khẳng định vai trò chủ đạo của các phản ứng trao đổi liên kết có hiệu quả trên dải nhiệt độ 2600 K đến 3500 K.
• Cung cấp nguồn dữ liệu vi cấu trúc chuẩn xác, mở ra hướng ứng dụng sâu rộng trong vật lý chất rắn, địa vật lý và công nghệ vật liệu silicat chịu nhiệt.

Các nhóm nghiên cứu có thể tiếp tục kế thừa mô hình này để mở rộng khảo sát cho các hệ vật liệu oxit phức hợp ba thành phần trong giai đoạn tiếp theo.

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

Chương 1 TỔNG QUAN Trong chương này, chúng tôi trình bày: tổng quan một số kết quả nghiên cứu cấu trúc và các tính chất về vật liệu SiO2 ở trạng thái lỏng bằng cả thực nghiệm, lý thuyết lẫn mô phỏng máy tính. Một số cơ chế khuếch tán trong vật liệu cũng được đề cập trong chương này. Cấu trúc và các tính chất về vật liệu SiO2 ở trạng thái lỏng Trong nhưng năm gần đây, nhiều nghiên cứu thực nghiệm và lý thuyết trong lĩnh vực khoa học vật liệu đã chứng minh rằng cách xắp xếp nguyên tử và quá trình hình thành cấu trúc có ảnh hưởng đến các tính chất vật lý của SiO2. Tuy nhiên, nghiên cứu sự chuyển pha cấu trúc của SiO2 rất khó đo được bằng thực nghiệm.

Sự hiểu biết biết về hiện tượng này vẫn chưa được thỏa đáng và vẫn đang còn nhiều vấn đề cần thảo luận, đặc biệt là sự thay đổi cấu trúc ở nhiệt độ 3200K trong một dải áp suất vẫn đang là một đề tài nóng. Trong lĩnh vực khoa học trái đất, SiO2 lỏng trong lòng trái đất lại tồn tại ở nhịêt độ 3200K nên được rất nhiều nhà khoa học tập trung nghiên cứu bằng nhiều kĩ thuật thực nghiệm và tính toán lý thuyết như nhiễu xạ tia X, nhiễu xạ nơtron, phổ Raman, phổ hấp thụ tia X, cộng hưởng từ hạt nhân. Còn đối với phương pháp mô phỏng, việc nghiên cứu sẽ đem lại nhiều kết quả có giá trị giúp so sánh dự báo các tính chất mới lạ của vật liệu. Trong công trình [23] trạng thái của SiO2 lỏng bao gồm các đơn vị cấu trúc cơ bản liên kết nhau trong một mạng liên tục trong không gian ba chiều hữu hạn và không có trật tự xa.

Mỗi đơn vị cấu trúc cơ bản là một khối tứ diện SiO4 với 4 nguyên tử ôxi(O) nằm ở đỉnh và tâm là nguyên tử silic(Si). Liên kết giữa hai tứ diện đòi hỏi một liên kết góc Si-O-Si và hai góc nhị diện. Sự biến đổi hai góc này được xem như nguyên nhân chính làm thay đổi cấu trúc của SiO2 lỏng. Trong công trình [20] bằng kĩ thuật nhiễu xạ tia X đã cho kết quả: trong mỗi đơn vị cấu trúc SiO4 các 4 thông tin cấu trúc được xác định bởi: số phối trí trung bình (SPTTB); độ dài liên kết Si-Si, O-O, Si-O; góc liên kết Si-O-Si và góc liên kết O-Si-O.

Trong thực nghiệm từ đường cong tán xạ tia X hoặc tán xạ nơtron ta có thể xác định được thông số quan trọng xác định cấu trúc vật liệu. Thừa số cấu trúc cho phép xác định số lượng trung bình các nguyên tử ở khoảng cách bất kì tính từ nguyên tử đang xét. Khi phân tích Phu-ri-ê thừa số cấu trúc ta còn thu được HPBXT, một thông số dùng để xác định trật tự gần của các vật liệu có cấu trúc mất trật tự. Thực nghiệm đã chứng tỏ rằng cấu trúc của SiO2 lỏng phụ thuộc khá mạnh vào áp suất và ít thay đổi theo nhiệt độ.

Phân tích thừa số cấu trúc nhiễu xạ tia X cho thấy khi tăng áp suất có sự thay đổi vị trí đỉnh và cường độ đỉnh nhiễu xạ thứ nhất. Khi áp suất tăng từ 0,1MPa đến 8GPa thì vị trí đỉnh nhiễu xạ thứ nhất dịch từ Q~ 1,55 đến 1,92 Å-1 trong khi cường độ của nó hấu như không thay đổi. Sự thay đổi của thừa số cấu trúc xảy ra mạnh nhất trong vùng áp suất từ 8 GPa đến 28 GPa, cường độ đỉnh nhiễu xạ thứ nhất giảm đi 50% trong khi vị trí đỉnh nhiễu xạ thứ nhất dịch từ Q~ 1,92 đến 2,29 Å-1 và xuất hiện thêm đỉnh nhiễu xạ mới ở vị trí 3,18 Å-1. Nếu tiếp tục tăng áp suất đến 42 GPa thì thừa số cấu trúc gần như không thay đổi nữa.