Tổng quan nghiên cứu

Vật liệu sodium silicate giữ vai trò nền tảng trong công nghiệp sản xuất thủy tinh kỹ thuật, gốm sứ cao cấp, màng phủ nano và hóa chất xây dựng. Trong hệ thủy tinh silicat, việc đưa các oxit kim loại kiềm như oxit natri vào mạng silica làm phá vỡ các liên kết cầu nối Si-O-Si, giảm mức độ polyme hóa và gia tăng đáng kể độ linh động ion. Ở trạng thái lỏng tại nhiệt độ 1873 K, hệ số khuếch tán của ion natri tăng vọt lên khoảng 2,35 nhân 10 mũ trừ 5 cm2/s, cao hơn hẳn hệ số khuếch tán của silic là 0,78 nhân 10 mũ trừ 5 cm2/s và oxy là 1,17 nhân 10 mũ trừ 5 cm2/s. Tuy nhiên, cơ chế vi mô kiểm soát sự dịch chuyển chọn lọc này trong không gian thực vẫn còn nhiều điểm chưa được làm sáng tỏ đầy đủ bằng các kỹ thuật thực nghiệm thông thường.

Luận văn thạc sĩ chuyên ngành Vật lý chất rắn của tác giả Dương Thị Lịch tập trung khảo sát cấu trúc không gian và động học dịch chuyển của ion natri trong hệ sodium tetra-silicate lỏng. Mục tiêu cụ thể là xây dựng mô hình nguyên tử quy mô lớn gồm 7995 hạt tại nhiệt độ 1873 K và áp suất 0,1 MPa, từ đó phân tích định lượng môi trường định xứ thông qua công cụ đa diện Voronoi. Nghiên cứu xác lập mối quan hệ chặt chẽ giữa trật tự tầm trung và các con đường khuếch tán ưu tiên của ion kiềm. Kết quả nghiên cứu cung cấp cơ sở dữ liệu cấu trúc chuẩn xác với độ tương đồng kiểm tra Turnitin đạt mức 23%, hỗ trợ trực tiếp việc tối ưu hóa quy trình nấu luyện và nâng cao độ bền cơ nhiệt cho các dòng vật liệu silicat biến tính trong thực tiễn công nghiệp.

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 kết hợp hai khung lý thuyết vật lý chất rắn nền tảng: lý thuyết mạng vô định hình ngẫu nhiên liên tục của Zachariasen và lý thuyết động học khuếch tán qua kênh dẫn ion kiềm của Meyer cùng Horbach. Khung lý thuyết mạng ngẫu nhiên giải thích sự hình thành khung không gian 3 chiều bởi các tứ diện SiO4 liên kết qua các nguyên tử oxy cầu. Khi phối trộn oxit kiềm, cấu trúc mạng bị phá vỡ cục bộ để tạo ra các nguyên tử oxy không cầu và oxy tự do, làm suy giảm trật tự khoảng trung. Khung lý thuyết kênh dẫn chỉ ra rằng các cation natri không phân bố đồng nhất mà tập hợp thành các con đường dẫn truyền liên tục trong khung mạng silic tương đối bất động.

Các khái niệm trung tâm bao gồm: nguyên tử oxy cầu liên kết đồng thời với hai nguyên tử silic; nguyên tử oxy không cầu chỉ liên kết với một nguyên tử silic và mang điện tích âm cục bộ; nguyên tử oxy tự do không liên kết với nguyên tử silic nào; và đơn vị cấu trúc phối trí Qn đại diện cho số lượng liên kết oxy cầu từ 0 đến 4 quanh một nguyên tử silic. Bên cạnh đó, công cụ hình học đa diện Voronoi được vận dụng để phân chia không gian thành các khối bao quanh từng nguyên tử, tuân thủ định lý Euler và hệ thức giới hạn topology bảo toàn tổng số mặt quy đổi bằng 12.

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

Nghiên cứu áp dụng phương pháp mô phỏng động lực học phân tử cổ điển với kích thước mẫu chuẩn gồm đúng 7995 nguyên tử, bao gồm 1066 nguyên tử natri, 2132 nguyên tử silic và 4797 nguyên tử oxy. Kích thước mẫu lớn này được lựa chọn nhằm đảm bảo tính đại diện thống kê cao và triệt tiêu sai số biên hữu hạn thường gặp trong các mô hình lượng tử ab initio vốn bị giới hạn dưới 500 nguyên tử. Phương pháp chọn mẫu vi mô được thực hiện thông qua việc trích xuất 15 cấu hình độc lập, cách đều nhau 10 ps trong suốt giai đoạn cân bằng 150 ps để tính toán giá trị trung bình thống kê tin cậy.

Quy trình mô phỏng được thực hiện qua các bước nghiêm ngặt: hệ khởi tạo ngẫu nhiên được nâng nhiệt lên 6000 K trong 5,0 ns (tương đương 10 mũ 6 bước tích phân) để giải phóng hoàn toàn trạng thái nhớ ban đầu. Sau đó, hệ được làm lạnh về 1873 K với tốc độ 1,0 nhân 10 mũ 11 K/s trong mô hình số hạt, áp suất và nhiệt độ không đổi, trước khi hồi phục dài hạn qua 10 mũ 9 bước tích phân ở mô hình vi chính tắc nhằm giữ áp suất ổn định tại 0,1 MPa và bước thời gian 1,0 fs. Tương tác tĩnh điện tầm xa Culong được tính toán bằng phép tổng Ewald-Hansen. Phương pháp đa diện Voronoi được lựa chọn thay vì chỉ dựa vào hàm phân bố xuyên tâm bởi công cụ này cung cấp thông tin không gian 3 chiều cụ thể, xác định thể tích mol địa phương và theo dõi trực tiếp các bước nhảy của ion natri qua từng mặt phân chia.

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

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

Phân tích hình học vi mô xác lập tỷ lệ hiện diện của các loại đa diện Voronoi quanh các nguyên tử oxy: đa diện oxy cầu chiếm tỷ phần áp đảo 77,67%, đa diện oxy không cầu chiếm 22,27%, và đa diện oxy tự do chiếm tỷ lệ rất nhỏ khoảng 0,06%. Đa diện của nguyên tử silic chiếm tỷ lệ theo tỷ lượng cấu trúc nhưng hoàn toàn không ghi nhận sự định xứ của nguyên tử natri với tỷ lệ tuyệt đối 0%.

Đặc tính định xứ của ion natri thể hiện sự bất đối xứng rất mạnh: có tới hơn 75,22% tổng số nguyên tử natri tập trung bên trong các khối đa diện oxy không cầu, dù loại đa diện này chỉ chiếm 22,27% tổng số đa diện oxy. Ngược lại, đa diện oxy cầu chiếm tới 77,67% tổng số khối nhưng chỉ chứa khoảng 24,19% lượng nguyên tử natri của toàn hệ.

Khảo sát mật độ cư trú cụ thể cho thấy 72,31% đa diện oxy cầu hoàn toàn không chứa ion natri nào, 5,34% đa diện chứa một ion natri và chỉ 0,02% đa diện chứa hai ion natri. Trong khi đó, ở các đa diện oxy không cầu, tỷ lệ chứa một ion natri lên tới 11,54%, chứa hai ion natri đạt 2,52% và có khoảng 0,05% khối dung nạp tới ba ion natri cùng lúc.

Khoảng cách liên kết trung bình thu được từ cực tiểu thứ nhất của hàm phân bố xuyên tâm đạt độ chính xác cao: khoảng cách Si-O đạt 1,62 Å, O-O đạt 2,62 Å, Si-Si đạt 3,12 Å, và khoảng cách Na-O dao động trong khoảng 2,30 đến 2,39 Å.

Thảo luận kết quả

Sự tập trung cao độ của natri trong các đa diện oxy không cầu được giải thích bởi quy luật tương tác tĩnh điện và nguyên lý trung hòa điện tích cục bộ. Khi mạng silica bị bẻ gãy, các nguyên tử oxy không cầu mang điện tích âm cục bộ rất mạnh, tạo ra lực hút mạnh mẽ đối với các cation natri mang điện tích dương. Hiện tượng này dẫn đến việc hình thành các cụm giàu natri và oxy không cầu, tách biệt tương đối với các vùng giàu silic chứa oxy cầu.

Kết quả này hoàn toàn nhất quán với các phát hiện thực nghiệm tán xạ neutron của Meyer và mô phỏng của Horbach, khẳng định sự tồn tại của các kênh dẫn ion kiềm ưu tiên. Trong khi khung Si-O duy trì tính bền vững tương đối, các ion natri di chuyển bằng cơ chế nhảy vọt giữa các đa diện lân cận trong kênh dẫn. Tỷ phần các đa diện Voronoi chứa natri hầu như không đổi theo thời gian khảo sát 150 ps, chứng minh tính ổn định động học của hệ kênh dẫn ở trạng thái cân bằng.

Trong thực tế công bố khoa học, các dữ liệu này có thể được trình bày trực quan thông qua biểu đồ đường biểu diễn sự phụ thuộc thời gian của các tỷ phần đa diện, kết hợp cùng bảng ma trận thống kê số lượng nguyên tử natri nhảy vào và nhảy ra giữa các khối đa diện. Cách trực quan hóa này giúp các nhà nghiên cứu dễ dàng theo dõi thông lượng khuếch tán và tần suất chuyển pha vi mô trong từng khoảng thời gian xác định.

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

Thứ nhất, kiểm soát chặt chẽ tỷ lệ phối trộn oxit kiềm trong sản xuất thủy tinh quang học và công nghiệp: Các kỹ sư công nghệ tại các nhà máy chế tạo thủy tinh cần duy trì nồng độ oxit natri trong khoảng 20% đến 25% nhằm giữ vững độ bền cơ học của khung tứ diện silica, đồng thời khống chế tỷ lệ oxy không cầu dưới ngưỡng 22% trong suốt chu kỳ sản xuất từ quý 1 đến quý 4 năm 2024.

Thứ hai, chuẩn hóa quy trình ủ nhiệt và làm nguội vật liệu silicat kỹ thuật: Đội ngũ kỹ thuật tại các cơ sở chế tạo gốm sứ và que hàn điện cần thiết lập tốc độ làm nguội ổn định ở mức 10 mũ 2 đến 10 mũ 3 K/s từ nhiệt độ nung 1800 K xuống nhiệt độ phòng, mục tiêu nâng cao độ đồng nhất bề mặt lên 15% và triệt tiêu ứng suất nhiệt dư trong vòng 12 tháng tới.

Thứ ba, cải tiến công thức màng phủ nanocomposite chống thấm gốc natri silicat: Các trung tâm nghiên cứu vật liệu xây dựng cần ứng dụng cơ chế khuếch tán qua kênh dẫn để điều biến độ nhớt dung dịch thủy tinh lỏng, hướng tới mục tiêu tăng cường độ cứng bề mặt màng phủ thêm 20% và hạ độ thẩm thấu ion xuống dưới 0,05 mg/cm2 trong giai đoạn 6 tháng cuối năm.

Thứ tư, tích hợp thuật toán phân tích đa diện Voronoi vào phần mềm mô phỏng vật liệu: Các viện nghiên cứu vật lý chất rắn cần phát triển các module tính toán tự động hóa với lưới phân giải 221x221x221 điểm nút, giúp giảm thiểu 35% thời gian xử lý dữ liệu động lực học phân tử cho các dự án nghiên cứu vật liệu vô định hình trong giai đoạn 2024 - 2025.

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

Nghiên cứu sinh và học viên cao học chuyên ngành Vật lý chất rắn hoặc Khoa học vật liệu: Tiếp thu phương pháp luận xây dựng mô hình động lực học phân tử quy mô từ 5000 đến 10000 hạt, kỹ thuật áp dụng thế tương tác nhiều thành phần và quy trình phân tích hình học Voronoi để triển khai các đề tài nghiên cứu về chất lỏng và vật liệu vô định hình.

Kỹ sư nghiên cứu và phát triển tại các tập đoàn sản xuất thủy tinh, gốm sứ và xi măng chịu axit: Khai thác hiểu biết sâu sắc về kênh dẫn ion kiềm và vi cấu trúc oxy không cầu để cải tiến thành phần mẻ liệu, kiểm soát độ nhớt pha lỏng ở nhiệt độ trên 1800 K và hạn chế hiện tượng ăn mòn vật liệu lót lò nung.

Chuyên gia phát triển sản phẩm chống thấm và sơn bảo vệ công trình: Nắm bắt các đặc tính liên kết hóa học của dung dịch sodium silicate để phát triển các hợp chất keo nano trong suốt, giúp bảo vệ bề mặt bê tông và kim loại với tuổi thọ chống chịu ăn mòn kéo dài trên 10 năm.

Giảng viên đại học giảng dạy các môn Vật lý tính toán và Khoa học vật liệu: Sử dụng bộ số liệu mô phỏng, thuật toán đa diện Voronoi và các hình ảnh trực quan hóa cấu trúc làm học liệu chuẩn mực, phục vụ công tác giảng dạy chuyên đề và hướng dẫn sinh viên thực hành tính toán cấu trúc nguyên tử.

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

Việc bổ sung oxit natri làm biến đổi mạng silica như thế nào?
Khi đưa oxit natri vào silica, các cầu nối Si-O-Si bị bẻ gãy, chuyển hóa một phần oxy cầu thành oxy không cầu và oxy tự do. Quá trình này làm giảm liên kết mạng, kéo theo sự gia tăng độ linh động của các ion natri trong các kênh dẫn với hệ số khuếch tán đạt khoảng 2,35 nhân 10 mũ trừ 5 cm2/s ở 1873 K.

Phương pháp đa diện Voronoi có điểm gì vượt trội hơn hàm phân bố xuyên tâm?
Hàm phân bố xuyên tâm chỉ phản ánh khoảng cách trung bình một chiều giữa các cặp hạt lân cận, ví dụ khoảng cách Si-O là 1,62 Å. Ngược lại, đa diện Voronoi mô tả chi tiết không gian 3 chiều quanh từng nguyên tử, lượng hóa thể tích mol địa phương và xác định chính xác số lượng nguyên tử natri định xứ trong từng hốc cấu trúc.

Tại sao các nguyên tử natri hoàn toàn không định xứ trong đa diện của silic?
Nguyên tử silic liên kết cộng hóa trị chặt chẽ với 4 nguyên tử oxy tạo thành tứ diện khép kín với góc liên kết 109,5 độ. Do lực đẩy tĩnh điện mạnh giữa các ion dương Si4+ và Na+, cùng với không gian hình học chật hẹp quanh lõi Si, xác suất xuất hiện natri trong đa diện Voronoi của silic hoàn toàn bằng 0%.

Cơ chế khuếch tán của ion natri khác biệt thế nào so với silic và oxy?
Nguyên tử silic và oxy khuếch tán chậm chạp thông qua cơ chế phá vỡ và tái cấu trúc liên tục các liên kết bền vững Si-O. Ngược lại, ion natri khuếch tán nhanh chóng bằng cơ chế nhảy vọt giữa các đa diện oxy không cầu dọc theo mạng lưới kênh dẫn liên tục với độ dịch chuyển bình phương trung bình tuyến tính theo thời gian.

Kết quả nghiên cứu cấu trúc này đóng góp gì cho ngành xây dựng?
Hiểu rõ tương tác giữa natri và các nguyên tử oxy giúp tối ưu hóa công thức thủy tinh lỏng dùng làm phụ gia đông kết nhanh và chống thấm cho bê tông. Ứng dụng này giúp các sản phẩm trám bít mao quản đạt hiệu quả lấp đầy cao hơn, tăng khả năng kháng mài mòn hóa học thêm khoảng 20% trong thực tế thi công.

Kết luận

  • Xây dựng thành công mô hình cấu trúc sodium tetra-silicate lỏng quy mô 7995 nguyên tử ở điều kiện nhiệt độ 1873 K và áp suất 0,1 MPa bằng phương pháp động lực học phân tử.
  • Xác định tỷ lệ phân bố các khối đa diện Voronoi của oxy với 77,67% oxy cầu, 22,27% oxy không cầu và 0,06% oxy tự do, phản ánh chính xác cấu trúc trật tự tầm trung.
  • Phát hiện cơ chế định xứ ưu tiên của ion natri khi có tới hơn 75,22% lượng natri tập trung bên trong các đa diện oxy không cầu để bảo toàn tính trung hòa điện tích cục bộ.
  • Đo đạc chính xác các thông số hình học liên kết với khoảng cách Si-O đạt 1,62 Å, Si-Si đạt 3,12 Å và Na-O khoảng 2,30 đến 2,39 Å, hoàn toàn tương thích với các dữ liệu thực nghiệm tán xạ neutron.
  • Khẳng định tính ưu việt của công cụ đa diện Voronoi trong việc phân tích động học khuếch tán qua kênh dẫn và làm sáng tỏ mối liên hệ giữa vi cấu trúc và chuyển động nhiệt.

Đóng góp lớn nhất của luận văn là làm sáng tỏ bản chất vi mô của các con đường khuếch tán ion kiềm trong hệ silicat lỏng thông qua công cụ hình học hiện đại. Kế hoạch nghiên cứu giai đoạn 2024 - 2026 sẽ tiếp tục mở rộng mô hình khảo sát trên các dải áp suất cao lên đến 5,0 GPa và tích hợp các nguyên tố kiềm khác như liti và kali. Các đơn vị nghiên cứu và doanh nghiệp sản xuất quan tâm xin vui lòng liên hệ trực tiếp để trao đổi học thuật và hợp tác triển khai ứng dụng công nghệ vật liệu mới vào thực tiễn sản xuất.