Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu cấu trúc và sự không đồng nhất động học trong vật liệu Silicát ba nguyên PbO.2SiO2 và Na2O.2SiO2 ở trạng thái lỏng và vô định hình" của nghiên cứu sinh Nguyễn Văn Yên, chuyên ngành Vật lý kỹ thuật (Mã số: 62520401), được thực hiện tại Trường Đại học Bách khoa Hà Nội dưới sự hướng dẫn khoa học của PGS. Lê Thế Vinh và PGS. Nguyễn Văn Hồng (năm 2017).

Vật liệu silicát đóng vai trò quan trọng trong nhiều ngành công nghiệp như vi điện tử, quang học, gốm sứ, vật liệu chịu lửa, công nghiệp cao su và che chắn phóng xạ hạt nhân (đặc biệt là chì-silicát). Cấu trúc của mạng silicát thường được mô tả bằng mô hình mạng ngẫu nhiên liên tục tạo bởi các đơn vị tứ diện $\text{SiO}_4$ kết nối qua nguyên tử oxy cầu. Khi pha tạp các ôxít kim loại kiềm, kiềm thổ hoặc kim loại chuyển tiếp vào mạng silica, các cation pha tạp đóng vai trò biến tính mạng hoặc tham gia tạo mạng, dẫn đến sự xuất hiện của oxy không cầu (NBO), oxy ba liên kết (tricluster O), các đơn vị phối trí cao ($\text{SiO}_5$, $\text{SiO}_6$, $\text{AlO}_5$, $\text{AlO}_6$) và tạo ra sự phân bố không đồng nhất về mặt không gian cũng như thành phần hóa học.

Khoảng trống nghiên cứu then chốt mà luận án xác định bao gồm:

  1. Thiếu các mô tả định lượng chi tiết về sự phân bố không gian của các đơn vị cấu trúc trong hệ silicát ba nguyên; chưa làm rõ liệu các cation phân bố đồng đều hay phân tách thành các cụm riêng biệt, kích thước và quy luật kết cụm của các mạng con (subnets) dưới tác động của áp suất và nồng độ pha tạp.
  2. Hiện tượng động học không đồng nhất (Dynamical Heterogeneity - DH) trong các hệ ôxít silicát lỏng còn rất ít được nghiên cứu định lượng. Trước luận án, các nghiên cứu (như công trình của K. Vargheese năm 2010 trên hệ nhôm-silicát) mới chỉ dừng ở mức định tính; chưa làm rõ kích thước, mật độ, số lượng liên kết của các cụm nguyên tử chuyển động nhanh và chuyển động chậm, cũng như quy luật suy giảm liên kết theo thời gian.
  3. Mối tương quan vi mô giữa cấu trúc hình học địa phương và các quá trình động học khuếch tán dị thường trong hệ silicát ba nguyên chưa được giải thích đầy đủ.

Mục tiêu nghiên cứu của luận án được xác định gồm 4 nội dung:

  1. Nghiên cứu vi cấu trúc của hệ vật liệu $\text{PbO.2SiO}_2$ và $\text{Al}_2\text{O}_3.2\text{SiO}_2$ (trạng thái lỏng và vô định hình) dưới ảnh hưởng của áp suất và nồng độ thông qua các hàm phân bố xuyên tâm, số phối trí và phân bố góc liên kết.
  2. Nghiên cứu cấu trúc không đồng nhất của hệ vật liệu $\text{Al}_2\text{O}_3.2\text{SiO}_2$ và $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng thông qua việc đề xuất và áp dụng hai phương pháp phân tích hình học mới: phương pháp quả cầu Simplex và phương pháp hạt Vỏ - Lõi (Shell-Core).
  3. Nghiên cứu động học không đồng nhất (DH) của hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$ và $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng qua việc phân tích định lượng chuyển động của nhóm 10% nguyên tử nhanh nhất và 10% nguyên tử chậm nhất so với nhóm nguyên tử ngẫu nhiên.
  4. Làm sáng tỏ mối tương quan giữa đặc trưng cấu trúc địa phương và các quá trình động học (hệ số khuếch tán, độ nhớt, khuếch tán dị thường, thời gian sống của các hạt siêu phân tử) trong các hệ silicát lỏng.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng: Các hệ silicát ba nguyên gồm chì-silicát ($\text{PbO.SiO}_2$), nhôm-silicát ($\text{Al}_2\text{O}_3.2\text{SiO}_2$, ký hiệu $\text{AS}_2$, và hệ $x\text{Al}_2\text{O}_3(1-x)\text{SiO}_2$) và natri-silicát ($\text{Na}_2\text{O}.2\text{SiO}_2$) ở trạng thái lỏng và vô định hình.
  • Phạm vi nghiên cứu: Khảo sát cấu trúc vi mô và động học phân tử trong không gian mô phỏng kích thước từ 3.000 đến 7.992 nguyên tử; dải áp suất từ 0 đến 60 GPa; thang nhiệt độ từ 600 K (vô định hình), 1.873 K, 3.200 K đến 3.500 K (trạng thái lỏng).

Tổng quan tài liệu và vị trí của luận án

Luận án đã hệ thống hóa các nghiên cứu thực nghiệm và mô phỏng lý thuyết về cấu trúc và động học của silica và ba nhóm vật liệu silicát tiêu biểu:

  1. Cấu trúc của silica ($\text{SiO}_2$): Wright (1988, 1994) xác định phân bố góc $\text{Si-O-Si}$ trong khoảng $120^\circ - 180^\circ$ với đỉnh cực đại tại $144^\circ$. Mauri và cộng sự (NMR) xác định góc phân bố $\text{Si-O-Si}$ có cực đại trong khoảng $150^\circ - 155^\circ$. Geissberger và Bray xác định góc $144^\circ$, Clark xác định $147^\circ$. Galeener và Gibbs nghiên cứu năng lượng mạch vòng và chỉ ra vòng 4 tứ diện $\text{SiO}_4$ có năng lượng cấu hình bền vững nhất. Meade và cộng sự, McMillan và Wolf khảo sát ảnh hưởng của quá trình nén áp suất cao (từ 0 đến 28 GPa và trên 40 GPa), chỉ ra sự gia tăng khoảng cách $\text{Si-O}$ từ 1,62 Å lên 1,66 Å ở 42 GPa và sự chuyển đổi số phối trí của Si từ 4 lên 5 và 6.

  2. Hệ ôxít nhôm-silicát ($\text{Al}_2\text{O}_3\text{-SiO}_2$): Winkler và cộng sự, Võ Văn Hoàng mô phỏng cấu trúc địa phương của $\text{Al}_2\text{O}_3.2\text{SiO}_2$, ghi nhận sự tồn tại của các đa diện $\text{AlO}_4$, $\text{AlO}_5$, $\text{AlO}_6$, các liên kết chung cạnh và nguyên tử oxy ba liên kết (tricluster O). Takahiro Takei và cộng sự khảo sát tỷ lệ các loại nguyên tử oxy liên kết với cation ($^{[2]}\text{O}$, $^{[3]}\text{O}$, $^{[4]}\text{O}$, $^{[5]}\text{O}$) theo hàm lượng $\text{Al}_2\text{O}_3$. Adama Tandia và cộng sự chỉ ra cơ chế khuếch tán gắn liền với sự đứt gãy và tái tạo liên kết giữa NBO và oxy tricluster. Yanhui Liu nghiên cứu độ nhớt của hệ xỉ thông qua tỷ số $\text{NBO/T}$ và mức độ polymer hóa $\text{DOP} = (\text{Q}^3+\text{Q}^4)/(\text{Q}^1+\text{Q}^2)$.

  3. Hệ ôxít chì-silicát ($\text{PbO-SiO}_2$): Bair (1936) thực hiện nghiên cứu thực nghiệm đầu tiên bằng tia X. Takaishi và cộng sự khảo sát hệ $\text{PbO.SiO}_2$ (25–90 mol% PbO), phát hiện khi tăng nồng độ PbO thì số phối trí của Pb giảm từ 6 xuống 3; ở nồng độ cao ($> 50%$), hình thành chuỗi polymer $(\text{PbO}_3)_n$ và Pb đóng vai trò nguyên tử tạo mạng. Kohara và cộng sự, Alderman và cộng sự (2013), Yoko và Fayon, Witkowska (EXAFS), Rybicki và cộng sự (MD và EXAFS) chỉ ra vai trò kép của Pb (vừa tạo mạng vừa làm biến tính mạng $\text{Si-O}$). B. Langanke và cộng sự, Mizuno khảo sát hệ số khuếch tán và độ nhớt theo hàm lượng $\text{SiO}_2$. Golubkov và cộng sự dự đoán sự tồn tại của vùng động học không đồng nhất kích thước $5 - 7\text{ \AA}$ ở nồng độ PbO trên 40%.

  4. Hệ ôxít natri-silicát ($\text{Na}_2\text{O-SiO}_2$): Warren và Biscoe (1938) chỉ ra sự hình thành oxy không cầu (NBO) khi đưa $\text{Na}_2\text{O}$ vào mạng silica, với khoảng cách $\text{Si-O} \approx 1,62\text{ \AA}$ và $\text{Na-O} \approx 2,35\text{ \AA}$. Misawa và cộng sự (tán xạ neutron), Henderson (EXAFS), Vedishcheva và cộng sự khẳng định sự phân bố không đồng nhất của ion $\text{Na}^+$ xung quanh các nguyên tử NBO tạo ra các kênh dẫn khuếch tán. Horbach và cộng sự phát hiện đỉnh tiền cấu trúc ở thừa số cấu trúc gắn liền với kênh dẫn ion.

  5. Động học không đồng nhất (DH): Antonio M. Puertas và Claudio Donati nghiên cứu DH trong hệ keo và chất lỏng nguội nhanh, chứng minh các hạt chuyển động nhanh có xu hướng kết cụm trong vùng mật độ thấp và các hạt chậm kết cụm trong vùng mật độ cao. K. Vargheese (2010) quan sát định tính DH trong nhôm-silicát lỏng.

Vị trí của luận án: Kế thừa các kết quả thực nghiệm và mô phỏng trên, luận án giải quyết khoảng trống về phân tích cấu trúc định lượng bằng cách phát triển hai công cụ hình học mới (Simplex và Shell-Core), đồng thời thực hiện khảo sát định lượng toàn diện về động học không đồng nhất (DH), phân lập các nhóm 10% hạt nhanh/chậm và giải thích cơ chế khuếch tán dị thường trong các hệ silicát ba nguyên lỏng và vô định hình dưới tác động của áp suất.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án sử dụng phương pháp mô phỏng Động lực học phân tử (Molecular Dynamics - MD) cổ điển, kết hợp các thuật toán phân tích hình học cấu trúc vi mô và kỹ thuật thống kê động học phân tử.

Khung lý thuyết và mô hình thế tương tác

  1. Thế tương tác Born-Mayer hai thành phần: Sử dụng cho hệ chì-silicát ($\text{PbO.SiO}2$) và nhôm-silicát ($\text{Al}2\text{O}3.2\text{SiO}2$). Biểu thức thế năng tương tác cặp $U(r{ij})$ có dạng: $$U(r{ij}) = \frac{q_i q_j}{r{ij}} + A{ij} \exp(-B_{ij} r_{ij})$$ Trong đó $q_i, q_j$ là điện tích hiệu dụng của các ion ($q_{\text{Si}} = +4$, $q_{\text{O}} = -2$, $q_{\text{Pb}} = +2$, $q_{\text{Al}} = +3$); $A_{ij}$ và $B_{ij}$ là các thông số thế tương tác đặc trưng cho từng cặp nguyên tử.
Hệ vật liệu Cặp tương tác $A_{ij}$ (eV) $B_{ij}$ ($\text{\AA}^{-1}$)
$\text{PbO.SiO}_2$ Si - Si 0,0000 0,0000
Si - O 1157,4483 3,4482
O - O 1003,3224 2,7778
Pb - Pb 0,0000 0,0000
Pb - O 895,3000 2,9412
Pb - Si 0,0000 0,0000
$\text{Al}_2\text{O}_3.2\text{SiO}_2$ Si - Si 0,0000 0,0000
Si - O 1844,0000 3,8460
O - O 1442,0000 3,4480
Al - Al 0,0000 0,0000
Al - O 1450,0000 3,6400
Al - Si 0,0000 0,0000
  1. Thế tương tác hai và ba thành phần cho hệ $\text{Na}_2\text{O}.2\text{SiO}_2$: Kết hợp thế hai hạt mở rộng (tương tác tĩnh điện, lực đẩy Born-Mayer, phân tán Van der Waals và liên kết cộng hóa trị) và thế ba hạt mô tả năng lượng uốn góc liên kết $\text{Si-O-Si}$: $$U_{ij} = \frac{q_i q_j}{r_{ij}} + f_0(b_i + b_j) \exp\left(\frac{a_i + a_j - r_{ij}}{b_i + b_j}\right) - \frac{c_i c_j}{r_{ij}^6} + D_{1ij}\exp(-\beta_{1ij} r_{ij}) - D_{2ij}\exp(-\beta_{2ij} r_{ij})$$ $$U_{ijk} = -f \left[\cos{2(\theta_{ijk} - \theta_0)} - 1\right] \exp\left[-g_r(r_{ij} - r_m)^2 - g_r(r_{ik} - r_m)^2\right]$$ với $f_0 = 41,865\text{ kJ}\cdot\text{nm}^{-1}\text{mol}^{-1}$, các thông số $a, b, c, D_1, D_2, \beta_1, \beta_2, f, \theta_0, r_m, g_r$ được nạp theo bộ thông số chuẩn của mô hình BKS biến tính.

Quy trình mô phỏng và phần mềm

  • Hệ $\text{PbO.SiO}_2$ và $\text{Al}_2\text{O}_3.2\text{SiO}_2$: Lập trình tính toán bằng ngôn ngữ C trên môi trường Linux, thuật toán tích phân Verlet, điều kiện biên tuần hoàn 3 chiều, bước thời gian mô phỏng $\Delta t = 1,0 - 1,5\text{ fs}$. Khởi tạo ngẫu nhiên ở 6.000 K trong 50 ps để xóa nhớ cấu hình, làm lạnh đẳng nhiệt - đẳng áp (NPT) qua $10^7$ bước mô phỏng về nhiệt độ khảo sát (3.200 K đối với $\text{PbO.SiO}_2$; 3.500 K và 600 K đối với $\text{AS}_2$), sau đó chuyển sang tập hợp vi chính tắc (NVE) để thu thập dữ liệu cân bằng.
  • Hệ $\text{Na}_2\text{O}.2\text{SiO}_2$: Tính toán bằng phần mềm Mxdorto (phát triển bởi Fumiya Noritake) tại Viện RIKEN (Nhật Bản) với 7.992 nguyên tử (1.776 Na, 1.776 Si, 4.440 O) ở 1.873 K tại 0 GPa và 8 GPa.

Phương pháp phân tích cấu trúc và động học

  1. Phân tích cấu trúc cổ điển:
  • Hàm phân bố xuyên tâm $g(r) = \frac{V}{N \cdot 4\pi r^2 dr} \langle dn(r) \rangle$.
  • Số phối trí trung bình tính qua tích phân đỉnh thứ nhất của $g(r)$ đến bán kính ngắt $r_{\text{cut}}$.
  • Hàm phân bố góc liên kết (BAD) cho các góc trong đơn vị ($\text{O-T-O}$) và góc nối giữa các đơn vị ($\text{T-O-T}$ với $\text{T} = \text{Si}, \text{Al}, \text{Pb}$).
  1. Phương pháp phân tích cấu trúc mới:
  • Phương pháp Simplex: Dựng các mặt cầu ngoại tiếp đi qua 4 nguyên tử bất kỳ trong mô hình. Phân loại thành: Void-simplex (VS - mặt cầu rỗng không chứa nguyên tử nào bên trong), Oxy-simplex (OS - chỉ chứa nguyên tử O bên trong), Cation-simplex (CS - chỉ chứa cation bên trong), và Cation-simplex-cluster (CSC - cụm các CS kết nối với nhau).
  • Phương pháp Shell-Core (SC): Xác định các hạt Shell-Core (SCP) gồm lõi (Core) chứa một hoặc nhiều cation và vỏ (Shell) chứa các nguyên tử oxy bao quanh; tập hợp các SCP liền kề tạo thành cụm Shell-Core-cluster (SCC).
  1. Phương pháp phân tích động học và Động học không đồng nhất (DH):
  • Dịch chuyển bình phương trung bình: $\langle r^2(t) \rangle = \frac{1}{N} \sum_{i=1}^N \langle |\vec{r}_i(t) - \vec{r}_i(0)|^2 \rangle$.
  • Hệ số khuếch tán $D$ xác định qua hệ thức Einstein: $D = \lim_{t \to \infty} \frac{\langle r^2(t) \rangle}{6t}$.
  • Độ nhớt $\eta$ tính toán qua hệ thức Eyring: $\eta = \frac{k_B T}{\lambda D}$ (với $\lambda = 2R$ là bước nhảy khuếch tán ion).
  • Đánh giá DH: Xác định bán kính ngắt liên kết $r_{lk}$ ($r_{lk} = 3,65\text{ \AA}$ cho $\text{Na}_2\text{O}.2\text{SiO}2$; $r{lk} = 4,00\text{ \AA}$ cho $\text{Al}_2\text{O}_3.2\text{SiO}_2$). Chọn 10% nguyên tử có dịch chuyển lớn nhất (nhanh nhất) và 10% nguyên tử có dịch chuyển nhỏ nhất (chậm nhất) của từng loại nguyên tử (O, Si, Al, Na), so sánh với 10% nguyên tử ngẫu nhiên qua các đại lượng:
    • Số lượng cụm $N_C$: $N_C(\text{nhanh/chậm}) < N_C(\text{ngẫu nhiên})$.
    • Số liên kết trung bình trên mỗi nguyên tử $\langle N_{LK} \rangle$: $\langle N_{LK}\rangle(\text{nhanh/chậm}) > \langle N_{LK}\rangle(\text{ngẫu nhiên})$.
    • Kích thước cụm trung bình $\langle S_C \rangle$: $\langle S_C \rangle(\text{nhanh/chậm}) > \langle S_C \rangle(\text{ngẫu nhiên})$.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 tổng thuật các cơ sở lý thuyết và các kết quả nghiên cứu trong và ngoài nước về cấu trúc mạng tứ diện của silica ($\text{SiO}_2$), hệ nhôm-silicát, hệ chì-silicát và hệ natri-silicát. Phân tích chi tiết cơ chế chuyển đổi giữa oxy cầu (BO), oxy không cầu (NBO) và tricluster O ảnh hưởng đến dòng chảy nhớt và quá trình khuếch tán. Trình bày tổng quan về hiện tượng động học không đồng nhất (DH) trong các hệ chất lỏng và nêu rõ khoảng trống nghiên cứu trong các hệ ôxít silicát lỏng ba nguyên.

Chương 2: Phương pháp tính toán

Chương 2 trình bày chi tiết cơ chế xây dựng mô hình mô phỏng MD cho 3 hệ vật liệu:

  • Hệ $\text{PbO.SiO}_2$ (5.000 hạt, $T = 3.200\text{ K}$, dải áp suất $0 - 35\text{ GPa}$, 8 mô hình).
  • Hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$ lỏng (5.500 hạt, $T = 3.500\text{ K}$, dải áp suất $0 - 20\text{ GPa}$, 5 mô hình).
  • Hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$ vô định hình (3.000 hạt, $T = 600\text{ K}$, dải áp suất $0 - 60\text{ GPa}$, 13 mô hình).
  • Hệ $x\text{Al}_2\text{O}_3(1-x)\text{SiO}_2$ lỏng với 4 tỷ phần mol: AS1 ($25%\text{ Al}_2\text{O}_3$, 4.997 hạt), AS2 ($33,3%$), AS3 ($50%$), AS4 ($60%$) ở 3.500 K và 0 GPa.
  • Hệ $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng (7.992 hạt, $T = 1.873\text{ K}$, tại 0,1 MPa và 8 GPa).

Trình bày phương pháp tính các thông số động học ($D, \eta$), phương pháp đánh giá định lượng DH qua nhóm 10% hạt và các thuật toán hình học phân tích cấu trúc không đồng nhất (Simplex và Shell-Core).

Hệ mô phỏng Số hạt ($N$) Nhiệt độ ($T$) Áp suất ($P$) Mật độ ($\text{g/cm}^3$) tại 0 GPa
$\text{PbO.SiO}_2$ (lỏng) 5.000 3.200 K 0 - 35 GPa 5,42 ($0\text{ GPa}$) $\to$ 8,96 ($35\text{ GPa}$)
$\text{AS}_2$ (lỏng) 5.500 3.500 K 0 - 20 GPa 2,55 ($0\text{ GPa}$) $\to$ 3,78 ($20\text{ GPa}$)
$\text{AS}_2$ (vô định hình) 3.000 600 K 0 - 60 GPa 2,68 ($0\text{ GPa}$) $\to$ 4,69 ($60\text{ GPa}$)
$\text{Na}_2\text{O}.2\text{SiO}_2$ (lỏng) 7.992 1.873 K 0 - 8 GPa 2,34 ($0\text{ GPa}$) $\to$ 3,12 ($8\text{ GPa}$)

Chương 3: Cấu trúc của hệ nhôm-silicát và chì-silicát

Chương 3 phân tích cấu trúc trật tự gần và trật tự khoảng trung của các hệ:

  1. Trật tự gần và sự chuyển pha phối trí:
  • Trong hệ $\text{PbO.SiO}_2$ lỏng, ở 0 GPa tồn tại chủ yếu các đơn vị $\text{SiO}_4$ và $\text{PbO}_3, \text{PbO}_4$. Khi nén áp suất lên 35 GPa, số phối trí của $\text{Si-O}$ tăng từ 4 lên 5 và 6 ($\text{SiO}_5, \text{SiO}6$), trong khi số phối trí của $\text{Pb-O}$ tăng từ 3,8 lên trên 6,2. Tỷ lệ các loại oxy: oxy liên kết cả Pb và Si ($\text{O}T$) chiếm ưu thế, bên cạnh các oxy chỉ liên kết Si ($\text{O}{\text{Si}}$) và oxy chỉ liên kết Pb ($\text{O}{\text{Pb}}$).
  • Trong hệ $\text{AS}_2$ lỏng và vô định hình, ở 0 GPa các đơn vị cơ bản là $\text{SiO}_4$ và $\text{AlO}_4$. Khi tăng áp suất lên 20 GPa (trạng thái lỏng) và 60 GPa (trạng thái vô định hình), diễn ra quá trình chuyển pha liên tục: $\text{TO}_4 \to \text{TO}_5 \to \text{TO}_6$ (với $\text{T} = \text{Si}, \text{Al}$). Tốc độ chuyển đổi số phối trí của Al diễn ra nhanh hơn và ở áp suất thấp hơn so với Si.
  1. Trật tự khoảng trung và kiểu liên kết đa diện:
  • Luận án thống kê chi tiết phân bố 3 kiểu liên kết không gian giữa hai khối đa diện $\text{TO}_x$ liền kề: liên kết chung góc ($N_c$), liên kết chung cạnh ($N_e$) và liên kết chung mặt ($N_f$). Ở áp suất thấp ($0\text{ GPa}$), liên kết chung góc $N_c$ chiếm ưu thế tuyệt đối ($> 90%$). Khi tăng áp suất, liên kết chung cạnh $N_e$ và chung mặt $N_f$ tăng mạnh (đặc biệt giữa các cặp $\text{Al-Al}$ và $\text{Al-Si}$).
  • Phân tích mạng con (subnet): Trong $\text{AS}2$ lỏng tồn tại mạng con $\text{Si-O}$ kích thước rất lớn ($S{mc}$ lớn, số lượng mạng con $N_{mc}$ nhỏ), bao phủ toàn bộ không gian mô hình; ngược lại, mạng con $\text{Al-O}$ gồm nhiều cụm kích thước nhỏ hơn ($S_{mc}$ nhỏ, $N_{mc}$ lớn) phân bố lồng vào mạng $\text{Si-O}$.

Chương 4: Phân tích cấu trúc bằng phương pháp Simplex và Shell-Core

Chương 4 áp dụng hai phương pháp hình học mới để làm rõ tính không đồng nhất về mặt vi cấu trúc:

  1. Phương pháp Simplex:
  • Trong hệ $\text{Na}_2\text{O}.2\text{SiO}_2$: Xác định các loại quả cầu Void-simplex (VS), Oxy-simplex (OS) và Cation-simplex (CS). Số lượng OS trên mỗi nguyên tử ($m_{\text{OS}}$) và bán kính trung bình $R_{\text{OS}}$ cho thấy sự tập trung của các nguyên tử oxy tạo thành các vùng giàu oxy. Các CS chứa các ion $\text{Na}^+$ và $\text{Si}^{4+}$ kết nối thành các dải cụm Cation-simplex-cluster (CSC).
  • Trong hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$: Khi tăng áp suất từ 0 lên 20 GPa, bán kính trung bình $R_{\text{OS}}$ giảm từ 1,85 Å xuống 1,52 Å; số lượng OS và CS tăng lên. Cụm CSC mở rộng kích thước $S_{\text{CSC}}$, tỷ lệ nguyên tử $\text{Al/Si}$ trong các CSC lớn phản ánh sự phân bố không đồng nhất của cation trong mạng.
  1. Phương pháp Shell-Core (SC) trong $\text{AS}_2$:
  • Phân loại các hạt SCP theo số cation trong lõi (1 cation hoặc 2 cation) và số oxy ở lớp vỏ (từ 3 đến 8 oxy).
  • Khi tăng áp suất từ 0 đến 20 GPa, bán kính lõi giảm, số lượng oxy trong vỏ tăng (phù hợp với sự tăng số phối trí).
  • Phân tích các cụm lớn SCC ($k > 6$ hạt SCP): Tỷ lệ nồng độ $C_{\text{O}}/(C_{\text{Al}} + C_{\text{Si}})$ và $C_{\text{Al}}/(C_{\text{Al}} + C_{\text{Si}})$ trong các SCC chứng minh sự tồn tại của các vùng giàu cation nhôm và các vùng giàu silica tách biệt trong cấu trúc không gian ba chiều.

Chương 5: Động học không đồng nhất của Natri và Nhôm Silicát

Chương 5 phân tích định lượng DH và mối quan hệ giữa cấu trúc và động học:

  1. Động học không đồng nhất trong $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng:
  • Hệ số khuếch tán của $\text{Na}^+$ ($D_{\text{Na}} \approx 1,2 \times 10^{-5}\text{ cm}^2/\text{s}$) lớn hơn hệ số khuếch tán của Si và O khoảng 100 lần ($D_{\text{Si}}, D_{\text{O}} \approx 10^{-7}\text{ cm}^2/\text{s}$).
  • Khảo sát nhóm 10% nguyên tử nhanh nhất và 10% chậm nhất: 10% nguyên tử Na nhanh nhất có số cụm $N_C$ nhỏ hơn rõ rệt và kích thước cụm $\langle S_C \rangle$, số liên kết $\langle N_{LK} \rangle$ lớn hơn nhóm ngẫu nhiên, chứng tỏ các nguyên tử Na nhanh kết cụm mạnh với nhau tạo thành các kênh dẫn khuếch tán dọc theo các nguyên tử NBO.
  • Trong hệ natri-silicát, vùng chuyển động nhanh là vùng có mật độ nguyên tử O thấp (giàu Na, giàu NBO), trong khi vùng chuyển động chậm là vùng mạng con $\text{Si-O}$ có mật độ cao và liên kết bền vững.
  1. Động học không đồng nhất và khuếch tán dị thường trong $\text{Al}_2\text{O}_3.2\text{SiO}_2$ lỏng:
  • Hệ số khuếch tán tuân theo thứ tự: ở áp suất thấp (0 GPa), $D_{\text{Al}} > D_{\text{O}} > D_{\text{Si}}$; ở áp suất cao (20 GPa), đảo trật tự thành $D_{\text{O}} > D_{\text{Al}} > D_{\text{Si}}$.
  • Dưới quá trình nén đẳng nhiệt (3.500 K), hệ số khuếch tán của tất cả các nguyên tử tăng lên, đạt cực đại tại vùng áp suất trung gian rồi sau đó giảm khi nén lên 20 GPa (hiện tượng khuếch tán dị thường). Hiện tượng này tương ứng trực tiếp với sự tăng và đạt cực đại của tỷ phần đơn vị phối trí 5 ($\text{SiO}_5$).
  • Khảo sát DH: 10% nguyên tử O nhanh nhất kết cụm với nhau tạo thành các vùng có mật độ nguyên tử O cao; 10% nguyên tử O chậm nhất nằm trong các vùng giàu Si có mật độ O thấp (ngược với hệ natri-silicát).
  • Tương quan cấu trúc - động học: Luận án xác định sự tồn tại của các hạt siêu phân tử (SM) và hạt cứng (SLP). Thời gian sống của các cấu trúc liên kết phụ thuộc vào độ bền liên kết: liên kết $\text{Si-O}$ bền vững quyết định cấu trúc khung chuyển động chậm, trong khi liên kết $\text{Al-O}$ kém bền hơn dễ bị đứt gãy và tái tạo, chi phối động học khuếch tán nhanh.
Hệ lỏng Nguyên tử $D \times 10^{-6}\text{ cm}^2/\text{s}$ ($0\text{ GPa}$) $D \times 10^{-6}\text{ cm}^2/\text{s}$ ($20\text{ GPa}$) Đặc trưng vùng chuyển động nhanh
$\text{PbO.SiO}_2$ (3.200 K) Si 1,34 0,21 Khuếch tán chi phối bởi Si ở mẫu giàu PbO
O 1,83 0,35 Chuyển động tương quan Si-O
Pb 2,45 0,48 Độ linh động cation Pb cao ở áp suất thấp
$\text{AS}_2$ (3.500 K) Si 0,82 0,45 Khuếch tán dị thường (cực đại ở áp suất trung gian)
O 1,15 1,12 Vùng nhanh có mật độ nguyên tử O cao
Al 1,46 0,98 Liên kết Al-O kém bền hơn Si-O
$\text{Na}_2\text{O}.2\text{SiO}_2$ (1.873 K) Si 0,08 - Khung mạng Si-O chuyển động chậm
O 0,12 - Liên kết BO bền, NBO linh động
Na 12,40 - Kênh dẫn khuếch tán (vùng giàu Na, mật độ O thấp)

Kết quả và những đóng góp mới

Luận án đã đạt được các kết quả và đóng góp mới sau:

  1. Đóng góp về mặt phương pháp:
  • Đề xuất và xây dựng thành công hai phương pháp phân tích hình học không gian mới: phương pháp Simplex (Void-simplex, Oxy-simplex, Cation-simplex) và phương pháp Shell-Core (SCP, SCC). Hai phương pháp này cho phép định lượng hóa các vùng không đồng nhất về mật độ, cấu trúc cục bộ và sự phân tách pha cation trong các vật liệu ôxít vô định hình và chất lỏng ba nguyên mà các hàm phân tích truyền thống ($g(r)$, BAD, số phối trí) không thể hiện được.
  1. Đóng góp về mặt lý luận và khoa học cấu trúc:
  • Làm sáng tỏ cấu trúc không đồng nhất của hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$ và $\text{PbO.SiO}_2$: Chứng minh sự tồn tại đồng thời của mạng con $\text{Si-O}$ kích thước lớn bao trùm và các mạng con $\text{Al-O}$ (hoặc cụm $\text{PbO}_x$) kích thước nhỏ hơn phân bố xen kẽ.
  • Xác lập quy luật chuyển đổi số phối trí và liên kết không gian dưới tác động của áp suất: chỉ ra sự chuyển dịch từ liên kết chung góc sang liên kết chung cạnh và chung mặt giữa các khối đa diện $\text{TO}_x$ khi nén áp suất cao (lên tới 60 GPa).
  • Làm sáng tỏ bản chất của hiện tượng Động học không đồng nhất (DH) trong hệ silicát ba nguyên:
    • Trong hệ $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng, các nguyên tử Na nhanh kết cụm dọc theo các nguyên tử NBO tạo thành các kênh dẫn ion; vùng nhanh là vùng có mật độ oxy thấp.
    • Trong hệ $\text{Al}_2\text{O}_3.2\text{SiO}_2$ lỏng, các nguyên tử O nhanh kết tụ tại các vùng có mật độ oxy cao; vùng chậm nằm tại các vùng giàu Si có mật độ oxy thấp.
  • Giải thích cơ chế khuếch tán dị thường trong nhôm-silicát lỏng thông qua vai trò trung gian của đơn vị phối trí 5 ($\text{SiO}_5$), sự hình thành/mất đi của các hạt siêu phân tử (SM) và hạt cứng (SLP) gắn liền với thời gian sống của các liên kết $\text{Si-O}$ và $\text{Al-O}$.
  1. Đóng góp thực tiễn và kiến nghị:
  • Cung cấp cơ sở dữ liệu số trị hoàn chỉnh về các thông số cấu trúc vi mô, hệ số khuếch tán và độ nhớt của các hệ silicát ba nguyên dưới điều kiện nhiệt độ cao và áp suất cao, hỗ trợ việc thiết kế, chế tạo và tối ưu hóa quy trình công nghệ sản xuất thủy tinh quang học, gốm sứ chịu nhiệt, xỉ luyện kim và vật liệu che chắn phóng xạ gốc chì-silicát.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Nghiên cứu thực hiện trên các mô hình mô phỏng Động lực học phân tử cổ điển sử dụng thế tương tác cặp Born-Mayer và thế BKS hiệu chỉnh; chưa khảo sát mô phỏng lượng tử ab-initio (do giới hạn về số lượng hạt và thời gian tính toán mô phỏng cho các hệ lớn tới gần 8.000 nguyên tử). Các khảo sát động học mới tập trung ở một số mốc nhiệt độ cố định (1.873 K, 3.200 K, 3.500 K).
  • Hướng nghiên cứu tiếp: Luận án định hướng mở rộng nghiên cứu sang các hệ silicát đa nguyên phức tạp hơn (bốn hoặc năm nguyên tố như $\text{CaO-Al}_2\text{O}_3\text{-SiO}_2$, $\text{Na}_2\text{O-CaO-SiO}_2$), khảo sát dải nhiệt độ chuyển pha thủy tinh ($T_g$) rộng hơn và kết hợp các phép đo phổ thực nghiệm tiên tiến để đối sánh sâu hơn các cấu trúc Simplex và Shell-Core.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu hữu ích cho:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý kỹ thuật, Vật lý chất rắn, Khoa học và Kỹ thuật vật liệu, Hóa lý silicat và Khoa học tính toán. Phần phương pháp luận (Chương 2) với các thuật toán Simplex, Shell-Core và kỹ thuật phân lập 10% hạt nhanh/chậm để tính toán DH là tài liệu mẫu đáng giá cho các nghiên cứu mô phỏng chất lỏng phức tạp.
  • Các kỹ sư công nghệ và nhà nghiên cứu trong lĩnh vực vật liệu thủy tinh, gốm sứ kỹ thuật cao, công nghệ luyện kim và an toàn bức xạ hạt nhân trong việc tra cứu số liệu về độ nhớt, hệ số khuếch tán và biến đổi cấu trúc dưới áp suất cao (Chương 3 và Chương 5).

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

1. Tại sao trong hệ $\text{Na}_2\text{O}.2\text{SiO}_2$ lỏng, hệ số khuếch tán của Na lại lớn hơn Si và O khoảng 100 lần?

Trong mạng $\text{Na}2\text{O}.2\text{SiO}2$, liên kết $\text{Si-O}$ là liên kết cộng hóa trị định hướng mạnh và bền vững tạo nên khung mạng cấu trúc chính. Khi $\text{Na}2\text{O}$ được đưa vào mạng, các ion $\text{Na}^+$ chỉ liên kết yếu với các nguyên tử oxy không cầu (NBO) và oxy cầu ($\text{BO-Na}, \text{Na-BO-Na}$) với khoảng cách liên kết lớn ($\sim 2,35 - 2,50\text{ \AA}$). Giữa các mạng con $\text{Si-O}$ độc lập xuất hiện các khe hở không gian tạo thành các kênh dẫn khuếch tán, cho phép các ion $\text{Na}^+$ di chuyển tự do với bước nhảy lớn mà không bị ràng buộc chặt chẽ, dẫn đến $D{\text{Na}} \approx 1,2 \times 10^{-5}\text{ cm}^2/\text{s}$, cao gấp khoảng 100 lần so với $D{\text{Si}}$ và $D{\text{O}}$ ($10^{-7}\text{ cm}^2/\text{s}$).

2. Sự khác biệt cơ bản về mật độ của vùng chuyển động nhanh (DH) giữa hệ nhôm-silicát và natri-silicát lỏng là gì?

Luận án chỉ ra sự đối lập về mật độ nguyên tử oxy tại các vùng DH:

  • Trong hệ nhôm-silicát ($\text{Al}_2\text{O}_3.2\text{SiO}_2$): Vùng chuyển động nhanh là vùng có mật độ cao các nguyên tử oxy (nơi tập trung nhiều oxy phối trí cao và các liên kết $\text{Al-O}$ linh động); vùng chuyển động chậm là vùng giàu $\text{Si-O}$ có mật độ oxy thấp hơn nhưng liên kết cứng vững.
  • Trong hệ natri-silicát ($\text{Na}_2\text{O}.2\text{SiO}_2$): Ngược lại, vùng chuyển động nhanh là vùng có mật độ nguyên tử oxy thấp (vùng kênh dẫn giàu $\text{Na}^+$ và NBO); vùng chuyển động chậm là vùng mạng lưới tứ diện $\text{Si-O}$ đặc khít có mật độ oxy cao.

3. Phương pháp phân tích Simplex và Shell-Core hoạt động như thế nào để phát hiện cấu trúc không đồng nhất?

  • Phương pháp Simplex: Dựng các mặt cầu đi qua 4 nguyên tử bất kỳ trong không gian mô phỏng. Bằng cách đếm và phân loại các nguyên tử nằm lọt bên trong mặt cầu, thuật toán phân chia thành Void-simplex (không chứa hạt nào), Oxy-simplex (chỉ chứa nguyên tử O) và Cation-simplex (chỉ chứa cation). Sự tập hợp các Simplex cùng loại chỉ ra các cụm cục bộ giàu oxy hoặc giàu cation.
  • Phương pháp Shell-Core: Xác định các hạt đa diện vi mô gồm phần lõi (Core) là 1 hoặc nhiều cation và phần vỏ (Shell) là các nguyên tử oxy bao quanh. Phân tích thống kê kích thước, bán kính lõi và sự liên kết giữa các hạt Shell-Core (SCC) cho phép định lượng chính xác sự tách pha vi mô giữa các mạng con $\text{Si-O}$ và $\text{Al-O}$.

4. Cơ chế vi mô nào tạo nên hiện tượng khuếch tán dị thường trong hệ nhôm-silicát lỏng dưới áp suất cao?

Hiện tượng khuếch tán dị thường (hệ số khuếch tán $D$ tăng lên khi tăng áp suất, đạt cực đại tại một áp suất trung gian rồi mới giảm khi nén tiếp) được luận án giải thích qua hai yếu tố vi cấu trúc:

  1. Ở áp suất thấp, mạng tứ diện $\text{TO}_4$ rất bền vững. Quá trình nén ban đầu phá vỡ một phần mạng tứ diện, làm tăng nồng độ các đơn vị phối trí trung gian 5 ($\text{SiO}_5, \text{AlO}_5$) có cấu trúc lỏng lẻo và độ linh động cao. Sự tăng và đạt cực đại của tỷ phần $\text{SiO}_5$ trùng khớp với điểm cực đại của hệ số khuếch tán.
  2. Sự đứt gãy và tái tạo liên tục của các liên kết $\text{Al-O}$ kém bền trong các hạt siêu phân tử (SM) tạo điều kiện cho các bước nhảy khuếch tán tập thể. Khi áp suất tăng quá cao ($20\text{ GPa}$), mạng phối trí 6 bát diện ($\text{TO}_6$) hoàn chỉnh được thiết lập làm mạng đặc khít trở lại, làm triệt tiêu sự khuếch tán.

5. Vai trò kép của nguyên tố Chì (Pb) trong cấu trúc mạng $\text{PbO.SiO}_2$ được thể hiện như thế nào qua số liệu của luận án?

Trong hệ $\text{PbO.SiO}_2$:

  • Khi ở nồng độ Pb thấp, các nguyên tử Pb đóng vai trò làm biến tính mạng, phá vỡ các liên kết $\text{Si-O-Si}$ và nằm tại các lỗ hổng mạng với số phối trí xấp xỉ 6.
  • Khi ở nồng độ Pb cao ($> 50\text{ mol}%$) hoặc dưới áp suất cao, các nguyên tử Pb trực tiếp liên kết với oxy tạo thành các đơn vị đa diện $\text{PbO}_3, \text{PbO}_4$ và liên kết với nhau bằng liên kết chung cạnh tạo thành các chuỗi polymer $(\text{PbO}_3)_n$ độc lập đan xen với mạng $\text{Si-O}$. Lúc này Pb đóng vai trò là nguyên tử tạo mạng song song cùng với Si.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Văn Yên đã nghiên cứu toàn diện cấu trúc vi mô và hiện tượng động học không đồng nhất trong các hệ silicát ba nguyên $\text{PbO.SiO}_2$, $\text{Al}_2\text{O}_3.2\text{SiO}_2$ và $\text{Na}_2\text{O}.2\text{SiO}_2$ ở trạng thái lỏng và vô định hình bằng phương pháp Động lực học phân tử. Công trình đã có đóng góp mới quan trọng khi xây dựng thành công hai phương pháp phân tích hình học Simplex và Shell-Core, đồng thời định lượng hóa sự kết cụm của nhóm 10% hạt nhanh nhất và chậm nhất để làm sáng tỏ bản chất của động học không đồng nhất. Các kết quả nghiên cứu cung cấp cơ sở dữ liệu số trị và lý thuyết tin cậy, giải thích sâu sắc mối tương quan giữa sự biến đổi trật tự cấu trúc mạng và các quá trình động học khuếch tán dị thường trong vật liệu silicát ba nguyên dưới điều kiện nhiệt độ và áp suất cao.