Tổng quan luận án
Tính cấp thiết và khoảng trống nghiên cứu
Vật liệu điôxít silíc ($\text{SiO}_2$) và manhê silicat ($\text{MgSiO}_3$) là hai thành phần khoáng vật chiếm hơn 80% khối lượng vỏ Trái Đất, đồng thời giữ vai trò then chốt trong nhiều ngành công nghệ cao như sản xuất gốm xốp chịu nhiệt, vật liệu nha khoa, vật liệu gạch chịu lửa và công nghệ thủy tinh hóa/ceramic hóa để cố định các đồng vị phóng xạ trong xử lý rác thải hạt nhân. Ở điều kiện nhiệt độ và áp suất cao trong lòng đất, hai hợp chất này tồn tại ở trạng thái nóng chảy (trạng thái lỏng). Việc hiểu rõ cấu trúc vi mô và cơ chế động học của chất lỏng có cấu trúc mạng (network-forming liquids) là nền tảng để giải thích các quá trình địa chất nhiệt động lực học và thiết kế vật liệu chuyên dụng.
Mặc dù đã có nhiều nghiên cứu thực nghiệm và mô phỏng trên thế giới, nhiều vấn đề cốt lõi vẫn chưa đạt được sự thống nhất hoặc còn bỏ ngỏ:
- Sự thay đổi cấu trúc vi mô và trật tự tầm trung của mạng silicat dưới tác động của áp suất cực cao chưa được mô tả đầy đủ.
- Bản chất cấu trúc của hiện tượng đa thù hình (polyamorphism) và tính xác thực của mô hình hai trạng thái (two-state model) trong pha lỏng vẫn thiếu các bằng chứng định lượng về cấu trúc hình học của từng pha mật độ.
- Cơ chế vi mô dẫn tới hiện tượng động học không đồng nhất (dynamical heterogeneity), hiện tượng thuyên giảm động học (dynamical slowdown) và khuếch tán dị thường (anomalous diffusion) chưa được liên hệ tường minh với quá trình tái sắp xếp cấu trúc mạng.
- Đối với hệ ôxít ba thành phần $\text{MgSiO}_3$ lỏng, các dữ liệu thực nghiệm ở nhiệt độ và áp suất cao còn rất hạn chế và tại Việt Nam chưa có công trình luận án nào khảo sát chuyên sâu về hệ vật liệu này.
Mục tiêu nghiên cứu
Luận án xác định mục tiêu tổng quát là nâng cao hiểu biết về cấu trúc và động học của các chất lỏng có cấu trúc mạng nói chung, làm rõ bản chất các hiện tượng vật lý vi mô trong hai hệ chất lỏng $\text{SiO}_2$ và $\text{MgSiO}_3$ ở điều kiện nhiệt độ và áp suất cao.
Các mục tiêu cụ thể được phân định theo hai hệ vật liệu:
- Đối với vật liệu $\text{SiO}_2$ lỏng: Khảo sát và giải thích cơ chế của hiện tượng đa thù hình dựa trên các đơn vị phối trí và mô hình hai trạng thái; làm rõ nguyên nhân vi mô của hiện tượng khuếch tán dị thường, hiện tượng thuyên giảm động học và hiện tượng động học không đồng nhất thông qua cơ chế chuyển đổi đơn vị phối trí và chuyển động tương quan tập thể.
- Đối với vật liệu $\text{MgSiO}_3$ lỏng: Khảo sát môi trường địa phương của nguyên tử silíc (Si) và manhê (Mg), sự phân bố các cầu nối liên kết, mức độ polymer hóa mạng silicat và kiểm chứng hiện tượng tách pha vi mô khi nén ép dưới áp suất cao.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Hai hệ vật liệu $\text{SiO}_2$ và $\text{MgSiO}_3$ ở trạng thái lỏng.
- Phạm vi không gian (kích thước mô hình):
- Mẫu $\text{SiO}_2$ lỏng gồm 1998 nguyên tử (666 nguyên tử Si và 1332 nguyên tử O).
- Mẫu $\text{MgSiO}_3$ lỏng gồm 5000 nguyên tử (1000 nguyên tử Si, 1000 nguyên tử Mg và 3000 nguyên tử O).
- Phạm vi nhiệt độ và áp suất:
- Hệ $\text{SiO}_2$ lỏng: Khảo sát trong dải nhiệt độ từ 2600 K đến 3500 K (các mốc 2600 K, 3000 K, 3200 K, 3500 K) và dải áp suất từ 0 GPa đến 30 GPa (các mốc 0, 5, 10, 15, 20, 25, 30 GPa).
- Hệ $\text{MgSiO}_3$ lỏng: Khảo sát ở nhiệt độ cố định 3500 K và dải áp suất từ 0 GPa đến 30 GPa (gồm 7 mẫu từ M1 đến M7 tương ứng với 0, 5, 10, 15, 20, 25, 30 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 đây
Luận án đã phân tích và tổng hợp 79 tài liệu tham khảo trong và ngoài nước, tập trung vào các nhánh nghiên cứu chính:
| Hướng nghiên cứu |
Tác giả, công trình tiêu biểu |
Phát hiện và phương pháp chính |
| Cấu trúc và chuyển pha của $\text{SiO}_2$ vô định hình và lỏng |
Sato và cộng sự (2010); Min Wu và cộng sự (2012); Mukherjee và cộng sự (2011) |
Sử dụng tán xạ tia X và mô phỏng chỉ ra sự tăng số phối trí của Si từ 4 lên 5 và 6 khi áp suất vượt 10–12 GPa; quan sát sự chuyển pha loại 1 từ pha mật độ thấp sang pha mật độ cao. |
| Hiện tượng đa thù hình và mô hình hai trạng thái |
Mishima và cộng sự (1984); Poole và cộng sự (1997); Saika-Voivod và cộng sự (2001, 2004) |
Phát hiện đa thù hình ở nước đá vô định hình; mở rộng giả thuyết mô hình hai trạng thái ("two-state", "two-domain") cho $\text{SiO}_2$, phân chia cấu trúc thành vùng mật độ thấp (mạng mở) và vùng mật độ cao (mạng nén). |
| Lý thuyết và mô hình động học chất lỏng quá nguội |
Adam & Gibbs (1965); Götze (lý thuyết Mode-Coupling); Cohen & Turnbull (thể tích tự do); Garrahan & Chandler (2002); Vogel và cộng sự (2004) |
Đưa ra khái niệm vùng tự sắp xếp (CRR); khảo sát tính không đồng nhất động học thông qua hàm tương quan đa điểm, chuyển động tập thể của các hạt linh động tạo thành đám kích thước lớn. |
| Khuếch tán dị thường và chuyển đổi đơn vị phối trí |
Horbach & Kob (1999); Shell, Debenedetti & Panagiotopoulos (2002); Angell và cộng sự |
Hệ số tự khuếch tán tăng khi tăng áp suất ở chất lỏng silicat; liên hệ khuếch tán với sự phá vỡ và tái lập liên kết thông qua cơ chế nhảy nguyên tử oxy giữa các quả cầu phối trí. |
| Cấu trúc $\text{MgSiO}_3$ vô định hình và lỏng |
Wilding và cộng sự (2004); Kohara và cộng sự (2004); Cormier & Cuello (2011); Oganov và cộng sự (2000); de Koker và cộng sự (2008) |
Thực nghiệm nhiễu xạ Neutron/tia X và mô phỏng lượng tử ab-initio; xác định khoảng cách liên kết $\text{Si-O} \approx 1.62\text{--}1.65\text{ \AA}$, $\text{Mg-O} \approx 2.0\text{--}2.1\text{ \AA}$; xuất hiện oxy không cầu (NBO) làm giảm mức độ polymer hóa. |
Khoảng trống nghiên cứu được lựa chọn giải quyết
Luận án tập trung giải quyết các hạn chế tồn đọng từ các nghiên cứu trước:
- Định lượng chi tiết cấu trúc hình học của 5 loại đơn vị phối trí cấu thành mạng $\text{SiO}_2$ lỏng ($\text{SiO}_4, \text{SiO}_5, \text{SiO}_6, \text{OSi}_2, \text{OSi}_3$), chứng minh rằng sự thay đổi mật độ vĩ mô là do thay đổi tỷ phần của các đơn vị phối trí chứ không phải do biến dạng hình học nội tại của từng đơn vị.
- Tách biệt và chứng minh vai trò của "chuyển đổi có ích" trong tổng số các biến cố chuyển đổi đơn vị phối trí, giải thích thỏa đáng cơ chế vi mô của hiện tượng khuếch tán dị thường và thuyên giảm động học.
- Đánh giá sự phân bố không gian và trật tự tầm trung của các cầu nối $\text{Si}_n\text{-O-Mg}_m$ trong $\text{MgSiO}_3$ lỏng ở điều kiện áp suất cao từ 0 đến 30 GPa, cung cấp bằng chứng định lượng về hiện tượng tách pha vi mô.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và khung phân tích
Luận án xây dựng khung phân tích dựa trên quan điểm: Cấu trúc của chất lỏng silicat nóng chảy là một mạng liên tục không gian được xây dựng từ các đơn vị phối trí (ĐVPT) cơ bản liên kết với nhau qua các nguyên tử ôxy dùng chung.
- Mô hình hai trạng thái (Two-state model): Hệ chất lỏng được xem như hỗn hợp nhiệt động của hai pha vi mô: pha mật độ thấp (Low-Density Liquid - LDL) đặc trưng bởi mạng tứ diện mở và pha mật độ cao (High-Density Liquid - HDL) đặc trưng bởi mạng đa diện nén chặt.
- Lý thuyết chuyển động tương quan và tạo đám: Động học vi mô được khảo sát dựa trên hành vi phân bố không gian của các nhóm nguyên tử có tính chất động học và cấu trúc đặc thù, theo dõi sự hình thành các đám (clusters) nguyên tử linh động hoặc kém linh động.
Phương pháp mô phỏng động lực học phân tử (ĐLHPT)
Mô phỏng sử dụng phương pháp ĐLHPT cổ điển, tích phân phương trình chuyển động Newton bằng thuật toán Verlet với điều kiện biên tuần hoàn 3 chiều nhằm triệt tiêu hiệu ứng bề mặt:
$$\mathbf{r}_i(t + \Delta t) = 2\mathbf{r}_i(t) - \mathbf{r}_i(t - \Delta t) + \frac{\mathbf{F}_i(t)}{m_i}\Delta t^2$$
Thế năng tương tác Cu-lông tầm xa được xử lý bằng kỹ thuật tính gần đúng Ewald-Hansen, phân tách thành thành phần không gian thực hội tụ nhanh, thành phần không gian ảo và số hạng hằng số:
$$U_{Coulomb} = U_{real} + U_{reciprocal} + U_{constant}$$
- Thế tương tác cho hệ $\text{SiO}_2$: Sử dụng hàm thế tương tác cặp Van Beets – Kramer – Van Santen (BKS), bắt nguồn từ tính toán ab-initio nguyên lý ban đầu:
$$U_{ij}(r_{ij}) = \frac{q_i q_j}{r_{ij}} + A_{ij}\exp(-B_{ij}r_{ij}) - \frac{C_{ij}}{r_{ij}^6}$$
Các thông số thế BKS:
-
$q_{\text{Si}} = +2.4e$, $q_{\text{O}} = -1.2e$.
-
Cặp $\text{O-O}$: $A_{ij} = 1388.7730\text{ eV}$, $B_{ij} = 2.76000\text{ \AA}^{-1}$, $C_{ij} = 175.0000\text{ eV}\cdot\text{\AA}^6$.
-
Cặp $\text{Si-O}$: $A_{ij} = 18003.7572\text{ eV}$, $B_{ij} = 4.87318\text{ \AA}^{-1}$, $C_{ij} = 133.5381\text{ eV}\cdot\text{\AA}^6$.
-
Cặp $\text{Si-Si}$: $A_{ij} = 0$, $B_{ij} = 0$, $C_{ij} = 0$.
-
Thế tương tác cho hệ $\text{MgSiO}_3$: Sử dụng hàm thế Oganov (OG), tích hợp các hiệu ứng thăng giáng nhiệt và chuyển động điểm không:
-
$q_{\text{Si}} = +2.9043e$, $q_{\text{Mg}} = +1.6049e$, $q_{\text{O}} = -1.5031e$.
-
Cặp $\text{Si-O}$: $A_{ij} = 1137.2\text{ eV}$, $B_{ij} = 4.2750\text{ \AA}^{-1}$, $C_{ij} = 0$.
-
Cặp $\text{Mg-O}$: $A_{ij} = 811.2\text{ eV}$, $B_{ij} = 3.6520\text{ \AA}^{-1}$, $C_{ij} = 0$.
-
Cặp $\text{O-O}$: $A_{ij} = 1018.6\text{ eV}$, $B_{ij} = 2.7600\text{ \AA}^{-1}$, $C_{ij} = 0$.
Quy trình tạo mẫu và kỹ thuật phân tích
- Quy trình tạo mẫu: Khởi tạo ngẫu nhiên vị trí các hạt $\rightarrow$ Nung nóng đẳng nhiệt ở 6000 K trong 50.000 bước $\rightarrow$ Làm nguội từng bước xuống 3500 K $\rightarrow$ Hồi phục hệ ở 3500 K và 0 GPa trong $10^6$ bước bằng mô hình NPT $\rightarrow$ Nén đẳng nhiệt tới các áp suất đích (0–30 GPa) $\rightarrow$ Cân bằng thống kê ở mô hình NVE trong $10^7$ bước. Bước thời gian mô phỏng là $\Delta t = 0.47\text{ fs}$ (hệ $\text{SiO}_2$) và $\Delta t = 0.40\text{ fs}$ (hệ $\text{MgSiO}_3$).
- Kỹ thuật phân tích vi cấu trúc:
- Hàm phân bố xuyên tâm cặp (PBXT) $g_{\alpha\beta}(r)$ và xác định bán kính ngắt $r_c$ tại cực tiểu thứ nhất ($r_c = 2.25\text{ \AA}$ cho $\text{Si-O}$ trong $\text{SiO}_2$; $r_c = 2.38\text{ \AA}$ cho $\text{Si-O}$ và $r_c = 3.32\text{ \AA}$ cho $\text{Mg-O}$ trong $\text{MgSiO}_3$).
- Tính toán số phối trí trung bình $Z_{\alpha\beta}$ và thống kê tỷ phần các đơn vị phối trí $\text{SiO}_x, \text{OSi}_y, \text{MgO}_y$.
- Tính toán hàm phân bố góc liên kết $P(\theta)$ cho các góc $\text{O-T-O}$ (trật tự gần) và $\text{T-O-T}$ (trật tự tầm trung, với $\text{T} = \text{Si, Mg}$).
- Kỹ thuật trực quan hóa đồ họa 3D để theo dõi sự phân bố không gian và kết cụm của các pha.
- Kỹ thuật phân tích động học:
- Hệ số tự khuếch tán $D$ xác định qua độ dịch chuyển bình phương trung bình (MSD) theo hệ thức Einstein:
$$D = \lim_{t \to \infty} \frac{\langle |\mathbf{r}_i(t) - \mathbf{r}_i(0)|^2 \rangle}{6t}$$
- Đếm số lượng biến cố chuyển đổi đơn vị phối trí theo thời gian và xác định tốc độ chuyển đổi $W$.
- Phân loại 5 nhóm nguyên tử đặc biệt: $A_1$ (ngẫu nhiên 50 Si, 100 O), $A_2$ (linh động nhất), $A_3$ (kém linh động nhất), $A_4$ (tần suất chuyển đổi cao nhất), $A_5$ (tần suất chuyển đổi thấp nhất).
- Đánh giá quá trình tạo đám thông qua số lượng đám $N_c$ và kích thước đám cực đại $S_{max}$.
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 1 hệ thống hóa các hiểu biết lý thuyết và thực nghiệm về cấu trúc, tính chất nhiệt động và động học của $\text{SiO}_2$ và $\text{MgSiO}_3$ ở các pha tinh thể, vô định hình và pha lỏng:
- Trình bày các dạng thù hình tinh thể của $\text{SiO}_2$ ($\alpha, \beta\text{-thạch anh}$, tridymite, cristobalite) và sự chuyển đổi trật tự tầm xa sang trật tự gần trong trạng thái vô định hình.
- Phân tích bản chất của nguyên tử ôxy cầu (BO - Bridging Oxygen) liên kết với 2 nguyên tử Si, hình thành mạng polymer hóa hoàn toàn, và nguyên tử ôxy không cầu (NBO - Non-Bridging Oxygen) xuất hiện khi pha tạp cation kim loại kiềm/kiềm thổ.
- Điểm luận hiện tượng đa thù hình ở chất lỏng tứ diện, các mô hình chuyển pha lỏng - lỏng loại 1, các đặc trưng dị thường (độ dốc âm của đường cong nóng chảy, hệ số giãn nở nhiệt âm, khuếch tán dị thường).
- Đánh giá các mô hình động học: lý thuyết Mode-Coupling, lý thuyết Adam-Gibbs về các vùng tự sắp xếp (CRR), lý thuyết thể tích tự do, mô hình thấm (percolation) và các kỹ thuật hàm tương quan 4 điểm $S_4(q,t)$.
- Tổng quan cấu trúc mạng $\text{MgSiO}_3$, sự phân bố các đơn vị $Q^n$ (với $n$ là số nguyên tử BO trong tứ diện $\text{SiO}_4$), vai trò phá vỡ mạng của ion $\text{Mg}^{2+}$ tạo ra các đơn vị $Q^2, Q^1, Q^3$, ôxy tự do và ảnh hưởng của mức độ polymer hóa lên độ nhớt của magma silicat.
Chương 2: Phương pháp tính toán
Chương 2 mô tả chi tiết phương pháp luận mô phỏng số và các giải thuật xử lý dữ liệu:
- Thiết lập phương trình chuyển động và thuật toán tích phân Verlet, kỹ thuật hiệu chỉnh vận tốc để kiểm soát nhiệt độ và biến đổi kích thước hộp mô phỏng để kiểm soát áp suất.
- Mô tả dạng toán học và các bộ tham số của thế tương tác BKS cho $\text{SiO}_2$ và thế tương tác OG cho $\text{MgSiO}_3$.
- Trình bày thuật toán xác định hàm PBXT $g(r)$, tích phân bán kính ngắt để tính số phối trí trung bình $Z$, và giải thuật xác định góc liên kết ba nguyên tử.
- Xây dựng thuật toán đếm chuyển đổi đơn vị phối trí $\text{SiO}x \rightarrow \text{SiO}{x'}$ và $\text{OSi}y \rightarrow \text{OSi}{y'}$.
- Thiết lập thuật toán nhận dạng đám (cluster analysis) dựa trên bán kính ngắt liên kết, đánh số thứ tự và gộp các nút mạng để tính toán kích thước đám $S$ và số đám $N_c$ theo thời gian mô phỏng.
Chương 3: Cấu trúc và động học của $\text{SiO}_2$ lỏng
Chương 3 trình bày các kết quả mô phỏng trên hệ $\text{SiO}_2$ lỏng:
- Đặc trưng cấu trúc theo nhiệt độ và áp suất: Vị trí đỉnh thứ nhất của hàm PBXT $r_{\text{Si-O}}$ tăng nhẹ từ $1.61\text{ \AA}$ lên $1.68\text{ \AA}$ khi áp suất tăng từ 0 lên 30 GPa (ở cả 3000 K và 3200 K), trong khi khoảng cách $r_{\text{O-O}}$ có xu hướng giảm.
| Đại lượng |
$P = 0\text{ GPa}$ |
$P = 10\text{ GPa}$ |
$P = 20\text{ GPa}$ |
$P = 30\text{ GPa}$ |
| Vị trí đỉnh $r_{\text{Si-O}}$ ($T = 3200\text{ K}$) |
$1.61\text{ \AA}$ |
$1.63\text{ \AA}$ |
$1.66\text{ \AA}$ |
$1.68\text{ \AA}$ |
| Tỷ phần $\text{SiO}_4$ |
$\sim 94%$ |
$\sim 48%$ |
$\sim 15%$ |
$\sim 3%$ |
| Tỷ phần $\text{SiO}_5$ |
$\sim 6%$ |
$\sim 44%$ (cực đại) |
$\sim 40%$ |
$\sim 22%$ |
| Tỷ phần $\text{SiO}_6$ |
$\sim 0%$ |
$\sim 8%$ |
$\sim 45%$ |
$\sim 75%$ |
| Tỷ phần $\text{OSi}_2$ |
$96%$ |
$\sim 70%$ |
$\sim 42%$ |
$24%$ |
| Tỷ phần $\text{OSi}_3$ |
$3%$ |
$\sim 28%$ |
$\sim 56%$ |
$68%$ |
- Bản chất hiện tượng đa thù hình: Luận án chứng minh mạng $\text{SiO}_2$ lỏng được cấu thành từ 5 đơn vị phối trí cơ bản. Cấu trúc hình học nội tại (góc liên kết $\text{O-Si-O}$, $\text{Si-O-Si}$ và khoảng cách liên kết) của từng loại ĐVPT riêng biệt không bị biến dạng khi tăng áp suất. Sự biến đổi mật độ toàn hệ thực chất là do sự chuyển dịch tỷ phần giữa hai pha:
- Pha mật độ thấp (LDL): Cấu thành từ các tứ diện $\text{SiO}_4$ kết nối qua các cầu nối $\text{OSi}_2$.
- Pha mật độ cao (HDL): Cấu thành từ các đa diện $\text{SiO}_5, \text{SiO}_6$ kết nối qua các liên kết $\text{OSi}_3$ (chia sẻ cạnh và mặt).
- Động học không đồng nhất: Khảo sát 5 nhóm nguyên tử $A_1 \div A_5$ cho thấy các nguyên tử linh động ($A_2$) và các nguyên tử có tần suất chuyển đổi cao ($A_4$) có xu hướng tự tập hợp thành các đám lớn có kích thước $S_{max}$ vượt trội so với nhóm ngẫu nhiên ($A_1$), chứng minh tính không đồng nhất về cấu trúc và động học trong không gian.
- Cơ chế khuếch tán: Khẳng định khuếch tán không đồng nhất với tổng số chuyển đổi đơn vị phối trí. Tác giả phân loại thành "chuyển đổi có ích" (chuyển đổi làm dịch chuyển khối tâm dẫn tới độ dời thực) và "chuyển đổi không có ích" (chuyển đổi dao động qua lại vị trí cũ). Tỷ số $\eta = M_{bt}/M_{bo}$ giữa chuyển đổi có ích và tổng số chuyển đổi quyết định hệ số khuếch tán $D$. Chuyển động tương quan tập thể của một nhóm nguyên tử là nguyên nhân trực tiếp dẫn tới hiện tượng khuếch tán dị thường và thuyên giảm động học.
Chương 4: Cấu trúc của $\text{MgSiO}_3$ lỏng
Chương 4 tập trung phân tích sự tiến triển cấu trúc của hệ $\text{MgSiO}_3$ lỏng ở 3500 K trong dải áp suất 0–30 GPa (các mẫu M1 đến M7):
- Môi trường địa phương của Si và Mg: Môi trường cục bộ xung quanh Si trong $\text{MgSiO}_3$ hoàn toàn tương đồng với trong $\text{SiO}_2$ lỏng. Ở 0 GPa, Si tồn tại chủ yếu ở dạng $\text{SiO}_4$; khi áp suất tăng lên 30 GPa, mạng $\text{Si-O}$ chuyển dịch từ cấu trúc tứ diện sang bát diện $\text{SiO}_6$ thông qua trạng thái trung gian $\text{SiO}_5$.
- Vai trò của cation $\text{Mg}^{2+}$: Các ion $\text{Mg}^{2+}$ định vị vào mạng làm bẻ gãy các liên kết $\text{Si-O-Si}$, phá vỡ mạng $\text{Si-O}$ liên tục thành các chuỗi và mạng con nhỏ hơn, làm xuất hiện các nguyên tử NBO và ôxy tự do.
- Sự tiến triển của các cầu nối theo áp suất:
- Ở áp suất thấp (0 GPa), nồng độ NBO và ôxy tự do chiếm tỷ lệ đáng kể.
- Khi áp suất tăng dần đến 30 GPa, nồng độ NBO và ôxy tự do giảm mạnh, trong khi tỷ phần nguyên tử ôxy cầu BO tăng lên rõ rệt. Đồng thời xuất hiện các liên kết cầu nối đa nhân $\text{OSi}_3$, $\text{OMg}_y$ và các cầu phức hợp $\text{Si}_n\text{-O-Mg}_m$ ($n, m = 1 \div 6$). Mức độ polymer hóa tổng thể của hệ tăng theo áp suất.
- Hiện tượng tách pha vi mô: Phân tích trực quan hóa không gian 3D và hàm phân bố các cầu nối khẳng định trong cấu trúc mạng $\text{MgSiO}_3$ lỏng luôn tồn tại sự phân tách không gian thành hai miền riêng biệt: các vùng giàu Si (tập trung các đám $\text{SiO}_x$) và các vùng giàu Mg (tập trung các đám $\text{MgO}_y$), minh chứng cho hiện tượng tách pha vi mô ở cấp độ phân tử.
Kết quả và những đóng góp mới
Đóng góp mới về lý luận và khoa học
- Làm sáng tỏ bản chất vi mô của hiện tượng đa thù hình: Luận án chỉ ra cấu trúc mạng của $\text{SiO}_2$ lỏng được hình thành từ 5 đơn vị phối trí cơ bản ($\text{SiO}_4, \text{SiO}_5, \text{SiO}_6, \text{OSi}_2, \text{OSi}_3$). Cấu trúc hình học của từng loại đơn vị phối trí không thay đổi theo áp suất; hiện tượng đa thù hình thực chất là sự thay đổi tỷ phần đóng góp giữa hai pha: 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$).
- Xác lập cơ chế vi mô của quá trình khuếch tán và động học không đồng nhất: Luận án chứng minh rằng khuếch tán trong chất lỏng có cấu trúc mạng không tỉ lệ thuận đơn thuần với tổng số chuyển đổi đơn vị phối trí, mà phụ thuộc quyết định vào số "chuyển đổi có ích". Hiện tượng khuếch tán dị thường và thuyên giảm động học được giải thích bởi sự xuất hiện của các chuyển động tương quan tập thể của các nhóm nguyên tử.
- Phát hiện tính tương đồng cấu trúc cục bộ và hiện tượng tách pha vi mô trong $\text{MgSiO}_3$ lỏng: Chứng minh môi trường địa phương của nguyên tử Si trong $\text{MgSiO}_3$ lỏng có sự tương đồng cao với trong $\text{SiO}_2$ lỏng khi bị nén ép. Dưới sự thay đổi của áp suất, vật liệu $\text{MgSiO}_3$ lỏng luôn biểu hiện hiện tượng tách pha vi mô thành các vùng giàu Si và các vùng giàu Mg đan xen trong mạng liên kết.
Ý nghĩa thực tiễn và định hướng ứng dụng
- Cung cấp các thông số nhiệt động, cấu trúc và hệ số truyền chất định lượng của silicat nóng chảy ở điều kiện sâu trong lòng đất, hỗ trợ các nhà địa vật lý xây dựng mô hình đối lưu manti và mô hình hóa quá trình tiến hóa magma.
- Cung cấp cơ sở khoa học vi mô cho công nghệ cố định chất thải phóng xạ: việc pha tạp các cation kim loại kiềm/kiềm thổ phá vỡ mạng $\text{SiO}_2$ tạo ra các vị trí mất cân bằng điện tích (NBO và ôxy tự do), giúp giam giữ các ion phóng xạ nguy hại trong cấu trúc thủy tinh hóa một cách bền vững.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của nghiên cứu
- Phạm vi khảo sát động học chi tiết (như đếm biến cố chuyển đổi ĐVPT, phân loại nhóm nguyên tử đặc biệt $A_1 \div A_5$, tính chuyển đổi có ích) mới được thực hiện tập trung trên hệ $\text{SiO}_2$ lỏng; hệ $\text{MgSiO}_3$ lỏng chủ yếu mới dừng lại ở phân tích đặc trưng vi cấu trúc tĩnh và trật tự tầm trung theo áp suất.
- Giới hạn kích thước mô hình mô phỏng (1998 nguyên tử cho $\text{SiO}_2$ và 5000 nguyên tử cho $\text{MgSiO}_3$) và thời gian hồi phục mô phỏng ($10^7$ bước) tuy đã khắc phục bằng điều kiện biên tuần hoàn nhưng vẫn chưa bao quát hết được các dao động thăng giáng cấu trúc ở thang độ dài lớn hơn (mesoscopic scale).
Hướng nghiên cứu tiếp
- Mở rộng thuật toán phân tích động học chuyển đổi có ích và chuyển động tương quan tập thể sang hệ ôxít ba thành phần $\text{MgSiO}_3$ lỏng và các hệ silicat đa thành phần khác ($\text{Na}_2\text{O-SiO}_2, \text{CaO-MgO-SiO}_2$).
- Khảo sát sâu hơn ảnh hưởng của trường nhiệt độ phối hợp với áp suất siêu cao (vượt trên 50 GPa) tương ứng với điều kiện manti dưới và lõi Trái Đất.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên khảo giá trị cho các đối tượng:
- Nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn, Vật lý tính toán, Khoa học vật liệu: Tham khảo phương pháp xây dựng mô hình ĐLHPT, kỹ thuật gần đúng Ewald-Hansen, thuật toán phân tích hàm PBXT, phân bố góc, phân tích đám và các kỹ thuật trực quan hóa dữ liệu mô phỏng.
- Giảng viên và nhà nghiên cứu Khoa học Trái Đất, Địa chất học: Sử dụng các dữ liệu về sự biến đổi pha silicat, mức độ polymer hóa, hành vi của các đơn vị phối trí Si và Mg dưới áp suất cao để giải thích tính chất lưu biến của magma trong vỏ và manti Trái Đất.
- Kỹ sư công nghệ vật liệu silicat và môi trường hạt nhân: Tham khảo cơ chế vi mô về sự phá vỡ mạng của ion kim loại và vai trò của ôxy không cầu (NBO) phục vụ phát triển vật liệu gốm chịu lửa và thủy tinh hóa chất thải phóng xạ.
Câu hỏi thường gặp
1. Mạng $\text{SiO}_2$ lỏng được cấu thành từ những đơn vị phối trí (ĐVPT) nào và tỷ phần của chúng biến đổi ra sao khi nén từ 0 đến 30 GPa?
Mạng $\text{SiO}_2$ lỏng được cấu thành từ 5 ĐVPT chính gồm: 3 ĐVPT đặc trưng cho trật tự gần là $\text{SiO}_4, \text{SiO}_5, \text{SiO}_6$ và 2 ĐVPT đặc trưng cho trật tự tầm trung là $\text{OSi}_2, \text{OSi}_3$. Khi áp suất tăng từ 0 lên 30 GPa (ở 3200 K), tỷ phần tứ diện $\text{SiO}_4$ giảm mạnh từ 94% xuống 3%, ĐVPT trung gian $\text{SiO}_5$ tăng và đạt cực đại tại khoảng 10–15 GPa ($\sim 44%$) rồi giảm về 22%, trong khi bát diện $\text{SiO}_6$ tăng từ 0% lên chiếm ưu thế tuyệt đối với 75%. Tương ứng, tỷ phần $\text{OSi}_2$ giảm từ 96% xuống 24% và $\text{OSi}_3$ tăng từ 3% lên 68%.
2. Luận án định nghĩa và chứng minh hiện tượng đa thù hình trong $\text{SiO}_2$ lỏng theo mô hình hai trạng thái như thế nào?
Luận án chỉ ra rằng tại mọi áp suất, $\text{SiO}_2$ lỏng luôn tồn tại đồng thời hai pha: pha mật độ thấp (cấu tạo từ các tứ diện $\text{SiO}_4$ liên kết qua $\text{OSi}_2$) và pha mật độ cao (cấu tạo từ các đa diện $\text{SiO}_5, \text{SiO}_6$ liên kết qua $\text{OSi}_3$). Cấu trúc hình học (khoảng cách liên kết và góc liên kết) của từng loại ĐVPT riêng biệt không đổi theo áp suất. Sự biến thiên mật độ vĩ mô khi nén ép thực chất là sự thay đổi tỷ phần đóng góp giữa hai pha này.
3. "Chuyển đổi có ích" giữa các ĐVPT là gì và liên hệ thế nào với hiện tượng khuếch tán trong $\text{SiO}_2$ lỏng?
Chuyển đổi có ích là các biến cố chuyển đổi giữa các ĐVPT ($\text{SiO}x \rightarrow \text{SiO}{x'}$) kèm theo sự dịch chuyển vị trí thực sự của khối tâm nguyên tử dẫn tới độ dời không gian, phân biệt với các chuyển đổi không có ích (chỉ là dao động liên kết quanh vị trí cân bằng cũ rồi hoàn nguyên). Luận án chứng minh hệ số khuếch tán $D$ không phụ thuộc vào tổng số chuyển đổi đơn thuần mà phụ thuộc vào tốc độ và tỷ số chuyển đổi có ích $\eta = M_{bt}/M_{bo}$.
4. Hiện tượng tách pha vi mô trong $\text{MgSiO}_3$ lỏng biểu hiện như thế nào qua kết quả mô phỏng?
Trong $\text{MgSiO}_3$ lỏng ở 3500 K, các ion $\text{Mg}^{2+}$ định vị vào mạng làm đứt gãy mạng $\text{Si-O}$, tạo ra các nguyên tử NBO và ôxy tự do. Phân tích hàm phân bố cầu nối $\text{Si}_n\text{-O-Mg}_m$ kết hợp trực quan hóa không gian 3D chứng minh các nguyên tử Mg và Si không phân bố trộn lẫn đồng nhất mà luôn tự tập hợp thành các miền cấu trúc riêng biệt: các vùng giàu Si (chứa các chuỗi/mạng con $\text{SiO}_x$) và các vùng giàu Mg (chứa các đám $\text{MgO}_y$), tạo nên hiện tượng tách pha vi mô.
Kết luận
Luận án của tác giả Luyện Thị San đã ứng dụng thành công phương pháp mô phỏng động lực học phân tử với các thế tương tác BKS và OG để khảo sát toàn diện cấu trúc vi mô và động học của hai hệ chất lỏng mạng $\text{SiO}_2$ và $\text{MgSiO}_3$ trong dải áp suất 0–30 GPa. Công trình đã làm sáng tỏ bản chất của hiện tượng đa thù hình thông qua 5 đơn vị phối trí và mô hình hai trạng thái, đồng thời chứng minh vai trò quyết định của chuyển đổi có ích và chuyển động tương quan tập thể đối với hiện tượng khuếch tán dị thường trong $\text{SiO}_2$ lỏng. Đối với hệ $\text{MgSiO}_3$ lỏng, luận án đã xác lập tính tương đồng của môi trường địa phương xung quanh nguyên tử silíc so với $\text{SiO}_2$, sự tiến triển của mức độ polymer hóa và cung cấp bằng chứng định lượng về hiện tượng tách pha vi mô giữa các miền giàu Si và giàu Mg. Các kết quả nghiên cứu đóng góp nền tảng lý thuyết quan trọng cho ngành vật lý chất lỏng cấu trúc mạng và cung cấp số liệu tham chiếu giá trị cho các nghiên cứu địa vật lý manti cũng như công nghệ xử lý rác thải hạt nhân.
