Tổng quan luận án

Đề tài "Nghiên cứu cấu trúc và một số tính chất của silica và sodium silicate" giải quyết các vấn đề cốt lõi về vi cấu trúc và động học của hệ vật liệu silica ($SiO_2$) và sodium silicate ($Na_2O.SiO_2$). Đây là các hệ vật liệu có ý nghĩa thực tiễn cao, xuất hiện phổ biến trong tự nhiên và được ứng dụng rộng rãi trong chế tạo linh kiện điện tử, thủy tinh, sợi quang, gốm sứ, xi măng và chất tẩy rửa.

Mặc dù các kỹ thuật thực nghiệm và mô phỏng trước đây đã xác định được các đơn vị cấu trúc cơ bản $SiO_x$ ($x = 4, 5, 6$) và $OSi_y$ ($y = 1, 2, 3$), nhiều khoảng trống khoa học quan trọng vẫn chưa được giải quyết đầy đủ:

  • Cấu trúc và sự chuyển đổi cấu trúc của silica lỏng theo nhiệt độ chưa được khảo sát chi tiết đến từng nguyên tử cụ thể.
  • Hiện tượng phân bố động học không đồng nhất và sự thuyên giảm động học trong silica lỏng chưa được làm rõ về mặt phân bố không gian - thời gian cũng như diễn biến kích thước các đám nguyên tử.
  • Cơ chế chuyển đổi từ cấu trúc tứ diện $SiO_4$ sang bát diện $SiO_6$ dẫn đến quá trình tinh thể hóa silica thủy tinh thành tinh thể stishovite dưới áp suất nén cao (lên đến 100 GPa) chưa có lời giải thích cụ thể về sự biến đổi độ dài và phân tách góc liên kết.
  • Trong hệ sodium silicate, vai trò biến đổi mạng của nguyên tử sodium, sự tồn tại và bản chất của kênh khuếch tán riêng (channel), diễn biến thay đổi cấu trúc theo thời gian và cơ chế khuếch tán tập thể kết hợp khuếch tán nhảy giữa các vị trí oxy cầu (BO) và oxy không cầu (NBO) cần được chứng minh chi tiết.

Mục tiêu nghiên cứu của luận án được xác định nhằm giải quyết ba nhóm nội dung cụ thể:

  1. Xác định cấu trúc và chuyển đổi cấu trúc theo nhiệt độ cụ thể tới từng nguyên tử trong mô hình silica lỏng; phân bố động học không đồng nhất trong silica lỏng theo không gian - thời gian, tương quan giữa cấu trúc và động học.
  2. Chuyển đổi cấu trúc silica thủy tinh khi bị nén dẫn đến tinh thể hóa.
  3. Cấu trúc và diễn biến thay đổi cấu trúc hệ sodium silicate, cơ chế khuếch tán và phân bố động học không đồng nhất trong mô hình.

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

  • Đối tượng nghiên cứu: Hệ silica lỏng và silica thủy tinh ($SiO_2$); hệ sodium silicate lỏng với hai thành phần $Na_2O.4SiO_2$ (ký hiệu NS4) và $Na_2O.3SiO_2$ (ký hiệu NS3).
  • Phạm vi điều kiện vật lý:
    • Silica lỏng: Khảo sát tại áp suất môi trường ($P = 0\text{ GPa}$) ở các mốc nhiệt độ 2600 K, 3000 K và 3500 K.
    • Silica thủy tinh: Khảo sát ở nhiệt độ $T = 500\text{ K}$ trong dải áp suất nén rộng từ 0 GPa đến 100 GPa.
    • Sodium silicate: Khảo sát các quá trình khuếch tán, chuyển đổi cấu trúc $SiO_4 \leftrightarrow SiO_3$ và $BO \leftrightarrow NBO$ theo thời gian mô phỏng (từ bước thời gian cỡ femtogiây đến các mốc thời gian $71,7\text{ ps}$ và $143,4\text{ ps}$).

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

Luận án tổng hợp toàn diện các kết quả nghiên cứu thực nghiệm và mô phỏng tính toán về silica và sodium silicate qua nhiều giai đoạn phát triển:

Các nghiên cứu thực nghiệm

  • Đặc trưng cấu trúc tinh thể và thủy tinh silica: Zachariasen (1932) đề xuất mô hình mạng ngẫu nhiên liên tục cho thủy tinh. Mozzi và Warren (1969) sử dụng nhiễu xạ tia X xác định khoảng cách $Si-O \approx 1,62\text{ \AA}$ và phân bố góc $Si-O-Si$ quanh $144^\circ$. Da Silva cùng cộng sự (1974) xác định đỉnh góc $Si-O-Si$ tại $152^\circ$. Kỹ thuật nhiễu xạ neutron của Wright (1988), Johnson (1983), Susman (1991) và cộng hưởng từ hạt nhân (NMR) của Dupree (1984), Pettifer (1988) tiếp tục khẳng định góc $Si-O-Si$ dao động trong khoảng $140^\circ - 155^\circ$.
  • Chuyển pha dưới áp suất cao: Nghiên cứu tán xạ tia X của Kingma (1993), Hemley (1988), Tsuchiya (1998) và nhiễu xạ Brillouin của Polian (1990) chứng minh sự chuyển đổi phối trí từ 4 sang 6, dẫn tới sự hình thành tinh thể stishovite ở áp suất cao (lên tới $92,1\text{ GPa}$).
  • Cấu trúc sodium silicate: Các công trình nhiễu xạ neutron của Misawa (1980), Wright (1988), Coombs (1995) và nhiễu xạ tia X của Waseda (1977), Imaoka (1983) cho thấy khi thêm $Na_2O$, số phối trí trung bình $O-Si$ giảm từ 2 xuống $1,5 - 1,7$, độ dài liên kết $Si-NBO$ ($1,587\text{ \AA}$) luôn ngắn hơn $Si-BO$ ($1,624 - 1,636\text{ \AA}$).
Dạng thù hình / Pha Khối lượng riêng $\rho$ ($\text{g/cm}^3$) Độ dài liên kết $Si-O$ ($\text{\AA}$) Góc liên kết $O-Si-O$ (độ) Góc liên kết $Si-O-Si$ (độ) Phương pháp / Tài liệu
Cristobalite $2,32 - 2,36$ $1,61$ $109,5$ $147 - 151$ Nhiễu xạ neutron [46, 47, 51]
Tridymite $2,37$ $1,57 - 1,62$ $109,5$ $140 - 173$ Nhiễu xạ tia X [41, 50, 51]
$\beta$-quartz $2,65$ $1,61$ $109,0$ $151$ Nhiễu xạ neutron [46, 47]
$\alpha$-quartz - $1,61$ $109,1 - 109,8$ $144 - 149$ Nhiễu xạ neutron [46]
Coesite $3,01$ - - - Thực nghiệm [51]
Stishovite $4,29 - 4,63$ $1,71$ - - Thực nghiệm [12, 51]
Silica thủy tinh (0 GPa) $2,20 - 2,50$ $1,59 - 1,62$ $109,0 - 109,8$ $144 - 152$ Nhiễu xạ tia X / neutron [8, 56, 58]

Các nghiên cứu mô phỏng động lực học phân tử

  • Bell và Dean (1972) xây dựng mô hình mạng silica thủy tinh liên kết qua oxy cầu.
  • Horbach và Kob (1999, 2008) mô phỏng silica lỏng bằng thế BKS và thế CHIK, chỉ ra sự thuyên giảm động học khi hạ nhiệt độ và hiện tượng khuếch tán dị thường ở vùng mật độ quanh $3,5\text{ g/cm}^3$ (áp suất $\approx 20\text{ GPa}$).
  • P. Hung cùng cộng sự (2006, 2007, 2013, 2018) khảo sát tương quan giữa phản ứng chuyển đổi $SiO_x \rightarrow SiO_{x\pm 1}$ và độ dịch chuyển nguyên tử, cũng như sự phân bố cấu trúc của hệ sodium silicate ở $1873\text{ K}$.
  • Badro (1997), Liping Huang (2006), Neng Li (2013) khảo sát silica nén ở áp suất cao, ghi nhận sự mở rộng phân bố góc $Si-O-Si$ thành hai đỉnh tách biệt ở $95^\circ$ và $130^\circ$.
  • Vogel và Glotzer (2004) phát hiện tính không đồng nhất động học thông qua các nhóm nguyên tử chuyển động nhanh và chuyển động chậm.

Vị trí và khoảng trống nghiên cứu luận án lựa chọn giải quyết

Luận án định vị trọng tâm vào việc vượt qua giới hạn của các phân tích cấu trúc trung bình bằng cách:

  1. Đưa ra phương pháp phân loại nút mạng chi tiết ($Z_{S_1S_2\dots S_Z}$) để phân biệt nút thường và nút khuyết tật tới từng nguyên tử trong silica lỏng.
  2. Thiết lập phương pháp phân tích động học không đồng nhất dựa trên ba tập hợp nguyên tử: nhanh nhất (SMA), chậm nhất (SIMA) và ngẫu nhiên (SRA), theo dõi sự tiến triển không gian - thời gian của các đám nguyên tử ($F_{link}(r,t)$) và xây dựng mô hình 2 miền.
  3. Giải thích cơ chế vi mô của hiện tượng dãn độ dài liên kết $Si-O$ và sự phân tách phân bố góc $Si-O-Si$ trong quá trình tinh thể hóa silica thủy tinh thành stishovite ở áp suất cao.
  4. Làm sáng tỏ động học khuếch tán của sodium thông qua phân tích các ô phối trí $NF_xB_y$, chứng minh cơ chế khuếch tán nhảy kết hợp khuếch tán tập thể đi kèm chuyển đổi $BO \leftrightarrow NBO$.

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

Thuật toán mô phỏng động lực học phân tử (MD)

Luận án ứng dụng phương pháp mô phỏng động lực học phân tử cổ điển. Hệ gồm $N$ nguyên tử được mô tả bằng phương trình chuyển động Newton: $$F_i = m_i a_i = m_i \frac{d^2 r_i}{dt^2} = -\sum_{j=1, j \neq i}^N \frac{\partial U_{ij}}{\partial r_{ij}}$$

Tọa độ và vận tốc tức thời của các nguyên tử được tích phân theo thời gian bằng thuật toán Verlet: $$r_i(t + \Delta t) = 2r_i(t) - r_i(t - \Delta t) + \frac{F_i(t)}{m_i} (\Delta t)^2$$ Bước thời gian mô phỏng được thiết lập $\Delta t \sim 10^{-15}\text{ s}$ (femtogiây), áp dụng điều kiện biên tuần hoàn (PBC) để loại bỏ hiệu ứng bề mặt của hệ kích thước hữu hạn.

Các hàm thế tương tác

  1. Thế tương tác BKS cho hệ silica ($SiO_2$): Do B. Van Santen đề xuất, gồm tương tác tĩnh điện Coulomb, tương tác đẩy Born-Mayer và tương tác hút van der Waals phân tán: $$U_{ij}(r_{ij}) = \frac{q_i q_j e^2}{r_{ij}} + A_{ij} \exp(-B_{ij} r_{ij}) - \frac{C_{ij}}{r_{ij}^6}$$
Cặp nguyên tử $A_{ij}$ ($\text{eV}$) $B_{ij}$ ($\text{\AA}^{-1}$) $C_{ij}$ ($\text{eV.\AA}^6$) Điện tích hiệu dụng ($e$)
$Si-O$ $18003,7572$ $4,87318$ $1133,5381$ $q_{Si} = +2,4$
$O-O$ $1388,7730$ $2,76000$ $175,0000$ $q_O = -1,2$
$Si-Si$ $0$ $0$ $0$ -
  1. Thế tương tác 2 và 3 thành phần cho hệ sodium silicate:
  • Thế hai thành phần: $$U_{ij}(r_{ij}) = \frac{q_i q_j e^2}{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})$$
  • Thế ba thành phần điều chỉnh góc liên kết $Si-O-Si$: $$U_{ijk}(\theta_{ijk}, r_{ij}, r_{ik}) = -f \left[\cos 2(\theta_{ijk} - \theta_0) - 1\right] k_{ij} k_{ik}$$ với $k_{ij} = \left[\exp\left(g_r (r_{ij} - r_m)\right) - 1\right]^{-1}$, $f_0 = 41,865\text{ kJ.nm}^{-1}\text{mol}^{-1}$, $f = 0,0006\text{ kJ.mol}^{-1}$, $\theta_0 = 147,0^\circ$, $r_m = 0,170\text{ nm}$, $g_r = 168,0\text{ nm}^{-1}$.

Quy trình xây dựng mô hình

  • Mô hình silica lỏng: Gồm 3000 nguyên tử (1000 nguyên tử $Si$ và 2000 nguyên tử $O$). Gieo ngẫu nhiên các nguyên tử trong ô mô phỏng ở áp suất môi trường, gia nhiệt lên nhiệt độ cao để phá vỡ trật tự ban đầu, làm nguội đẳng nhiệt - đẳng áp về các nhiệt độ $2600\text{ K}$, $3000\text{ K}$, $3500\text{ K}$, sau đó hồi phục cấu trúc trong hệ vi chính tắc đẳng tích - đẳng năng (NVE).
  • Mô hình silica thủy tinh: Gồm 4998 nguyên tử (1666 nguyên tử $Si$ và 3332 nguyên tử $O$). Tạo cấu hình lỏng, hạ nhiệt độ xuống $500\text{ K}$, nén đẳng nhiệt - đẳng áp (NPT) ở các mức áp suất từ $0\text{ GPa}$ đến $100\text{ GPa}$.
  • Mô hình sodium silicate: Hai hệ $Na_2O.4SiO_2$ (NS4) và $Na_2O.3SiO_2$ (NS3) được xây dựng tại Viện RIKEN (Nhật Bản).

Kỹ thuật phân tích vi cấu trúc và động học

  • Hàm phân bố xuyên tâm (PBXT $g(r)$) và bán kính ngắt ($r_c$): Vị trí cực tiểu thứ nhất của $g_{Si-O}(r)$ tại $r_c = 2,35\text{ \AA}$ được dùng làm bán kính ngắt để xác định liên kết $Si-O$ và số phối trí.
  • Phân loại nút mạng: Gán bộ chỉ số $Z_{S_1S_2\dots S_Z}$ cho từng nguyên tử ($Z$ là số phối trí nguyên tử trung tâm, $S_k$ là số phối trí của nguyên tử lân cận thứ $k$). Nút thường gồm $Si$ loại $4_{2222}$ và $O$ loại $2_{44}$; các dạng còn lại là nút khuyết tật (như $Si$ loại $4_{3222}$, $O$ loại $2_{54}$).
  • Cấu trúc Subnet: Xác định subnet thường, subnet khuyết tật (SDS: nhỏ, LDS: lớn, SNS: thường nhỏ) và các mạng con $SiO_4, SiO_5, SiO_6$.
  • Khảo sát động học: Tính độ dịch chuyển bình phương trung bình (MSD): $$\langle r^2(t) \rangle = \frac{1}{N} \sum_{i=1}^N \langle |r_i(t) - r_i(0)|^2 \rangle$$ Hệ số khuếch tán $D$ xác định theo phương trình Einstein: $$D = \lim_{t \to \infty} \frac{\langle r^2(t) \rangle}{6t}$$
  • Phân loại tập hợp động học: Tách các nguyên tử có độ dịch chuyển lớn nhất thành tập hợp nhanh nhất (SMA), nhỏ nhất thành tập hợp chậm nhất (SIMA) và tập hợp đối chứng ngẫu nhiên (SRA). Phân tích hàm phân bố đám $F_{link}(r, t)$ và $F_{inb}(r)$ để đánh giá mức độ kết đám trong không gian - thời gian.

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

Chương 1: Tổng quan cấu trúc và động học hệ silica và sodium silicate

Chương này tổng kết toàn bộ bức tranh lý thuyết và thực nghiệm về silica ở các pha lỏng, vô định hình/thủy tinh, tinh thể ($\alpha$-quartz, $\beta$-quartz, tridymite, cristobalite, coesite, stishovite) và hệ thủy tinh kiềm sodium silicate. Tác giả chỉ ra các quy luật biến đổi cấu trúc phụ thuộc vào áp suất và nhiệt độ: cấu trúc silica biến đổi rất ít theo nhiệt độ nhưng thay đổi đột biến dưới tác dụng của áp suất nén. Đối với sodium silicate, sự có mặt của ion $Na^+$ phá vỡ mạng lưới liên kết $Si-O-Si$, làm giảm số phối trí $O-Si$ từ $2,0$ xuống $1,6 - 1,7$, hình thành các oxy không cầu ($NBO$) và làm xuất hiện các kênh khuếch tán riêng biệt cho sodium. Chương này làm rõ các mâu thuẫn và giới hạn của các công trình trước đó, thiết lập cơ sở khoa học cho việc lựa chọn đề tài.

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

Trình bày chi tiết cơ sở toán - lý của phương pháp mô phỏng động lực học phân tử (MD), thuật toán Verlet và việc áp dụng điều kiện biên tuần hoàn. Thiết lập hệ thông số thế tương tác tĩnh điện và đẩy/hút BKS cho silica và hàm thế 2-3 thành phần cho sodium silicate. Định nghĩa tường minh các đại lượng vi cấu trúc: hàm PBXT $g(r)$, số phối trí $Z_{\alpha\beta}$, phân bố góc liên kết $\theta_{O-Si-O}$ và $\theta_{Si-O-Si}$, phân loại liên kết đa diện (liên kết góc, liên kết cạnh, liên kết mặt). Đưa ra phương pháp gán chỉ số nhận diện nút thường ($4_{2222}, 2_{44}$) và nút khuyết tật ($4_{3222}, 2_{54}, SiO_3, SiO_5, SiO_6$), thuật toán gom cụm xây dựng mạng con (subnet), cùng phương pháp phân tích động học không đồng nhất thông qua việc chia tách các nhóm SMA, SIMA, SRA và xác định các ô $NF_xB_y$.

Chương 3: Cấu trúc và động học hệ silica

Tập trung phân tích các kết quả mô phỏng trên hệ silica lỏng (2600 K, 3000 K, 3500 K) và quá trình nén silica thủy tinh (500 K, $0 - 100\text{ GPa}$):

  • Cấu trúc nút mạng và Subnet trong silica lỏng: Ở điều kiện áp suất môi trường, nút thường ($Si$ loại $4_{2222}$ và $O$ loại $2_{44}$) chiếm ưu thế tuyệt đối ở nhiệt độ thấp nhưng tỉ phần giảm dần khi nhiệt độ tăng. Khi nhiệt độ tăng từ 2600 K lên 3500 K, số lượng và kích thước của các subnet khuyết tật (LDS, SDS) tăng lên, phản ánh sự phá vỡ mạng ngẫu nhiên liên tục ở cấp độ nguyên tử. Tỉ phần nút khuyết tật tăng tuyến tính theo nhiệt độ.
  • Động học không đồng nhất và mô hình 2 miền: Phân tích các tập hợp SMA, SIMA và SRA tại các thời điểm $71,7\text{ ps}$ và $143,4\text{ ps}$ cho thấy các nguyên tử chuyển động nhanh nhất (SMA) có xu hướng kết tụ không gian thành các đám lớn vượt trội so với phân bố ngẫu nhiên SRA. Ngược lại, các nguyên tử chậm nhất (SIMA) tạo thành các miền cấu trúc bền vững. Tính không đồng nhất động học giảm khi nhiệt độ tăng và biến đổi chậm theo thời gian. Mô hình 2 miền không gian được thiết lập gồm miền linh động (chứa nhiều nút khuyết tật và các nguyên tử nhanh SMA) và miền trơ/đóng băng (chứa mạng nút thường và các nguyên tử chậm SIMA).
  • Chuyển đổi cấu trúc và tinh thể hóa silica thủy tinh dưới áp suất cao: Trong dải áp suất từ 0 đến 100 GPa ở 500 K:
    • Tỉ phần tứ diện $SiO_4$ giảm mạnh và chuyển dần sang $SiO_5$, đạt trạng thái chiếm ưu thế của bát diện $SiO_6$ tại áp suất trên $40\text{ GPa}$.
    • Tỉ phần liên kết cạnh và liên kết mặt tăng mạnh ở áp suất cao, thay thế dần liên kết góc đơn thuần.
    • Phân bố góc $Si-O-Si$ từ một đỉnh duy nhất tại $148^\circ$ ở 0 GPa bị dẹt xuống và tách rõ rệt thành hai đỉnh đặc trưng tại $95^\circ$ và $130^\circ$ khi áp suất đạt $60 - 100\text{ GPa}$.
    • Đồ thị hàm PBXT và phân bố không gian ở $100\text{ GPa}$ chứng minh cấu trúc silica thủy tinh vô định hình đã chuyển pha trật tự hóa thành tinh thể stishovite. Luận án giải thích hiện tượng độ dài liên kết $Si-O$ tăng từ $1,62\text{ \AA}$ lên $1,72 - 1,75\text{ \AA}$ trong khoảng $15 - 30\text{ GPa}$ là do sự hình thành phối trí 6 ($SiO_6$), sau đó giảm nhẹ trở lại ở áp suất cực cao do lực nén cơ học trực tiếp lên mạng tinh thể.
Áp suất nén ($P$) Độ dài liên kết $Si-O$ ($\text{\AA}$) Phân bố góc $Si-O-Si$ (độ) Phân bố góc $O-Si-O$ (độ) Hiện tượng vi cấu trúc đặc trưng
0 GPa $1,60 - 1,65$ Đỉnh đơn $\sim 147 - 153$ $109,0 - 109,5$ Cấu trúc tứ diện $SiO_4$, liên kết góc
10 GPa $1,63 - 1,70$ $137,4 - 140,0$ $107,0 - 108,0$ Xuất hiện $SiO_5$, bắt đầu nén góc
15 GPa $1,72$ $133,9$ $106,9$ Bắt đầu chuyển pha phối trí
20 GPa - $135,0$ $109,0$ Dịch chuyển đỉnh PBXT
30 GPa $1,72 - 1,75$ $120,0$ $103,0$ Tỉ phần $SiO_6$ tăng nhanh, xuất hiện liên kết cạnh
40 - 50 GPa $1,69$ Tách 2 đỉnh: $100$ và $130$ $90,0$ Bát diện $SiO_6$ chiếm ưu thế
60 - 100 GPa $1,66 - 1,69$ Tách 2 đỉnh rõ: $95$ và $130$ - Tinh thể hóa hoàn chỉnh pha stishovite

Chương 4: Cấu trúc và động học hệ sodium silicate

Khảo sát chuyên sâu cấu trúc mạng và động học khuếch tán trong hai hệ $Na_2O.4SiO_2$ và $Na_2O.3SiO_2$:

  • Phân bố vi cấu trúc: Bổ sung $Na_2O$ dẫn đến sự phân cắt mạng lưới $Si-O$, làm xuất hiện các nguyên tử oxy không cầu ($NBO$) và oxy tự do ($FO$), trong đó $BO$ giảm và $NBO$ tăng mạnh. Cấu trúc tứ diện $Q^4$ giảm sâu, thay thế bằng các đơn vị $Q^3, Q^2, Q^1$. Các nguyên tử $Na$ có xu hướng tập trung bao quanh các nguyên tử $NBO$, tạo ra các vùng giàu sodium đan xen với các vùng giàu silicon.
  • Kênh khuếch tán và phân tích ô $NF_xB_y$: Sự phân bố không gian của sodium tạo thành các kênh khuếch tán (channel) riêng biệt. Phân tích các ô $NF_xB_y$ (trong đó $x$ là số $FO/NBO$, $y$ là số $BO$) làm sáng tỏ hình ảnh các kênh này. Mật độ sodium trong channel cao hơn rõ rệt so với ngoài channel ở nhiệt độ thấp.
  • Cơ chế khuếch tán của Sodium: Luận án chứng minh sodium khuếch tán thông qua hai cơ chế đồng thời:
    1. Cơ chế nhảy (jump diffusion): Sodium nhảy giữa các vị trí lân cận $NBO$ dọc theo các channel rỗng.
    2. Cơ chế khuếch tán tập thể (collective diffusion): Xảy ra đồng thời với quá trình chuyển đổi phối trí liên tục giữa $BO \leftrightarrow NBO$ và $SiO_4 \leftrightarrow SiO_3$ theo thời gian. Sự di chuyển của sodium đi kèm với sự tái cấu trúc cục bộ của mạng lưới silicat, giải thích tính linh động vượt trội của sodium so với khung mạng $Si-O$.

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

Đóng góp về mặt khoa học

  1. Xác định cấu trúc silica lỏng ở cấp độ từng nút nguyên tử: Luận án là công trình đầu tiên thiết lập phương pháp phân loại và định lượng các loại nút mạng silicon và oxy ($Z_{S_1S_2\dots S_Z}$), phân biệt tường minh các nút thường ($4_{2222}, 2_{44}$) và các loại nút khuyết tật ($4_{3222}, 2_{54}$, v.v.), cũng như sự tiến triển kích thước của các mạng con (subnet thường và subnet khuyết tật) theo nhiệt độ (2600 K, 3000 K, 3500 K).
  2. Làm rõ bản chất không gian - thời gian của động học không đồng nhất: Thông qua việc phân tích động học ba tập hợp SMA, SIMA và SRA kết hợp hàm phân bố đám $F_{link}(r,t)$, luận án chứng minh tính không đồng nhất động học trong silica lỏng giảm dần theo nhiệt độ và biến đổi chậm theo thời gian; đồng thời đề xuất mô hình 2 miền không gian giải thích mối tương quan giữa vi cấu trúc khuyết tật và tính linh động nguyên tử.
  3. Lý giải cơ chế chuyển đổi cấu trúc và tinh thể hóa silica thủy tinh dưới áp suất cao: Chứng minh quá trình nén từ 0 đến 100 GPa ở 500 K dẫn đến sự chuyển đổi từ tứ diện $SiO_4$ sang bát diện $SiO_6$, làm xuất hiện liên kết cạnh và làm phân tách phân bố góc $Si-O-Si$ thành hai đỉnh ($95^\circ$ và $130^\circ$). Xác định sự hình thành pha tinh thể stishovite ở áp suất cao và giải thích bản chất hiện tượng biến thiên độ dài liên kết $Si-O$.
  4. Xác lập cơ chế khuếch tán kép của sodium trong sodium silicate: Khẳng định các quá trình chuyển đổi $SiO_4 \leftrightarrow SiO_3$ và $BO \leftrightarrow NBO$ luôn diễn ra động học theo thời gian; chứng minh sodium khuếch tán kết hợp giữa cơ chế nhảy cục bộ quanh $NBO$ trong các ô $NF_xB_y$ và cơ chế khuếch tán tập thể phối hợp với sự chuyển đổi liên kết của mạng nền.

Đóng góp về mặt thực tiễn

  • Các dữ liệu định lượng về vi cấu trúc, độ bền nhiệt động, sự biến dạng góc liên kết và hệ số khuếch tán cung cấp cơ sở dữ liệu vật lý quan trọng cho việc mô phỏng, tính toán và tối ưu hóa quy trình công nghệ chế tạo vật liệu thủy tinh chịu lực, gốm sứ kỹ thuật cao, sợi quang học và vật liệu silicat kiềm.

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

  • Hạn chế: Các kết quả của luận án thu được từ phương pháp mô phỏng động lực học phân tử cổ điển với kích thước hệ từ 3000 đến 4998 nguyên tử và thời gian mô phỏng trong thang picogiây đến nanogiây. Mô hình thế BKS và thế 2-3 thành phần tuy đã được tối ưu hóa nhưng vẫn là thế thực nghiệm bán kinh nghiệm, chưa mô tả trực tiếp sự tái phân bố mật độ điện tử cơ học lượng tử trong các trạng thái chuyển tiếp phức tạp.
  • Hướng nghiên cứu tiếp:
    • Mở rộng khảo sát với các mô hình kích thước lớn hơn và thang thời gian dài hơn để khảo sát sâu hơn động thái của các channel ở nhiệt độ trung gian.
    • Áp dụng mô phỏng động lực học phân tử lượng tử (ab initio MD) tại các vùng chuyển pha áp suất cao để kiểm chứng vi cấu trúc điện tử.
    • Nghiên cứu thêm các hệ silicat kiềm hỗn hợp (như $Na_2O-K_2O-SiO_2$, $Li_2O-SiO_2$) dưới các điều kiện tải trọng và nhiệt độ khắc nghiệt.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị cao đối với:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý chất rắn, Vật lý tính toán, Khoa học vật liệu và Hóa học vô cơ: Cung cấp phương pháp luận chuẩn xác về thuật toán mô phỏng MD, kỹ thuật phân tích hàm phân bố xuyên tâm, phân loại nút mạng topo, và phương pháp xử lý số liệu phân bố đám động học không đồng nhất.
  • Các kỹ sư công nghệ vật liệu silicat, thủy tinh và gốm sứ: Cung cấp các thông số định lượng chính xác về vi cấu trúc, sự chuyển đổi phối trí, cơ chế di chuyển của ion kiềm và ngưỡng áp suất chuyển pha vật liệu.

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

1. Luận án sử dụng tiêu chuẩn nào để phân biệt nút thường và nút khuyết tật trong silica lỏng?

Luận án sử dụng bộ chỉ số đặc trưng $Z_{S_1S_2\dots S_Z}$, trong đó $Z$ là số phối trí của nguyên tử trung tâm và $S_k$ là số phối trí của từng nguyên tử lân cận (với bán kính ngắt xác định liên kết $Si-O$ là $r_c = 2,35\text{ \AA}$). Nút thường được định nghĩa duy nhất là các nút $Si$ loại $4_{2222}$ ($Si$ liên kết với 4 nguyên tử $O$, mỗi $O$ đều liên kết với 2 $Si$) và nút $O$ loại $2_{44}$ ($O$ liên kết với 2 nguyên tử $Si$, mỗi $Si$ đều liên kết với 4 $O$). Tất cả các cấu hình phối trí khác (như $Si$ loại $4_{3222}$, $SiO_3, SiO_5, SiO_6$, $O$ loại $2_{54}, OSi_1, OSi_3$) đều được xác định là nút khuyết tật.

2. Quá trình nén silica thủy tinh ở 500 K đến 100 GPa làm thay đổi phân bố góc liên kết $Si-O-Si$ như thế nào?

Ở áp suất 0 GPa, phân bố góc $Si-O-Si$ có một đỉnh phân bố đơn tập trung quanh $147^\circ - 153^\circ$. Khi nén tăng áp suất, góc liên kết bị ép hẹp lại và dẹt dần. Khi áp suất đạt mức cao ($60 - 100\text{ GPa}$), phân bố góc $Si-O-Si$ tách rõ rệt thành hai đỉnh riêng biệt tại khoảng $95^\circ$ và $130^\circ$, phản ánh sự hình thành mạng lưới cấu trúc bát diện $SiO_6$ đặc trưng cho pha tinh thể stishovite.

3. Tại sao độ dài liên kết $Si-O$ trong silica thủy tinh lại tăng lên khi bị nén trong dải áp suất 15 - 30 GPa?

Trong khoảng áp suất $15 - 30\text{ GPa}$, hệ xảy ra sự chuyển đổi phối trí mạnh mẽ từ tứ diện $SiO_4$ sang bát diện $SiO_6$ thông qua cấu trúc trung gian $SiO_5$. Do độ dài liên kết trung bình $Si-O$ trong khối đa diện phối trí 6 ($SiO_6$) về bản chất lớn hơn ($\approx 1,72 - 1,75\text{ \AA}$) so với trong khối tứ diện phối trí 4 ($SiO_4$, $\approx 1,60 - 1,62\text{ \AA}$), sự chiếm ưu thế của cấu trúc $SiO_6$ đã làm độ dài liên kết trung bình $Si-O$ toàn hệ tăng lên trước khi bị nén ngắn lại dưới lực ép ở các áp suất cao hơn ($> 40\text{ GPa}$).

4. Sodium trong hệ sodium silicate khuếch tán theo những cơ chế cụ thể nào?

Sodium khuếch tán theo hai cơ chế song hành:

  1. Cơ chế nhảy (jump diffusion): Các ion $Na^+$ dịch chuyển nhảy bước giữa các vị trí lân cận nguyên tử oxy không cầu ($NBO$) trong không gian các kênh (channel) khuếch tán.
  2. Cơ chế khuếch tán tập thể (collective diffusion): Xảy ra đồng thời với quá trình phá vỡ và tái tạo liên kết giữa mạng lưới $Si-O$ ($BO \leftrightarrow NBO$ và $SiO_4 \leftrightarrow SiO_3$), tạo điều kiện cho các ion sodium di chuyển phối hợp theo từng nhóm nguyên tử.

Kết luận

Luận án "Nghiên cứu cấu trúc và một số tính chất của silica và sodium silicate" đã giải quyết thành công các vấn đề vật lý then chốt về vi cấu trúc và động học của silica và sodium silicate bằng phương pháp mô phỏng động lực học phân tử. Công trình đã làm sáng tỏ cấu trúc mạng nút thường - nút khuyết tật của silica lỏng, bản chất không gian - thời gian của hiện tượng động học không đồng nhất và cơ chế tinh thể hóa silica thủy tinh thành stishovite dưới áp suất cực cao lên tới 100 GPa. Đồng thời, luận án đã chứng minh quy luật biến đổi cấu trúc theo thời gian và cơ chế khuếch tán kép của ion sodium trong hệ sodium silicate, đóng góp những dữ liệu lý thuyết và thực nghiệm tính toán có độ tin cậy cao cho ngành vật lý chất rắn và công nghệ vật liệu.