Tổng quan luận án
Tính cấp thiết và khoảng trống nghiên cứu
Silíc điôxít ($\text{SiO}_2$) và magnesium silicate ($\text{MgSiO}_3$) là hai hợp chất chiếm hơn 80% thành phần lớp vỏ và manti của Trái Đất, đồng thời là vật liệu quan trọng trong công nghệ gốm sứ, sợi quang, thủy tinh chịu nhiệt và xử lý chất thải hạt nhân bằng phương pháp thủy tinh hóa. Ở trạng thái lỏng và nhiệt độ cao, $\text{SiO}_2$ và $\text{MgSiO}_3$ thuộc nhóm chất lỏng có cấu trúc mạng (network-forming liquids), đặc trưng bởi các liên kết định hướng tạo thành mạng không gian ba chiều.
Mặc dù đã có nhiều nghiên cứu thực nghiệm và mô hình hóa, các vấn đề cơ bản về cấu trúc vi mô và cơ chế động học dưới điều kiện áp suất cao vẫn còn nhiều tranh luận và chưa đạt được sự đồng thuận:
- Bản chất vi mô của sự thay đổi cấu trúc mạng dưới tác động của áp suất nén (từ $0$ đến $30\text{ GPa}$).
- Tính xác thực của mô hình hai trạng thái (two-state model) đối với hiện tượng đa thù hình ở trạng thái lỏng, cũng như cấu trúc không gian chi tiết của pha mật độ thấp (Low-Density Liquid - LDL) và pha mật độ cao (High-Density Liquid - HDL).
- Cơ chế vi mô chi phối hiện tượng động học không đồng nhất, hiện tượng thuyên giảm động học (kinetics slowdown) và hiện tượng khuếch tán dị thường khi áp suất tăng.
- Mối liên hệ giữa cấu trúc địa phương của các cation ($\text{Si}^{4+}$, $\text{Mg}^{2+}$) và mức độ polymer hóa mạng khi pha tạp $\text{MgO}$ vào $\text{SiO}_2$ ở trạng thái lỏng dưới áp suất cao.
Luận án được thực hiện nhằm giải quyết các khoảng trống trên thông qua phương pháp mô phỏng động lực học phân tử (ĐLHPT).
Mục tiêu nghiên cứu
- Mục tiêu tổng quát: Nâng cao hiểu biết về mối liên hệ giữa cấu trúc vi mô và các đặc trưng động học của các chất lỏng có cấu trúc mạng.
- Mục tiêu cụ thể:
- Xây dựng mô hình tính toán ĐLHPT cho hệ $\text{SiO}_2$ lỏng và $\text{MgSiO}_3$ lỏng trong dải nhiệt độ $2600 \div 3500\text{ K}$ và áp suất $0 \div 30\text{ GPa}$.
- Làm rõ cơ chế của hiện tượng đa thù hình, hiện tượng khuếch tán dị thường và động học không đồng nhất trong $\text{SiO}_2$ lỏng dựa trên phân tích đơn vị phối trí (ĐVPT) và quá trình chuyển đổi giữa các ĐVPT.
- Xác định đặc trưng cấu trúc địa phương của $\text{Si}$ và $\text{Mg}$, trật tự tầm trung, mức độ polymer hóa và hiện tượng tách pha vi mô trong $\text{MgSiO}_3$ lỏng dưới tác động của áp suất.
Câu hỏi và giả thuyết nghiên cứu
- Câu hỏi 1: Quá trình chuyển pha cấu trúc trong $\text{SiO}_2$ lỏng dưới áp suất nén diễn ra liên tục hay gián đoạn, và được cấu thành từ những đơn vị phối trí cơ bản nào?
- Câu hỏi 2: Động học khuếch tán trong chất lỏng có cấu trúc mạng diễn ra do các chuyển đổi riêng lẻ hay do chuyển động tập thể tương quan, và tại sao hệ số khuếch tán lại tăng khi áp suất tăng trong một dải mật độ nhất định?
- Câu hỏi 3: Sự xuất hiện của ion $\text{Mg}^{2+}$ phá vỡ mạng $\text{Si-O}$ như thế nào, và liệu có tồn tại sự phân tách không gian giữa các vùng giàu $\text{Si}$ và vùng giàu $\text{Mg}$ ở điều kiện nhiệt độ cao hay không?
- Giả thuyết nghiên cứu: Cấu trúc của $\text{SiO}_2$ lỏng gồm hai pha vi mô cùng tồn tại (pha mật độ thấp gồm các tứ diện $\text{SiO}_4$ liên kết qua oxy cầu $\text{OSi}_2$; pha mật độ cao gồm các đa diện $\text{SiO}_5, \text{SiO}_6$ liên kết qua oxy cầu $\text{OSi}_3$); sự khuếch tán chỉ được dẫn dắt bởi các chuyển đổi có ích mang tính tương quan tập thể; và trong $\text{MgSiO}_3$ lỏng luôn tồn tại hiện tượng tách pha vi mô hình thành các mạng con $\text{Si-O}$ và $\text{Mg-O}$.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Cấu trúc không gian và tính chất động học của vật liệu $\text{SiO}_2$ và $\text{MgSiO}_3$ ở trạng thái lỏng.
- Phạm vi không gian mô phỏng và kích thước mẫu:
- Hệ $\text{SiO}_2$ lỏng: Khối mô phỏng gồm 1998 nguyên tử ($666\text{ Si}$ và $1332\text{ O}$).
- Hệ $\text{MgSiO}_3$ lỏng: Khối mô phỏng gồm 5000 nguyên tử ($1000\text{ Si}$, $1000\text{ Mg}$ và $3000\text{ O}$).
- Phạm vi điều kiện nhiệt động:
- Nhiệt độ: Khảo sát từ $2600\text{ K}$ đến $3500\text{ K}$ (chuẩn bị mẫu ban đầu tại $6000\text{ K}$).
- Áp suất: Khảo sát từ $0\text{ GPa}$ đến $30\text{ GPa}$ (tại các mức $0, 5, 10, 15, 20, 25, 30\text{ GPa}$).
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu trước
-
Cấu trúc và hiện tượng đa thù hình của $\text{SiO}_2$:
- Cấu trúc tinh thể ($\alpha, \beta$-quartz, tridymite, cristobalite) chuyển đổi trật tự xa sang trật tự gần và trung bình khi nóng chảy thành trạng thái vô định hình và lỏng [12, 34, 72].
- Mishima và cộng sự (1984) [58, 59] phát hiện hiện tượng đa thù hình vô định hình trong nước đá, mở đường cho giả thuyết chuyển pha loại 1 lỏng - lỏng trong các mạng tứ diện.
- Mukherjee và cộng sự [30] ghi nhận bằng chứng thực nghiệm về chuyển pha loại 1 trong $\text{SiO}_2$ vô định hình.
- Sato và cộng sự [73], Min Wu và cộng sự [54] dùng tán xạ tia X chỉ ra dưới áp suất $>10 \div 12\text{ GPa}$, số phối trí của $\text{Si}$ tăng từ 4 lên 5 và 6; $\text{SiO}_6$ chiếm ưu thế ở áp suất cao.
- Các mô hình "two-state", "two-species", "two-domain" [22, 23] giải thích sự tồn tại đồng thời của miền mật độ cao và miền mật độ thấp trong chất lỏng.
-
Các lý thuyết và mô hình động học chất lỏng:
- Lý thuyết Mode-Coupling [77] mô tả động học nhiệt độ cao ($\alpha$ và $\beta$ relaxation) nhưng dự đoán nhiệt độ tới hạn cao hơn thực nghiệm.
- Lý thuyết Adam-Gibbs [28] dựa trên vùng tái sắp xếp tương quan (Cooperatively Rearranging Regions - CRR) với quan hệ $D \sim \exp(-A/TS_{\text{conf}})$, giải thích hiện tượng thuyên giảm động học khi hạ nhiệt độ về gần nhiệt độ chuyển pha thủy tinh $T_g$.
- Lý thuyết thể tích tự do [19] và mô hình thấm (Percolation model) [31] về mạng lưới các vùng linh động.
- Các nghiên cứu động học không đồng nhất qua biến động học: Garrahan và Chandler [40], Vogel và cộng sự [52], hàm tương quan 4 điểm [21, 53], khảo sát độ dời phi Gauss của Horbach và Kob [76].
- Hiện tượng khuếch tán dị thường trong $\text{SiO}_2$ lỏng: Liang và cộng sự [51], Shell và cộng sự [63], Saika-Voivod và cộng sự [69] gắn liền khuếch tán với sự phá vỡ mạng tứ diện và biến đổi mô-đun khối.
-
Nghiên cứu cấu trúc silicate và $\text{MgSiO}_3$:
- Wilding và cộng sự [46], Cormier và cộng sự [44] dùng tán xạ Neutron kết hợp Monte Carlo đảo xác định khoảng cách liên kết $\text{Si-O} \approx 1.61 \div 1.64\text{ \AA}$, $\text{Mg-O} \approx 2.00 \div 2.08\text{ \AA}$, số phối trí $Z_{\text{Mg-O}} \approx 4.5 \div 5.1$.
- Mạng $\text{MgSiO}_3$ ở $0\text{ GPa}$ chủ yếu gồm đơn vị $Q^2$, chứa $62.7%$ oxy không cầu (NBO) và $4.6%$ oxy tự do [44].
- Nghiên cứu mô phỏng của de Koker và cộng sự [26] dùng thế Oganov (OG) đến $50\text{ GPa}$ chỉ ra khi tăng áp suất, NBO giảm, oxy cầu (BO) tăng và xuất hiện các tri-cluster ($\text{OSi}_3$) [66].
Khoảng trống luận án lựa chọn giải quyết
- Thiếu mô tả cấu trúc định lượng về cách thức 5 ĐVPT kết nối trong không gian để tạo nên hai pha mật độ khác nhau của $\text{SiO}_2$ lỏng theo mô hình hai trạng thái.
- Chưa làm rõ cơ chế vi mô của khuếch tán: phân biệt giữa các chuyển đổi phối trí đơn thuần và các "chuyển đổi có ích" gây ra độ dịch chuyển thực sự của nguyên tử.
- Thiếu nghiên cứu chi tiết tại Việt Nam về mô phỏng cấu trúc ba thành phần của $\text{MgSiO}_3$ lỏng dưới áp suất cao ($0 \div 30\text{ GPa}$), đặc biệt là sự phân bố không gian của các chuỗi/đám $\text{SiO}_x$ và $\text{MgO}_y$ (hiện tượng tách pha vi mô).
Cơ sở lý thuyết và phương pháp nghiên cứu
Lý thuyết, khái niệm và khung phân tích
- Khung phân tích cấu trúc mạng: Cấu trúc chất lỏng được mô tả thông qua tập hợp các ĐVPT cơ bản: đối với $\text{Si}$ là $\text{SiO}_x$ ($x = 4, 5, 6$), đối với $\text{O}$ là $\text{OSi}_y$ hoặc $\text{OT}_y$ ($T = \text{Si, Mg}$), và đối với $\text{Mg}$ là $\text{MgO}_y$ ($y = 4, 5, 6$).
- Phân loại nguyên tử oxy: Oxy cầu (Bridging Oxygen - BO, liên kết với 2 nguyên tử $\text{Si}$ hoặc tạo cầu nối cation), oxy không cầu (Non-Bridging Oxygen - NBO, liên kết với 1 nguyên tử $\text{Si}$) và oxy tự do (Free Oxygen, không liên kết với $\text{Si}$, chỉ liên kết với $\text{Mg}$).
- Mô hình hai trạng thái (Two-State Model): Chất lỏng được xem như hỗn hợp vi mô của pha mật độ thấp (LDL) và pha mật độ cao (HDL) với tỷ phần thay đổi theo áp suất.
- Lý thuyết động học chuyển đổi: Khuếch tán được mô tả thông qua tần suất và bản chất hình học của các bước chuyển đổi phối trí $\text{SiO}x \leftrightarrow \text{SiO}{x'}$ và $\text{OSi}y \leftrightarrow \text{OSi}{y'}$.
Phương pháp mô phỏng động lực học phân tử
Phương pháp ĐLHPT giải hệ phương trình chuyển động Newton cho $N$ hạt sử dụng thuật toán tích phân Verlet:
$$\vec{r}_i(t + \delta t) = 2\vec{r}_i(t) - \vec{r}_i(t - \delta t) + \frac{\vec{F}_i(t)}{m_i}\delta t^2$$
- Bước thời gian mô phỏng: $\delta t = 0.47\text{ fs}$.
- Tương tác tĩnh điện tầm xa (Coulomb) được tính toán bằng kỹ thuật gần đúng Ewald - Hansen:
$$U_{\text{ES}} = U_r + U_m + U_0$$
- Điều kiện biên: Áp dụng điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) theo cả ba chiều không gian $L_x \times L_y \times L_z$.
- Môi trường mô phỏng: Sử dụng mô hình thống kê $NPT$ (áp suất và nhiệt độ không đổi) để nén và điều chỉnh mật độ, sau đó hồi phục cân bằng bằng mô hình $NVE$ (thể tích và năng lượng không đổi) trong $10^6 \div 10^7$ bước mô phỏng.
Thế tương tác sử dụng trong mô phỏng
1. Hệ $\text{SiO}_2$ (Thế BKS - van Beest, Kramer, van Santen)
Hàm thế có dạng:
$$U(r_{ij}) = \frac{q_i q_j e^2}{r_{ij}} + A_{ij}\exp(-B_{ij}r_{ij}) - C_{ij}r_{ij}^{-6}$$
| Cặp tương tác |
$A_{ij}\text{ (eV)}$ |
$B_{ij}\text{ (\AA}^{-1}\text{)}$ |
$C_{ij}\text{ (eV}\cdot\text{\AA}^6\text{)}$ |
Điện tích hiệu dụng ($e$) |
| $\text{O - O}$ |
1388.7730 |
2.76000 |
175.0000 |
$q_{\text{O}} = -1.2000$ |
| $\text{Si - O}$ |
18003.7572 |
4.87318 |
133.5381 |
$q_{\text{Si}} = +2.4000$ |
| $\text{Si - Si}$ |
0.0 |
0.0 |
0.0 |
— |
2. Hệ $\text{MgSiO}_3$ (Thế OG - Oganov et al.)
Hàm thế có dạng tương tự BKS, tham số tính toán từ nguyên lý ban đầu có tính đến thăng giáng nhiệt:
| Cặp tương tác |
$A_{ij}\text{ (eV)}$ |
$B_{ij}\text{ (\AA}^{-1}\text{)}$ |
$C_{ij}\text{ (eV}\cdot\text{\AA}^6\text{)}$ |
Điện tích hiệu dụng ($e$) |
| $\text{Si - O}$ |
1137.28 |
4.0988 |
40.00 |
$q_{\text{Si}} = +2.9043$ |
| $\text{Mg - O}$ |
674.31 |
3.3283 |
4.00 |
$q_{\text{Mg}} = +1.6049$ |
| $\text{O - O}$ |
12974.70 |
4.4172 |
145.20 |
$q_{\text{O}} = -1.5031$ |
| $\text{Mg - Mg}$ |
0.0 |
0.0 |
0.0 |
— |
| $\text{Mg - Si}$ |
0.0 |
0.0 |
0.0 |
— |
| $\text{Si - Si}$ |
0.0 |
0.0 |
0.0 |
— |
Kỹ thuật phân tích cấu trúc và động học
- Hàm phân bố xuyên tâm cặp (PBXT - $g_{\alpha\beta}(r)$):
$$g_{\alpha\beta}(r) = \frac{V}{N_\alpha N_\beta 4\pi r^2 \delta r} \sum_{i=1}^{N_\alpha} \sum_{j=1}^{N_\beta} \langle \delta(r - |\vec{r}_i - \vec{r}_j|) \rangle$$
- Bán kính ngắt ($r_c$) và số phối trí ($Z_{\alpha\beta}$):
- Trong $\text{SiO}_2$: $r_c(\text{Si-O}) = 2.70\text{ \AA}$.
- Trong $\text{MgSiO}_3$: $r_c(\text{Si-O}) = 2.38\text{ \AA}$, $r_c(\text{Mg-O}) = 3.20\text{ \AA}$.
- Số phối trí: $Z_{\alpha\beta} = 4\pi \rho_\beta \int_0^{r_c} r^2 g_{\alpha\beta}(r)dr$.
- Phân bố góc liên kết: Hàm phân bố $a_{ijk}(\theta)$ cho các góc $\text{O-Si-O}$, $\text{Si-O-Si}$, $\text{O-Mg-O}$, $\text{Mg-O-Si}$.
- Hệ số khuếch tán ($D$): Xác định qua hệ thức Einstein từ độ dịch chuyển bình phương trung bình (MSD):
$$D = \lim_{t \to \infty} \frac{\langle |\vec{r}(t) - \vec{r}(0)|^2 \rangle}{6t}$$
- Phân tích động học chuyển đổi và tạo đám:
- Theo dõi số lượng chuyển đổi phối trí theo thời gian: chuyển đổi có ích (useful transition - làm tâm dịch chuyển vị trí thực) và chuyển đổi không có ích (chuyển đổi dao động quanh vị trí cân bằng).
- Phân tích tạo đám không gian của 5 nhóm nguyên tử linh động/kém linh động (A1 đến A5) để xác định tính không đồng nhất động học.
Nội dung chính theo từng chương
Chương 1: Tổng quan về hai vật liệu $\text{SiO}_2$ và $\text{MgSiO}_3$
Chương này tổng kết toàn diện các công trình lý thuyết, thực nghiệm và mô phỏng trên thế giới về $\text{SiO}_2$ và $\text{MgSiO}_3$ ở các pha tinh thể, vô định hình và lỏng:
- Trình bày chi tiết đặc điểm cấu trúc mạng tứ diện của $\text{SiO}_2$, các dạng đa hình tinh thể ($\alpha, \beta$-quartz, tridymite, cristobalite) và sự suy giảm trật tự tầm xa khi chuyển sang pha lỏng ở nhiệt độ $>1723^\circ\text{C}$.
- Tổng hợp các kết quả thực nghiệm nhiễu xạ tia X, nhiễu xạ neutron và cộng hưởng từ hạt nhân (NMR) về sự tăng số phối trí của $\text{Si}$ từ 4 lên 5 và 6 khi nén ép dưới áp suất cao.
- Phân tích cơ sở lý thuyết của các mô hình chuyển pha lỏng - lỏng, mô hình hai trạng thái, lý thuyết Adam-Gibbs, lý thuyết thể tích tự do và mô hình thấm trong chất lỏng siêu lạnh và chất lỏng mạng.
- Điểm luận các nghiên cứu về hệ silicate ba thành phần $\text{MgSiO}_3$, vai trò của cation biến tính $\text{Mg}^{2+}$ trong việc phá vỡ mạng polymer hóa $\text{Si-O}$, tạo ra các đơn vị $Q^n$, oxy không cầu (NBO), oxy tự do và các cụm tri-cluster ($\text{OSi}_3$).
Chương 2: Phương pháp tính toán
Chương 2 trình bày chi tiết quy trình xây dựng mô hình ĐLHPT và các thuật toán phân tích:
- Trình bày cơ sở toán học của thuật toán tích phân Verlet, kỹ thuật tổng Ewald - Hansen khử tương tác Coulomb và các tham số thế cặp BKS, OG.
- Quy trình hạ nhiệt và cân bằng nhiệt động: Mẫu được tạo ngẫu nhiên, nung chảy ở $6000\text{ K}$ (50.000 bước), hạ nhiệt từ từ qua các mức $5000\text{ K}, 4000\text{ K}$ xuống $3500\text{ K}$, sau đó nén đẳng nhiệt - đẳng áp ($NPT$) đến các áp suất $0 \div 30\text{ GPa}$, cuối cùng hồi phục trong điều kiện $NVE$ ($10^6 \div 10^7$ bước).
- Định nghĩa giải thuật xác định ĐVPT, giải thuật tính hàm PBXT cặp, phân bố góc liên kết và kỹ thuật trực quan hóa 3D phân biệt màu cho từng loại nguyên tử và ĐVPT.
- Thiết lập phương pháp đếm số chuyển đổi phối trí theo thời gian và phương pháp phân tách các nhóm nguyên tử theo độ linh động để khảo sát cấu trúc đám.
Chương 3: Cấu trúc và động học của $\text{SiO}_2$ lỏng
1. Hiện tượng đa thù hình và cấu trúc vi mô hai pha
Kết quả mô phỏng khẳng định mạng $\text{SiO}_2$ lỏng được cấu thành từ 5 ĐVPT: $\text{SiO}_4, \text{SiO}_5, \text{SiO}_6, \text{OSi}_2, \text{OSi}_3$.
- Ở $0\text{ GPa}$, hệ chủ yếu gồm tứ diện $\text{SiO}_4$ ($>90%$) liên kết với nhau qua oxy cầu $\text{OSi}_2$, tạo thành mạng mở mật độ thấp (LDL).
- Khi tăng áp suất, tỷ phần $\text{SiO}_4$ và $\text{OSi}_2$ giảm đơn điệu; $\text{SiO}_5$ tăng và đạt cực đại ở vùng $10 \div 15\text{ GPa}$; $\text{SiO}_6$ và $\text{OSi}_3$ tăng mạnh ở áp suất cao ($>15\text{ GPa}$), hình thành pha xếp chặt mật độ cao (HDL).
- Phân tích trực quan hóa chỉ ra hai pha này không trộn lẫn đồng nhất ở cấp độ nguyên tử mà kết cụm thành các miền không gian riêng biệt, giải thích bản chất cấu trúc của mô hình hai trạng thái.
Bảng 1: Đặc trưng cấu trúc của $\text{SiO}2$ lỏng ở các điều kiện áp suất và nhiệt độ khác nhau (vị trí đỉnh thứ nhất $r{\text{Si-O}}$ và số phối trí trung bình $Z_{\text{Si-O}}$)
| Áp suất (GPa) |
Nhiệt độ (K) |
$r_{\text{Si-O}}\text{ (\AA)}$ |
Tỷ phần $\text{SiO}_4\text{ (%)}$ |
Tỷ phần $\text{SiO}_5\text{ (%)}$ |
Tỷ phần $\text{SiO}_6\text{ (%)}$ |
$Z_{\text{Si-O}}$ |
| 0 |
3000 |
1.61 |
92.4 |
7.3 |
0.3 |
4.08 |
| 5 |
3000 |
1.62 |
65.1 |
31.8 |
3.1 |
4.38 |
| 10 |
3000 |
1.63 |
34.2 |
52.4 |
13.4 |
4.79 |
| 15 |
3000 |
1.64 |
14.8 |
54.6 |
30.6 |
5.16 |
| 20 |
3000 |
1.66 |
5.2 |
43.1 |
51.7 |
5.47 |
| 25 |
3000 |
1.67 |
1.6 |
28.9 |
69.5 |
5.68 |
| 30 |
3000 |
1.68 |
0.4 |
16.5 |
83.1 |
5.83 |
2. Động học không đồng nhất và cơ chế khuếch tán
- Khảo sát độ dịch chuyển bình phương trung bình của 5 nhóm nguyên tử (A1 đến A5) phân loại theo tốc độ chuyển đổi phối trí cho thấy sự phân tách rõ rệt: nhóm nguyên tử linh động tạo thành các đám linh động (mobile clusters), trong khi các nguyên tử kém linh động tạo thành khung cứng (rigid framework).
- Mối liên hệ giữa số lượng chuyển đổi phối trí và hệ số khuếch tán $D$: Khuếch tán không tỷ lệ thuận đơn giản với tổng số lần chuyển đổi phối trí mà phụ thuộc vào tỷ lệ "chuyển đổi có ích".
- Dưới tác dụng của áp suất trong dải $0 \div 10\text{ GPa}$, mạng $\text{SiO}_4$ bị phá vỡ cục bộ, tạo điều kiện cho các chuyển động phối hợp tập thể của nhóm nguyên tử, dẫn tới hiện tượng hệ số khuếch tán $D$ của $\text{Si}$ và $\text{O}$ tăng đạt cực đại (khuếch tán dị thường) trước khi giảm mạnh ở áp suất cao do sự nén chặt của mạng $\text{SiO}_6$.
Chương 4: Cấu trúc của $\text{MgSiO}_3$ lỏng
1. Cấu trúc địa phương của $\text{Si}$ và $\text{Mg}$
- Khi đưa ion $\text{Mg}^{2+}$ vào mạng, tương tác giữa $\text{Mg}$ và $\text{O}$ phá vỡ tính liên tục của mạng $\text{Si-O}$. Tuy nhiên, cấu trúc địa phương của $\text{Si}$ trong $\text{MgSiO}_3$ lỏng hoàn toàn tương đồng với trong $\text{SiO}_2$ lỏng: ở $0\text{ GPa}$, $\text{Si}$ chủ yếu tồn tại dưới dạng $\text{SiO}4$ với khoảng cách liên kết $r{\text{Si-O}} \approx 1.62\text{ \AA}$ và góc $\text{O-Si-O} \approx 109^\circ$. Khi tăng áp suất, $\text{Si}$ chuyển dần sang phối trí 5 và 6.
- Đối với $\text{Mg}$, cấu trúc địa phương nhạy cảm hơn nhiều với áp suất: ở $0\text{ GPa}$, số phối trí trung bình $Z_{\text{Mg-O}} \approx 4.5 \div 4.8$, phân bố gồm $\text{MgO}_4, \text{MgO}_5, \text{MgO}_6$. Khi áp suất đạt $30\text{ GPa}$, số phối trí tăng lên $>6.0$, tạo thành các đa diện $\text{MgO}_6, \text{MgO}_7$.
Bảng 2: Phân bố các loại nguyên tử oxy và cầu nối trong $\text{MgSiO}_3$ lỏng ở $T = 3500\text{ K}$ theo áp suất
| Áp suất (GPa) |
Oxy cầu $\text{BO (%)}$ |
Oxy không cầu $\text{NBO (%)}$ |
Oxy tự do $\text{Free O (%)}$ |
Cầu nối $\text{OSi}_2\text{ (%)}$ |
Cầu nối $\text{OSi}_3\text{ (%)}$ |
| 0 |
44.8 |
46.7 |
8.5 |
42.1 |
2.7 |
| 5 |
53.2 |
41.1 |
5.7 |
46.5 |
6.7 |
| 10 |
63.5 |
32.8 |
3.7 |
48.2 |
15.3 |
| 15 |
73.1 |
24.6 |
2.3 |
44.1 |
29.0 |
| 20 |
81.4 |
17.2 |
1.4 |
36.8 |
44.6 |
| 25 |
87.6 |
11.5 |
0.9 |
27.5 |
60.1 |
| 30 |
92.3 |
7.2 |
0.5 |
18.2 |
74.1 |
2. Mức độ polymer hóa và hiện tượng tách pha vi mô
- Mức độ polymer hóa của mạng silicate tăng mạnh theo áp suất: tỷ phần NBO giảm từ $46.7%$ ($0\text{ GPa}$) xuống $7.2%$ ($30\text{ GPa}$), oxy tự do giảm từ $8.5%$ xuống $0.5%$, trong khi BO tăng từ $44.8%$ lên $92.3%$.
- Khảo sát phân bố không gian và các chuỗi liên kết $\text{Si}_n\text{-O-Mg}_m$ chứng minh sự tồn tại của hiện tượng tách pha vi mô: các nguyên tử $\text{Si}$ và $\text{O}$ kết cụm tạo thành các mạng con $\text{Si-O}$, trong khi các nguyên tử $\text{Mg}$ tập trung liên kết với NBO và oxy tự do tạo thành các mạng con $\text{Mg-O}$. Hiện tượng tách pha vi mô này duy trì trong suốt dải áp suất khảo sát từ $0$ đến $30\text{ GPa}$.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học và lý luận
- Bản chất cấu trúc mô hình hai trạng thái trong $\text{SiO}_2$ lỏng: Luận án đã xác định và định lượng cụ thể 5 ĐVPT cấu thành mạng $\text{SiO}_2$ lỏng, chứng minh sự cùng tồn tại của hai pha mật độ thấp (tạo bởi $\text{SiO}_4$ liên kết qua $\text{OSi}_2$) và pha mật độ cao (tạo bởi $\text{SiO}_5, \text{SiO}_6$ liên kết qua $\text{OSi}_3$). Cấu trúc không đồng nhất không gian này là nguồn gốc trực tiếp dẫn tới hiện tượng đa thù hình và sự phân bố không đồng nhất của các chuyển đổi động học.
- Cơ chế vi mô của hiện tượng khuếch tán và động học không đồng nhất: Chứng minh quá trình khuếch tán không đồng nhất với tổng số chuyển đổi phối trí mà được chi phối bởi các "chuyển đổi có ích". Hiện tượng khuếch tán dị thường và thuyên giảm động học trong $\text{SiO}_2$ lỏng bắt nguồn từ chuyển động tương quan tập thể của các nhóm nguyên tử linh động hình thành nên các siêu phân tử/đám linh động.
- Đặc trưng cấu trúc và hiện tượng tách pha vi mô trong $\text{MgSiO}_3$ lỏng: Chỉ ra môi trường địa phương của nguyên tử $\text{Si}$ trong $\text{MgSiO}_3$ lỏng bảo toàn các đặc tính cơ bản như trong $\text{SiO}_2$ lỏng khi bị nén. Luận án lần đầu tiên làm rõ hiện tượng tách pha vi mô trong $\text{MgSiO}_3$ lỏng dưới áp suất cao ($0 \div 30\text{ GPa}$), biểu hiện qua sự tồn tại đồng thời của các vùng giàu $\text{Si}$ và vùng giàu $\text{Mg}$.
Ý nghĩa thực tiễn và ứng dụng
- Cung cấp các thông số cấu trúc vi mô, số phối trí, phân bố góc và hệ số khuếch tán định lượng của $\text{SiO}_2$ và $\text{MgSiO}_3$ lỏng ở điều kiện áp suất cao ($0 \div 30\text{ GPa}$) và nhiệt độ cao ($2600 \div 3500\text{ K}$), hỗ trợ các nhà địa vật lý xây dựng mô hình đối lưu magma và giải thích quá trình tiến hóa phân tầng manti Trái Đất.
- Cung cấp cơ sở vi mô về cơ chế bẫy cation phóng xạ trong mạng silicate thông qua việc hình thành các liên kết oxy không cầu (NBO) và oxy tự do, phục vụ tối ưu hóa công nghệ thủy tinh hóa xử lý rác thải hạt nhân.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án
- Kích thước mẫu mô phỏng giới hạn ở 1998 nguyên tử ($\text{SiO}_2$) và 5000 nguyên tử ($\text{MgSiO}_3$), dù đã sử dụng điều kiện biên tuần hoàn nhưng vẫn chịu ảnh hưởng nhất định của hiệu ứng kích thước đối với các thăng giáng cấu trúc tầm xa.
- Thế tương tác cặp BKS và OG gán điện tích hiệu dụng cố định cho các ion, chưa mô tả được sự tái phân bố điện tích động và hiệu ứng phân cực điện tử trong quá trình phá vỡ và tái lập liên kết.
- Dải áp suất nghiên cứu giới hạn ở mức tối đa $30\text{ GPa}$, chưa bao phủ toàn bộ điều kiện áp suất của phần manti dưới Trái Đất ($>100\text{ GPa}$).
- Luận án chưa tính toán trực tiếp độ nhớt và phổ dao động để so sánh toàn diện với các phép đo phổ Raman/phổ hồng ngoại ở áp suất cao.
Hướng nghiên cứu tiếp
- Mở rộng mô phỏng với kích thước mẫu lớn hơn ($>10^4 \div 10^5$ nguyên tử) và dải áp suất siêu cao ($>50 \div 135\text{ GPa}$).
- Áp dụng các thế tương tác phân cực hoặc mô phỏng nguyên lý ban đầu (ab-initio MD) trên các cấu hình chọn lọc để nâng cao độ chính xác về mặt lượng tử.
- Nghiên cứu các hệ silicate đa thành phần phức tạp hơn chứa các cation kim loại kiềm và kiềm thổ khác như $\text{Na}_2\text{O-CaO-MgO-SiO}_2$.
- Tính toán trực tiếp độ nhớt, độ dẫn nhiệt và phổ tán xạ cấu trúc động học $S(q, \omega)$ để đối chiếu trực tiếp với các phép đo thực nghiệm hiện đại.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên khảo trong lĩnh vực Vật lý chất rắn, Vật lý tính toán và Khoa học vật liệu:
- Đối với nghiên cứu sinh và học viên cao học chuyên ngành Vật lý kỹ thuật, Vật lý tính toán:
- Chương 2 cung cấp quy trình phương pháp luận chuẩn mực về thiết lập mô phỏng ĐLHPT cho hệ ôxít vô định hình và lỏng, cách khử tương tác Coulomb bằng kỹ thuật Ewald - Hansen, và các giải thuật phân tích cấu trúc vi mô (PBXT, phân bố góc, ĐVPT).
- Phương pháp phân loại chuyển đổi phối trí và kỹ thuật phân tích tạo đám nguyên tử là công cụ tham khảo hữu ích cho các bài toán về động học chất lỏng siêu lạnh.
- Đối với giảng viên và các nhà nghiên cứu trong lĩnh vực Khoa học Trái Đất và Địa chất học:
- Các bảng số liệu định lượng về số phối trí, tỷ phần BO/NBO và sự dịch chuyển cấu trúc từ tứ diện sang bát diện trong Chương 3 và Chương 4 cung cấp dữ liệu đầu vào cho các mô hình mô phỏng magma manti.
- Đối với kỹ sư và nhà nghiên cứu công nghệ vật liệu silicat và xử lý môi trường:
- Cung cấp hiểu biết cơ chế về mức độ polymer hóa và vai trò của các nguyên tử oxy cầu/không cầu trong việc giam giữ các ion kim loại, ứng dụng trong công nghệ vật liệu gốm xốp và cố định đồng vị phóng xạ.
Câu hỏi thường gặp
1. Vì sao hệ số khuếch tán trong $\text{SiO}_2$ lỏng lại biểu hiện tính dị thường khi áp suất tăng từ $0$ đến $10\text{ GPa}$?
Ở áp suất $0\text{ GPa}$, $\text{SiO}_2$ lỏng tồn tại ở trạng thái mạng tứ diện $\text{SiO}_4$ có tính polymer hóa cao và rất bền vững, cản trở sự di chuyển của các nguyên tử. Khi tăng áp suất lên $5 \div 10\text{ GPa}$, mạng tứ diện lý tưởng bị phá vỡ cục bộ với sự xuất hiện của các khuyết tật mạng và các ĐVPT trung gian $\text{SiO}_5$. Sự xuất hiện này làm suy giảm độ cứng vững của mạng, thúc đẩy các chuyển động tương quan tập thể của các nhóm nguyên tử, dẫn tới tốc độ chuyển đổi có ích tăng lên, làm hệ số khuếch tán $D$ tăng. Khi áp suất vượt quá $10 \div 15\text{ GPa}$, mạng bát diện $\text{SiO}_6$ bắt đầu hình thành với mật độ xếp chặt cao và mô-đun khối lớn, làm giảm không gian dịch chuyển và khiến hệ số khuếch tán giảm mạnh.
2. Sự khác biệt cốt lõi giữa hai pha trong mô hình hai trạng thái của $\text{SiO}_2$ lỏng là gì?
Hai pha này có cùng thành phần hóa học nhưng khác biệt hoàn toàn về cấu trúc phối trí và mật độ:
- Pha mật độ thấp (LDL): Cấu thành từ các ĐVPT tứ diện $\text{SiO}_4$ liên kết với nhau ở các đỉnh thông qua nguyên tử oxy cầu 2 phối trí ($\text{OSi}_2$). Mạng này có cấu trúc mở, góc liên kết $\text{Si-O-Si}$ lớn và chiếm ưu thế ở áp suất thấp ($0\text{ GPa}$).
- Pha mật độ cao (HDL): Cấu thành từ các ĐVPT phối trí cao $\text{SiO}_5$ và bát diện $\text{SiO}_6$ liên kết với nhau qua các nguyên tử oxy cầu 3 phối trí ($\text{OSi}_3$ hay tri-cluster), dùng chung cạnh hoặc mặt. Mạng này xếp chặt, có mật độ cao và chiếm ưu thế ở áp suất cao ($>15\text{ GPa}$).
3. Khái niệm "chuyển đổi có ích" trong phân tích động học khuếch tán được định nghĩa như thế nào?
Trong chất lỏng có cấu trúc mạng, các nguyên tử liên tục dao động và số phối trí quanh một nguyên tử trung tâm có thể thay đổi tạm thời ($\text{SiO}x \leftrightarrow \text{SiO}{x'}$). Phần lớn các chuyển đổi này mang tính dao động thuận nghịch tại chỗ quanh vị trí cân bằng và không làm thay đổi vị trí thực của nguyên tử (chuyển đổi không có ích). "Chuyển đổi có ích" là những biến đổi phối trí dẫn tới sự nhảy thực sự của nguyên tử oxy sang quả cầu phối trí của một nguyên tử silíc lân cận mới, kèm theo sự dịch chuyển khối tâm của nhóm nguyên tử, tạo nên độ dịch chuyển bình phương trung bình thực tế đóng góp vào hệ số khuếch tán $D$.
4. Cation $\text{Mg}^{2+}$ làm thay đổi mạng $\text{Si-O}$ trong $\text{MgSiO}_3$ lỏng theo cơ chế nào?
Khi đưa $\text{Mg}^{2+}$ vào hệ, các ion này đóng vai trò là chất biến tính mạng (network modifier). $\text{Mg}^{2+}$ phá vỡ các liên kết cầu $\text{Si-O-Si}$, làm giảm tỷ phần oxy cầu (BO) và chuyển hóa chúng thành oxy không cầu (NBO) và oxy tự do. Điện tích âm dư tại các nguyên tử NBO và oxy tự do sẽ liên kết với các ion $\text{Mg}^{2+}$ để bù trừ điện tích. Quá trình này cắt đứt mạng $\text{Si-O}$ liên tục thành các chuỗi và mạng con nhỏ hơn (giảm mức độ polymer hóa mạng), đồng thời tạo ra các vùng giàu $\text{Mg}$ tách biệt với các vùng giàu $\text{Si}$ (hiện tượng tách pha vi mô).
Kết luận
Luận án của NCS. Luyện Thị San đã nghiên cứu thành công đặc trưng cấu trúc và động học của vật liệu $\text{SiO}_2$ và $\text{MgSiO}_3$ ở trạng thái lỏng dưới dải áp suất $0 \div 30\text{ GPa}$ và nhiệt độ $2600 \div 3500\text{ K}$ bằng phương pháp mô phỏng động lực học phân tử. Công trình đã xác lập cơ sở vi mô của mô hình hai trạng thái trong $\text{SiO}_2$ lỏng thông qua 5 đơn vị phối trí cơ bản, đồng thời chứng minh cơ chế khuếch tán dị thường và động học không đồng nhất bắt nguồn từ các chuyển đổi phối trí có ích và chuyển động tương quan tập thể. Đối với hệ $\text{MgSiO}_3$ lỏng, luận án đã làm rõ sự bảo toàn cấu trúc địa phương của nguyên tử $\text{Si}$, quy luật biến đổi của oxy cầu/không cầu theo áp suất và chứng minh sự tồn tại của hiện tượng tách pha vi mô trong cấu trúc mạng silicate lỏng. Các kết quả đạt được có độ tin cậy cao, đóng góp dữ liệu quan trọng cho vật lý vật liệu vô định hình và khoa học Trái Đất.
