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.
