Tổng quan về luận án

Vật liệu silica ($SiO_2$) và các hệ thủy tinh, silicat nóng chảy đa nguyên ($Na_2O-SiO_2$, $Al_2O_3-SiO_2$, $PbO-SiO_2$) giữ vai trò trung tâm trong khoa học vật liệu hiện đại, chi phối từ công nghệ vi điện tử, siêu tụ điện, pin natri trạng thái rắn (SSNBs) cho đến công nghệ gốm sứ chịu nhiệt và vật liệu chắn bức xạ hạt nhân. Tuy nhiên, việc khảo sát thực nghiệm cấu trúc nguyên tử và động học vi mô của silicat ở điều kiện nhiệt độ cực cao (lên đến 3500 K) hoặc áp suất lớn (đạt 30 GPa) gặp phải những rào cản kỹ thuật nghiêm trọng. Luận án tiến sĩ khoa học vật liệu của nghiên cứu sinh Mai Văn Dũng với tiêu đề "Mô phỏng cấu trúc và động học của vật liệu silicát - Na/Al/Pb" (chuyên ngành Vật liệu điện tử, mã số: 9440123, thực hiện dưới sự hướng dẫn của PGS. TS. Lê Thế Vinh và PGS. TS. Nguyễn Mạnh Tuấn tại Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam, 2022) đại diện cho một nghiên cứu tiên phong áp dụng phương pháp Động lực học phân tử (Molecular Dynamics - MD) cổ điển kết hợp các công cụ topo học tiên tiến nhằm giải mã bản chất trật tự gần (SRO), trật tự tầm trung (MRO) và hiện tượng động học không đồng nhất (Dynamical Heterogeneity - DH).

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án hướng tới giải quyết bắt nguồn từ những tranh luận học thuật kéo dài về hiện tượng biến đổi đa hình vô định hình (polyamorphism) và sự chuyển pha cấu trúc dưới áp suất cao của $SiO_2$, cơ chế hình thành các cụm $SiO_x$ ($x=4, 5, 6$), vai trò lưỡng tính (network former vs. network modifier) của cation kim loại nặng $Pb^{2+}$, và đặc biệt là sự thiếu vắng một bức tranh định lượng hoàn chỉnh về cơ chế khuếch tán không đồng nhất của ion $Na^+$ trong môi trường mạng silicat bị phân tách pha vi mô.

Nghiên cứu được xây dựng trên hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) chặt chẽ:

  • RQ1: Sự biến đổi mật độ từ 2,499 g/cm³ đến 3,994 g/cm³ ảnh hưởng như thế nào đến tỷ lệ phối trí $SiO_4 \rightarrow SiO_5 \rightarrow SiO_6$ và hiện tượng tách đỉnh trong hàm phân bố xuyên tâm (HPBXT) của cặp Si-Si trong silica vô định hình?
  • H1: Quá trình nén ép silica dẫn đến sự tái sắp xếp topo từ liên kết chung đỉnh sang liên kết chung cạnh và chung mặt giữa các đa diện $SiO_x$, tạo nên các cụm cấu trúc không đồng nhất.
  • RQ2: Dưới tác động đồng thời của áp suất (0–30 GPa) và nhiệt độ (2100–3500 K), cấu trúc topo của hệ nhôm-silicat $Al_2O_3 \cdot 2SiO_2$ (AS2) biến đổi ra sao thông qua các miền thể tích tự do và cụm cation?
  • H2: Sự nén ép của pha lỏng AS2 không diễn ra đồng nhất mà bắt nguồn từ sự co rút thể tích của các Void-simplex cùng sự hình thành mạng lưới oxy ba phối trí (O-triclusters).
  • RQ3: Nồng độ $PbO$ (từ $x = 0,05$ đến $0,75$) chi phối sự chuyển đổi chức năng cấu trúc của chì trong thủy tinh $xPbO \cdot (1-x)SiO_2$ như thế nào?
  • H3: Cation $Pb$ chuyển hóa từ vai trò phá vỡ mạng (modifier) ở nồng độ thấp sang vai trò tạo mạng (former) ở nồng độ $x > 0,50$ thông qua sự hình thành các kênh ống liên kết đa diện $PbO_y$.
  • RQ4: Động học khuếch tán của ion $Na^+$ trong $Na_2O \cdot 2SiO_2$ (NS2) và $Na_2O \cdot 4SiO_2$ (NS4) diễn ra theo cơ chế đơn lẻ hay tập thể, và tương quan thế nào với các đa diện Voronoi của nguyên tử oxy?
  • H4: Khuếch tán của $Na^+$ bị chi phối bởi hiện tượng động học không đồng nhất, trong đó ion $Na$ tập trung áp đảo trong các đa diện chứa oxy không cầu nối (NBO) và oxy tự do (FO), tạo ra sự phân tách thành các kênh dẫn truyền ion tốc độ cao.

Khung lý thuyết của luận án tích hợp Lý thuyết Mạng ngẫu nhiên liên tục (Continuous Random Network - CRN) của Zachariasen (1932), Lý thuyết Mạng ngẫu nhiên biến tính (Modified Random Network - MRN) của Greaves (1985), kết hợp các mô hình thế tương tác BKS (van Beest, Kramer, van Santen, 1990) và Born–Mayer. Tác động đột phá của công trình được định lượng qua việc xác định chính xác tỷ lệ phân bố của hơn 3000 nguyên tử trong các hộp mô phỏng tuần hoàn, theo dõi động học trong thang thời gian 10–150 ps, làm sáng tỏ các thông số cấu trúc vi mô với sai số dưới 3% so với dữ liệu thực nghiệm nhiễu xạ tia-X (XRD) và phổ cộng hưởng từ hạt nhân (NMR).

Literature Review và Positioning

Nghiên cứu cấu trúc silica và silicat khởi đầu từ công trình kinh điển của Zachariasen (1932) và thực nghiệm nhiễu xạ tia-X của Mozzi & Warren (1969) [20], xác định góc liên kết Si-O-Si trung bình trong silica thủy tinh đạt $144^\circ$ cùng độ dài liên kết $Si-O = 1,62\text{ Å}$. Tiếp đó, Bell & Dean (1972) [57] thiết lập mô hình quả cầu - thanh nối đưa ra góc liên kết $153^\circ$, phù hợp với các quan sát tán xạ neutron hiện đại. Tuy nhiên, các kỹ thuật thực nghiệm truyền thống bộc lộ nhiều hạn chế khi tiếp cận trạng thái vô định hình chịu nén. Nghiên cứu thực nghiệm XRD áp suất siêu cao của Prescher et al. (2017) [26] trên $SiO_2$ đến 172 GPa chỉ ra số phối trí của Si tăng từ 4 lên 6 trong khoảng 14–40 GPa và đạt 6,8 ở 172 GPa, nhưng không thể phân tách chi tiết sự phân bố không gian cục bộ của từng loại đa diện $SiO_4, SiO_5, SiO_6$.

Trong hệ nhôm-silicat, Morikawa et al. (1982) [92] và Wilding et al. (2010) [94] qua phổ NMR $^{27}Al$ đã chứng minh sự tồn tại đồng thời của $AlO_4, AlO_5, AlO_6$, nhưng cấu trúc trật tự tầm trung và sự hình thành các liên kết ba oxy (O-triclusters) vẫn tồn tại nhiều mâu thuẫn giữa các kết quả mô phỏng ab initio của Kubicki et al. (1993) [96] và Benoit et al. (2001) [104]. Đối với hệ chì-silicat, nghiên cứu của Bair (1936) [126] và Warren & Loring (1935) [127] khẳng định cấu trúc mạng định hướng ngẫu nhiên, trong khi Mizuno et al. (2005) [135] qua nhiễu xạ tia-X phát hiện cấu trúc dạng đường ống độc lập khi $PbO < 35\text{ mol}%$, chuyển thành mạng không gian 3 chiều gồm các đường ống lồng nhau khi $PbO$ đạt 35–50 mol%.

TIẾN TRÌNH NGHIÊN CỨU CẤU TRÚC SILICAT TRONG LỊCH SỬ
                             Đa diện Voronoi Oxy giải mã Động học không đồng nhất
Tiêu chí so sánh Nghiên cứu của Prescher et al. (2017) [26] & Mizuno et al. (2005) [135] Luận án Mai Văn Dũng (2022)
Đối tượng & Hệ mẫu $SiO_2$ thủy tinh nén đơn thuần; $xPbO-(100-x)SiO_2$ ở nhiệt độ phòng Khảo sát đồng thời $SiO_2$ (2,499–3,994 g/cm³), $Al_2O_3 \cdot 2SiO_2$ (0–30 GPa, 2100–3500 K), $xPbO \cdot (1-x)SiO_2$ và $Na_2O \cdot (2,4)SiO_2$
Phương pháp tiếp cận Thực nghiệm nhiễu xạ XRD, EXAFS, phổ NMR gián đoạn Mô phỏng ĐLHPT quy mô lớn, tính tổng Ewald, phân tích vi cấu trúc theo thời gian thực
Độ phân giải topo Chỉ đo được số phối trí trung bình và khoảng cách cặp tổng thể Phân tách chi tiết cụm $SiO_x$, phân tích Void-simplex, Shell-core, và đa diện Voronoi oxy
Phân tích động học Giới hạn ở trạng thái tĩnh, không theo dõi được quỹ đạo ion Làm sáng tỏ cơ chế khuếch tán $Na^+$ theo bước nhảy đơn lẻ và dịch chuyển tập thể dạng chuỗi

Luận án định vị chính xác vị thế học thuật bằng việc khắc phục triệt để điểm mù của các phương pháp thực nghiệm: lượng hóa cấu trúc vi mô không đồng nhất và xác lập mối liên hệ nhân quả trực tiếp giữa cấu trúc topo cục bộ với các quá trình động học ion ở điều kiện khắc nghiệt.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng sâu sắc Lý thuyết Mạng ngẫu nhiên liên tục (CRN) của Zachariasen và Lý thuyết Mạng ngẫu nhiên biến tính (MRN) của Greaves thông qua các phát hiện mang tính chuyển đổi nhận thức (paradigm shift):

  1. Mô hình hóa sự chuyển pha đa hình vô định hình: Chứng minh sự chuyển dịch từ cấu trúc tứ diện $SiO_4$ sang bát diện $SiO_6$ dưới áp suất nén không xảy ra đột ngột mà thông qua trạng thái trung gian $SiO_5$. Sự nén ép làm xuất hiện liên kết chung cạnh và chung mặt giữa các đa diện, bẻ gãy giả định truyền thống của Zachariasen rằng các tứ diện chỉ chia sẻ đỉnh oxy.
  2. Khung lý thuyết lưỡng tính cation đa hóa trị: Mở rộng lý thuyết oxide thủy tinh đối với $Pb^{2+}$, xác định điểm chuyển tiếp topo (topological crossover) tại nồng độ $x = 0,50$, nơi chì chính thức chuyển hóa từ vai trò phá vỡ mạng sang vai trò trung tâm tạo mạng thông qua các phức đa diện $PbO_3$ và $PbO_4$.
  3. Lý thuyết phân tách pha vi mô và động học phi liên tục: Phát triển mô hình động học không đồng nhất dựa trên sự phân chia không gian thành "miền mạng con $Si-O$" tĩnh tại và "miền cụm $O-O$" linh động. Luận án bác bỏ quan điểm khuếch tán đồng nhất đẳng hướng, xác lập nguyên lý tập trung ion dẫn đường: ion $Na^+$ bị giam giữ ưu tiên trong các đa diện oxy không cầu nối ($NBO$) và oxy tự do ($FO$).
KHUNG LÝ THUYẾT ĐỘNG HỌC KHÔNG ĐỒNG NHẤT TRONG HỆ NATRI-SILICAT

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp đồng thời ba cấu trúc hình học tính toán cao cấp:

  • Phương pháp Simplex: Chia không gian thành các đơn vị tứ diện cơ bản cấu thành từ 4 hạt lân cận, phân loại thành Void-simplex (VS) (miền thể tích tự do rỗng), Oxy-simplex (OS) (miền chỉ chứa nguyên tử O), và T-simplex/Cation-simplex (CS) (miền chỉ chứa các cation tạo mạng như Si, Al). Phương pháp này cho phép đo kích thước thể tích trống cục bộ mà không cần phụ thuộc vào bán kính cắt ($r_{cut}$).
  • Phương pháp Shell-Core: Xác định các vi hạt cấu trúc gồm lõi (core) đặc và vỏ (shell) liên kết, từ đó khảo sát sự kết tụ của các đơn vị thành cụm SC-cluster, định lượng mức độ phân tách pha vi mô ở nhiệt độ 2100 K và 3500 K.
  • Phương pháp Đa diện Voronoi Oxy: Xây dựng các đa diện Voronoi bao quanh từng nguyên tử oxy trong không gian 3 chiều, theo dõi số lượng cation đi vào ($\langle n_{VYP} \rangle$) hoặc hiện diện ($\langle n_{AYP} \rangle$) trong đa diện qua từng bước thời gian $\Delta t$, cho phép tách biệt chính xác hai trạng thái oxy cầu nối (BO) và oxy không cầu nối (NBO).

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên lập trường triết học thực chứng (positivism) thông qua mô phỏng nguyên tử số trị tất định (deterministic atomistic simulation). Thiết kế mô phỏng đa mức độ (multi-level) kết nối chặt chẽ từ tương tác cơ bản giữa các cặp ion đến các biến đổi nhiệt động học vĩ mô.

Quy mô mẫu và thông số mô phỏng được chuẩn hóa chính xác:

  • Hệ $SiO_2$ vô định hình: Hộp mô phỏng gồm 1998 nguyên tử (666 Si, 1332 O) được khảo sát ở ba mức khối lượng riêng chuẩn: 2,499 g/cm³, 3,524 g/cm³ và 3,994 g/cm³.
  • Hệ nhôm-silicat ($Al_2O_3 \cdot 2SiO_2$ - AS2): Gồm 3000 nguyên tử thiết lập ở các mức áp suất $P = 0, 2, 5, 10, 15, 20, 30\text{ GPa}$ và khoảng nhiệt độ $T = 2100\text{ K} \div 3500\text{ K}$.
  • Hệ chì-silicat ($xPbO \cdot (1-x)SiO_2$): Hệ 3000 nguyên tử tại $T = 3000\text{ K}$, khảo sát 8 dải nồng độ $x = 0,05; 0,10; 0,15; 0,20; 0,33; 0,50; 0,67; 0,75$.
  • Hệ natri-silicat ($Na_2O \cdot 2SiO_2$ và $Na_2O \cdot 4SiO_2$): Hệ gồm 3000 nguyên tử lỏng tại $T = 1500\text{ K} \div 3000\text{ K}$, ghi nhận dữ liệu liên tục trong thang thời gian 10–150 ps.
QUY TRÌNH TÍNH TOÁN VÀ PHÂN TÍCH ĐỘNG LỰC HỌC PHÂN TỬ

Quy trình nghiên cứu và cơ chế tính toán

  1. Thuật toán tích phân Verlet: Giải phương trình vi phân chuyển động Newton cho từng nguyên tử: $$r_i(t + \Delta t) = r_i(t) + v_i(t)\Delta t + \frac{1}{2}a_i(t)\Delta t^2$$ $$v_i(t) = \frac{1}{2\Delta t}[r_i(t + \Delta t) - r_i(t - \Delta t)]$$
  2. Phương pháp tính tổng Ewald cho thế tầm xa Coulomb: Triệt tiêu sự phân kỳ chậm của tương tác tĩnh điện trong điều kiện biên tuần hoàn bằng cách phân rã thế năng Coulomb thành tổng trong không gian thực ($U_t$), không gian đảo ($U_d$) và số hạng tự tương tác ($U_h$): $$U_t = \frac{1}{2}\sum_{\mathbf{n}}!!^{'}\sum_{i=1}^N \sum_{j=1}^N \frac{q_i q_j}{|r_i - r_j| + \mathbf{n}},\text{erfc}\big(\alpha(|r_i - r_j| + \mathbf{n})\big)$$ $$U_d = \frac{1}{2\pi V_0}\sum_{i=1}^N \sum_{j=1}^N q_i q_j \sum_{\mathbf{m} \neq 0} \frac{\exp\big(-\frac{\pi^2 \mathbf{m}^2}{\alpha^2} + 2\pi i \mathbf{m}\cdot(r_i - r_j)\big)}{\mathbf{m}^2}$$
  3. Thế tương tác BKS và Born–Mayer: Đối với $SiO_2$, thế BKS có dạng: $$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}$$ với $q_{Si} = +2,4e$, $q_O = -1,2e$, $A_{Si-O} = 18003,7572\text{ eV}$, $B_{Si-O} = 4,87318\text{ \AA}^{-1}$, $C_{Si-O} = 1133,5381\text{ eV}\cdot\text{\AA}^6$. Đối với hệ $Na_2O-SiO_2$, thế góc ba phần tử Si-O-Si được bổ sung để đảm bảo tính định hướng liên kết cộng hóa trị: $$U_{ijk} = -f\big[\cos{2(\theta_{kij} - \theta_0)} - 1\big]\sqrt{k_{ij} k_{ik}}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án mang lại 4 phát hiện mang tính bước ngoặt, được hỗ trợ bởi các bằng chứng dữ liệu chi tiết trích xuất trực tiếp từ mô hình:

BIẾN ĐỔI PHỐI TRÍ VÀ CẤU TRÚC TOPO THEO MẬT ĐỘ TRONG SILICA

  1. Hiện tượng tách đỉnh hàm phân bố xuyên tâm Si-Si trong silica nén cao: Khi khối lượng riêng của silica vô định hình tăng từ 2,499 g/cm³ lên 3,994 g/cm³, hàm phân bố xuyên tâm $g_{Si-Si}(r)$ xuất hiện sự tách đỉnh thứ nhất rõ rệt thành hai đỉnh phụ tại $r_1 \approx 2,75\text{ \AA}$ và $r_2 \approx 3,12\text{ \AA}$. Luận án chứng minh nguyên nhân trực tiếp là sự hình thành các liên kết đa diện chung cạnh (edge-sharing, tương ứng với đỉnh $r_1$) và chung mặt (face-sharing) bên cạnh liên kết chung đỉnh truyền thống (corner-sharing, tương ứng với $r_2$). Tỷ lệ các khối $SiO_4$ giảm dần, nhường chỗ cho các cấu trúc bát diện $SiO_6$ thông qua mắt xích trung gian $SiO_5$.

  2. Cơ chế nén ép pha lỏng nhôm-silicat qua tái cấu trúc Simplex: Trong hệ $Al_2O_3 \cdot 2SiO_2$ lỏng dưới áp suất 0–30 GPa, sự cô đặc thể tích diễn ra chủ yếu do sự co rút mạnh mẽ của các Void-simplex (VS) (bán kính trung bình $R_V$ giảm từ $1,35\text{ \AA}$ xuống dưới $0,95\text{ \AA}$) kết hợp với sự gia tăng của các nguyên tử oxy ba phối trí ([3]O, chiếm từ 20% đến 31%). Các nguyên tử Al và Si tạo thành hai phân mạng con đan xen nhưng không hòa tan lý tưởng, hình thành sự phân tách pha vi mô rõ rệt ở cấp độ nano.

  3. Sự chuyển đổi chức năng cấu trúc của ion $Pb^{2+}$ theo nồng độ: Trong hệ nóng chảy $xPbO \cdot (1-x)SiO_2$ tại 3000 K:

    • Ở vùng $x < 0,35$: Chì đóng vai trò modifier làm đứt gãy mạng Si-O, tạo ra các đường ống $Pb-O$ biệt lập.
    • Ở vùng $0,35 \le x \le 0,50$: Hình thành mạng không gian ba chiều gồm các kênh ống đan lồng nhau.
    • Ở vùng $x > 0,50$: Đa diện $PbO_y$ ($y=3, 4$) đóng vai trò trung tâm tạo mạng (former), chiếm ưu thế hơn cả $SiO_4$, trong khi hàm lượng oxy cầu nối biến đổi tỷ lệ nghịch với sự gia tăng của $NBO$ và $FO$.

  4. Bằng chứng định lượng về Động học không đồng nhất trong Natri-Silicat: Phân tích đa diện Voronoi oxy trong $Na_2O \cdot 2SiO_2$ và $Na_2O \cdot 4SiO_2$ làm sáng tỏ cơ chế khuếch tán ion:

    "Các nguyên tử Na phân bố không đồng nhất trong các đa diện quanh nguyên tử O, nhưng chúng chủ yếu nằm trong các đa diện NBO và FO với khoảng 75.86 % tổng số các nguyên tử Na."

    Đồng thời, luận án phát hiện:

    "Mật độ các nguyên tử Na trong cụm O-O lớn hơn mười lần so với miền Si-O."

    Khuếch tán của $Na^+$ diễn ra theo hai cơ chế song song: (i) bước nhảy đơn lẻ giữa hai đa diện Voronoi oxy liền kề và (ii) sự dịch chuyển tập thể (collective motion) của một chuỗi 4–5 ion $Na$ trượt dọc theo đường rãnh hình thành bởi cụm oxy không cầu nối.

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập cơ sở vi mô vững chắc cho lý thuyết dẫn truyền ion phi cổ điển trong môi trường thủy tinh; bổ sung cơ chế chuyển pha vô định hình dưới áp suất cao cho ngành khoáng vật học và địa vật lý magma.
  • Ý nghĩa phương pháp luận: Chuẩn hóa quy trình phân tích hình học topo (Simplex kết hợp Shell-Core và Voronoi theo thời gian), tạo tiền đề nghiên cứu cho các hệ vật liệu phức tạp khác như geopolymer, thủy tinh hoạt tính sinh học và pin thể rắn.
  • Ý nghĩa thực tiễn & công nghiệp:
    • Tối ưu hóa thành phần thủy tinh chứa chì ($PbO > 50\text{ mol}%$) phục vụ chế tạo kính chắn bức xạ gamma và tia-X có độ bền cơ học cao.
    • Cung cấp định hướng thiết kế cấu trúc mạng vi mô cho chất điện phân rắn trong pin Natri thể rắn (SSNBs), tận dụng các kênh dẫn liên kết $NBO-FO$ nhằm đạt độ dẫn ion tối ưu ở nhiệt độ phòng ($D_{Na} > 10^{-5}\text{ cm}^2/\text{s}$).

Limitations và Future Research

Luận án thừa nhận một số giới hạn phương pháp luận khách quan:

  1. Giới hạn của thế tương tác cổ điển: Các hàm thế BKS và Born–Mayer dạng cặp hữu hiệu chưa phản ánh đầy đủ hiệu ứng tái phân bố đám mây điện tử (polarizability) và liên kết cộng hóa trị định hướng mạnh ở các trạng thái biến dạng cực hạn của cặp $Pb-O$ và $Al-O$.
  2. Thang thời gian mô phỏng (Time-scale limit): Giới hạn mô phỏng trong khoảng 10–150 ps chưa đủ dài để quan sát hiện tượng hồi phục cấu trúc chậm (alpha-relaxation) khi vật liệu tiến sát nhiệt độ chuyển dịch thủy tinh $T_g$.
  3. Kích thước hộp mô phỏng (Finite-size effect): Quy mô từ 1998 đến 3025 nguyên tử dù đại diện tốt cho trật tự gần và trật tự tầm trung nhưng chưa phản ánh trọn vẹn các thăng giáng mật độ tầm xa ($q < 0,1\text{ \AA}^{-1}$).

Chương trình nghiên cứu tiếp nối trong tương lai:

  • Hướng 1: Ứng dụng thế tương tác máy học (Machine Learning Interatomic Potentials - MLIPs) huấn luyện từ tính toán ab initio DFT để mô phỏng chính xác hơn hiệu ứng phân cực của cặp electron tự do trên ion $Pb^{2+}$.
  • Hướng 2: Mở rộng nghiên cứu sang hệ bốn nguyên $Na_2O-Al_2O_3-SiO_2$ và các hệ borosilicat nhằm khảo sát "hiệu ứng kiềm hỗn hợp" (mixed-alkali effect).
  • Hướng 3: Áp dụng các thuật toán tăng tốc lấy mẫu thời gian (Accelerated Molecular Dynamics) để nghiên cứu quá trình kết tinh và chuyển dịch thủy tinh ở thang thời gian micro-giây.

Tác động và ảnh hưởng

  • Ảnh hưởng học thuật: Các bài báo công bố từ luận án trên các tạp chí chuyên ngành vật lý và khoa học vật liệu quốc tế cung cấp bộ dữ liệu chuẩn (benchmark data) về hàm phân bố xuyên tâm, phân bố góc và phổ động học, dự kiến thu hút trích dẫn cao trong cộng đồng mô phỏng silicat và địa hóa học.
  • Chuyển đổi công nghiệp: Thúc đẩy nghiên cứu chế tạo pin natri trạng thái rắn (SSNBs) nhằm thay thế pin Lithium-ion trong lưu trữ năng lượng tái tạo quy mô lưới, giải quyết bài toán phụ thuộc nguyên liệu hiếm.
  • An toàn hạt nhân và môi trường: Làm sáng tỏ cơ chế giữ các ion kim loại nặng nguy hiểm trong mạng thủy tinh chì-silicat và nhôm-silicat, đóng góp trực tiếp vào công nghệ thủy tinh hóa chất thải phóng xạ hạt nhân (nuclear waste vitrification).

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận phân tích Simplex, Shell-Core và thuật toán Voronoi nguyên tử oxy độc đáo để áp dụng cho các đề tài mô phỏng chất rắn vô định hình.
  • Các nhà khoa học vật lý lý thuyết và tính toán: Sử dụng bộ tham số thế và cơ chế động học không đồng nhất để phát triển các mô hình giải tích về độ nhớt và độ dẫn ion.
  • Kỹ sư R&D trong ngành công nghiệp bán dẫn và pin: Khai thác các quy luật biến tính mạng thủy tinh để tạo ra màng mỏng cách điện silica chất lượng cao và chất điện phân rắn dẫn ion natri siêu nhanh.
  • Các nhà hoạch định chính sách công nghệ: Nắm bắt cơ sở khoa học để xây dựng tiêu chuẩn kỹ thuật về vật liệu an toàn bức xạ và lộ trình phát triển chuỗi cung ứng pin lưu trữ thế hệ mới.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Mạng ngẫu nhiên biến tính (MRN) của Greaves sang không gian động học thời gian thực, chứng minh rằng sự dẫn ion trong silicat kiềm không phải là quá trình khuếch tán ngẫu nhiên đẳng hướng mà là sự vận chuyển có định hướng trong các kênh cấu trúc nano tạo bởi các đa diện chứa oxy không cầu nối ($NBO$) và oxy tự do ($FO$). Mật độ ion $Na^+$ trong các cụm $O-O$ cao hơn gấp 10 lần so với miền mạng con $Si-O$.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây là gì?

So với các nghiên cứu của Prescher et al. (2017) [26] và Wilding et al. (2010) [94] chỉ dừng lại ở các hàm phân bố bán kính cắt truyền thống, luận án đã lần đầu tiên kết hợp đồng thời bộ ba công cụ topo: Phân tích Simplex (tách biệt Void/Oxy/Cation-simplex), Mô hình Shell-Core (đo kích thước cụm phân tách pha) và Phân tích Đa diện Voronoi theo thời gian (theo dõi dòng cation $\langle n_{VYP} \rangle$ dịch chuyển qua các cửa sổ đa diện).

3. Phát hiện bất ngờ nhất (counter-intuitive) thu được từ dữ liệu là gì?

Phát hiện bất ngờ nhất là hiện tượng tách đôi đỉnh thứ nhất của hàm phân bố xuyên tâm Si-Si ở mật độ cao ($3,994\text{ g/cm}^3$) trong silica vô định hình do sự xuất hiện đột ngột của liên kết chia sẻ cạnh và mặt giữa các khối bát diện $SiO_6$, trái ngược hoàn toàn với quan niệm truyền thống cho rằng mạng silica chỉ duy trì liên kết chia sẻ đỉnh oxy.

4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả mô phỏng không?

Luận án cung cấp chi tiết toàn bộ tham số: thế BKS cho $SiO_2$, bộ hệ số Born–Mayer cho $Al_2O_3-SiO_2$ và $PbO-SiO_2$, các hằng số lực góc ba phần tử Si-O-Si, quy trình tích phân Verlet với bước thời gian $\Delta t = 1,0\text{ fs}$, các phương trình tính tổng Ewald khử sai số Coulomb và quy trình ủ nhiệt – cân bằng nhiệt động học NPT/NVT.

5. Định hướng nghiên cứu 10 năm tiếp theo từ kết quả luận án?

Lộ trình 10 năm hướng tới: (1) Tích hợp học sâu (Deep Learning) để xây dựng thế nguyên tử lượng tử Neural Network Potentials; (2) Thiết kế tối ưu hóa trên máy tính các màng điện phân rắn natri-silicat pha tạp đa nguyên nhằm đạt độ dẫn ion ngang bằng chất điện phân lỏng ở nhiệt độ phòng; (3) Khảo sát các pha cấu trúc silicat sâu trong lớp manti Trái đất ở điều kiện siêu áp suất trên 100 GPa.

Kết luận

  1. Xây dựng thành công hệ thống mô hình mô phỏng ĐLHPT quy mô lớn cho 4 nhóm vật liệu silicat trọng điểm: $SiO_2$, $Al_2O_3 \cdot 2SiO_2$, $xPbO \cdot (1-x)SiO_2$ và $Na_2O \cdot (2,4)SiO_2$, đạt độ tương thích chuẩn xác với thực nghiệm XRD, NMR và tính toán ab initio.
  2. Khám phá cơ chế vi mô của sự chuyển pha nén trong silica vô định hình: quá trình chuyển dịch cấu trúc $SiO_4 \rightarrow SiO_5 \rightarrow SiO_6$ kéo theo sự xuất hiện của liên kết chung cạnh/mặt làm tách đỉnh $g_{Si-Si}(r)$.
  3. Giải mã sự co đặc thể tích và phân tách pha vi mô trong pha lỏng nhôm-silicat (AS2) dưới áp suất 0–30 GPa thông qua sự thu hẹp bán kính Void-simplex và sự tái tổ hợp của các O-triclusters.
  4. Xác lập ranh giới chuyển tiếp chức năng cấu trúc của $Pb^{2+}$ tại mốc 50 mol%: từ yếu tố phá vỡ mạng (modifier) thành yếu tố tạo mạng (former) dạng ống liên kết 3D.
  5. Cung cấp bằng chứng định lượng vững chắc về hiện tượng động học không đồng nhất trong $Na_2O-SiO_2$: ion $Na^+$ khuếch tán ưu tiên dọc theo các chuỗi đa diện chứa $NBO/FO$ với mật độ cao gấp 10 lần so với miền $Si-O$ nền.
  6. Đặt nền móng phương pháp luận topo tính toán vững chắc cho các nghiên cứu thiết kế vật liệu thủy tinh kỹ thuật cao, pin thể rắn và gốm hạt nhân trong tương lai.