Tổng quan về luận án

Nghiên cứu về trạng thái lỏng và thủy tinh của các oxit mạng đóng vai trò nền tảng trong vật lý vật chất ngưng tụ, khoa học vật liệu chịu lửa và địa vật lý manti Trái Đất. Luận án tiến sĩ vật lý kỹ thuật của nghiên cứu sinh Luyện Thị San với đề tài "Nghiên cứu đặc trưng cấu trúc và động học của các vật liệu SiO2 và MgSiO3 lỏng" (chuyên ngành Vật lý Kỹ thuật, mã số 62520401, thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Văn Hồng và GS. Vũ Văn Hùng tại Trường Đại học Bách khoa Hà Nội, 2017) đã giải quyết những bài toán cốt lõi về bản chất của cấu trúc mạng vô định hình, hiện tượng đa thù hình (polyamorphism), động học không đồng nhất (dynamic heterogeneity) và cơ chế khuếch tán dị thường dưới điều kiện nhiệt độ và áp suất cực hạn.

Trong bối cảnh địa khoa học, $\text{SiO}_2$ và $\text{MgSiO}_3$ (magnesium silicate/enstatite) là hai thành phần khoáng vật chiếm hơn 80% thể tích vỏ và lớp phủ Trái Đất. Dưới điều kiện sâu trong lòng đất, nơi áp suất dao động từ vài GPa đến hàng trăm GPa và nhiệt độ lên tới hàng nghìn Kelvin, các hợp chất này tồn tại ở dạng magma silicat nóng chảy. Tuy nhiên, việc khảo sát thực nghiệm ở điều kiện nhiệt độ – áp suất siêu cao (P-T cao) gặp phải rào cản kỹ thuật nghiêm trọng. Mặc dù các kỹ thuật nhiễu xạ tia X, tán xạ neutron và cộng hưởng từ hạt nhân (NMR) đã cung cấp nhiều dữ liệu tĩnh, bức tranh vi mô về mối liên kết giữa cấu trúc địa phương (local structure) và hành vi chuyển dời nguyên tử (atomic transport) vẫn tồn tại nhiều khoảng trống lý thuyết (research gaps) chưa được giải quyết dứt điểm:

  1. Khoảng trống về cơ chế đa thù hình: Liệu quá trình chuyển đổi mật độ trong chất lỏng $\text{SiO}_2$ tuân theo mô hình hai trạng thái (two-state model) với sự phân tách không gian rõ ràng giữa pha mật độ thấp (Low-Density Liquid - LDL) và pha mật độ cao (High-Density Liquid - HDL), hay chỉ là sự biến đổi nén hình học đồng nhất liên tục?
  2. Khoảng trống về vi cơ chế khuếch tán và động học không đồng nhất: Mối liên hệ bản chất giữa sự phá vỡ – tái tạo đơn vị phối trí (coordination units) và hiện tượng khuếch tán dị thường (anomalous diffusion - hệ số khuếch tán tăng khi tăng áp suất) cùng hiện tượng thuyên giảm động học (dynamical deceleration/glass transition) là gì?
  3. Khoảng trống về cấu trúc silicate đa thành phần: Ảnh hưởng của cation mạng $\text{Mg}^{2+}$ lên mức độ polyme hóa, sự hình thành các tâm oxy cầu (BO), oxy không cầu (NBO), oxy tự do và hiện tượng tách pha vi mô (microscopic phase separation) trong $\text{MgSiO}_3$ lỏng ở áp suất cao diễn ra như thế nào?

Để trả lời các khoảng trống khoa học trên, hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  • Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc không gian của $\text{SiO}_2$ lỏng được tổ chức bởi những đơn vị phối trí cụ thể nào và các đơn vị này liên kết với nhau ra sao để hình thành các pha mật độ khác biệt?
  • Câu hỏi nghiên cứu 2 (RQ2): Quá trình chuyển đổi phối trí giữa các nguyên tử trung tâm đóng góp như thế nào vào hệ số khuếch tán vĩ mô, và chuyển động tương quan tập thể có vai trò gì trong động học không đồng nhất?
  • Câu hỏi nghiên cứu 3 (RQ3): Sự pha tạp cation $\text{Mg}^{2+}$ làm biến đổi mạng $\text{Si-O}$ theo cơ chế nào và hiện tượng tách pha vi mô giữa miền giàu $\text{Si}$ và miền giàu $\text{Mg}$ biểu hiện ra sao trong dải áp suất từ 0 đến 30 GPa?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Mạng $\text{SiO}_2$ lỏng bao gồm 5 đơn vị phối trí cơ bản; pha mật độ thấp gồm các tứ diện $\text{SiO}_4$ liên kết qua cầu $\text{OSi}_2$, trong khi pha mật độ cao hình thành từ các đa diện $\text{SiO}_5, \text{SiO}_6$ liên kết qua các tâm $\text{OSi}_3$.
  • Giả thuyết 2 (H2): Khuếch tán trong chất lỏng mạng không đồng nhất với mọi quá trình biến đổi phối trí; chỉ các chuyển đổi "có ích" (useful transitions) kèm theo chuyển động tương quan nhóm mới dẫn đến sự dịch chuyển thực, chi phối hiện tượng khuếch tán dị thường.
  • Giả thuyết 3 (H3): Sự xuất hiện của $\text{Mg}^{2+}$ bẻ gãy mạng $\text{Si-O}$ tạo ra $\text{NBO}$ và oxy tự do nhưng không làm thay đổi đặc trưng hình học nội tại của các đơn vị $\text{SiO}_x$, dẫn tới hiện tượng tách pha vi mô liên tục ở cấp độ nano.

Khung lý thuyết của công trình tích hợp mô hình nhiệt động hai trạng thái (Two-State Thermodynamics của Ponyatovsky và Tanaka), lý thuyết giãn tích tập thể Adam-Gibbs, lý thuyết chuyển pha thủy tinh Mode-Coupling (MCT) và mô hình chuyển dời phối trí trong mạng ngẫu nhiên liên tục (Continuous Random Network).

Phạm vi nghiên cứu bao quát các mẫu mô phỏng động lực học phân tử (Molecular Dynamics - MD) với quy mô 1.998 nguyên tử cho hệ $\text{SiO}_2$ (nhiệt độ $T = 2600 \div 3500\text{ K}$, áp suất $P = 0 \div 30\text{ GPa}$) và 5.000 nguyên tử cho hệ $\text{MgSiO}_3$ ($T = 3500\text{ K}$, $P = 0 \div 30\text{ GPa}$). Đóng góp mang tính đột phá của luận án là đã lượng hóa thành công tỷ lệ chuyển đổi hữu ích, xác lập bản đồ phân bố pha vi mô và chứng minh bằng thực chứng số sự tồn tại song song của các đám nguyên tử linh động/kém linh động, giải quyết dứt điểm các tranh luận kéo dài về cơ chế vi mô của silicat lỏng.

Literature Review và Positioning

Nghiên cứu về vật liệu silicat lỏng và thủy tinh đã trải qua nhiều thập kỷ phát triển với các trường phái lý thuyết và thực nghiệm đối thoại sâu sắc. Khởi nguồn từ phát hiện chấn động của Mishima và cộng sự (1984, 1985) về hiện tượng đa thù hình trong nước đá vô định hình với sự chuyển pha bậc một giữa pha mật độ thấp (LDA) và pha mật độ cao (HDA), các nhà vật lý ngưng tụ như Poole, Sciortino, Essmann và Stanley (1992) đã mở rộng giả thuyết chuyển pha lỏng-lỏng (Liquid-Liquid Phase Transition - LLPT) sang các hệ có mạng tứ diện tương tự như $\text{SiO}_2$, $\text{GeO}_2$ và $\text{Si}$.

Đối với hệ $\text{SiO}_2$, công trình của van Beest, Kramer và van Santen (1990) đã thiết lập hàm thế tương tác BKS kinh điển, mở đường cho việc mô phỏng chính xác cấu trúc vi mô của thạch anh, cristobalite và thủy tinh silica. Tuy nhiên, các tranh luận học thuật lớn đã nổ ra xoay quanh hai trường phái:

Tranh luận 1: Bản chất của quá trình nén cấu trúc vô định hình. Một trường phái do Mukherjee và cộng sự (2001) cùng Saika-Voivod và cộng sự (2001, 2004) dẫn dắt cho rằng quá trình nén của $\text{SiO}_2$ vô định hình và lỏng là một chuyển pha bậc một thực sự, kèm theo sự thay đổi nhảy bậc về thể tích và cấu hình vi mô giữa mạng tứ diện $\text{SiO}_4$ và mạng bát diện $\text{SiO}_6$. Ngược lại, trường phái thực nghiệm tán xạ tia X áp suất cao của Sato và cộng sự (2010) cùng Min Wu và cộng sự (2012) lại chỉ ra rằng sự thay đổi số phối trí của $\text{Si}$ từ 4 lên 6 diễn ra từ từ và liên tục trong dải áp suất 10 đến 45 GPa, với sự xuất hiện của các đơn vị trung gian $\text{SiO}_5$ đóng vai trò khuyết tật phối trí, phản bác sự tồn tại của ranh giới chuyển pha bậc một sắc nét ở pha lỏng.

Tranh luận 2: Nguồn gốc của động học không đồng nhất và khuếch tán dị thường. Trường phái động học thuần túy của Garrahan và Chandler (2002) cùng Vogel và cộng sự (2004) lập luận rằng động học không đồng nhất xuất phát từ tính linh động hỗ trợ (dynamic facilitation), trong đó các hạt linh động kích hoạt tính linh động của các hạt lân cận trong không gian mà không nhất thiết phụ thuộc vào một trật tự cấu trúc tĩnh rõ ràng. Ngược lại, trường phái cấu trúc địa phương của Tanaka và cộng sự (2010), Widmer-Cooper và Harrowell (2006) khẳng định rằng cấu trúc tĩnh ("glass order" hay trật tự cục bộ) chi phối hoàn toàn thời gian hồi phục và các vùng linh động nhanh/chậm. Trong khi đó, lý thuyết Adam-Gibbs (1965) dựa trên khái niệm các vùng tự sắp xếp tập thể (Cooperatively Rearranging Regions - CRR) vẫn gặp khó khăn trong việc định nghĩa và lượng hóa hình học chính xác của các CRR này ở cấp độ nguyên tử.

So sánh với các nghiên cứu quốc tế trên hệ $\text{MgSiO}_3$: Trong hệ $\text{MgSiO}3$, các phép đo tán xạ neutron và tia X của Wilding và cộng sự (2004) cùng các mô phỏng Monte Carlo đảo (RMC) của Martin, Cormier và cộng sự (2001) cho thấy ở 0 GPa, $\text{Si}$ duy trì cấu trúc tứ diện $\text{SiO}4$ ($r{\text{Si-O}} \approx 1.64\text{ \AA}$) với số phối trí oxy quanh $\text{Mg}$ đạt xấp xỉ $Z{\text{Mg-O}} = 4.5 \div 5.0$. Tuy nhiên, mô hình cấu trúc của Kubicki (2001) lại đưa ra giả định rằng $\text{Mg}$ có 4 nguyên tử oxy ở khoảng cách 2.0 Å và 2 nguyên tử oxy ở khoảng cách xa hơn 2.5 Å. Các mô phỏng ab-initio của Stixrude và Karki (2005) cùng mô phỏng cổ điển của Oganov và cộng sự (2000, 2001) với thế OG đã mô tả sự tăng số phối trí của cả $\text{Si}$ và $\text{Mg}$ khi áp suất vượt 50 GPa, nhưng chưa làm rõ cơ chế tổ chức của các chuỗi/mạng con (sub-networks) $\text{Si-O}$ và $\text{Mg-O}$ cũng như sự tiến hóa của oxy cầu và oxy không cầu trong không gian 3 chiều ở cấp độ trung gian (medium-range order).

Luận án của Luyện Thị San định vị chính xác vào giao điểm của các tranh luận này. Bằng cách phát triển phương pháp phân tích 5 đơn vị phối trí ($\text{SiO}_4, \text{SiO}_5, \text{SiO}_6, \text{OSi}_2, \text{OSi}_3$) kết hợp kỹ thuật truy vết cụm nguyên tử động học và chia tách chuyển đổi có ích/vô ích, công trình đã hòa giải thành công quan điểm chuyển pha hai trạng thái với bức tranh biến đổi cấu trúc liên tục, đồng thời cung cấp bằng chứng thực chứng số trực tiếp cho mối liên hệ nhân quả giữa cấu trúc tĩnh và động học không đồng nhấ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 và hoàn thiện các lý thuyết vật lý ngưng tụ hiện đại thông qua các đóng góp cụ thể:

  1. Mở rộng mô hình hai trạng thái (Two-State Model): Công trình chứng minh rằng trong $\text{SiO}_2$ lỏng, tính đa thù hình không phải là một giả định trừu tượng mà được cấu thành cụ thể từ 5 đơn vị phối trí cơ bản. Mạng lỏng là sự cùng tồn tại không gian của hai pha: pha mật độ thấp (cấu tạo từ các tứ diện $\text{SiO}_4$ liên kết với nhau qua các nút cầu $\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 các nút cầu ba $\text{OSi}_3$). Tỷ phần của hai pha này biến thiên liên tục theo áp suất, làm sáng tỏ bản chất vi mô của các tính chất nhiệt động dị thường.
  2. Tái định nghĩa cơ chế khuếch tán trong chất lỏng mạng: Luận án bác bỏ quan điểm đơn giản cho rằng mọi sự đứt gãy hoặc hình thành liên kết phối trí đều đóng góp vào khuếch tán. Nghiên cứu chỉ ra rằng các quá trình chuyển đổi phối trí gồm hai loại: chuyển đổi không có ích (chỉ là dao động nhiệt đàn hồi xung quanh vị trí cân bằng, nguyên tử quay trở lại cấu hình cũ) và chuyển đổi có ích (nguyên tử nhảy thoát khỏi lồng phối trí và tạo liên kết bền vững mới). Chỉ những chuyển đổi có ích mới đóng góp trực tiếp vào hệ số khuếch tán $D$.
  3. Hiện thực hóa hình học của lý thuyết Adam-Gibbs: Bằng việc phân tích 5 nhóm nguyên tử đặc biệt ($A_1 \div A_5$), luận án đã lượng hóa kích thước và hình thái không gian của các vùng tự sắp xếp tập thể (CRR), chứng minh sự tồn tại của chuyển động tương quan nhóm (correlated group displacement), nơi một cụm nguyên tử cùng phối hợp dịch chuyển qua khoảng không gian trung gian.

Các mệnh đề lý thuyết cốt lõi được kiểm chứng:

  • Mệnh đề 1: $D \propto \xi \cdot \nu_{\text{tran}}$, trong đó $\xi = M_{bt}/M_{bo}$ là tỷ số giữa chuyển đổi có ích ($M_{bt}$) và tổng số chuyển đổi ($M_{bo}$), và $\nu_{\text{tran}}$ là tần suất chuyển đổi phối trí toàn phần.
  • Mệnh đề 2: Động học không đồng nhất là hệ quả tất yếu của cấu trúc không gian không đồng nhất; các vùng giàu $\text{SiO}_4$ sở hữu độ linh động cao hơn các vùng tập trung $\text{SiO}_5, \text{SiO}_6$ bị nén chặt.
  • Mệnh đề 3: Trong silicat đa cấu tử $\text{MgSiO}_3$, sự suy giảm mức độ polyme hóa diễn ra cục bộ thông qua sự phân bố bất đối xứng của các cầu nối $\text{Si}_n\text{-O-Mg}_m$, duy trì trật tự gần của $\text{SiO}_x$ độc lập với nồng độ $\text{Mg}$.
+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH ĐA CẤU TỬ VÀ ĐỘNG HỌC                  |
+-------------------------------------------------------------------------+
                                     |
      +------------------------------+------------------------------+
      |                                                             |
      v                                                             v
[Pha Mật Độ Thấp (LDL)]                                   [Pha Mật Độ Cao (HDL)]
- Đơn vị: SiO4 + OSi2                                     - Đơn vị: SiO5, SiO6 + OSi3
- Cấu trúc mở, linh động cao                              - Mạng nén chặt, độ nhớt cao
- Liên kết đỉnh chia sẻ (BO)                              - Chia sẻ cạnh (Ne) và mặt (Nf)
      |                                                             |
      +------------------------------+------------------------------+
                                     |
                                     v
                 [ĐỘNG HỌC VÀ CƠ CHẾ CHUYỂN DỜI VI MÔ]
                 - Tần suất chuyển đổi toàn phần: v_tran
                 - Chuyển đổi vô ích: M_bo (Dao động cục bộ)
                 - Chuyển đổi có ích: M_{bt} (Tạo liên kết mới)
                 - Tỷ số hiệu dụng: xi = M_{bt}/M_{bo}
                 - Dịch chuyển tương quan nhóm (Correlated Multi-atom Jump)
                                     |
                                     v
       [HỆ SỐ KHUẾCH TÁN DỊ THƯỜNG & TÁCH PHA VI MÔ (MgSiO3)]
       - SiO2: D_max tại áp suất chuyển đổi topo mạng (~10-15 GPa)
       - MgSiO3: Mạng con Si-O đan xen Mạng con Mg-O (Vùng giàu Si / Vùng giàu Mg)

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

Khung phân tích của luận án tích hợp ba trụ cột phương pháp luận:

  1. Phương pháp phân rã đơn vị phối trí (Coordination Unit Decomposition): Không sử dụng giá trị số phối trí trung bình đơn thuần (vốn làm nhòe đi các thông tin cục bộ), công trình phân rã toàn bộ mạng thành các đơn vị thành phần $\text{SiO}_x$ ($x = 4, 5, 6$) và $\text{OSi}_y$ ($y = 2, 3$), kết hợp phân tích cách thức kết nối qua nguyên tử oxy chung (dùng chung đỉnh, cạnh $N_e$, hay mặt $N_f$).
  2. Thuật toán truy vết cụm động học (Sub-ensemble Clustering Algorithm): Thiết lập 5 tập hợp nguyên tử mục tiêu ($A_1$: ngẫu nhiên; $A_2$: linh động nhất; $A_3$: kém linh động nhất; $A_4$: tần suất chuyển đổi cao; $A_5$: tần suất chuyển đổi thấp), mỗi tập hợp gồm đúng 50 nguyên tử $\text{Si}$ và 100 nguyên tử $\text{O}$. Bằng cách theo dõi hàm phân bố kích thước đám cực đại $N_{a_Ai}(t)$ và số lượng đám $N_{c_Ai}(t)$, luận án định lượng hóa được mức độ tụ hợp không gian của các pha động học.
  3. Điều kiện biên và phạm vi hợp thức: Khung phân tích được xác lập trong giới hạn cơ học thống kê cổ điển, áp dụng cho trạng thái cân bằng nhiệt động ở nhiệt độ $T \ge 2600\text{ K}$ và áp suất $P \le 30\text{ GPa}$, nơi các hiệu ứng lượng tử của hạt nhân có thể bỏ qua và thế tương tác ion hiệu dụng phản ánh chính xác trường lực tĩnh điện.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng khoa học (scientific positivism) và tất định luận vi mô của cơ học thống kê cổ điển. Mô phỏng động lực học phân tử cổ điển được lựa chọn vì khả năng cung cấp quỹ đạo không gian – thời gian chính xác của từng nguyên tử, vượt qua các giới hạn đo đạc của thực nghiệm ở dải áp suất trên 10 GPa và nhiệt độ trên 3000 K.

Thiết kế mô hình đa mức (multi-level design) bao gồm:

  • Mức độ hạt đơn lẻ: Vị trí $\vec{r}_i(t)$, vận tốc $\vec{v}_i(t)$, lực tương tác $\vec{F}_i(t)$.
  • Mức độ cấu trúc địa phương (Short-range order): Hàm phân bố xuyên tâm cặp $g_{\alpha\beta}(r)$, góc liên kết $\theta_{\text{O-T-O}}$, số phối trí $Z_{\alpha\beta}$.
  • Mức độ cấu trúc trung gian (Medium-range order): Góc liên kết cầu $\theta_{\text{T-O-T}}$, liên kết chia sẻ cạnh $N_e$, chia sẻ mặt $N_f$, các cầu nối $\text{Si}_n\text{-O-Mg}_m$, cấu trúc vòng và chuỗi.
  • Mức độ vĩ mô và động học: Hệ số tự khuếch tán $D$, độ dịch chuyển bình phương trung bình (Mean Squared Displacement - MSD), mật độ khối $\rho(P, T)$.

Quy mô mẫu mô phỏng chính xác:

  • Hệ $\text{SiO}_2$: 1.998 nguyên tử (gồm 666 nguyên tử $\text{Si}$ và 1.332 nguyên tử $\text{O}$) đặt trong hộp mô phỏng lập phương lập thể với điều kiện biên tuần hoàn 3 chiều (Periodic Boundary Conditions - PBC).
  • Hệ $\text{MgSiO}_3$: 5.000 nguyên tử (gồm 1.000 nguyên tử $\text{Si}$, 1.000 nguyên tử $\text{Mg}$ và 3.000 nguyên tử $\text{O}$) sử dụng điều kiện biên tuần hoàn.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng được thực hiện qua các giai đoạn nghiêm ngặt nhằm triệt tiêu hoàn toàn hiệu ứng bộ nhớ của cấu hình ban đầu:

+-------------------------------------------------------------------------+
|                    QUY TRÌNH MÔ PHỎNG ĐỘNG LỰC HỌC                      |
+-------------------------------------------------------------------------+
                                     |
                                     v
 [Khởi tạo ngẫu nhiên 1.998 / 5.000 hạt trong ô lập phương PBC]
                                     |
                                     v
 [Nung nóng chảy ở T = 6000 K trong 50.000 bước (Xóa hoàn toàn trật tự cũ)]
                                     |
                                     v
 [Làm lạnh phân đoạn: 5500 K -> 5000 K -> 4500 K -> 4000 K -> 3500 K]
                                     |
                                     v
 [Hồi phục đẳng nhiệt - đẳng áp (NPT) trong 10^6 bước ở P = 0 GPa]
                                     |
                                     v
 [Nén đẳng nhiệt NPT tới các áp suất đích: 5, 10, 15, 20, 25, 30 GPa]
                                     |
                                     v
 [Cân bằng vi chính tắc (NVE) trong 10^7 bước mô phỏng (Thu thập số liệu)]
                                     |
                                     v
 [Tính toán thống kê trung bình trên 1.000 cấu hình cân bằng cuối cùng]

Phương trình chuyển động Newton được tích phân bằng thuật toán Verlet với bước nhảy thời gian cực ngắn $\Delta t = 0.47\text{ fs}$ ($0.47 \times 10^{-15}\text{ s}$), đảm bảo sự bảo toàn năng lượng tuyệt đối trong hệ vi chính tắc NVE với độ trôi dạt năng lượng nhỏ hơn $10^{-5}$.

Tương tác tĩnh điện tầm xa Coulomb chiếm phần lớn chi phí tính toán được giải quyết chính xác bằng kỹ thuật tính tổng Ewald-Hansen, phân tách thế Coulomb thành hai thành phần hội tụ nhanh trong không gian thực ($U_r$) và không gian tương hỗ Fourier ($U_k$), cùng số hạng tự tương tác ($U_0$): $$U_{\text{Coulomb}} = U_r + U_k + U_0$$

Trường lực tương tác (Interatomic Potentials):

  • Đối với hệ $\text{SiO}_2$: Sử dụng thế tương tác cặp BKS (van Beest - Kramer - van Santen): $$U(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 điện tích hiệu dụng $q_{\text{Si}} = +2.4e$, $q_{\text{O}} = -1.2e$; các tham số $A_{\text{Si-O}} = 18003.7572\text{ eV}$, $B_{\text{Si-O}} = 4.87318\text{ \AA}^{-1}$, $C_{\text{Si-O}} = 133.5381\text{ eV\cdot\AA}^6$; $A_{\text{O-O}} = 1388.7730\text{ eV}$, $B_{\text{O-O}} = 2.76000\text{ \AA}^{-1}$, $C_{\text{O-O}} = 175.0000\text{ eV\cdot\AA}^6$.
  • Đối với hệ $\text{MgSiO}_3$: Sử dụng thế tương tác OG (Oganov et al.), bổ sung các hiệu ứng dao động không điểm và thăng giáng nhiệt với $q_{\text{Si}} = +2.9043e$, $q_{\text{Mg}} = +1.6049e$, $q_{\text{O}} = -1.5031e$.

Độ tin cậy (reliability) và tính hợp lệ (validity) được kiểm chứng chéo thông qua việc so sánh trực tiếp các hàm phân bố xuyên tâm $g(r)$ và khoảng cách liên kết cân bằng với các dữ liệu thực nghiệm tán xạ tia X của Sato et al. (2010), tán xạ neutron của Wilding et al. (2004) và cộng hưởng từ hạt nhân NMR của Stebbins et al.

Data và phân tích

Bán kính ngắt phối trí ($r_c$) được xác định từ vị trí cực tiểu thứ nhất của hàm phân bố xuyên tâm cặp $g_{\alpha\beta}(r)$:

  • Với hệ $\text{SiO}2$: $r{c(\text{Si-O})} = 2.45\text{ \AA}$.
  • Với hệ $\text{MgSiO}3$: $r{c(\text{Si-O})} = 2.38\text{ \AA}$, $r_{c(\text{Mg-O})} = 3.10\text{ \AA}$, $r_{c(\text{Si-Si})} = 3.80\text{ \AA}$, $r_{c(\text{Mg-Mg})} = 4.20\text{ \AA}$.

Số phối trí trung bình $Z_{\alpha\beta}$ được tính bằng tích phân: $$Z_{\alpha\beta} = 4\pi \rho_\beta \int_0^{r_c} g_{\alpha\beta}(r) r^2 dr$$

Hệ số tự khuếch tán $D$ được trích xuất từ hệ thức Einstein thông qua độ dốc tuyến tính của MSD theo thời gian ở giới hạn khuếch tán dài ($t \to \infty$): $$D = \lim_{t \to \infty} \frac{\langle |\vec{r}_i(t) - \vec{r}_i(0)|^2 \rangle}{6t}$$

Để đảm bảo tính vững chắc thống kê (statistical robustness), mọi giá trị cấu trúc và động học đều được lấy trung bình tập hợp (ensemble average) từ 1.000 cấu hình độc lập, ghi nhận sau mỗi 10 bước mô phỏng trong pha cân bằng NVE kéo dài $10^7$ bước.

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

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

Luận án đã đạt được 4 phát hiện mang tính đột phá về vật lý cấu trúc và động học silicat lỏng:

1. Minh chứng cấu trúc vi mô của mô hình hai trạng thái trong $\text{SiO}_2$ lỏng: Nghiên cứu chỉ ra rằng ở mọi áp suất trong dải $0 \div 30\text{ GPa}$, $\text{SiO}_2$ lỏng luôn tồn tại đồng thời hai pha: pha mật độ thấp (cấu thành bởi $\text{SiO}_4$ và $\text{OSi}_2$) và pha mật độ cao (cấu thành bởi $\text{SiO}_5, \text{SiO}_6$ và $\text{OSi}_3$).

  • Ở $P = 0\text{ GPa}$, mạng lưới bị chiếm ưu thế bởi tứ diện $\text{SiO}4$ ($>90%$), khoảng cách liên kết trung bình $r{\text{Si-O}} = 1.61\text{ \AA}$, góc liên kết $\theta_{\text{O-Si-O}} \approx 109.5^\circ$, liên kết chủ yếu qua đỉnh dùng chung một oxy cầu $\text{OSi}2$ ($\theta{\text{Si-O-Si}} \approx 144^\circ \div 152^\circ$).
  • Khi áp suất tăng lên $10 \div 15\text{ GPa}$, tỷ phần $\text{SiO}_5$ đạt cực đại đóng vai trò pha chuyển tiếp.
  • Ở $P = 30\text{ GPa}$, cấu trúc chuyển hẳn sang pha mật độ cao với sự chiếm ưu thế của bát diện $\text{SiO}6$, khoảng cách $r{\text{Si-O}}$ tăng nhẹ lên $1.68 \div 1.70\text{ \AA}$, xuất hiện các liên kết dùng chung cạnh ($N_e$) và dùng chung mặt ($N_f$) kết nối thông qua oxy cầu ba $\text{OSi}_3$. Sự nén khối diễn ra nhờ sự tái cấu trúc topo không gian chứ không phải do sự co ngắn liên kết cộng hóa trị $\text{Si-O}$.

2. Giải mã hiện tượng động học không đồng nhất qua 5 nhóm nguyên tử đặc biệt: Phân tích quá trình tạo đám của các nhóm nguyên tử $A_1 \div A_5$ đã cung cấp bằng chứng trực tiếp về sự phân tách pha động học:

  • Nhóm nguyên tử linh động nhất ($A_2$) và nhóm có tần suất chuyển đổi cao ($A_4$) có xu hướng kết tụ không gian mạnh mẽ, tạo thành các đám có kích thước cực đại $N_{a_A2}$ lớn hơn đáng kể so với nhóm chọn ngẫu nhiên ($A_1$).
  • Ngược lại, nhóm nguyên tử kém linh động ($A_3$) và ít chuyển đổi ($A_5$) tạo thành các "đám cứng" (rigid cages) giam hãm các hạt xung quanh.
  • Tại thời điểm trễ động học $t_d = 47.8\text{ ps}$ ($P = 0\text{ GPa}$), kích thước đám linh động đạt cực đại, phản ánh trạng thái kích hoạt tập thể (dynamic facilitation).

3. Bản chất của khuếch tán dị thường và vai trò của chuyển đổi có ích: Nghiên cứu giải thích hiện tượng nghịch lý: hệ số khuếch tán $D$ của $\text{Si}$ và $\text{O}$ tăng khi áp suất tăng từ 0 lên khoảng $10 \div 15\text{ GPa}$ trước khi suy giảm ở áp suất cao hơn. Luận án chứng minh rằng:

  • Tốc độ chuyển đổi phối trí toàn phần $\nu_{\text{tran}}$ tăng tuyến tính theo áp suất do mật độ va chạm tăng.
  • Tuy nhiên, tỷ số chuyển đổi có ích $\xi = M_{bt}/M_{bo}$ đạt cực đại ở vùng áp suất trung gian $10 \div 15\text{ GPa}$. Ở vùng này, năng lượng kích hoạt chuyển đổi giữa $\text{SiO}_4$ và $\text{SiO}_6$ là thấp nhất do sự tương thích mô-đun khối, tạo điều kiện cho các chuyển động tương quan nhóm (chuyển động của một cụm nguyên tử qua kênh thể tích tự do) diễn ra thuận lợi nhất.
Bảng 1: Tiến hóa cấu trúc và động học của SiO2 lỏng theo áp suất (T = 3000 K)
+--------------+---------------+---------------+---------------+--------------------+
| Áp suất (GPa)| Tỷ phần SiO4  | Tỷ phần SiO5  | Tỷ phần SiO6  | D_Si (10^-5 cm2/s) |
+--------------+---------------+---------------+---------------+--------------------+
|      0       |     92.4%     |      7.3%     |      0.3%     |        1.85        |
|      5       |     68.1%     |     27.2%     |      4.7%     |        2.40        |
|     10       |     35.6%     |     46.8%     |     17.6%     |        3.15 (Đỉnh) |
|     15       |     12.3%     |     43.2%     |     44.5%     |        2.90        |
|     20       |      3.1%     |     28.5%     |     68.4%     |        1.75        |
|     30       |      0.2%     |      8.9%     |     90.9%     |        0.65        |
+--------------+---------------+---------------+---------------+--------------------+

4. Hiện tượng tách pha vi mô và mức độ polyme hóa trong $\text{MgSiO}_3$ lỏng: Trong hệ $\text{MgSiO}_3$ lỏng ở $T = 3500\text{ K}$:

  • Ở $P = 0\text{ GPa}$, sự hiện diện của cation $\text{Mg}^{2+}$ phá vỡ tính liên tục của mạng $\text{Si-O}$, làm tỷ phần oxy không cầu ($\text{NBO}$) chiếm tới $45%$, oxy cầu ($\text{BO}$) chiếm $45%$ và oxy tự do (chỉ liên kết với $\text{Mg}$) chiếm khoảng $8 \div 10%$.
  • Cấu trúc địa phương của tứ diện $\text{SiO}4$ ($r{\text{Si-O}} = 1.62\text{ \AA}$) hoàn toàn được bảo toàn giống như trong $\text{SiO}_2$ nguyên chất.
  • Kỹ thuật trực quan hóa 3D phát hiện hiện tượng tách pha vi mô rõ rệt: các đơn vị $\text{SiO}_x$ liên kết với nhau tạo thành các chuỗi/mạng con giàu $\text{Si}$, trong khi các đơn vị $\text{MgO}_y$ ($y = 4, 5, 6$) liên kết tạo thành các miền giàu $\text{Mg}$. Hai mạng con này đan xen trong không gian thông qua các cầu nối $\text{Si-O-Mg}$.
  • Khi áp suất tăng lên $30\text{ GPa}$, tỷ phần $\text{NBO}$ và oxy tự do giảm mạnh, tỷ phần $\text{BO}$ tăng lên $>80%$, số phối trí trung bình $Z_{\text{Si-O}}$ tăng từ 4.0 lên 5.8 và $Z_{\text{Mg-O}}$ tăng từ 4.6 lên 6.5, chứng tỏ mức độ polyme hóa toàn mạng tăng mạnh dưới tác dụng của áp suất nén.
Bảng 2: Phân bố các loại nguyên tử Oxy trong MgSiO3 lỏng theo áp suất (T = 3500 K)
+--------------+------------------+----------------------+--------------------+
| Áp suất (GPa)| Oxy cầu (BO, %)  | Oxy không cầu (NBO,%)| Oxy tự do (Free O%)|
+--------------+------------------+----------------------+--------------------+
|      0       |      45.2%       |        46.5%         |        8.3%        |
|      5       |      54.8%       |        39.1%         |        6.1%        |
|     10       |      66.3%       |        29.8%         |        3.9%        |
|     20       |      78.5%       |        19.7%         |        1.8%        |
|     30       |      86.2%       |        13.1%         |        0.7%        |
+--------------+------------------+----------------------+--------------------+

Implications đa chiều

  • Về mặt lý thuyết: Thống nhất mô hình nhiệt động hai trạng thái với động học hồi phục bất đối xứng, cung cấp một khung giải thích toàn diện cho các chất lỏng tạo mạng tứ diện ($\text{H}_2\text{O}, \text{GeO}_2, \text{BeF}_2$).
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn phân tích cụm vi mô và phân loại chuyển đổi phối trí hữu ích, có khả năng áp dụng rộng rãi cho các hệ vật liệu vô định hình phức tạp khác như kim loại thủy tinh (bulk metallic glasses) hay chất điện ly rắn.
  • Về mặt địa vật lý và địa hóa học: Cung cấp các tham số truyền chất chính xác (hệ số khuếch tán, độ nhớt suy biến) của dung nham silicat trong điều kiện manti sâu, giúp mô hình hóa chính xác sự đối lưu manti, quá trình phân tách magma trong đại dương magma nguyên thủy và sự tiến hóa nhiệt của Trái Đất sơ khai.
  • Về mặt công nghệ vật liệu: Cung cấp cơ sở khoa học để thiết kế thủy tinh hoạt tính chịu lực, màng gốm xốp chịu nhiệt và đặc biệt là tối ưu hóa ma trận thủy tinh hóa (vitrification) đa thành phần để giam giữ an toàn các chất thải phóng xạ chu kỳ bán rã dài.

Limitations và Future Research

Luận án thừa nhận một số giới hạn phương pháp luận xuất phát từ điều kiện tính toán:

  1. Mô hình thế tương tác cặp ion: Thế tương tác BKS và OG là thế cặp hiệu dụng với điện tích cố định, chưa mô tả được hiệu ứng phân cực điện tử phi tuyến (electronic polarizability) và sự chuyển dịch điện tích động (charge transfer) khi các liên kết hóa học bị biến dạng mạnh dưới áp suất siêu cao ($>30\text{ GPa}$).
  2. Hiệu ứng kích thước và thời gian mô phỏng: Mẫu mô phỏng quy mô 1.998 và 5.000 nguyên tử với thời gian quan sát hàng chục nano-giây tuy đủ để khảo sát trạng thái lỏng cân bằng ở $T \ge 2600\text{ K}$, nhưng chưa đủ để tiếp cận trực tiếp động học chậm chạp sát nhiệt độ chuyển pha thủy tinh thực nghiệm $T_g \approx 1400\text{ K}$.
  3. Thành phần hóa học lý tưởng hóa: Nghiên cứu tập trung vào hai hệ tinh khiết $\text{SiO}_2$ và $\text{MgSiO}_3$, chưa xét tới sự hiện diện của các oxit phụ và chất dễ bay hơi như $\text{Al}_2\text{O}_3, \text{CaO}, \text{FeO}, \text{H}_2\text{O}$ và $\text{CO}_2$ vốn luôn tồn tại trong magma tự nhiên.

Chương trình nghiên cứu tiếp nối trong giai đoạn 5–10 năm tới bao gồm:

  • Ứng dụng thế tương tác học máy chính xác lượng tử (Machine Learning Interatomic Potentials - MLIPs như DeepMD, GAP) dựa trên dữ liệu phiếm hàm mật độ (DFT) để mô tả chính xác liên kết cộng hóa trị phân cực.
  • Mở rộng mô phỏng sang các hệ đa cấu tử phức tạp của lớp phủ Trái Đất ($\text{MgO-FeO-SiO}_2\text{-Al}_2\text{O}_3\text{-CaO}$).
  • Khảo sát đáp ứng lưu biến phi tuyến (non-linear rheology) và độ nhớt cắt dưới trường biến dạng ngoài bằng động lực học phân tử không cân bằng (NEMD).
  • Kết hợp mô phỏng với các phép đo thực nghiệm tại nguồn bức xạ synchrotron thế hệ mới và tế bào đe kim cương nung laser (Laser-Heated Diamond Anvil Cell).

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

Công trình luận án của TS. Luyện Thị San đã tạo ra dấu ấn học thuật rõ nét với 06 công trình khoa học công bố, trong đó có 04 bài báo trên các tạp chí quốc tế uy tín thuộc danh mục ISI:

  1. Journal of Non-Crystalline Solids / Modern Physics B 30, pp. 1650059 (2016).
  2. European Physical Journal B 89, pp. 73 (2016).
  3. High Pressure Research 36, pp. 187-197 (2016).
  4. Physica B: Condensed Matter 501, pp. 27-33 (2016). Cùng 02 công trình trên tạp chí chuyên ngành trong nước và kỷ yếu hội nghị quốc tế (Journal of Science of HNUE, Journal of Physics: Conference Series).

Tác động lan tỏa của luận án thể hiện trên ba phương diện:

  • Học thuật: Cung cấp bộ dữ liệu chuẩn mực về cấu trúc và động học silicat lỏng, được trích dẫn và sử dụng bởi cộng đồng mô phỏng vật liệu vô định hình quốc tế.
  • Công nghiệp: Đóng góp dữ liệu nhiệt động học cho công nghệ chế tạo gốm kỹ thuật cao, sợi quang silica tinh khiết và công nghệ cố định hóa chất thải hạt nhân bằng thủy tinh silicat.
  • Địa chất: Hoàn thiện mô hình động lực học magma, nâng cao độ chính xác của các mô phỏng máy tính dự báo phun trào núi lửa và chu trình kiến tạo mảng toàn cầu.

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

  • Nghiên cứu sinh và học giả chuyên ngành Vật lý Vật chất ngưng tụ: Tiếp cận phương pháp phân tích cấu trúc phối trí và giải thuật truy vết động học cụm nguyên tử.
  • Các nhà Địa vật lý và Khoáng vật học: Sử dụng các số liệu về cấu trúc $g(r)$, độ nén ép thể tích, hệ số khuếch tán của $\text{SiO}_2$ và $\text{MgSiO}_3$ ở điều kiện manti để xây dựng mô hình tiến hóa hành tinh.
  • Kỹ sư R&D ngành Vật liệu Thủy tinh – Gốm sứ: Ứng dụng các quy luật biến tính mạng và tỷ lệ oxy không cầu ($\text{NBO}$) để nâng cao độ bền cơ nhiệt của vật liệu chịu lửa và kính chuyên dụng.
  • Cơ quan Quản lý An toàn Hạt nhân: Có thêm luận cứ khoa học thực chứng về độ bền vững cấu trúc vi mô của khối thủy tinh hóa chứa đồng vị phóng xạ.

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ì? Trả lời: Đóng góp độc đáo nhất là việc mở rộng Mô hình hai trạng thái (Two-State Model) bằng cách lượng hóa cụ thể cấu trúc vi mô thông qua 5 đơn vị phối trí ($\text{SiO}_4, \text{SiO}_5, \text{SiO}_6, \text{OSi}_2, \text{OSi}_3$). Luận án chứng minh rằng tính đa thù hình trong $\text{SiO}_2$ lỏng là sự cùng tồn tại không gian của hai pha LDL (mật độ thấp) và HDL (mật độ cao), hòa giải mâu thuẫn giữa mô hình chuyển pha nhiệt động học và các quan sát thực nghiệm biến đổi liên tục.

2. Đổi mới phương pháp luận nổi bật nhất so với các nghiên cứu trước đây là gì? Trả lời: So với các nghiên cứu trước đây (như Sato et al., 2010 hay Kubicki, 2001 chỉ dùng số phối trí trung bình), luận án đã phát triển phương pháp phân loại chuyển đổi phối trí thành "chuyển đổi có ích" ($M_{bt}$) và "chuyển đổi vô ích" ($M_{bo}$), kết hợp thuật toán truy vết cụm động học trên 5 tập hợp nguyên tử mục tiêu ($A_1 \div A_5$). Phương pháp này cho phép liên kết trực tiếp một đại lượng cấu trúc vi mô với hệ số khuếch tán vĩ mô $D$.

3. Phát hiện bất ngờ nhất từ dữ liệu mô phỏng là gì? Trả lời: Đó là cơ chế vi mô của hiện tượng khuếch tán dị thường: sự gia tăng của hệ số khuếch tán $D$ ở dải áp suất $10 \div 15\text{ GPa}$ không xuất phát từ việc tăng tần suất va chạm đơn thuần, mà do tỷ số chuyển đổi có ích $\xi = M_{bt}/M_{bo}$ đạt cực đại nhờ sự tương thích mô-đun đàn hồi giữa các đơn vị $\text{SiO}_4$ và $\text{SiO}_6$, kích hoạt các chuyển động tương quan nhóm.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Trả lời: Có. Luận án cung cấp đầy đủ các thông số hàm thế (BKS và OG), điện tích hạt, kích thước mẫu (1.998 và 5.000 nguyên tử), điều kiện biên tuần hoàn PBC, bước nhảy tích phân $\Delta t = 0.47\text{ fs}$, quy trình nung - ủ - nén chi tiết qua các ensemble NPT và NVE, cùng bán kính ngắt phối trí $r_c$ chuẩn xác, cho phép tái lập 100% kết quả trên mọi phần mềm MD tiêu chuẩn (như LAMMPS hay DL_POLY).

5. Định hướng nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Trả lời: Tập trung vào việc tích hợp thế học máy nơ-ron (Machine Learning Potentials) huấn luyện từ tính toán ab-initio DFT, mở rộng khảo sát các hệ magma silicat đa cấu tử chứa $\text{FeO}, \text{CaO}, \text{Al}_2\text{O}_3, \text{H}_2\text{O}$, và nghiên cứu độ nhớt dưới biến dạng cắt phi tuyến ở điều kiện manti dưới.

Kết luận

Luận án tiến sĩ của tác giả Luyện Thị San là một công trình nghiên cứu công phu, chuẩn mực và xuất sắc trong lĩnh vực vật lý chất ngưng tụ và vật liệu vô định hình. Tóm tắt 5 đóng góp cốt lõi:

  1. Xác lập bức tranh vi mô hoàn chỉnh về tính đa thù hình của $\text{SiO}_2$ lỏng dựa trên 5 đơn vị phối trí và mô hình hai trạng thái (pha mật độ thấp $\text{SiO}_4\text{-OSi}_2$ và pha mật độ cao $\text{SiO}_5/\text{SiO}_6\text{-OSi}_3$).
  2. Chứng minh thực chứng số vai trò quyết định của "chuyển đổi phối trí có ích" và chuyển động tương quan nhóm đối với cơ chế khuếch tán và hiện tượng khuếch tán dị thường.
  3. Cung cấp bằng chứng định lượng trực tiếp về hiện tượng động học không đồng nhất thông qua sự phân tách không gian giữa các đám nguyên tử linh động và đám nguyên tử kém linh động.
  4. Làm sáng tỏ cơ chế biến tính mạng $\text{Si-O}$ bởi cation $\text{Mg}^{2+}$, phát hiện hiện tượng tách pha vi mô giữa miền giàu $\text{Si}$ và miền giàu $\text{Mg}$ trong $\text{MgSiO}_3$ lỏng ở áp suất cao.
  5. Định lượng hóa sự tiến hóa của oxy cầu ($\text{BO}$), oxy không cầu ($\text{NBO}$) và oxy tự do dưới tác dụng của áp suất nén từ 0 đến $30\text{ GPa}$, đóng góp cơ sở dữ liệu quan trọng cho khoa học Trái Đất và công nghệ vật liệu silicat.