Tổng quan về luận án

Nghiên cứu cấu trúc vi mô và cơ chế chuyển dịch pha ở cấp độ nguyên tử đóng vai trò then chốt trong vật lý chất rắn và khoa học vật liệu hiện đại. Luận án tiến sĩ mang tên "Mô phỏng cấu trúc và quá trình chuyển pha của các vật liệu Fe, FeB và SiO2" (Chuyên ngành: Vật lý kỹ thuật, Mã số: 9520401) do nghiên cứu sinh Giáp Thị Thùy Trang thực hiện dưới sự hướng dẫn của PGS. Phạm Khắc Hùng và PGS. Phạm Hữu Kiên tại Trường Đại học Bách khoa Hà Nội (2020) đã giải quyết những bài toán cốt lõi về trạng thái vô định hình, động học kết tinh và cơ chế đậm đặc hóa dưới áp suất cao.

                  +---------------------------------------------------------+
                  |  HỆ THỐNG MÔ PHỎNG ĐỘNG LỰC HỌC PHÂN TỬ (MD SIMULATION)  |
                  +---------------------------------------------------------+
                                               |
                       +-----------------------+-----------------------+
                       |                                               |
                       v                                               v
        +------------------------------+                +------------------------------+
        |   HỆ HẠT NANO KIM LOẠI/HỢP KIM|                |   HỆ VẬT LIỆU MẠNG OXIT       |
        |   - Fe (5.000 & 10.000 nguyên tử) |            |   - SiO2 lỏng (20.000 nguyên tử)  |
        |   - FexB100-x (x=90, 95)     |                |   - SiO2 VĐH (5.000 nguyên tử)    |
        +------------------------------+                +------------------------------+
                       |                                               |
          [Thế tương tác Pak-Doyama]                        [Thế tương tác cặp BKS]
                       |                                               |
                       v                                               v
        +------------------------------+                +------------------------------+
        | ĐẶC TRƯNG CHUYỂN PHA VĐH-TT  |                | ĐỘNG HỌC & ĐẬM ĐẶC HÓA (P, T)|
        | - Cấu trúc Lõi - Vỏ          |                | - Vùng vi mô tinh khiết      |
        | - Tạo mầm trung gian bcc/ih  |                | - Động học không đồng nhất   |
        | - Bo kìm hãm tại biên pha    |                | - Chuyển đổi D4 -> D5 -> D6  |
        +------------------------------+                +------------------------------+

Khoảng trống nghiên cứu (research gap) lớn nhất được xác định là: Phần lớn các mô phỏng trước đây đối với hạt nano Fe và FeB chỉ tiếp cận ở quy mô cụm nguyên tử rất nhỏ ($N \le 331 - 686$ nguyên tử) hoặc dựa hoàn toàn vào lý thuyết tạo mầm cổ điển (CNT) mà bỏ qua các trạng thái cấu trúc trung gian đa thù hình tại vùng biên bề mặt. Đối với $\text{SiO}_2$, dù mô hình mạng ngẫu nhiên của Zachariasen [57] đã định hình nền tảng cấu trúc tứ diện $\text{SiO}_4$, nhưng sự tồn tại của các vùng vi mô tinh khiết (pure micro-regions), cơ chế động học không đồng nhất (dynamical heterogeneity) và bức tranh phân rã - hợp nhất của các trường miền cấu trúc (domain $D_x$) dưới áp suất siêu cao từ 0 đến 100 GPa vẫn chưa có lời giải chi tiết ở quy mô mô hình lớn.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Cấu trúc đa thù hình của hạt nano vô định hình (VĐH) Fe và $\text{Fe}x\text{B}{100-x}$ tiến hóa như thế nào qua các giai đoạn ủ nhiệt?
    • H1: Hạt nano VĐH phân tách thành cấu trúc lõi - vỏ; quá trình tinh thể hóa tuân theo cơ chế mọc mầm đa bước qua các pha trung gian giả bền (bcc, icosahedral).
  • RQ2: Nguyên tử á kim Bo (B) ảnh hưởng như thế nào đến động học tạo mầm và sự phát triển tinh thể trong hạt nano FeB?
    • H2: Nguyên tử B có xu hướng dịch chuyển ra vùng biên pha và lớp vỏ ngoài, đóng vai trò rào cản cản trở sự tái sắp xếp mạng tinh thể của Fe.
  • RQ3: Bản chất của động học không đồng nhất và sự mất trật tự hóa học trong $\text{SiO}_2$ lỏng ở quy mô lớn là gì?
    • H3: Trong $\text{SiO}_2$ lỏng tồn tại các vùng vi mô tinh khiết chỉ chứa riêng biệt nguyên tử Si hoặc O cùng các mạng con bền vững (SSLA) có tốc độ hồi phục khác biệt.
  • RQ4: Cơ chế đậm đặc hóa không gian của $\text{SiO}_2$ lỏng và vô định hình dưới tác động của áp suất cực hạn diễn ra theo quy luật cấu trúc nào?
    • H4: Sự tăng áp suất điều khiển quá trình chuyển đổi phối trí liên tục từ tứ diện $\text{SiO}_4$ ($D_4$) sang bát diện $\text{SiO}_6$ ($D_6$) thông qua trạng thái trung gian $\text{SiO}_5$ ($D_5$), đi kèm sự co cụm thể tích Voronoi và tái cấu trúc mạng domain.

Khung lý thuyết của luận án được xây dựng trên nền tảng nhiệt động học chuyển pha Gibbs ($dG = VdP - SdT$), lý thuyết động học chuyển pha Johnson-Mehl-Avrami (JMA) [$X = 1 - \exp(-kt^\eta)$], lý thuyết tạo mầm cổ điển (Classical Nucleation Theory - CNT), quy tắc bước Ostwald (Ostwald's step rule) và mô hình mạng ngẫu nhiên liên tục (Continuous Random Network - CRN).

Quy mô khảo sát bao gồm: Hạt nano Fe và $\text{Fe}x\text{B}{100-x}$ ($x=90, 95$) quy mô $5.000 - 10.000$ nguyên tử tại nhiệt độ $300 - 900\text{ K}$; hệ $\text{SiO}_2$ lỏng quy mô lớn $20.000$ nguyên tử (6.666 Si + 13.332 O) ở $3.000\text{ K}$ ($P=0$) và $3.500\text{ K}$ ($P = 0 - 45\text{ GPa}$); hệ $\text{SiO}_2$ vô định hình $5.000$ nguyên tử ở $500\text{ K}$ ($P = 0 - 100\text{ GPa}$). Kết quả công trình tạo đột phá lý thuyết với 08 bài báo khoa học, trong đó có 06 công trình trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Journal of Non-Crystalline Solids, The European Physical Journal B, Materials Research Express, AIP Advances).

Literature Review và Positioning

Lý thuyết cấu trúc mạng oxit vô định hình khởi nguồn từ công trình kinh điển của Zachariasen (1932) [57], xác lập mô hình các tứ diện $\text{SiO}4$ liên kết qua cầu oxy $\text{Si-O-Si}$ tạo mạng không gian ba chiều không có trật tự tầm xa. Mô hình này sau đó được Mozzi và Warren (1969) [58] chứng minh thực nghiệm bằng nhiễu xạ tia X, xác định góc liên kết $\theta{\text{Si-O-Si}} \approx 144^\circ$ và khoảng cách liên kết $r_{\text{Si-O}} = 1,62\text{ \AA}$. Khi vật liệu chịu áp suất nén cao, các công trình của Sato et al. (2010) [84] và Min Wu et al. (2012) [83] ghi nhận sự gia tăng số phối trí của Si từ 4 lên 6 trong khoảng áp suất từ $12\text{ GPa}$ đến trên $40\text{ GPa}$. Tuy nhiên, tranh luận học thuật sâu sắc vẫn tồn tại giữa hai quan điểm: (1) Quá trình tăng số phối trí diễn ra đồng nhất trên toàn mạng tinh thể hay (2) Diễn ra cục bộ thông qua sự phân chia các pha vi mô (domain) không đồng nhất.

                             TIẾN TRÌNH NGHIÊN CỨU VÀ ĐỊNH VỊ HỌC THUẬT
====================================================================================================
Giai đoạn / Tác giả               Phát hiện & Hạn chế                      Định vị Luận án này
----------------------------------------------------------------------------------------------------
Zachariasen (1932) [57]            Mô hình mạng ngẫu nhiên CRN,             Mở rộng từ mạng tứ diện đơn
Mozzi & Warren (1969) [58]         tứ diện SiO4, trật tự gần.              thuần sang cấu trúc 5 domain
                                   (Chưa xét áp suất cao)                  (D4 -> D5 -> D6) dưới nén.
----------------------------------------------------------------------------------------------------
Belashchenko (1980s) [24]          Mô phỏng Fe lỏng/VĐH quy mô nhỏ         Nâng quy mô lên 10.000 nguyên
Evteev et al. (2003) [25]          (686 nguyên tử); tính chất mầm           tử; giải mã chi tiết cơ chế
Xiaohua Li et al. (2014) [19]      chưa phản ánh hiệu ứng bề mặt.          lõi - vỏ và mầm giả bền bcc/ih.
----------------------------------------------------------------------------------------------------
van Beest et al. (BKS, 1990) [96]  Thế tương tác BKS mô phỏng              Ứng dụng BKS trên hệ 20.000
M. Scott Shell et al. (2002) [53]  chính xác SiO2 lỏng/rắn.               nguyên tử; phát hiện vùng vi mô
Sato et al. (2010) [84]            (Chưa khảo sát trường domain)            tinh khiết và động học SSLA.
====================================================================================================

Đối với kim loại và hợp kim vô định hình, nghiên cứu tiên phong của Pak và Doyama (1969) [95] đã thiết lập thế tương tác cặp thực nghiệm biểu diễn chính xác năng lượng tương tác giữa các nguyên tử Fe. Yamamoto et al. (1978) [23] và Belashchenko (1989) [24] đã áp dụng thế Pak-Doyama để mô phỏng Fe lỏng và vô định hình, tái tạo thành công đặc trưng tách đỉnh thứ hai của hàm phân bố xuyên tâm $g(r)$ – dấu hiệu đặc trưng của trạng thái thủy tinh kim loại. Gần đây hơn, Xiaohua Li et al. (2014) [19] khảo sát hạt nano Fe 331 nguyên tử đường kính $2\text{ nm}$, kết hợp CNT và lý thuyết bề mặt khuếch tán để giải thích tốc độ tạo mầm tại $750 - 850\text{ K}$. Tuy nhiên, mô hình 331 nguyên tử quá nhỏ để mô tả sự phát triển độc lập của mầm tinh thể bên trong lõi mà không bị chi phối áp đảo bởi hiệu ứng giam cầm bề mặt.

Luận án này định vị sự đột phá học thuật thông qua việc so sánh đối chuẩn quốc tế:

  1. So với nghiên cứu của Xiaohua Li et al. [19]: Luận án mở rộng quy mô hạt nano Fe lên $5.000 - 10.000$ nguyên tử ($R = 28 - 34\text{ \AA}$), cho phép tách biệt rõ ràng động học tinh thể hóa tại vùng lõi và vùng bề mặt, phát hiện sự cạnh tranh năng lượng giữa cấu trúc icosahedral và bcc theo quy tắc bước Ostwald.
  2. So với nghiên cứu của Sato et al. [84] và Min Wu et al. [83]: Luận án không chỉ dừng lại ở việc tính toán tỷ phần phối trí $z_{\text{Si}}$ trung bình, mà lần đầu tiên thiết lập bức tranh cấu trúc 5 domain không gian ($D_4, D_5, D_6$, domain biên $\text{BD}$), giải thích tường tận cơ chế biến thiên thể tích Voronoi và sự thay đổi góc liên kết $\theta_{\text{Si-O-Si}}$ từ $146^\circ$ ($0\text{ GPa}$) xuống $105^\circ$ ($45\text{ GPa}$).

Đó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à bổ sung sâu sắc cho các học thuyết nhiệt động học và động học chất rắn:

  1. Bổ sung Lý thuyết Tạo mầm Cổ điển (CNT): Chứng minh rằng quá trình tinh thể hóa hạt nano kim loại vô định hình không diễn ra theo một bước nhảy trực tiếp vượt qua rào cản thế năng $\Delta G^*$, mà tuân thủ quy tắc bước Ostwald thông qua các trạng thái cấu trúc trung gian. Mầm kết tinh ban đầu xuất hiện dưới dạng các cụm icosahedral (đối xứng bậc năm) hoặc bcc giả bền, sau đó tái sắp xếp cấu trúc để đạt cực tiểu thế năng.
  2. Lý thuyết về vai trò pha á kim trong hợp kim vô định hình: Luận án chỉ ra rằng trong hệ $\text{Fe}x\text{B}{100-x}$, các nguyên tử Bo không phân bố ngẫu nhiên đồng nhất trong quá trình ủ nhiệt mà bị đẩy dồn về vùng biên giữa pha vô định hình và pha tinh thể, hình thành nên rào cản nhiệt động làm tăng hoạt năng khuếch tán bề mặt $Q_d$, ức chế mạnh mẽ tốc độ phát triển mầm tinh thể của Fe.
  3. Mô hình Mạng vi mô không đồng nhất trong chất lỏng oxit: Thách thức quan điểm truyền thống cho rằng chất lỏng $\text{SiO}_2$ hoàn toàn đồng nhất về mặt hóa học. Luận án xác lập sự tồn tại của các "vùng vi mô tinh khiết" (chỉ chứa riêng biệt ion Si hoặc ion O) và các mạng con bền vững ($SSLA$), đặt nền móng lý thuyết mới cho hiện tượng động học không đồng nhất.
                              KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN
+--------------------------------------------------------------------------------------------------+
| PHƯƠNG PHÁP PHÂN TÍCH HÌNH HỌC VÀ ĐỘNG HỌC MỨC NGUYÊN TỬ                                         |
+--------------------------------------------------------------------------------------------------+
  |---> 1. Phân tích Lân cận chung (CNA): Xác định mạng tinh thể bcc, fcc, hcp, icosahedral
  |---> 2. Tần số cấu trúc động học <fx>: Đánh giá mức độ ổn định thời gian của các cụm tinh thể
  |---> 3. Phân tích Hạt Lõi - Vỏ (Core-Shell): Tách biệt mật độ DACP (lõi) và DASP (vỏ)
  |---> 4. Phân tích Không gian Domain (Dx): Phân loại trường cấu trúc tứ diện D4, tam giác đôi D5, bát diện D6
  |---> 5. Đa diện Voronoi & SPFA: Đo lường thể tích tự do và kích thước lỗ rỗng cấu trúc
+--------------------------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp đồng thời 5 công cụ tính toán hình học không gian và động học tiên tiến:

  • Tần số cấu trúc động học $\langle f_x \rangle$: Đo lường xác suất một nguyên tử duy trì trạng thái cấu trúc tinh thể $x$ qua các bước thời gian liên tiếp, cho phép phân biệt chính xác giữa các dao động nhiệt ngẫu nhiên và mầm kết tinh thực thụ.
  • Phương pháp phân tích Hạt lõi - vỏ (Core-Shell Analysis): Thiết lập các thông số định lượng mật độ nguyên tử trong lõi ($\text{DACP}$) và trong vỏ ($\text{DASP}$), làm sáng tỏ sự khác biệt về năng lượng bề mặt $G_S$ của hạt nano.
  • Phương pháp phân tích Miền cấu trúc (Domain Analysis): Phân chia không gian vật liệu thành các miền $D_4, D_5, D_6$ (vùng mật độ cao) và miền biên $\text{BD}$ (vùng mật độ thấp), định lượng hóa tương tác giữa các mạng con.
  • Hình cầu đi qua 4 nguyên tử (SPFA - Sphere Passing Four Atoms) và Đa diện Voronoi: Xác định chính xác sự biến đổi thể tích tự do cục bộ và kích thước các lỗ xốp trong cấu trúc vô định hình.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý thực chứng tính toán (Computational Positivism), kết hợp mô phỏng động lực học phân tử cổ điển với phân tích thống kê cấu trúc đa cấp độ. Thiết kế nghiên cứu bao phủ từ cấp độ vi mô (nguyên tử riêng lẻ, đa diện phối trí) đến cấp độ trung mô (đám hạt, mạng con, domain) và cấp độ vĩ mô (mật độ khối, áp suất, thế năng toàn phần).

                      THIẾT KẾ CÁC MÔ HÌNH MÔ PHỎNG ĐỘNG LỰC HỌC PHÂN TỬ
====================================================================================================
Hệ vật liệu     Quy mô mẫu       Thế tương tác      Bước thời gian (dt)  Nhiệt độ (T)    Áp suất (P)
----------------------------------------------------------------------------------------------------
Hạt nano Fe     5.000 - 10.000   Pak-Doyama         4,6 x 10^-14 s       300 - 900 K     Chân không (tự do)
Hạt nano FeB    5.000 nguyên tử  Pak-Doyama mở rộng 4,6 x 10^-14 s       300 - 900 K     Chân không (tự do)
SiO2 lỏng       20.000 nguyên tử BKS                1,0 fs (10^-15 s)    3000 - 3500 K   0 - 45 GPa
SiO2 VĐH        5.000 nguyên tử  BKS                0,48 fs              500 K           0 - 100 GPa
====================================================================================================

Quy trình nghiên cứu rigorous

Quy trình giải phương trình chuyển động Newton được thực hiện thông qua thuật toán Verlet vận tốc, đảm bảo bảo toàn năng lượng và độ ổn định tích phân quỹ đạo: $$\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$$ $$\vec{v}_i(t) = \frac{\vec{r}_i(t + \Delta t) - \vec{r}_i(t - \Delta t)}{2\Delta t}$$

  1. Đối với hệ hạt nano Fe và FeB: Sử dụng thế tương tác cặp Pak-Doyama có dạng đa thức bậc bốn cắt cụt tại bán kính $r_{cut}$: $$U(r) = a(r - b)^4 + c(r - d)^2 + e \quad (0 \le r \le r_{cut})$$ Các tham số thế cho cặp $\text{Fe-Fe}$ gồm: $a = -0,18892\text{ eV/\AA}^4$, $b = -1,82709\text{ \AA}$, $c = 1,70192\text{ eV/\AA}^2$, $d = -2,50849\text{ \AA}$, $e = -0,19829\text{ eV}$, $r_{cut} = 3,14\text{ \AA}$. Mẫu được thiết lập với điều kiện biên tự do trong quả cầu bán kính $R = 28\text{ \AA}$ ($5.000$ hạt) và $R = 34\text{ \AA}$ ($10.000$ hạt), hồi phục nhiệt qua $2 \times 10^7$ bước MD, sau đó ủ nhiệt đẳng nhiệt tới $8 \times 10^7 - 9 \times 10^7$ bước.

  2. Đối với hệ $\text{SiO}_2$: Áp dụng thế cặp BKS (van Beest - Kramer - van Santen) kết hợp tương tác Coulomb tầm xa và thế Buckingham tầm ngắn: $$U_{ij}(r_{ij}) = \frac{q_i q_j e^2}{r_{ij}} + A_{ij}\exp(-B_{ij}r_{ij}) - \frac{C_{ij}}{r_{ij}^6}$$ Điện tích hiệu dụng: $q_{\text{Si}} = +2,4e$, $q_{\text{O}} = -1,2e$. Các tham số: cặp $\text{Si-O}$ có $A = 18003,7572\text{ eV}$, $B = 4,87318\text{ \AA}^{-1}$, $C = 1133,5381\text{ eV}\cdot\text{\AA}^6$; cặp $\text{O-O}$ có $A = 1388,7730\text{ eV}$, $B = 2,76000\text{ \AA}^{-1}$, $C = 175,0000\text{ eV}\cdot\text{\AA}^6$. Tương tác tầm xa Coulomb được xử lý chính xác bằng phương pháp tổng Ewald trong điều kiện biên tuần hoàn 3 chiều.

Data và phân tích

Độ tin cậy của mô hình được đối chuẩn nghiêm ngặt với các dữ liệu thực nghiệm nhiễu xạ tia X và neutron:

  • Mô hình $\text{SiO}_2$ lỏng tại $P = 0$: Tính toán đạt $r_{\text{Si-Si}} = 3,12\text{ \AA}$, $r_{\text{Si-O}} = 1,60\text{ \AA}$, $r_{\text{O-O}} = 2,60\text{ \AA}$, số phối trí $z_{\text{Si}} = 4,03$, $z_{\text{O}} = 2,01$, phân bố góc $\theta_{\text{O-Si-O}} = 106,0^\circ$, $\theta_{\text{Si-O-Si}} = 146,0^\circ$, khớp hoàn hảo với dữ liệu thực nghiệm của Mozzi & Warren [58] ($r_{\text{Si-O}} = 1,62\text{ \AA}$, $\theta_{\text{Si-O-Si}} = 144^\circ$) và Meade et al. [63].
  • Động học hạt: Đánh giá thông qua độ dịch chuyển bình phương trung bình $\text{MSDA} = \langle |\vec{r}_i(t) - \vec{r}_i(0)|^2 \rangle$ và số lượng nguyên tử từng là phối trí ($\text{NAVBU}$), phân tách rõ hai nhóm nguyên tử linh động ($\text{MS}$) và không linh động ($\text{IMS}$).

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

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

[!IMPORTANT] Trích dẫn minh chứng từ luận án:

  1. "Mẫu hạt nano Fe tinh thể hóa hoàn toàn bao gồm: phần lõi với cấu trúc tinh thể và phần bề mặt với cấu trúc vô định hình xốp." (Chương 3).
  2. "Phân tích các hạt lõi/vỏ cho thấy trong $\text{SiO}_2$ lỏng tồn tại các vùng vi mô tinh khiết, ở đó chỉ chứa hoặc nguyên tử O hoặc Si." (Chương 4).
  3. "Khi áp suất tăng, $\text{SiO}_2$ lỏng (VĐH) xảy ra chuyển đổi cấu trúc từ domain $D_4$ sang domain $D_6$ thông qua domain $D_5$, kèm theo hai quá trình chia tách và hợp nhất các domain." (Chương 4).
                      CƠ CHẾ TIẾN HÓA VÀ CHUYỂN PHA THEO ÁP SUẤT CỦA SiO2
+--------------------------------------------------------------------------------------------------+
| Dải áp suất         Cấu trúc Domain chủ đạo               Đặc trưng phối trí SiOx / Thể tích     |
+--------------------------------------------------------------------------------------------------+
| 0 - 10 GPa          1 Domain D4 duy nhất                  Tứ diện SiO4 chiếm ưu thế (~100%)       |
| Tại điểm 10 GPa     Tồn tại đồng thời D4 & D5             Bắt đầu phân tách domain, r_Si-O tăng   |
| 10 - 30 GPa         1 Domain D5 chiếm ưu thế              Trạng thái trung gian phối trí 5        |
| Tại điểm 30 GPa     Tồn tại đồng thời D5 & D6             Chuyển hóa mạnh mẽ sang bát diện SiO6   |
| 30 - 45 GPa         1 Domain D6 hoàn chỉnh                Bát diện SiO6 chiếm ưu thế (~100%)      |
+--------------------------------------------------------------------------------------------------+

Bốn phát hiện khoa học mang tính đột phá của luận án bao gồm:

  1. Cấu trúc Lõi - Vỏ và Cơ chế Tinh thể hóa Đa bước của Hạt Nano Fe: Hạt nano Fe vô định hình phân hóa rõ rệt thành hai vùng: lõi đặc khít có cấu trúc tương tự vật liệu khối và vỏ ngoài xốp có độ dày khoảng $8\text{ \AA}$. Quá trình ủ nhiệt tại $750 - 900\text{ K}$ cho thấy sự xuất hiện của các mầm tinh thể bcc tại vùng lõi. Tốc độ phát triển đám tinh thể tuân theo hàm mũ theo thời gian; hình thái đám tinh thể có xu hướng tiến dần đến hình cầu nhằm cực tiểu hóa năng lượng tự do bề mặt $\Delta G_S$.
  2. Vai trò Kìm hãm Kết tinh của Nguyên tử Bo trong $\text{Fe}x\text{B}{100-x}$: Khi ủ nhiệt mẫu $\text{Fe}{95}\text{B}5$ và $\text{Fe}{90}\text{B}{10}$ tại $900\text{ K}$, các nguyên tử Bo bị đẩy ra lớp vỏ cầu bề mặt và các ranh giới biên giữa pha vô định hình ($\text{Am}$) và pha tinh thể ($\text{Cr}$). Nồng độ Bo tại bề mặt tăng cao đóng vai trò như một lớp màng ngăn khuếch tán, làm tăng thời gian sống của mầm tới hạn $\langle t_{nucleus} \rangle$ và ức chế tốc độ kết tinh toàn phần của hạt nano.
  3. Phát hiện Vùng Vi mô Tinh khiết và Động học Không đồng nhất trong $\text{SiO}_2$ Lỏng: Lần đầu tiên chứng minh bằng mô phỏng quy mô $20.000$ nguyên tử rằng chất lỏng $\text{SiO}2$ ở $3.000\text{ K}$ không đồng nhất về mặt hóa học mà tồn tại các cụm hạt lõi - vỏ thuần nhất chỉ chứa riêng biệt nguyên tử Si hoặc O. Hệ phân rã thành hai tập hợp nguyên tử: nhóm linh động ($\text{MS}$) và nhóm không linh động ($\text{IMS}$). Nhóm $\text{IMS}$ liên kết thành mạng con bền vững ($SSLA$), duy trì liên kết trong thời gian quan sát dài ($t{obs} > 0,4\text{ ns}$), làm suy giảm tốc độ khuếch tán tổng cộng khi nhiệt độ hạ thấp.
  4. Cơ chế Chuyển pha Domain $D_4 \to D_5 \to D_6$ và Quy luật Đậm đặc hóa Dưới Áp suất Cao: Khảo sát hệ $\text{SiO}2$ lỏng ($3.500\text{ K}$, $0 - 45\text{ GPa}$) và vô định hình ($500\text{ K}$, $0 - 100\text{ GPa}$) phát hiện cấu trúc 5 domain. Dưới tác dụng của áp suất nén, góc liên kết $\theta{\text{Si-O-Si}}$ giảm mạnh từ $146^\circ$ ($0\text{ GPa}$) xuống $105^\circ$ ($45\text{ GPa}$) đi kèm sự tăng độ dài liên kết trung bình $r_{\text{Si-O}}$ từ $1,60\text{ \AA}$ lên $1,72\text{ \AA}$ do lực đẩy Coulomb giữa các ion O lân cận tăng lên khi số phối trí tăng từ 4 lên 6. Quá trình đậm đặc hóa không gian thực chất là sự thu hẹp thể tích đa diện Voronoi của các miền biên ($\text{BD}$) và sự sáp nhập các miền $D_x$.

Implications đa chiều

  • Về mặt Lý thuyết: Cung cấp bức tranh toàn diện về chuyển dịch cấu trúc pha phi cổ điển, hoàn thiện lý thuyết mạng ngẫu nhiên và giải thích thỏa đáng các dị thường động học chất lỏng silicat.
  • Về Khoa học Địa cầu và Địa vật lý: Giải mã trạng thái cấu trúc của macma silicat và biến chất khoáng vật $\text{SiO}_2$ (từ thạch anh $\to$ coesite $\to$ stishovite) tại điều kiện áp suất - nhiệt độ cực hạn sâu trong lớp manti Trái Đất, cung cấp thông số đầu vào chuẩn xác cho các mô hình địa chấn.
  • Về Ứng dụng Công nghệ Vật liệu:
    • Tối ưu hóa quy trình ủ nhiệt kiểm soát kích thước hạt nano Fe nhằm nâng cao hoạt tính xúc tác xử lý môi trường (khử độc nước ngầm nhiễm kim loại nặng $\text{Pb}^{2+}, \text{Ni}^{2+}, \text{Hg}^{2+}, \text{U}$).
    • Định hướng công nghệ chế tạo hợp kim từ mềm FeB vô định hình - tinh thể dùng cho lõi biến áp hiệu suất cao, cảm biến từ và vật liệu lưu trữ hydro.

Limitations và Future Research

Luận án thừa nhận một cách khách quan 3 hạn chế nội tại:

  1. Giới hạn của thế tương tác cổ điển: Thế cặp thực nghiệm Pak-Doyama và BKS dù tái tạo xuất sắc cấu trúc hình học và năng lượng cơ học nhưng không mô tả được trạng thái kích thích điện tử, liên kết hóa trị định hướng phi cổ điển hay hiện tượng phân cực ion bậc cao.
  2. Thang thời gian mô phỏng (Time-scale limit): Giới hạn tính toán của phương pháp MD cổ điển nằm trong khoảng hàng chục nanogiây ($ns$), chưa thể bao quát trọn vẹn quá trình ủ nhiệt thực nghiệm kéo dài nhiều giờ ở nhiệt độ thấp gần điểm chuyển pha thủy tinh $T_g$.
  3. Ảnh hưởng của môi trường hóa học bề mặt: Mô hình hạt nano Fe và FeB được khảo sát trong điều kiện chân không lý tưởng, chưa tính đến ảnh hưởng của lớp màng oxit thụ động hóa hoặc các phân tử dung môi phân cực trong dung dịch thực tế.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Hướng 1: Ứng dụng mô phỏng lượng tử Ab-initio Molecular Dynamics (AIMD) hoặc thế tương tác máy học (Machine Learning Potentials - MLP/DeepMD) để mô tả chính xác hơn các trạng thái liên kết lai hóa tại vùng chuyển tiếp $D_5$.
  • Hướng 2: Nghiên cứu hệ hợp kim từ tính ba thành phần $\text{Fe-B-Si}$ và $\text{Fe-Co-B}$ với kích thước hạt mở rộng trên $50.000$ nguyên tử ($D > 10\text{ nm}$).
  • Hướng 3: Khảo sát quá trình chuyển pha của $\text{SiO}_2$ trong điều kiện kết hợp đồng thời áp suất siêu cao ($P > 130\text{ GPa}$) và nhiệt độ siêu cao ($T > 4.000\text{ K}$) mô phỏng ranh giới lõi - manti Trái Đất.
  • Hướng 4: Mô phỏng động học solvat hóa và khả năng hấp phụ kim loại nặng của hạt nano Fe trong môi trường nước ngầm có độ pH thay đổi.

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

Các kết quả nghiên cứu của luận án đã đóng góp trực tiếp vào kho tàng tri thức vật lý tính toán quốc tế thông qua 08 công trình công bố, bao gồm các bài báo trên các tạp chí chuyên ngành hàng đầu:

  • Journal of Non-Crystalline Solids (IF ~ 3.5 - 4.0, Q1): Khẳng định cấu trúc domain và cơ chế đậm đặc hóa của $\text{SiO}_2$.
  • The European Physical Journal B (Springer): Công bố động học không đồng nhất và mạng con bền vững $SSLA$.
  • Materials Research Express & AIP Advances: Động học kết tinh hạt nano Fe và hợp kim FeB.
  • Journal of Physics: Conference Series (IOP Publishing): 02 bài báo hội nghị quốc tế.

Về mặt công nghệ và xã hội, các thông số mô phỏng từ luận án cung cấp cơ sở dữ liệu vật lý quan trọng cho các ngành công nghiệp luyện kim bột nano, công nghệ sản xuất sợi quang học cáp viễn thông, gốm kỹ thuật điện tử và vật liệu lưu trữ năng lượng tương thích sinh học.

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

                                  ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP
+-------------------------------+------------------------------------------------------------------+
| Nhóm đối tượng                | Giá trị thụ hưởng cụ thể                                         |
+-------------------------------+------------------------------------------------------------------+
| Nghiên cứu sinh & Nhà khoa học| Khung phân tích domain Dx, chỉ số <fx>, phương pháp hạt lõi-vỏ   |
| (Vật lý chất rắn / Tính toán) | để giải mã các bài toán cấu trúc vô định hình phức tạp.          |
+-------------------------------+------------------------------------------------------------------+
| Kỹ sư R&D Công nghệ Nano      | Điều kiện nhiệt độ ủ tối ưu (750 - 900 K) để tạo cấu trúc hạt    |
| & Luyện kim bột               | nano Fe/FeB có từ tính kháng từ cao hoặc hoạt tính xúc tác mạnh. |
+-------------------------------+------------------------------------------------------------------+
| Nhà nghiên cứu Địa chất       | Dữ liệu chuyển đổi số phối trí và thể tích Voronoi của SiO2 phục |
| & Khoa học Trái Đất           | vụ mô hình hóa macma và giải mã truyền sóng địa chấn manti sâu.  |
+-------------------------------+------------------------------------------------------------------+
| Kỹ sư Vật liệu Quang điện tử  | Hiểu biết về khuyết tật nút mạng và mạch vòng SiO2 giúp nâng     |
| & Gốm kỹ thuật                | cao độ bền cơ học và độ truyền dẫn quang của thủy tinh thạch anh.|
+-------------------------------+------------------------------------------------------------------+

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 học thuyết nào?

Đóng góp lý thuyết độc đáo nhất là Mô hình 5 Domain cấu trúc ($D_4, D_5, D_6$) và Khái niệm Vùng vi mô tinh khiết trong chất lỏng có cấu trúc mạng, mở rộng trực tiếp Mô hình Mạng ngẫu nhiên liên tục (CRN) của Zachariasen [57]. Luận án chứng minh rằng sự chuyển đổi số phối trí của Si dưới áp suất không diễn ra đơn lẻ mà theo cơ chế phân tách và hợp nhất các trường không gian domain, trong đó các domain biên đóng vai trò vùng đệm mật độ thấp.

2. Phương pháp nghiên cứu có cải tiến gì so với các nghiên cứu quốc tế trước đây?

So với các mô phỏng của Xiaohua Li et al. (2014) [19] ($N = 331$ nguyên tử) và Belashchenko (1989) [24] ($N = 686$ nguyên tử), luận án đã nâng quy mô tính toán lên $5.000 - 10.000$ nguyên tử (đối với hạt nano Fe/FeB) và $20.000$ nguyên tử (đối với $\text{SiO}_2$ lỏng). Đồng thời, luận án phát triển bộ công cụ phân tích cấu trúc phối hợp độc quyền gồm: Tần số cấu trúc động học $\langle f_x \rangle$, phân tích đa diện Voronoi, thuật toán SPFA và phân tích hạt lõi - vỏ.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu định lượng tương ứng?

Phát hiện bất ngờ nhất là sự tồn tại của các vùng vi mô tinh khiết (chỉ chứa riêng biệt nguyên tử Si hoặc O) trong $\text{SiO}_2$ lỏng ở $3.000\text{ K}$. Trước đây, $\text{SiO}2$ lỏng luôn được giả định là dung dịch đồng nhất lý tưởng ở mức phân tử. Dữ liệu phân tích cụm hạt lõi - vỏ chỉ ra các mạng con Si-O bền vững ($SSLA$) có thời gian sống vượt trội, và khi nén áp suất lên $30\text{ GPa}$, khoảng cách liên kết $r{\text{Si-O}}$ không hề co lại mà bị kéo dài nghịch lý từ $1,60\text{ \AA}$ lên $1,72\text{ \AA}$ do lực đẩy tĩnh điện giữa các anion oxy trong khối bát diện $\text{SiO}_6$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) chi tiết không?

Có. Toàn bộ thông số thế tương tác thực nghiệm (bảng tham số Pak-Doyama cho Fe-Fe, Fe-B, B-B; bảng hằng số BKS cho Si-O, O-O), các bước tích phân thời gian Verlet ($dt = 4,6 \times 10^{-14}\text{ s}$ và $dt = 1,0\text{ fs}$), quy trình ủ nhiệt đẳng nhiệt NPT/NVT và tiêu chuẩn phân loại mạng tinh thể qua CNA đều được mô tả chi tiết và công bố minh bạch trong Chương 2 của luận án.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình nghiên cứu 10 năm tập trung vào: (1) Ứng dụng trí tuệ nhân tạo xây dựng thế tương tác Machine Learning Potentials mô phỏng hệ oxit phức tạp ($\text{MgSiO}_3, \text{Al}_2\text{O}_3-\text{SiO}_2$); (2) Khảo sát thực nghiệm in-situ kết hợp mô phỏng đa quy mô cho hạt nano hợp kim từ tính sinh học Fe-B-Si; (3) Mở rộng thang mô phỏng lên hàng triệu nguyên tử nhằm giải mã các ranh giới hạt và vết nứt tế vi trong vật liệu gốm áp suất cực hạn.

Kết luận

Luận án tiến sĩ của NCS. Giáp Thị Thùy Trang đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 kết luận cốt lõi:

  1. Xác lập mô hình cấu trúc lõi - vỏ hoàn chỉnh cho hạt nano kim loại Fe và FeB vô định hình, chứng minh quá trình kết tinh khởi phát từ lõi đặc và tiến hóa qua các cấu trúc trung gian giả bền bcc và icosahedral.
  2. Khám phá cơ chế cản trở tinh thể hóa của nguyên tử Bo: Bo tập trung tại bề mặt và biên pha vô định hình - tinh thể, làm gia tăng năng lượng hoạt hóa khuếch tán và kìm hãm sự phát triển của mầm tinh thể Fe.
  3. Phát hiện tính không đồng nhất hóa học và động học trong $\text{SiO}_2$ lỏng: Tồn tại các vùng vi mô tinh khiết và mạng con nguyên tử bền vững ($SSLA$), làm thay đổi bản chất của các mô hình lý thuyết về chất lỏng silicat.
  4. Giải mã trọn vẹn cơ chế chuyển pha 5 domain $D_4 \to D_5 \to D_6$ dưới áp suất cao: Xác định quy luật biến đổi góc liên kết $\theta_{\text{Si-O-Si}}$, độ dài liên kết $r_{\text{Si-O}}$ và thể tích Voronoi từ áp suất khí quyển đến $100\text{ GPa}$.
  5. Đóng góp hệ thống phương pháp luận tính toán tiên tiến: Ứng dụng thành công tần số cấu trúc động học $\langle f_x \rangle$, phân tích domain và thuật toán hạt lõi - vỏ, mở ra hướng tiếp cận chuẩn mực cho vật lý chất rắn và khoa học vật liệu tính toán tại Việt Nam và quốc tế.