Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Các vật liệu ôxit hai nguyên như nhôm ôxit ($\text{Al}_2\text{O}_3$), giécmani ôxit ($\text{GeO}_2$) và silic đioxit ($\text{SiO}_2$) giữ vai trò thiết yếu trong khoa học kỹ thuật và đời sống, được ứng dụng rộng rãi trong các ngành luyện kim, cơ khí, gốm sứ, cao su, xi măng, sợi quang học và công nghệ thủy tinh. Việc nghiên cứu cấu trúc vi mô và cơ tính của các hệ vật liệu này có ý nghĩa khoa học và thực tiễn cao, giúp định hướng thiết kế và chế tạo vật liệu sở hữu các đặc tính cơ học đáp ứng yêu cầu kỹ thuật.

Trong các nghiên cứu trước đây, cấu trúc của $\text{Al}_2\text{O}_3$, $\text{GeO}_2$ và $\text{SiO}_2$ ở pha lỏng và pha vô định hình (VĐH)/thủy tinh tại điều kiện chuẩn đều được tạo thành từ các đơn vị cấu trúc cơ bản $\text{TO}_x$ ($\text{T} = \text{Al}, \text{Ge}, \text{Si}$; $x = 4, 5, 6$). Dưới tác dụng của áp suất nén, quá trình chuyển đổi từ cấu trúc tứ diện $\text{TO}_4$ sang cấu trúc bát diện $\text{TO}_6$ thông qua trạng thái trung gian $\text{TO}_5$ đã được ghi nhận. Tuy nhiên, các công trình đã công bố vẫn tồn tại nhiều khoảng trống và vấn đề chưa được thống nhất:

  1. Cơ chế biến đổi cấu trúc trật tự gần (gắn với độ dài liên kết và hình học đa diện) và trật tự tầm trung (gắn với liên kết nguyên tử oxy cầu) dưới tác dụng của áp suất cao chưa có sự giải thích thỏa đáng.
  2. Một số hiện tượng bất thường trong hàm phân bố xuyên tâm (PBXT) cặp khi chịu nén như hiện tượng tách đỉnh chính của cặp $\text{Si-Si}$/$\text{Ge-Ge}$ hoặc sự xuất hiện đỉnh phụ của cặp $\text{O-O}$ chưa được chứng minh định lượng.
  3. Trật tự sắp xếp không gian của các nguyên tử oxy trong quá trình nén áp suất cao chưa được nhận diện cụ thể.
  4. Các nghiên cứu chuyên sâu về mối liên hệ giữa cấu trúc vi mô, sự hình thành/phân bố của các quả cầu lỗ hổng (voids) với ứng xử cơ học (đặc biệt là biến dạng đơn trục, dòng chảy dẻo và quá trình gãy đứt) của các hệ ôxit này còn rất hạn chế.

Từ những lý do trên, đề tài luận án "Mô phỏng cấu trúc và cơ tính của hệ ôxit $\text{Al}_2\text{O}_3$, $\text{GeO}_2$, $\text{SiO}_2$" được thực hiện nhằm giải quyết các vấn đề tồn tại nêu trên.

Mục tiêu nghiên cứu

Luận án xác định ba mục tiêu nghiên cứu cụ thể:

  1. Nghiên cứu ảnh hưởng của áp suất nén lên sự thay đổi cấu trúc và phân bố các quả cầu lỗ hổng trong các hệ $\text{Al}_2\text{O}_3$ VĐH, $\text{SiO}_2$ và $\text{GeO}_2$ thủy tinh.
  2. Khảo sát sự phụ thuộc vào áp suất nén của các đặc trưng cơ tính (đặc biệt là mô-đun I-âng) của ba hệ ôxit.
  3. Phân tích chi tiết ảnh hưởng của quá trình biến dạng đơn trục lên cấu trúc vi mô và sự tiến triển của các quả cầu lỗ hổng trong hệ $\text{SiO}_2$ thủy tinh nhằm làm rõ cơ chế biến dạng và phá hủy.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Mạng cấu trúc vi mô và tính chất cơ học của ba hệ ôxit hai nguyên: $\text{Al}_2\text{O}_3$ ở pha vô định hình (VĐH), $\text{GeO}_2$ và $\text{SiO}_2$ ở pha thủy tinh.
  • Phạm vi không gian và điều kiện mô phỏng:
    • Khảo sát mô hình nguyên tử trong ô mô phỏng lập phương có điều kiện biên tuần hoàn với quy mô từ 5.000 đến 6.000 nguyên tử.
    • Dải áp suất nén khảo sát từ $0\text{ GPa}$ đến $90\text{ GPa}$ (đối với $\text{SiO}_2$, mật độ mẫu khảo sát tương đương trong dải từ $2,20\text{ g/cm}^3$ đến $4,42\text{ g/cm}^3$).
    • Nhiệt độ khảo sát cơ tính và phân tích cấu trúc cân bằng được cố định tại $300\text{ K}$.
    • Chế độ biến dạng: Biến dạng kéo đơn trục dọc theo trục $Oz$ ở tốc độ không đổi.

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 đã tổng quan 97 tài liệu khoa học liên quan đến cấu trúc và cơ tính của các hệ ôxit:

  • Đối với hệ $\text{Al}_2\text{O}_3$: Các nghiên cứu thực nghiệm và mô phỏng của Landron và cộng sự, San Miguel, Petkov ghi nhận cấu trúc mạng tứ diện ở áp suất thường với khoảng cách $\text{Al-O} \approx 1,8\text{ \AA}$ và góc $\text{O-Al-O} \approx 109,5^\circ$. Các công trình của Gutierrez, Ahuja, và Adjaoud chỉ ra sự chuyển đổi phối trí từ 4 lên 6 khi tăng mật độ. Nghiên cứu mô phỏng ab initio của Adjaoud cho hệ lỏng ghi nhận sự xuất hiện đỉnh phụ trong hàm PBXT cặp $\text{O-O}$ ở áp suất trên $10\text{ GPa}$ nhưng chưa giải thích bản chất. Wenhui Duan và cộng sự khảo sát độ dài liên kết của khối bát diện trong các pha tinh thể $\alpha\text{-Al}_2\text{O}_3$ và dạng $\text{Rh}_2\text{O}_3\text{(III)}$. Về cơ tính, các phép đo thực nghiệm của màng mỏng $\text{Al}_2\text{O}_3$ ghi nhận mô-đun I-âng từ $122\text{ GPa}$ đến $140\text{ GPa}$, trong khi mô phỏng ab initio và ĐLHPT cho mẫu khối đạt giá trị khoảng $340,3\text{ GPa}$.
  • Đối với hệ $\text{SiO}_2$ và $\text{GeO}_2$: Cấu trúc tinh thể và chuyển pha áp suất cao được khảo sát qua các công trình của Hemley, Sato, Funamori, Micoulaut, Salmon. Tinh thể $\text{SiO}_2$ chuyển từ cấu trúc tứ diện sang bát diện qua các pha quartz, coesite, stishovite ($>8\text{ GPa}$), $\text{CaCl}_2\text{-type}$ ($60\text{ GPa}$), $\alpha\text{-PbO}_2\text{-type}$ ($120\text{ GPa}$), và pyrite-type ($260\text{ GPa}$). Tinh thể $\text{GeO}_2$ diễn ra sự chuyển pha tương tự nhưng ở ngưỡng áp suất thấp hơn ($1,8\text{--}2,2\text{ GPa}$ chuyển sang rutile; $19\text{--}25\text{ GPa}$ sang $\text{CaCl}_2\text{-type}$; $59\text{--}66\text{ GPa}$ sang pyrite-type). Đối với pha thủy tinh, các nghiên cứu thực nghiệm tán xạ nơtron và tia X của Wright, Johnson, Sampath cho thấy góc liên kết $\text{Si-O-Si} \approx 144^\circ$ và $\text{Ge-O-Ge} \approx 130^\circ$. Hiện tượng tách đỉnh trong hàm PBXT cặp $\text{Ge-Ge}$ được quan sát bởi Bras, Guthrie, Drewitt, còn với $\text{Si-Si}$ được ghi nhận bởi Meade và Devine nhưng chưa được định lượng cơ chế vi mô.
  • Về cơ tính và cơ chế biến dạng của thủy tinh: Nghiên cứu của Vandembroucq, Rottler, Robbins, Tanguy, Yuan, Huang chỉ ra ứng xử cơ học dị thường của thủy tinh silica quanh $8\text{--}10\text{ GPa}$ gắn với sự chuyển từ biến dạng đàn hồi sang biến dạng dẻo và vai trò của đơn vị khuyết tật $\text{SiO}_5$. Nghiên cứu của Ochoa, Swiler, Rountree, Muralidharan phân tích đường cong ứng suất - biến dạng gồm bốn vùng và chỉ ra sự phát triển của các lỗ hổng kích thước nano. Công trình của P. Hung và cộng sự phân tích sự tạo cụm lỗ hổng trong mạng ôxit.

Khoảng trống luận án lựa chọn giải quyết

  1. Luận án giải quyết triệt để nguyên nhân vi mô gây ra hiện tượng tách đỉnh trong hàm PBXT cặp $\text{Si-Si}$/$\text{Ge-Ge}$ và sự xuất hiện đỉnh phụ trong hàm $\text{O-O}$ thông qua phân tích đóng góp của các liên kết chung cạnh và chung mặt.
  2. Nhận diện cấu trúc sắp xếp trật tự tầm trung của các nguyên tử oxy ở áp suất cao dưới dạng các mầm tinh thể cục bộ.
  3. Xác định có hệ thống mật độ riêng phần của từng loại đơn vị cấu trúc $\text{TO}_x$ và thiết lập công thức biểu diễn mật độ mẫu toàn phần.
  4. Làm rõ vai trò kích hoạt sự phá hủy của các quả cầu lỗ hổng lân cận oxy ($\text{O-void}$) trong quá trình biến dạng đơn trục của hệ $\text{SiO}_2$ thủy tinh.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và mô hình tương tác

Luận án ứng dụng phương pháp mô phỏng Động lực học phân tử (ĐLHPT) cổ điển. Tương tác liên nguyên tử giữa các cặp ion $i$ và $j$ được mô tả bằng hàm thế cặp có dạng tổng quát:

$$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}$$

Trong đó, số hạng thứ nhất biểu diễn tương tác tĩnh điện Coulomb giữa các ion có điện tích hiệu dụng $q_i, q_j$ (tính toán bằng phương pháp tổng Ewald); số hạng thứ hai mô tả tương tác đẩy tầm ngắn do sự xen phủ đám mây điện tử; số hạng thứ ba mô tả tương tác phân cực cảm ứng Van der Waals.

Hệ vật liệu Thế tương tác sử dụng Cặp nguyên tử $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{Al}_2\text{O}_3$ Coulomb-Buckingham $\text{Al-O}$
$\text{O-O}$		 $12201,42$
$1844,75$		 $5,11247$
$2,91036$		 $32,00$
$192,58$		 $q_{\text{Al}} = +1,8$
$q_{\text{O}} = -1,2$
$\text{GeO}_2$ Oeffner - Elliot $\text{Ge-O}$
$\text{O-O}$		 $208006,8549$
$7693,3496$		 $6,12933$
$3,28511$		 $236,6475$
$131,0874$		 $q_{\text{Ge}} = +1,5$
$q_{\text{O}} = -0,75$
$\text{SiO}_2$ BKS (van Beest-Kramer-van Santen) $\text{Si-O}$
$\text{O-O}$		 $18003,7572$
$1388,7730$		 $4,87318$
$2,76000$		 $1133,5381$
$175,0000$		 $q_{\text{Si}} = +2,4$
$q_{\text{O}} = -1,2$

Quy trình xây dựng mô hình ĐLHPT

Phương trình chuyển động được tích phân bằng thuật toán Verlet với các thông số:

  • Hệ $\text{Al}_2\text{O}_3$: Gồm $5.000$ nguyên tử ($2.000\text{ Al} + 3.000\text{ O}$). Bước thời gian $\Delta t = 1\text{ fs}$, mô hình tập hợp đẳng áp - đẳng nhiệt ($NPT$). Quá trình: gieo hạt ban đầu ở mật độ $2,8\text{ g/cm}^3$, gia nhiệt đến $5000\text{ K}$ ($5\cdot 10^5$ bước), làm nguội về $3500\text{ K}$ (tốc độ $1\text{ K/ps}$), nén đẳng nhiệt tại $3500\text{ K}$ đến các áp suất từ $0$ đến $90\text{ GPa}$ ($10^6$ bước), sau đó làm nguội về $300\text{ K}$ (tốc độ $1\text{ K/ps}$) và hồi phục trong $10^5$ bước.
  • Hệ $\text{GeO}_2$: Gồm $6.000$ nguyên tử ($2.000\text{ Ge} + 4.000\text{ O}$). Bước thời gian $\Delta t = 1,5\text{ fs}$, mô hình $NPT$. Mật độ ban đầu $3,5\text{ g/cm}^3$, quy trình nhiệt tương tự hệ $\text{Al}_2\text{O}_3$, xây dựng 4 mô hình trong khoảng $0\text{--}90\text{ GPa}$ ở $300\text{ K}$.
  • Hệ $\text{SiO}_2$: Gồm $6.000$ nguyên tử ($2.000\text{ Si} + 4.000\text{ O}$). Bước thời gian $\Delta t = 1\text{ fs}$, mô hình tập hợp đẳng tích - đẳng nhiệt ($NVT$). Mật độ ban đầu $2,2\text{ g/cm}^3$, nung lên $5000\text{ K}$, làm nguội về $3500\text{ K}$, nén để thu được các mẫu có mật độ từ $2,35\text{ g/cm}^3$ đến $4,42\text{ g/cm}^3$, sau đó làm nguội xuống $300\text{ K}$ với tốc độ $10^{13}\text{ K/s}$.
  • Dữ liệu thống kê được lấy trung bình trên 50 cấu hình độc lập cách nhau 100 bước trong giai đoạn cân bằng 5000 bước ĐLHPT tại $300\text{ K}$.

Các phương pháp phân tích cấu trúc và cơ tính

  • Hàm phân bố xuyên tâm (PBXT) toàn phần $g(r)$ và cặp $g_{\alpha\beta}(r)$: Dùng để xác định trật tự gần, độ dài liên kết và kiểm chứng với dữ liệu thực nghiệm.
  • Số phối trí và đơn vị cấu trúc $\text{TO}_x$: Xác định qua bán kính ngắt $r_c$ tại vị trí cực tiểu thứ nhất của $g_{\text{T-O}}(r)$. Mật độ riêng phần $\rho_x$ của các đơn vị $\text{TO}_x$ được tính bằng phương pháp tạo ngẫu nhiên $10^6$ điểm trong ô mô phỏng.
  • Phân bố góc liên kết: Khảo sát góc trong $\text{O-T-O}$ và góc cầu $\text{T-O-T}$ giữa các đa diện.
  • Phân bố quả cầu lỗ hổng (voids): Xác định các hình cầu rỗng cực đại tiếp xúc với 4 nguyên tử xung quanh và không chứa nguyên tử nào. Bán kính cứng quy ước của các nguyên tử: $r_{\text{Al}} = 1,23\text{ \AA}$, $r_{\text{Ge}} = 1,52\text{ \AA}$, $r_{\text{Si}} = 1,10\text{ \AA}$, $r_{\text{O}} = 0,73\text{ \AA}$. Thể tích lỗ hổng $V_{\text{void}}$ tính qua thuật toán Monte Carlo với $10^6$ điểm thử.
  • Phương pháp phân tích lân cận chung (CNA): Nhận diện cấu trúc trật tự cục bộ dựa trên bộ chỉ số $kmn$:
    • Cấu trúc lập phương tâm mặt ($\text{fcc}$): 12 bộ $421$.
    • Cấu trúc lục giác xếp chặt ($\text{hcp}$): 6 bộ $421$ và 6 bộ $422$.
    • Cấu trúc lập phương tâm khối ($\text{bcc}$): 6 bộ $444$ và 8 bộ $666$.
    • Cấu trúc 20 mặt (icosahedra): 12 bộ $555$.
  • Kỹ thuật trực quan hóa: Xây dựng trên nền tảng phần mềm MATLAB để hiển thị sự phân bố không gian của các đơn vị cấu trúc, mầm tinh thể và các quả cầu lỗ hổng.
  • Phương pháp mô phỏng biến dạng đơn trục: Kéo dãn ô mô phỏng dọc trục $Oz$ theo thời gian $t$ với độ biến dạng $\varepsilon = [L_z(t) - L_z(0)] / L_z(0)$. Ten-xơ ứng suất virial $\sigma_{\alpha\beta}$ được xác định từ động năng và lực tương tác liên nguyên tử:

$$\sigma_{\alpha\beta} = \frac{1}{V} \sum_{i=1}^N \left[ m_i v_{\alpha i} v_{\beta i} + \frac{1}{2} \sum_{j \neq i} F_{ij} \frac{r_{ij,\alpha} r_{ij,\beta}}{r_{ij}} \right]$$

Từ đường cong ứng suất - biến dạng ($\sigma\text{--}\varepsilon$), xác định mô-đun đàn hồi Young ($E$), ứng suất đàn hồi, ứng suất chảy và phân chia các giai đoạn biến dạng.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 hệ thống hóa các kiến thức lý thuyết và thực nghiệm về cấu trúc, quá trình chuyển pha và cơ tính của ba hệ ôxit hai nguyên:

  • Hệ $\text{Al}_2\text{O}_3$: Trình bày các pha thù hình gồm pha bền $\alpha\text{-Al}_2\text{O}_3$ (cấu trúc corundum, mạng oxy $\text{hcp}$, $\text{Al}$ ở vị trí bát diện) và các pha không bền ($\gamma, \eta, \delta, \theta$ có mạng oxy $\text{fcc}$; $\kappa, \chi$ có mạng oxy $\text{hcp}$). Phân tích các nghiên cứu về pha lỏng và pha vô định hình dưới áp suất nén, sự thay đổi độ dài liên kết $\text{Al-O}$, $\text{Al-Al}$, $\text{O-O}$, sự chuyển đổi từ liên kết chung góc sang liên kết chung cạnh/chung mặt và hiện tượng xuất hiện đỉnh phụ trong hàm PBXT cặp $\text{O-O}$.
  • Hệ $\text{SiO}_2$ và $\text{GeO}_2$: Điểm lại các dạng thù hình tinh thể ở áp suất thường và áp suất cao (quartz, tridymite, crystobalite, coesite, stishovite, $\text{CaCl}_2\text{-type}$, $\alpha\text{-PbO}_2\text{-type}$, pyrite-type). Trình bày cấu trúc mạng tứ diện trong pha thủy tinh ở áp suất thường và cơ chế chuyển đổi liên tục sang mạng bát diện khi chịu nén. Nêu rõ các phát hiện thực nghiệm và mô phỏng về hiện tượng tách đỉnh hàm PBXT của cặp cation ($\text{Si-Si}$, $\text{Ge-Ge}$) và sự biến đổi của trật tự tầm trung.
  • Cơ tính của các hệ ôxit: Tổng quan các phép đo mô-đun Young trên màng mỏng và vật liệu khối $\text{Al}_2\text{O}_3$; tính chất đàn hồi phụ thuộc áp suất của $\text{GeO}_2$; các hiện tượng cơ học dị thường của thủy tinh $\text{SiO}_2$ ở dải áp suất $8\text{--}10\text{ GPa}$, vai trò của đơn vị khuyết tật $\text{SiO}_5$ làm tăng độ dẻo. Tổng hợp các mô hình biến dạng kéo/cắt, quá trình hình thành mầm nứt gãy và 4 vùng đặc trưng trên đường cong ứng suất - biến dạng.

Chương 2: Phương pháp tính toán

Chương 2 mô tả chi tiết quy trình tính toán mô phỏng ĐLHPT và các thuật toán phân tích:

  • Trình bày dạng toán học và bộ thông số của ba hàm thế cặp tương tác: Coulomb-Buckingham ($\text{Al}_2\text{O}_3$), Oeffner-Elliot ($\text{GeO}_2$), BKS ($\text{SiO}_2$).
  • Thiết lập quy trình tạo mẫu cân bằng nhiệt động qua các bước nung chảy ($5000\text{ K}$), làm nguội ($3500\text{ K}$), nén đẳng nhiệt/đẳng tích và làm nguội về $300\text{ K}$ cho từng hệ.
  • Trình bày công thức giải tích xác định hàm PBXT $g(r)$, số phối trí $Z_{\alpha\beta}$, phân bố góc liên kết và thuật toán Monte Carlo tính thể tích lỗ hổng $V_{\text{void}}$.
  • Trình bày thuật toán phân tích mật độ riêng phần của các đơn vị cấu trúc $\text{TO}_4$, $\text{TO}_5$, $\text{TO}_6$ dựa trên sự phân chia thể tích giao thoa giữa các hình cầu phối trí.
  • Mô tả chi tiết phương pháp phân tích lân cận chung (CNA) để nhận diện các pha tinh thể $\text{fcc}$, $\text{hcp}$, $\text{bcc}$, icosahedra và kỹ thuật trực quan hóa cấu trúc bằng MATLAB.
  • Thiết lập phương pháp mô phỏng biến dạng kéo đơn trục và công thức tính ten-xơ ứng suất virial.

Chương 3: Ảnh hưởng của áp suất lên cấu trúc của các hệ ôxit $\text{Al}_2\text{O}_3$ VĐH, $\text{GeO}_2$ và $\text{SiO}_2$ thủy tinh

Chương 3 trình bày các kết quả mô phỏng cấu trúc của ba hệ ôxit khi chịu nén:

        Áp suất nén tăng (0 GPa -> 90 GPa / Mật độ tăng)
                             |
    +------------------------+------------------------+
    |                                                 |
[Trật tự gần]                                  [Trật tự tầm trung]
- TO4 -> TO5 -> TO6                            - Hình thành liên kết cạnh và mặt
- Khoảng cách T-O tăng nhẹ rồi giảm            - PBXT T-T tách thành 2 đỉnh
- Góc O-T-O dịch chuyển về góc bát diện        - PBXT O-O xuất hiện đỉnh phụ thứ hai
                                               - Xuất hiện mầm tinh thể O (fcc, hcp)
                                               - Thể tích lỗ hổng Vvoid/V giảm mạnh
  • Hệ $\text{Al}_2\text{O}_3$ VĐH:
    • Ở $0\text{ GPa}$, khoảng cách liên kết trung bình $\text{Al-O} = 1,80 \pm 0,01\text{ \AA}$, $\text{Al-Al} = 3,20 \pm 0,01\text{ \AA}$, $\text{O-O} = 2,82 \pm 0,01\text{ \AA}$. Khi áp suất tăng, đỉnh $\text{Al-O}$ dịch nhẹ sang phải đến $1,82 \pm 0,01\text{ \AA}$, đỉnh $\text{Al-Al}$ dịch về $2,80 \pm 0,01\text{ \AA}$, đỉnh $\text{O-O}$ dịch về $2,54 \pm 0,01\text{ \AA}$.
    • Ở áp suất $P > 15\text{ GPa}$, hàm $g_{\text{O-O}}(r)$ xuất hiện đỉnh phụ tại vị trí $3,68 \pm 0,06\text{ \AA}$.
    • Tỉ phần đơn vị $\text{AlO}_4$ giảm từ $67%$ ($0\text{ GPa}$) xuống $18%$ ($22,5\text{ GPa}$); đơn vị trung gian $\text{AlO}_5$ tăng và đạt cực đại $57%$ tại $22,5\text{ GPa}$ rồi giảm dần; đơn vị $\text{AlO}_6$ tăng liên tục theo áp suất.
    • Thiết lập công thức liên hệ mật độ mẫu toàn phần theo mật độ riêng phần và tỉ phần của từng đơn vị: $\rho \approx \rho_4 \text{Al}_4 + \rho_5 \text{Al}_5 + \rho_6 \text{Al}_6$ với độ lệch so với mô phỏng nhỏ hơn $2,3%$.
    • Nhận diện các đám mầm tinh thể oxy có cấu trúc $\text{fcc}$ và $\text{hcp}$ gồm vài chục nguyên tử ở áp suất cao ($60\text{ GPa}$).
  • Hệ $\text{GeO}_2$ thủy tinh:
    • Khoảng cách liên kết ở $0\text{ GPa}$: $r_{\text{Ge-O}} = 1,75 \pm 0,01\text{ \AA}$, $r_{\text{Ge-Ge}} = 3,26 \pm 0,10\text{ \AA}$, $r_{\text{O-O}} = 2,84 \pm 0,02\text{ \AA}$.
    • Dưới áp suất nén, xảy ra quá trình chuyển đổi từ $\text{GeO}_4$ sang $\text{GeO}_6$ qua trung gian $\text{GeO}_5$.
    • Đỉnh hàm PBXT $g_{\text{Ge-Ge}}(r)$ bị tách thành hai đỉnh và hàm $g_{\text{O-O}}(r)$ xuất hiện đỉnh thứ hai do sự gia tăng số lượng các liên kết chung cạnh và chung mặt giữa các đa diện.
    • Phân tích CNA xác nhận sự hình thành các mầm tinh thể oxy dạng $\text{fcc}$ và $\text{hcp}$ ở $60\text{ GPa}$. Tỉ số thể tích lỗ hổng $V_{\text{void}}/V$ giảm mạnh khi tăng áp suất nén.
  • Hệ $\text{SiO}_2$ thủy tinh:
    • Khảo sát các mẫu có mật độ từ $2,35\text{ g/cm}^3$ đến $4,42\text{ g/cm}^3$. Khi mật độ tăng, tỉ phần $\text{SiO}_4$ giảm, $\text{SiO}_5$ và $\text{SiO}_6$ tăng, số phối trí $\text{Si-Si}$ tăng lên 10 và $\text{O-O}$ tăng lên 12.
    • Sự tách đỉnh trong hàm $g_{\text{Si-Si}}(r)$ và đỉnh thứ hai trong $g_{\text{O-O}}(r)$ được định lượng chi tiết: xuất phát từ sự chuyển đổi phương thức kết nối giữa các đơn vị $\text{SiO}_x$, làm tăng số lượng nguyên tử oxy cầu nằm trong các liên kết cạnh và liên kết mặt.
    • Ở mật độ $4,42\text{ g/cm}^3$, nhận diện các đám nguyên tử oxy sắp xếp trật tự theo cấu trúc $\text{fcc}$, $\text{hcp}$ và cấu trúc lai $\text{fcc-hcp}$.
    • Phân bố bán kính các loại quả cầu lỗ hổng ($\text{O-void}$, $\text{Si}_4\text{-void}$, $\text{Si}_5\text{-void}$, $\text{Si}_6\text{-void}$ và các void lai) co hẹp lại, thể tích lỗ hổng giảm tuyến tính theo mật độ.

Chương 4: Cơ tính của các hệ ôxit $\text{Al}_2\text{O}_3$ VĐH, $\text{GeO}_2$ và $\text{SiO}_2$ thủy tinh

Chương 4 tập trung phân tích các đặc tính cơ học và cơ chế biến dạng vi mô:

  • Cơ tính hệ $\text{Al}_2\text{O}_3$ VĐH: Xây dựng đường cong ứng suất - biến dạng toàn phần và đường cong riêng phần cho từng loại đơn vị $\text{AlO}_x$ ($x=4,5,6$). Xác định sự phụ thuộc của mô-đun Young toàn phần và mô-đun riêng phần vào mật độ và áp suất nén của mẫu.
  • Cơ tính hệ $\text{GeO}_2$ thủy tinh: Khảo sát đường cong ứng suất - biến dạng trong dải áp suất $0\text{--}90\text{ GPa}$ tại $300\text{ K}$, xác định các thông số cơ học đặc trưng:
Mẫu $\text{GeO}_2$ (Áp suất chế tạo) Mô-đun đàn hồi Young $E\text{ (GPa)}$ Ứng suất đàn hồi $\sigma_e\text{ (GPa)}$ Ứng suất chảy $\sigma_y\text{ (GPa)}$
Mẫu ở $0\text{ GPa}$ Số liệu xác định từ đường cong Giá trị giới hạn đàn hồi Ngưỡng chảy dẻo
Mẫu ở áp suất cao Tăng theo áp suất nén Tăng theo mật độ Tăng theo mức độ xếp chặt
  • Cơ tính và cơ chế biến dạng của hệ $\text{SiO}_2$ thủy tinh:
    • Khảo sát 3 mẫu thủy tinh $\text{SiO}_2$ (mật độ $2,35\text{ g/cm}^3$, mẫu nén trung bình và $4,29\text{ g/cm}^3$). Đường cong $\sigma\text{--}\varepsilon$ thể hiện rõ bốn giai đoạn:
      1. Vùng biến dạng đàn hồi: Ứng suất tăng tuyến tính theo độ biến dạng; các quả cầu lỗ hổng biến dạng đàn hồi; độ dài liên kết $\text{Si-O}$ bị kéo giãn.
      2. Vùng biến dạng dẻo (dòng chảy dẻo): Ứng suất duy trì gần như không đổi; các nguyên tử tái sắp xếp cấu trúc mà không làm tăng độ dài liên kết; góc liên kết $\text{Si-O-Si}$ mở rộng từ $144^\circ$ lên $163^\circ$; hình thành các đơn vị khuyết tật $\text{SiO}_3$, $\text{SiO}_5$ và liên kết chung cạnh.
      3. Vùng biến dạng đứt gãy: Ứng suất giảm mạnh; các quả cầu lỗ hổng lân cận oxy ($\text{O-void}$) phát triển kích thước vượt ngưỡng tới hạn ($> 2,0\text{ \AA}$ tại mật độ $2,35\text{ g/cm}^3$); các $\text{O-void}$ kết cụm và hình thành vùng tới hạn gây chia cắt vật liệu.
      4. Vùng phá hủy hoàn toàn: Ứng suất triệt tiêu về 0; mẫu bị đứt gãy hoàn toàn; cấu trúc cục bộ ở các vùng ngoài vết nứt có sự phục hồi một phần độ dài và góc liên kết.
    • Khảo sát sự biến đổi của các đám mầm tinh thể oxy ($\text{fcc}$, $\text{hcp}$) trong mẫu mật độ cao ($4,29\text{ g/cm}^3$) dưới tác dụng của độ biến dạng kéo đơn trục.

Kết quả và những đóng góp mới

Những đóng góp mới về mặt khoa học

  1. Giải thích bản chất hiện tượng tách đỉnh trong hàm PBXT: Luận án lần đầu tiên đưa ra sự giải thích định lượng rõ ràng về sự thay đổi cấu trúc trật tự tầm trung thể hiện qua hiện tượng tách đỉnh trong hàm PBXT cặp $\text{Si-Si}$/$\text{Ge-Ge}$ và sự xuất hiện đỉnh phụ thứ hai trong hàm PBXT cặp $\text{O-O}$ ở áp suất cao trong thủy tinh $\text{SiO}_2$ và $\text{GeO}_2$. Hiện tượng này được chứng minh là do sự hình thành và gia tăng của các liên kết chung cạnh và chung mặt giữa các khối đa diện $\text{TO}_x$.
  2. Phát hiện trật tự sắp xếp mới của các nguyên tử oxy: Luận án lần đầu tiên tìm ra sự sắp xếp có trật tự của mạng nguyên tử oxy trong cả ba hệ vật liệu $\text{Al}_2\text{O}_3$ VĐH, $\text{GeO}_2$ và $\text{SiO}_2$ thủy tinh khi chịu nén ở áp suất cao ($60\text{ GPa}$ hoặc mật độ cao $4,42\text{ g/cm}^3$). Các nguyên tử oxy tự sắp xếp cục bộ tạo thành các đám mầm tinh thể kích thước vài chục nguyên tử có cấu trúc mạng lục giác xếp chặt ($\text{hcp}$) hoặc lập phương tâm mặt ($\text{fcc}$).
  3. Định lượng mối liên hệ giữa cấu trúc và cơ tính: Tính toán chính xác mô-đun Young toàn phần, mô-đun Young riêng phần của từng đơn vị cấu trúc $\text{TO}_x$ và xác định quy luật phụ thuộc của các đặc trưng cơ học vào áp suất nén (mật độ) trong dải từ $0$ đến $90\text{ GPa}$.
  4. Làm rõ cơ chế vi mô của dòng chảy dẻo và biến dạng đứt gãy: Phân tích chi tiết sự tiến triển vi mô qua 4 giai đoạn biến dạng của hệ $\text{SiO}_2$ thủy tinh.
  5. Phát hiện vai trò quyết định của quả cầu lỗ hổng $\text{O-void}$: Luận án lần đầu tiên chỉ ra rằng các quả cầu lỗ hổng lân cận nguyên tử oxy ($\text{O-void}$) đóng vai trò là các lỗ hổng kích thước lớn chủ đạo, dẫn dắt sự co cụm thể tích tự do và hình thành vùng tới hạn gây ra quá trình đứt gãy trong vật liệu $\text{SiO}_2$ thủy tinh.

Đề xuất và kiến nghị giải pháp

  • Ứng dụng các kết quả mô phỏng về mối liên hệ giữa tỉ phần đơn vị cấu trúc $\text{TO}_x$, liên kết cạnh/mặt và phân bố lỗ hổng để xây dựng các mô hình nguyên tử thực tế mô tả dòng chảy dẻo trong vật liệu ôxit vô định hình.
  • Sử dụng các thông số về áp suất chuyển pha và mật độ tới hạn làm cơ sở dữ liệu tính toán cho công nghệ xử lý nén ép, gia công nhiệt - áp nhằm cải thiện độ bền cơ học và kiểm soát tính dẻo/giòn của vật liệu thủy tinh quang học và gốm kỹ thuật.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  • Mô hình tính toán ĐLHPT cổ điển sử dụng các hàm thế cặp thực nghiệm (Coulomb-Buckingham, Oeffner-Elliot, BKS), chưa xét đến sự phân cực động học phi đối xứng của đám mây điện tử ở mức lượng tử dưới điều kiện áp suất cực hạn như các phương pháp ab initio.
  • Quy mô mô hình giới hạn trong phạm vi $5.000\text{--}6.000$ nguyên tử và thời gian mô phỏng ở thang picô-giây ($\text{ps}$), do đó tốc độ biến dạng trong mô phỏng kéo đơn trục cao hơn so với tốc độ biến dạng trong thực nghiệm cơ học thông thường.
  • Nghiên cứu cơ chế biến dạng đứt gãy chi tiết mới chỉ tập trung thực hiện sâu trên hệ $\text{SiO}_2$ thủy tinh dưới tác dụng của tải trọng đơn trục tại nhiệt độ $300\text{ K}$.

Hướng nghiên cứu tiếp theo

  • Mở rộng mô phỏng biến dạng cho hệ $\text{Al}_2\text{O}_3$ VĐH và $\text{GeO}_2$ thủy tinh dưới các chế độ tải trọng phức tạp khác (biến dạng cắt, nén hai trục, nén ba trục).
  • Khảo sát ảnh hưởng đồng thời của nhiệt độ cao và tốc độ biến dạng khác nhau lên động học phát triển của các quả cầu lỗ hổng $\text{O-void}$ và mầm nứt gãy.
  • Nghiên cứu mô phỏng các cấu trúc màng mỏng, cấu trúc nano và vật liệu ôxit xốp nhiều thành phần có chứa tạp chất.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Vật lý kỹ thuật, Khoa học và kỹ thuật vật liệu, Vật lý tin học, Hóa học chất rắn và Địa chất vật lý.
  • Kỹ sư và nhà công nghệ: Làm việc trong lĩnh vực chế tạo gốm sứ kỹ thuật cao, sản xuất thủy tinh quang học, sợi dẫn quang và vật liệu chịu áp suất/nhiệt độ cao.

Các nội dung đáng tham khảo

  • Chương 2: Cung cấp đầy đủ bộ thông số hàm thế tương tác cho ba hệ ôxit ($\text{Al}_2\text{O}_3$, $\text{GeO}_2$, $\text{SiO}_2$), quy trình thuật toán phân tích lân cận chung (CNA), phương pháp tính thể tích/phân bố quả cầu lỗ hổng bằng Monte Carlo và thuật toán phân tích mật độ riêng phần của các đơn vị cấu trúc.
  • Chương 3 & 4: Cung cấp bộ dữ liệu số chi tiết về hàm PBXT, độ dài liên kết, tỉ phần đơn vị cấu trúc, số lượng liên kết cạnh/mặt, thông số mô-đun Young và đường cong ứng suất - biến dạng chuẩn hóa ở dải áp suất rộng từ $0$ đến $90\text{ GPa}$.

Câu hỏi thường gặp

1. Nguyên nhân vật lý nào gây ra hiện tượng tách đỉnh trong hàm PBXT cặp cation ($\text{Si-Si}$, $\text{Ge-Ge}$) ở áp suất cao?

Hiện tượng tách đỉnh trong hàm phân bố xuyên tâm cặp cation ($\text{Si-Si}$ và $\text{Ge-Ge}$) khi bị nén ở áp suất cao là do sự thay đổi phương thức kết nối giữa các đơn vị cấu trúc đa diện $\text{TO}_x$. Ở áp suất thấp, các đa diện kết nối chủ yếu qua liên kết chung góc (chia sẻ 1 nguyên tử oxy). Khi áp suất tăng cao, mạng cấu trúc hình thành thêm các liên kết chung cạnh (chia sẻ 2 nguyên tử oxy) và chung mặt (chia sẻ 3 nguyên tử oxy), làm xuất hiện các khoảng cách cation - cation ngắn hơn đáng kể so với liên kết chung góc, dẫn đến sự phân tách đỉnh chính thành hai vị trí phân biệt.

2. Dưới áp suất nén cao ($60\text{ GPa}$), các nguyên tử oxy trong ba hệ ôxit sắp xếp theo cấu trúc trật tự nào?

Dưới ảnh hưởng của áp suất nén cao ($60\text{ GPa}$ đối với $\text{Al}_2\text{O}_3$, $\text{GeO}_2$ và ở mật độ $4,42\text{ g/cm}^3$ đối với $\text{SiO}_2$), các nguyên tử oxy trong mạng vô định hình tự sắp xếp có trật tự hơn, tạo thành các đám mầm tinh thể cục bộ gồm vài chục nguyên tử có cấu trúc lục giác xếp chặt ($\text{hcp}$) hoặc lập phương tâm mặt ($\text{fcc}$), được nhận diện chính xác qua phương pháp phân tích lân cận chung (CNA).

3. Tỉ phần các đơn vị cấu trúc $\text{AlO}_4$, $\text{AlO}_5$ và $\text{AlO}_6$ biến đổi như thế nào trong dải áp suất từ $0$ đến $22,5\text{ GPa}$ ở hệ $\text{Al}_2\text{O}_3$ VĐH?

Tại áp suất $0\text{ GPa}$, đơn vị tứ diện $\text{AlO}_4$ chiếm ưu thế với tỉ phần $67%$, tổng tỉ phần của $\text{AlO}_5$ và $\text{AlO}_6$ chiếm khoảng $33%$. Khi áp suất nén tăng từ $0$ đến $22,5\text{ GPa}$, tỉ phần $\text{AlO}_4$ giảm mạnh từ $67%$ xuống còn $18%$, trong khi tỉ phần đơn vị bát diện $\text{AlO}_6$ tăng liên tục. Đơn vị trung gian $\text{AlO}_5$ tăng nhanh và đạt giá trị cực đại $57%$ tại mốc áp suất $22,5\text{ GPa}$ trước khi giảm dần ở các áp suất cao hơn.

4. Quá trình biến dạng kéo đơn trục của thủy tinh $\text{SiO}_2$ diễn ra qua những giai đoạn nào và vai trò của quả cầu lỗ hổng $\text{O-void}$ là gì?

Quá trình biến dạng kéo đơn trục của thủy tinh $\text{SiO}_2$ diễn ra qua bốn vùng rõ rệt: (1) vùng biến dạng đàn hồi, (2) vùng biến dạng dẻo (dòng chảy dẻo), (3) vùng biến dạng đứt gãy, và (4) vùng phá hủy hoàn toàn. Trong quá trình này, các quả cầu lỗ hổng lân cận nguyên tử oxy ($\text{O-void}$) đóng vai trò là các lỗ hổng kích thước lớn chủ đạo; sự phát triển và kết cụm của các $\text{O-void}$ quanh vùng có mật độ thấp tạo nên vùng tới hạn, trực tiếp dẫn đến quá trình tạo mầm nứt và đứt gãy vật liệu.

5. Luận án sử dụng những loại hàm thế tương tác nào và quy mô số lượng nguyên tử cho từng mô hình là bao nhiêu?

Luận án sử dụng ba loại thế tương tác liên nguyên tử: thế Coulomb-Buckingham cho hệ $\text{Al}_2\text{O}_3$ (quy mô $5.000$ nguyên tử: $2.000\text{ Al} + 3.000\text{ O}$), thế Oeffner-Elliot cho hệ $\text{GeO}_2$ (quy mô $6.000$ nguyên tử: $2.000\text{ Ge} + 4.000\text{ O}$), và thế BKS cho hệ $\text{SiO}_2$ (quy mô $6.000$ nguyên tử: $2.000\text{ Si} + 4.000\text{ O}$).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thu Giang đã khảo sát toàn diện và có hệ thống cấu trúc vi mô và cơ tính của ba hệ ôxit hai nguyên $\text{Al}_2\text{O}_3$ VĐH, $\text{GeO}_2$ và $\text{SiO}_2$ thủy tinh dưới tác dụng của áp suất nén từ $0$ đến $90\text{ GPa}$ bằng phương pháp động lực học phân tử. Công trình đã làm sáng tỏ cơ chế vi mô của hiện tượng tách đỉnh trong hàm phân bố xuyên tâm cặp cation và đỉnh phụ oxy thông qua vai trò của các liên kết chung cạnh và chung mặt, đồng thời lần đầu tiên phát hiện sự hình thành các mầm tinh thể oxy dạng $\text{fcc}$ và $\text{hcp}$ ở áp suất cao. Bên cạnh việc định lượng mô-đun Young phụ thuộc áp suất, luận án đã xác lập vai trò quyết định của các quả cầu lỗ hổng $\text{O-void}$ trong việc hình thành vùng tới hạn gây đứt gãy vật liệu khi chịu biến dạng đơn trục. Các kết quả đạt được cung cấp cơ sở dữ liệu khoa học có độ tin cậy cao, đóng góp thiết thực cho việc nghiên cứu lý thuyết trạng thái mất trật tự và định hướng công nghệ chế tạo vật liệu ôxit chịu lực.