Tổng quan nghiên cứu

Silica và silicat là nhóm vật liệu chiếm hơn 12% khối lượng lớp vỏ Trái Đất, giữ vai trò nền tảng trong các ngành công nghiệp mũi nhọn như vi điện tử, quang học, gốm sứ và sản xuất thủy tinh kỹ thuật cao. Dưới tác động của các điều kiện khắc nghiệt về nhiệt độ và áp suất, các vật liệu alkali silicat thể hiện nhiều đặc tính vật lý phức tạp như chuyển pha đa thù hình, khuếch tán dị thường và bất đồng nhất động học. Tuy nhiên, các kỹ thuật thực nghiệm truyền thống gặp nhiều rào cản khi theo dõi chi tiết quá trình tái cấu trúc ở cấp độ nguyên tử.

Luận văn thạc sĩ chuyên ngành Vật lí lí thuyết và vật lí toán (mã số: 8440130.01) của học viên Nguyễn Hoàng Mạnh Vượng, dưới sự hướng dẫn khoa học của Tiến sĩ Nguyễn Thị Thanh Hà và Giáo sư Nguyễn Quang Báu tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, đã tập trung giải quyết bài toán trên. Mục tiêu chính của đề tài là khảo sát tường minh vi cấu trúc và các đặc trưng động học của chất lỏng natri silicat dưới các điều kiện nén ép khác nhau.

Phạm vi nghiên cứu được thực hiện trên hai mô hình tính toán chuẩn mực: hệ natri silicat gồm 3000 nguyên tử ở nhiệt độ 3500 K với áp suất biến thiên từ 0 đến 60 GPa, và hệ 7995 nguyên tử ở nhiệt độ 1873 K tại áp suất môi trường 0,1 GPa. Kết quả nghiên cứu cung cấp bộ số liệu toàn diện về sự biến đổi trật tự tầm gần, trật tự tầm trung và nguồn gốc vi mô của hiện tượng phân tách kênh khuếch tán ion natri. Đây là cơ sở khoa học quan trọng giúp nâng cao độ chính xác trong việc kiểm soát các chỉ số cơ lý và định hướng công nghệ chế tạo vật liệu chịu áp lực 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 dựa trên nền tảng lý thuyết mạng ngẫu nhiên liên tục của Zachariasen kết hợp với mô hình biến tính mạng thủy tinh của Warren. Trong hệ vật liệu natri silicat, các ion silic (Si) đóng vai trò là chất tạo mạng, trong khi các cation natri (Na) hoạt động như chất điều chỉnh mạng, làm đứt gãy các liên kết cầu nối Si-O-Si để tạo ra các tâm oxy không cầu.

Khung lý thuyết của luận văn tập trung vào 4 khái niệm và mô hình cấu trúc then chốt:

  • Oxy cầu (Bridging Oxygen - BO) và Oxy không cầu (Non-bridging Oxygen - NBO): Đại diện cho mức độ polyme hóa và tính liên kết không gian của mạng silica.
  • Đơn vị phối trí SiOx (với x biến thiên từ 4, 5 đến 6) và OSiy (với y biến thiên từ 2 đến 4): Phản ánh trật tự tầm gần và trật tự tầm trung của vật liệu.
  • Phương pháp hình học Simplex: Định nghĩa các quả cầu rỗng đi qua 4 nguyên tử bất kỳ trong không gian để xác định kích thước lỗ xốp và vị trí không gian cư trú của các hạt kiềm.
  • Hàm phân bố đám Flk(r,t) và lý thuyết Mode-Coupling: Dùng để định lượng sự kết cụm không gian và phân tích sự chuyển tiếp động học giữa trạng thái chất lỏng bền vững và chất lỏng dễ vỡ.

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ử (Molecular Dynamics - MD) cổ điển, tích phân phương trình chuyển động Newton thông qua thuật toán Verlet với bước thời gian tích phân siêu nhỏ là 1,0 femtosecond. Phương pháp này cho phép theo dõi chính xác quỹ đạo không gian và vận tốc của từng hạt theo thời gian thực.

Cỡ mẫu và quy trình chọn mẫu được phân tách thành hai mô hình chuyên biệt:

  • Mô hình nén áp suất cao (3000 nguyên tử): Bao gồm 900 nguyên tử Si, 1900 nguyên tử O và 200 nguyên tử Na. Tương tác giữa các hạt được mô tả thông qua hàm thế Buckingham với bộ tham số Teter chuẩn hóa. Hệ được nung chảy ngẫu nhiên ở 5000 K trong 50.000 bước mô phỏng để xóa bỏ trạng thái ghi nhớ ban đầu, làm nguội nhanh về 3500 K trong 100.000 bước ở trạng thái đẳng nhiệt - đẳng áp (NPT), sau đó tiến hành nén đẳng nhiệt từ 0 đến 60 GPa trong 1.000.000 bước ở trạng thái thể tích - nhiệt độ không đổi (NVT).
  • Mô hình động học vi mô (7995 nguyên tử): Gồm 2132 nguyên tử Si, 4797 nguyên tử O và 1066 nguyên tử Na ở nhiệt độ 1873 K và áp suất 0,1 GPa. Mô hình sử dụng hàm thế cặp kết hợp thế ba phần Kawamura để mô tả chính xác sự phân cực góc O-Si-O và lực liên kết cộng hóa trị.

Để triệt tiêu các sai số ngẫu nhiên và đảm bảo độ hội tụ thống kê, toàn bộ dữ liệu về số phối trí, hàm phân bố xuyên tâm và phân bố góc liên kết đều được tính toán lấy trung bình từ 1000 cấu hình cân bằng cuối cùng, với chu kỳ trích xuất 50 bước mô phỏng một lần.

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

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

Quá trình phân tích dữ liệu mô phỏng đã đưa ra 4 phát hiện quan trọng về cấu trúc và động học của hệ natri silicat:

  • Chuyển pha cấu trúc liên tục theo áp suất: Tại áp suất 0 GPa, mạng lưới chủ yếu cấu thành từ các tứ diện SiO4 chiếm tỷ trọng 89%, trong khi SiO5 chỉ chiếm 4%. Khi tăng áp suất lên 30 GPa, tỷ phần SiO4 giảm mạnh xuống 49%, nhường chỗ cho sự gia tăng của SiO5 lên mức 41% và bắt đầu xuất hiện SiO6 với 9%. Tại 40 GPa, đơn vị phối trí trung gian SiO5 đạt đỉnh cực đại 47%. Đến áp suất cực hạn 60 GPa, cấu trúc bát diện SiO6 vươn lên chiếm ưu thế tuyệt đối với 47%, SiO5 duy trì ở mức 38%, còn SiO4 suy giảm xuống chỉ còn 8%.
  • Tái tổ chức trật tự tầm trung qua cầu nối oxy: Ở điều kiện thường 0 GPa, đơn vị phối trí OSi2 chiếm ưu thế tuyệt đối. Dưới lực nén 60 GPa, tỷ phần OSi2 suy giảm mạnh về mức 43%, trong khi đơn vị OSi3 tăng vọt từ dưới 10% lên 49,36%, trở thành thành phần liên kết chủ đạo trong chất lỏng nén ép.
  • Phân bố không đồng nhất của cation kiềm: Kỹ thuật phân tích Simplex chứng minh có tới 70% tổng số nguyên tử natri tập trung bên trong các khối simplex O4 (đi qua 4 nguyên tử oxy), trong khi khối SO3 chỉ chứa khoảng 30%. Đặc biệt, mật độ natri trong các cụm oxy không cầu (CLNBO) đạt giá trị trung bình 4,88 ion/cụm, cao hơn gấp gần 3 lần so với cụm Si (khoảng 1,63 ion/cụm) và cụm BO (khoảng 2,45 ion/cụm).
  • Hiện tượng bất đồng nhất động học rõ nét: Nhóm 10% nguyên tử di động nhanh nhất (tương đương 700 nguyên tử trong mô hình 7995 hạt) có độ dịch chuyển bình phương trung bình tăng trưởng tuyến tính vượt trội theo thời gian. Ngược lại, nhóm 10% nguyên tử bất động chỉ dao động quanh vị trí cân bằng với giá trị hàm phân bố đám Flk(r,t) tại khoảng cách 1,95 Å suy giảm đột ngột từ 700 xuống 382.

Thảo luận kết quả

Các kết quả trên có thể được minh họa trực quan thông qua biểu đồ đường phân bố hàm xuyên tâm g(r) và đồ thị tỷ phần cấu trúc SiOx theo dải áp suất từ 0 đến 60 GPa. Đồ thị cho thấy sự giao thoa rõ nét giữa ba đường cong SiO4, SiO5 và SiO6 trong khoảng áp suất trung gian từ 25 đến 40 GPa, khẳng định cơ chế chuyển tiếp cấu trúc từ tứ diện sang bát diện bắt buộc phải đi qua trạng thái trung gian năm phối trí.

Nguyên nhân của quá trình biến đổi cấu trúc bắt nguồn từ sự nén ép thể tích tự do, buộc các nguyên tử Si và O phải gia tăng số phối trí cục bộ để đạt cấu hình đóng gói tối ưu. Số phối trí trung bình của Na-O tăng đều đặn từ 5,61 ở 0 GPa lên hơn 10,2 ở 60 GPa. Đỉnh phân bố xuyên tâm của cặp Na-Na ghi nhận mức co rút khoảng cách lên tới 46,7%, phản ánh sự thu hẹp không gian tự do giữa các hạt kiềm.

Hiện tượng các ion natri tập trung cao độ tại các vùng NBO được giải thích bởi lực tương tác tĩnh điện Coulomb mạnh mẽ giữa cation Na+ và các anion NBO mang điện tích âm cục bộ. Sự tách biệt giữa các vi vùng giàu natri (gắn với NBO) và vi vùng nghèo natri (gắn với Si và BO) đã hình thành nên các kênh dẫn truyền ion ưu tiên. Kết quả này hoàn toàn tương thích với các dữ liệu thực nghiệm tán xạ neutron và cộng hưởng từ hạt nhân đa chiều, giải thích triệt để nguồn gốc của hiệu ứng hỗn hợp kiềm và tính linh động dị thường của ion natri trong chất lỏng silicat.

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

Dựa trên các phát hiện chuyên sâu về vi cấu trúc và động học, luận văn đề xuất 4 khuyến nghị thực tiễn nhằm tối ưu hóa quy trình sản xuất và nghiên cứu phát triển vật liệu mới:

  • Tối ưu hóa chế độ nhiệt áp trong sản xuất thủy tinh kỹ thuật: Các cơ sở sản xuất cần duy trì áp lực tạo hình khuôn dập dưới 5 GPa và nhiệt độ làm việc trong dải 1400 - 1873 K để kiểm soát tỷ trọng đơn vị tứ diện SiO4 luôn đạt trên 85%. Việc này giúp hạn chế tối đa sự hình thành sớm của các pha trung gian SiO5 kém bền, từ đó nâng cao 15 - 20% độ bền uốn và độ trong suốt quang học của sản phẩm.
  • Kiểm soát chính xác hàm lượng kiềm Na2O phối liệu: Doanh nghiệp sản xuất vật liệu cách điện và cáp quang cần điều chỉnh tỷ lệ natri oxit ở ngưỡng 15 - 20% mol. Tỷ lệ này giúp hạn chế sự kết nối thành chuỗi của các kênh oxy không cầu NBO, ngăn ngừa hiện tượng suy giảm hệ số đàn hồi cơ học C11 và giảm thiểu dòng rò ion không mong muốn.
  • Ứng dụng thiết kế vật liệu màng điện giải rắn cho pin ion kiềm: Các viện nghiên cứu công nghệ năng lượng nên tận dụng cơ chế tạo kênh dẫn truyền ion trong cấu trúc giàu NBO để phát triển màng điện ly rắn thế hệ mới. Mục tiêu hướng tới việc nâng cao độ dẫn ion natri lên trên 25% ở nhiệt độ phòng trong lộ trình nghiên cứu từ năm 2026 đến 2028.
  • Nâng cấp mô hình mô phỏng tính toán vật liệu đa tỉ lệ: Các nhóm nghiên cứu lý thuyết cần tích hợp thuật toán máy học với động lực học phân tử để mở rộng quy mô mô phỏng lên trên 50.000 nguyên tử, giúp khảo sát sâu hơn động học hồi phục cấu trúc tại điểm chuyển pha thủy tinh trong dải nhiệt độ thấp dưới 1000 K.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chính sau:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý lý thuyết, Vật lý chất rắn: Cung cấp phương pháp luận chặt chẽ trong việc xây dựng thế tương tác, giải thuật toán tích phân chuyển động và khai thác công cụ hình học Simplex để phân tích cấu trúc vô định hình.
  • Kỹ sư công nghệ vật liệu silicat, thủy tinh và gốm sứ cao cấp: Cung cấp các số liệu định lượng chính xác về sự phụ thuộc của số phối trí và góc liên kết vào áp suất nén, giúp tối ưu hóa công thức phối liệu và thông số gia công nhiệt áp.
  • Chuyên gia nghiên cứu và phát triển pin thế hệ mới: Giúp nắm bắt tường minh cơ chế khuếch tán theo kênh của ion natri, phục vụ việc thiết kế và chế tạo các chất điện phân rắn có độ linh động ion tối ưu.
  • Giảng viên và sinh viên các trường đại học khối kỹ thuật, vật lý và khoa học vật liệu: Sử dụng làm tài liệu giảng dạy và học tập chuyên đề về mô phỏng động lực học phân tử và nhiệt động học chất lỏng phức tạp.

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

  • Câu hỏi 1: Cấu trúc của hệ natri silicat biến đổi như thế nào khi áp suất tăng từ 0 lên 60 GPa? Khi áp suất tăng từ 0 lên 60 GPa, hệ chuyển đổi từ cấu trúc tứ diện SiO4 (giảm từ 89% xuống 8%) sang cấu trúc bát diện SiO6 (tăng từ 0% lên 47%) thông qua đơn vị trung gian SiO5 đạt đỉnh 47% tại 40 GPa.

  • Câu hỏi 2: Nguyên nhân cốt lõi gây ra hiện tượng bất đồng nhất cấu trúc trong natri silicat lỏng là gì? Hiện tượng này bắt nguồn từ lực hút tĩnh điện mạnh giữa cation natri và các nguyên tử oxy không cầu NBO mang điện tích âm, dẫn tới việc 70% hạt natri tập trung cục bộ tại các khối simplex O4 thay vì phân bố đều trong mạng lưới.

  • Câu hỏi 3: Phương pháp phân tích Simplex đem lại ưu thế gì so với các phép tính cấu trúc truyền thống? Phương pháp Simplex cho phép phân chia không gian rỗng thành các quả cầu hình học xác định qua 4 nguyên tử, giúp định lượng chính xác kích thước lỗ xốp và nhận diện rõ nét mạng lưới kênh khuếch tán vi mô của ion kiềm.

  • Câu hỏi 4: Số phối trí trung bình của ion Natri và Silic biến động ra sao dưới áp suất cực hạn? Dưới áp lực nén 60 GPa, số phối trí trung bình của Na-O tăng mạnh từ 5,61 lên hơn 10,2, trong khi số phối trí của Si-O chuyển dịch từ mức 4 liên kết tứ diện lên mức 6 liên kết bát diện bền vững.

  • Câu hỏi 5: Kết quả nghiên cứu về kênh khuếch tán có ý nghĩa gì trong công nghệ chế tạo pin năng lượng? Việc chứng minh ion natri di chuyển ưu tiên qua các cụm NBO giúp các nhà khoa học định hướng cấu trúc mạng silicat nhằm tối ưu hóa độ dẫn ion, mở ra hướng chế tạo màng điện giải rắn hiệu suất cao cho pin natri-ion.

Kết luận

Công trình nghiên cứu đã mang lại những đóng góp khoa học nổi bật và định hình hướng phát triển tương lai:

  • Xác lập quy luật chuyển pha cấu trúc từ tứ diện SiO4 (89% tại 0 GPa) sang bát diện SiO6 (47% tại 60 GPa) qua pha trung gian SiO5.
  • Làm rõ sự chuyển dịch căn bản của trật tự tầm trung từ đơn vị OSi2 sang OSi3 (chiếm 49,36% tại 60 GPa) dưới điều kiện nén ép cao.
  • Chứng minh bằng số liệu định lượng về sự cư trú ưu tiên của 70% ion natri trong các khối simplex O4 giàu oxy không cầu.
  • Làm sáng tỏ cơ chế bất đồng nhất động học giữa nhóm nguyên tử di động nhanh và nhóm bất động thông qua phân tích hàm phân bố đám Flk(r,t).
  • Cung cấp bộ thông số thế tương tác chuẩn xác cho mô hình 3000 và 7995 nguyên tử, làm giàu cơ sở dữ liệu tính toán cho ngành vật lý chất lỏng và vô định hình.

Trong giai đoạn tiếp theo, nhóm tác giả định hướng mở rộng nghiên cứu sang các hệ kiềm silicat đa nguyên tố dưới dải nhiệt độ rộng. Các nhà khoa học và kỹ sư vật liệu có thể khai thác trực tiếp các bộ dữ liệu trong luận văn để đẩy mạnh ứng dụng trong công nghệ chế tạo thủy tinh chịu lực và pin năng lượng mới.