Tổng quan luận án

Luận án tiến sĩ Vật lí với đề tài "Nghiên cứu tính chất nhiệt động của hợp kim xen kẽ nhị nguyên và tam nguyên có khuyết tật với các cấu trúc lập phương tâm diện và lập phương tâm khối" được thực hiện bởi nghiên cứu sinh Lê Hồng Việt, dưới sự hướng dẫn khoa học của TS Nguyễn Quang Học và TS Phạm Thị Minh Hạnh, thuộc chuyên ngành Vật lí lý thuyết và Vật lí toán (mã số: 9 44 01 03) tại Trường Đại học Sư phạm Hà Nội 2 vào năm 2022. Toàn văn luận án gồm 141 trang, tích hợp 25 bảng số liệu, 71 hình vẽ và đồ thị, cùng 129 tài liệu tham khảo.

Tính cấp thiết của đề tài xuất phát từ vai trò của vật liệu hợp kim xen kẽ (IA) trong các lĩnh vực quân sự, khoa học kỹ thuật và công nghệ vật liệu. Các hợp kim nhẹ, bền, có khả năng chịu nhiệt độ và áp suất cao được sử dụng trong chế tạo vũ khí, tên lửa, tàu vũ trụ, máy bay và ô tô. Các hợp kim có độ bền hóa học và cơ học cao phục vụ chế tạo thiết bị ngành dầu mỏ và công nghiệp hóa chất; hợp kim không gỉ ứng dụng trong dụng cụ y tế và dân dụng; hợp kim của vàng với đồng, silic ứng dụng trong dây siêu dẫn, trang sức, ổ trục, đúc, hàn bước và che chắn bức xạ. Nghiên cứu hợp kim chứa sắt (FeC, FeSi, FeH) cung cấp dữ liệu về thành phần và sự tiến hóa của lõi Trái Đất cùng các thiên thể; việc bổ sung silic vào thép làm tăng điện trở suất và độ từ thẩm, giúp giảm tổn hao do dòng xoáy trong máy biến áp và động cơ.

Khoảng trống nghiên cứu được tác giả xác định là: các nghiên cứu trước đây chủ yếu tập trung vào kim loại nguyên chất, hợp kim thay thế (SA) có khuyết tật, hoặc hợp kim xen kẽ lý tưởng (không có khuyết tật mạng). Trong thực tế ở vùng nhiệt độ cao, các khuyết tật điểm, đặc biệt là nút khuyết cân bằng, xuất hiện với nồng độ lớn (bậc độ lớn khoảng $10^{-3}$ ở kim loại có điểm nóng chảy thấp và $10^{-2}$ ở kim loại có điểm nóng chảy cao), gây ra các ảnh hưởng phi tuyến mạnh đến nhiệt dung và độ dãn nở nhiệt. Ở dải áp suất cao trên 100 GPa, các sai hỏng mạng đóng góp trực tiếp vào nhiệt độ nóng chảy của tinh thể. Tuy nhiên, việc nghiên cứu định lượng ảnh hưởng của khuyết tật lên tính chất nhiệt động, nóng chảy và chuyển pha cấu trúc của hợp kim xen kẽ nhị nguyên và tam nguyên cấu trúc lập phương tâm diện (FCC) và lập phương tâm khối (BCC) bằng phương pháp lý thuyết giải tích vẫn còn ít, nhỏ lẻ và là vấn đề bỏ ngỏ.

Mục tiêu nghiên cứu của luận án:

  • Áp dụng phương pháp thống kê mômen (SMM) để nghiên cứu các tính chất nhiệt động, nóng chảy và chuyển pha cấu trúc của các hợp kim xen kẽ nhị nguyên và hợp kim xen kẽ tam nguyên có khuyết tật với cấu trúc lập phương, có tính đến ảnh hưởng đồng thời của nhiệt độ, áp suất, nồng độ nguyên tử thay thế, nồng độ nguyên tử xen kẽ và nồng độ nút khuyết cân bằng.

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

  • Đối tượng nghiên cứu: Các tính chất nhiệt động, đặc tính nóng chảy và chuyển pha cấu trúc của các hợp kim xen kẽ nhị nguyên (TaSi, WSi, FeC, AuSi, PtSi) và hợp kim xen kẽ tam nguyên (FeCrSi, VWSi, AuCuSi, PtCuSi).
  • Phạm vi nghiên cứu:
    • Nhiệt độ: Từ 0 K đến nhiệt độ nóng chảy của kim loại chính trong hợp kim.
    • Áp suất: Từ áp suất không ($P = 0$) đến dưới 100 GPa.
    • Nồng độ nguyên tử thay thế: Từ 0 đến 10%.
    • Nồng độ nguyên tử xen kẽ: Từ 0 đến 5%.
    • Nồng độ nút khuyết cân bằng: Khảo sát tương ứng theo điều kiện nhiệt độ và áp suất thực nghiệm.

Tổng quan tài liệu và vị trí của luận án

Theo phân loại của Goldschmidt, trong hợp kim xen kẽ, kim loại chính chiếm tỷ lệ từ 90% trở lên và quyết định kiểu cấu trúc tinh thể (thường là FCC, BCC hoặc lục giác xếp chặt HCP), trong khi các nguyên tố phi kim như H, C, Si đóng vai trò là tác nhân tạo hợp kim với nồng độ thường dưới 1%.

Đối với hợp kim nhị nguyên, giản đồ pha của hệ Fe và FeC đã được xây dựng thực nghiệm bởi Okamoto, Chipman. Nghiên cứu hệ FeH ở áp suất cao nhằm tìm hiểu cấu trúc lõi Trái Đất được thực hiện bởi Hirao và cộng sự, Mao và cộng sự, Pépin và cộng sự, Shibazaki và cộng sự, Terasaki và cộng sự. Kuopanportti và cộng sự nghiên cứu thế tương tác của FeH trong mô phỏng chiếu xạ; Lee và Jang khảo sát thế tương tác nhúng biến dạng cho FeH. Về lý thuyết nóng chảy FeC, Melnykov và Davidchack sử dụng mô phỏng động lực học phân tử (MD) ở áp suất không với ba loại thế: thế EAM của Lau và cộng sự, thế bậc liên kết giải tích của Henriksson và cộng sự, thế EAM của Hepburn và cộng sự, nhưng các mô hình này bộc lộ sai lệch về vùng ổn định của pha $\gamma$-FeC hoặc làm tăng quá cao nhiệt độ nóng chảy. Phương pháp tính toán giản đồ pha (CALPHAD) của Fei và Brosh khắc phục được một số nhược điểm định tính nhưng đòi hỏi lượng dữ liệu thực nghiệm lớn để tinh chỉnh hàm nhiệt động Gibbs và chỉ áp dụng tại một số áp suất cố định. Hệ Fe-Cr được tổng kết thực nghiệm và mô phỏng bởi Xiong và cộng sự cùng Okamoto. Đối với hệ TaSi và WSi tại $P = 0$, Guo và cộng sự tính toán giản đồ pha bằng CALPHAD; đường cong nóng chảy áp suất cao của Ta và W được nghiên cứu bởi Strachan và cộng sự (phương pháp qEAM FF), Taioli và cộng sự (ab initio), Wang và cộng sự (định luật Lindemann), Errandonea và cộng sự (buồng kim cương DAC), Liu và cộng sự (MD-PDEAM và MD-EFS), Burakovsky và cộng sự. Maier và cộng sự xác định nồng độ nút khuyết bằng kỹ thuật hủy positron; Kraftmakher phân tích mẫu làm nguội. Các hệ AuSi, CuSi được Hoc và cộng sự nghiên cứu bằng SMM ở trạng thái lý tưởng.

Đối với hợp kim tam nguyên, Zhang và cộng sự nghiên cứu ảnh hưởng của Si lên nhiệt độ Debye và hệ số dãn nở nhiệt của FeCr; Liang và cộng sự đo nhiệt độ nóng chảy của AlSi và AlCu; Hallstedt và cộng sự mô tả hệ AlCuSi bằng CALPHAD; Nguyễn Thành Cường và cộng sự cùng Nguyễn Quang Học và cộng sự tính toán các tính chất nhiệt động, nóng chảy của AlCuSi (dải áp suất 0–100 GPa), FeCrSi, FeCrC, AuCuSi và AuCuLi lý tưởng bằng SMM. Về khuyết tật, Røyset và cộng sự nghiên cứu nút khuyết cân bằng trong hợp kim AlMgSi; Li và cộng sự khảo sát khuyết tật điểm nâng cao hiệu suất vật liệu chiếu xạ.

Về lý thuyết khuyết tật điểm, Frenkel dự đoán sự hình thành cặp Frenkel; Wagner và Schottky thiết lập mô hình khuyết tật Schottky. Huntington và Seitz tính toán thông số tạo khuyết tật; Carpenter và Pochapsky, Vineyard và Dienes, Rasor và cộng sự, Cezairliyan và cộng sự, Kramer và Nolting, Miiller và Cezairliyan, Kraftmakher phân tích đóng góp của nút khuyết vào nhiệt dung và độ dãn nở nhiệt. Krivoglaz và Smirnov cùng Smirnov xây dựng lý thuyết xác định nồng độ nút khuyết cân bằng phụ thuộc nhiệt độ và nồng độ xen kẽ cho mạng FCC và BCC.

Khoảng trống luận án lựa chọn giải quyết: Các phương pháp mô phỏng như ab initio bị hạn chế ở kích thước hệ nhỏ, nhiệt độ thấp và áp suất thấp; CALPHAD phụ thuộc tham số thực nghiệm; MD phụ thuộc vào việc lựa chọn thế tương tác. Trong khi đó, các nghiên cứu SMM trước đây mới áp dụng cho hợp kim thay thế có khuyết tật hoặc hợp kim xen kẽ lý tưởng. Luận án giải quyết việc xây dựng biểu thức giải tích SMM và tính toán định lượng cho hợp kim xen kẽ nhị nguyên và tam nguyên có khuyết tật mạng dưới tác động đồng thời của nhiệt độ và áp suất cao.

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

Luận án sử dụng phương pháp thống kê mômen (SMM) do Nguyễn Tăng và Terletsky đề xuất từ những năm 1960. SMM phát triển từ cơ học thống kê, cho phép nghiên cứu dao động mạng tinh thể có tính đến các số hạng phi điều hòa trong biểu thức Hamiltonian: $$\hat{H} = \hat{H}0 - \sum_i a_i \hat{Q}i$$ Phương pháp thiết lập mối liên hệ giữa toán tử bất kỳ $\hat{F}$ và tọa độ suy rộng $\hat{Q}k$ thông qua chuỗi các hệ số Bernoulli $B{2m}$ và các đạo hàm theo tham số ngoại lực $a_k$. Từ đó, phương pháp dẫn ra công thức truy chứng cho các mômen cấp cao và xác định năng lượng tự do Helmholtz $\psi = - \theta \ln Z$ (với $\theta = k{Bo}T$). Năng lượng tự do của nguyên tử $X$ ($X = A, A_1, A_2, C$) được khai triển tường minh chứa các thông số mạng $k_X$ và các tham số phi điều hòa $\gamma{1X}, \gamma_{2X}, \gamma_X$.

Về lý thuyết nút khuyết, sự hình thành nút khuyết làm biến đổi thế nhiệt động Gibbs của hệ. Nồng độ nút khuyết cân bằng $n_v$ thỏa mãn điều kiện cực tiểu hóa thế Gibbs $(\partial \Delta G / \partial n){P,T} = 0$: $$n_v = \exp\left(-\frac{g_v^f}{k{Bo}T}\right)$$ Trong hợp kim xen kẽ $AC$, độ biến thiên thế Gibbs tạo nút khuyết được xác định qua phân bố nồng độ: $$g_v^f(AC) = \sum_{X} c_X g_v^f(X)$$ với $g_v^f(X) = n_1(\psi_X^{(1)} - \psi_X) + (B_X - 1)\psi_X$, trong đó $n_1$ là số phối vị trên quả cầu thứ nhất quanh nút khuyết, $\psi_X^{(1)}$ là năng lượng tự do của nguyên tử lân cận nút khuyết, và $B_X \approx 1 + u_{0X} / (4\psi_0X)$. Nồng độ nút khuyết cân bằng trong hợp kim xen kẽ rút ra dạng: $$n_v^{AC} = n_v^A \exp\left(-\frac{c_C g_v^f(C)}{k_{Bo}T}\right)$$

Phương pháp quả cầu phối vị được sử dụng để xác định năng lượng liên kết $u_{0X}$ và các thông số tinh thể $k_X, \gamma_{1X}, \gamma_{2X}$:

  • Mạng BCC: Nguyên tử xen kẽ C được tính trong phép gần đúng hai quả cầu phối vị ($r_{1C}, r_{2C} = \sqrt{2}r_{1C}$); các nguyên tử kim loại chính $A_1, A_2$ được tính trong phép gần đúng ba quả cầu phối vị.
  • Mạng FCC: Nguyên tử C tính trong phép gần đúng hai quả cầu phối vị ($r_{1C}, r_{2C} = \sqrt{3}r_{1C}$); các nguyên tử kim loại chính $A_1, A_2$ tính trong phép gần đúng ba quả cầu phối vị.

Nguồn dữ liệu và công cụ tính toán:

  • Tương tác giữa các nguyên tử được mô tả bằng thế cặp Mie-Lennard-Jones (MLJ n-m): $$\phi(r) = \frac{D}{n - m} \left[ m \left(\frac{r_0}{r}\right)^n - n \left(\frac{r_0}{r}\right)^m \right]$$ với các tham số độ sâu thế $D$ và khoảng cách cân bằng $r_0$ trích xuất từ dữ liệu thực nghiệm kim loại sạch và hợp kim.
  • Khoảng cách lân cận gần nhất $r_1(P,T)$ và độ dời nguyên tử $y(P,T)$ được xác định từ phương trình trạng thái áp suất $P \cdot v = -r_1 \frac{1}{6} \frac{\partial u_0}{\partial r_1} + \frac{\theta x \coth x}{2k} \frac{\partial k}{\partial r_1}$, với thể tích nguyên tử $v = 4r_1^3 / (3\sqrt{3})$ cho mạng BCC và $v = r_1^3 / \sqrt{2}$ cho mạng FCC.
  • Hệ phương trình trạng thái phi tuyến được giải số bằng phương pháp lặp gần đúng trên phần mềm Maple.

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

Chương 1: Tổng quan về đối tượng nghiên cứu và phương pháp nghiên cứu

Chương này tổng kết cơ sở lý luận về cấu trúc hợp kim xen kẽ, phân loại theo mạng cơ sở FCC và BCC, vai trò của các nguyên tử xen kẽ phi kim (C, Si, H) và các nguyên tố thay thế (Cu, Cr, W). Tác giả phân tích lý thuyết tạo thành nút khuyết theo cơ chế Frenkel và Schottky, sự biến thiên của entanpi $H_F = E_F + pV_F$, entrôpi dao động $S_F$ và entrôpi cấu hình $S_c$. Chương 1 đánh giá chi tiết các phương pháp lý thuyết hiện đại: phương pháp động lực mạng tinh thể (sai lệch ở nhiệt độ cao), trường phonon tự hợp (phức tạp khi tính bậc hai), mô phỏng Monte Carlo (MCS), ab initio (DFT), động lực học phân tử (MD) và CALPHAD. Trên cơ sở đó, tác giả luận giải việc lựa chọn SMM làm công cụ nghiên cứu nhờ tính giải tích tường minh và khả năng bao hàm tính phi điều hòa.

Chương 2: Tính chất nhiệt động của hợp kim xen kẽ nhị nguyên có khuyết tật với cấu trúc lập phương

Chương 2 mô hình hóa cấu trúc hợp kim xen kẽ nhị nguyên $AC$:

  • Mạng FCC: Nguyên tử xen kẽ C nằm ở tâm khối; nguyên tử kim loại chính $A_1$ ở các tâm mặt; nguyên tử $A_2$ ở các đỉnh. Tỷ lệ nồng độ hình học: $c_A = 1 - 15c_C$, $c_{A1} = 6c_C$, $c_{A2} = 8c_C$.
  • Mạng BCC: Nguyên tử xen kẽ C nằm ở các tâm mặt; nguyên tử $A_1$ ở tâm khối; nguyên tử $A_2$ ở các đỉnh. Tỷ lệ nồng độ hình học: $c_A = 1 - 7c_C$, $c_{A1} = 2c_C$, $c_{A2} = 4c_C$.

Tác giả thiết lập biểu thức giải tích của năng lượng tự do Helmholtz $\psi_{AC}^R$, độ dời nguyên tử trung bình $y(P,T)$ và nồng độ nút khuyết cân bằng $n_v^{AC}$. Từ các đạo hàm nhiệt động của $\psi_{AC}^R$, chương này rút ra các biểu thức giải tích cho 8 đại lượng nhiệt động: hệ số nén đẳng nhiệt $\chi_T$, môđun đàn hồi đẳng nhiệt $B_T$, hệ số dãn nở nhiệt $\alpha_T$, nhiệt dung đẳng tích $C_V$, nhiệt dung đẳng áp $C_P$, hệ số nén đoạn nhiệt $\chi_S$, môđun đàn hồi đoạn nhiệt $B_S$ và thông số Grüneisen $\gamma_G$.

Kết quả tính số được thực hiện cho các hợp kim FCC gồm AuSi và PtSi có khuyết tật:

  • So sánh các đại lượng nhiệt động của Au, Pt sạch với hợp kim AuSi, PtSi tại áp suất $P = 0$ và nhiệt độ cao.
  • Tính toán nhiệt dung đẳng áp $C_P(T, c_{Si})$ của PtSi tại áp suất $P = 80\text{ GPa}$, chỉ ra sự gia tăng phi tuyến rõ rệt của nhiệt dung khi tính đến đóng góp của nút khuyết nhiệt so với mô hình hợp kim lý tưởng.

Chương 3: Tính chất nhiệt động của hợp kim xen kẽ tam nguyên có khuyết tật với cấu trúc lập phương

Chương này mở rộng hệ thống biểu thức giải tích cho hợp kim xen kẽ tam nguyên $ABC$ có khuyết tật (gồm kim loại cơ sở A, nguyên tố thay thế B và nguyên tố xen kẽ C) cho cả hai cấu trúc FCC và BCC. Tác giả thiết lập công thức năng lượng tự do $\psi_{ABC}^R$, độ dời mạng $y(P,T)$, nồng độ nút khuyết cân bằng $n_v^{ABC}$ và các đại lượng nhiệt động $\chi_T, B_T, \alpha_T, C_V, C_P, \chi_S, B_S, \gamma_G$.

Tính số được áp dụng trên 4 hệ vật liệu đại diện:

  • Mạng FCC: Hợp kim AuCuSi (tính tại $P = 8\text{ GPa}, 12\text{ GPa}, 80\text{ GPa}$; $c_{Cu} = 10%$, $c_{Si} = 1%$) và PtCuSi (tại $P = 80\text{ GPa}$; $c_{Cu} = 10%$, $c_{Si} = 1%$). Khảo sát sự biến thiên của $n_v$, $\alpha_T$ và $C_P$ theo nhiệt độ đến 1300 K và 2000 K.
  • Mạng BCC: Hợp kim FeCrSi (tính tại $P = 5\text{ GPa}, 10\text{ GPa}$; $c_{Cr} = 10%$, $c_{Si} = 5%$) và VWSi (tại $P = 80\text{ GPa}$; $c_W = 10%$, $c_{Si} = 1%$). Đánh giá độ dãn nở nhiệt và nhiệt dung đẳng áp phụ thuộc áp suất và nhiệt độ đến 2100 K.

Bảng tổng hợp các hệ hợp kim và điều kiện khảo sát tính số trong Chương 2 và Chương 3:

Hệ hợp kim Cấu trúc mạng Loại hợp kim Khoảng áp suất ($P$) Nồng độ thành phần Các đại lượng tính số chính
AuSi FCC Xen kẽ nhị nguyên $P = 0$ $c_{Si}$ biến thiên $\chi_T, B_T, \alpha_T, C_V, C_P, n_v$
PtSi FCC Xen kẽ nhị nguyên $P = 0$ đến 80 GPa $c_{Si}$ biến thiên $\alpha_T, C_P, n_v$
AuCuSi FCC Xen kẽ tam nguyên 8, 12, 80 GPa $c_{Cu} = 10%$, $c_{Si} = 1%$ $n_v, \alpha_T, C_P$ theo $T$ và $P$
PtCuSi FCC Xen kẽ tam nguyên 80 GPa $c_{Cu} = 10%$, $c_{Si} = 1%$ $C_P$ theo $T$ và $P$
FeCrSi BCC Xen kẽ tam nguyên 5, 10 GPa $c_{Cr} = 10%$, $c_{Si} = 5%$ $n_v, \alpha_T, C_P$ theo $T$ và $P$
VWSi BCC Xen kẽ tam nguyên 80 GPa $c_W = 10%$, $c_{Si} = 1%$ $n_v, \alpha_T, C_P$ theo $T$ và $P$

Chương 4: Nóng chảy và chuyển pha cấu trúc của hợp kim xen kẽ nhị nguyên và tam nguyên có khuyết tật với cấu trúc lập phương

Chương 4 xây dựng lý thuyết nóng chảy và chuyển pha cấu trúc cho hợp kim xen kẽ nhị nguyên và tam nguyên có khuyết tật bằng SMM. Tiêu chuẩn nóng chảy được thiết lập dựa trên điều kiện cân bằng pha nhiệt động học và tiêu chuẩn mất ổn định cơ học của mạng khi nồng độ nút khuyết đạt giá trị tới hạn tại biên rắn - lỏng (kết hợp quan hệ thể tích tạo nút khuyết $v_f \approx 8v_m$ theo Bollman).

Kết quả tính số và đối chứng thực nghiệm:

  • Tính đường cong nóng chảy $T_m(P)$ cho Ta và W ở áp suất đến 80 GPa, so sánh với dữ liệu thực nghiệm buồng kim cương DAC của Errandonea và cộng sự, tính toán qEAM FF của Strachan và cộng sự, ab initio của Taioli và cộng sự, và mô phỏng MD của Liu và cộng sự.
  • Tính toán nhiệt độ nóng chảy phụ thuộc nồng độ Si đối với TaSi và WSi tại $P = 0$, so sánh phù hợp với giản đồ pha CALPHAD của Guo và cộng sự; xác định độ dốc đường cong nóng chảy của TaSi và WSi tại $P = 0$.
  • Khảo sát ảnh hưởng của nguyên tử xen kẽ C lên nhiệt độ nóng chảy của hợp kim $\gamma$-FeC tại các nồng độ $c_C = 0, 3, 6, 9%$, giải thích quy luật giảm điểm nóng chảy khi tăng hàm lượng cacbon.
  • Khảo sát sự phụ thuộc nhiệt độ của thế nhiệt động Gibbs đối với pha $\alpha$-Fe và pha $\gamma$-Fe, từ đó xây dựng đường chuyển pha cấu trúc $\alpha - \gamma$ của sắt theo áp suất và nhiệt độ, cho kết quả phù hợp với thực nghiệm giản đồ pha của sắt.

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

Luận án tự nêu các đóng góp mới về mặt lý luận và thực tiễn như sau:

  • Về lý luận: Đã xây dựng được hệ thống biểu thức giải tích hoàn chỉnh cho các đại lượng nhiệt động ($\chi_T, B_T, \alpha_T, C_V, C_P, \chi_S, B_S, \gamma_G$), nồng độ nút khuyết cân bằng, nhiệt độ nóng chảy và nhiệt độ chuyển pha cấu trúc của các hợp kim xen kẽ nhị nguyên ($AC$) và tam nguyên ($ABC$) có khuyết tật với các cấu trúc FCC và BCC bằng phương pháp thống kê mômen (SMM).
  • Về thực tiễn và tính số: Áp dụng các công thức giải tích để tính số chi tiết cho các kim loại và hợp kim xen kẽ gồm AuSi, PtSi, AuCuSi, PtCuSi (cấu trúc FCC) và TaSi, WSi, FeC, FeCrSi, VWSi cùng sự chuyển pha của Fe (cấu trúc BCC). Kết quả tính số có sự phù hợp với số liệu thực nghiệm đo đạc và các kết quả tính bằng phương pháp khác như CALPHAD, ab initio, MD.
  • Khả năng dự báo: Kết quả tính số về các đại lượng nhiệt động, nhiệt độ nóng chảy và chuyển pha cấu trúc dưới áp suất cao (lên tới 80–100 GPa) và nồng độ khuyết tật cân bằng cung cấp dữ liệu định hướng cho các nghiên cứu thực nghiệm trong tương lai ở điều kiện áp suất cực đoan, đồng thời bổ sung dữ liệu lý thuyết về nút khuyết trong hợp kim xen kẽ cấu trúc lập phương.
  • Kiến nghị đề xuất: Đề xuất sử dụng các biểu thức giải tích SMM kết hợp thế tương tác MLJ như một phương pháp tính toán nhanh, hỗ trợ dự đoán tính chất vật liệu ở dải áp suất cao mà thực nghiệm buồng kim cương hay sóng xung kích gặp khó khăn kỹ thuật.

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

Hạn chế của nghiên cứu được xác định dựa trên phạm vi và mô hình của luận án:

  • Mô hình lý thuyết mới chỉ áp dụng cho hai dạng cấu trúc lập phương (FCC và BCC), chưa mở rộng cho cấu trúc lục giác xếp chặt (HCP) hoặc các mạng tinh thể phức tạp khác.
  • Khuyết tật mạng mới chỉ xét đến khuyết tật điểm là nút khuyết cân bằng đơn, chưa tính đến các dạng khuyết tật phức tạp như nút khuyết kép, cụm nút khuyết hay lệch mạng vốn có ảnh hưởng lớn ở áp suất trên 100 GPa.
  • Tương tác giữa các nguyên tử dựa trên thế cặp hiệu dụng Mie-Lennard-Jones, chưa xét đến các hiệu ứng tương tác nhiều hạt phi cặp phức tạp.
  • Giới hạn nồng độ pha tạp khảo sát ở mức nhỏ ($c_C \le 5%$ đối với nguyên tử xen kẽ, $c_B \le 10%$ đối với nguyên tử thay thế).

Hướng nghiên cứu tiếp theo (trình bày tại trang 111 của văn bản):

  • Mở rộng phương pháp thống kê mômen để nghiên cứu tính chất nhiệt động, nóng chảy và chuyển pha của hợp kim xen kẽ có cấu trúc lục giác xếp chặt (HCP).
  • Nghiên cứu ảnh hưởng của lệch mạng và các loại khuyết tật khác lên tính chất nhiệt động và độ bền của hợp kim ở dải áp suất cực cao trên 100 GPa.
  • Khảo sát các hệ hợp kim đa nguyên phức tạp hơn và mở rộng sang các nguyên tử xen kẽ khác như H, N, B.

Giá trị tham khảo

Luận án có giá trị tham khảo cho các đối tượng và lĩnh vực chuyên môn cụ thể:

  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Vật lí lý thuyết, Vật lí toán và Khoa học vật liệu: Tài liệu tham khảo về phương pháp luận áp dụng cơ học thống kê lượng tử (SMM) để giải quyết các bài toán phi điều hòa và khuyết tật mạng trong chất rắn kết tinh.
  • Các nhà nghiên cứu công nghệ vật liệu và luyện kim: Dữ liệu tính toán nhiệt động, hệ số dãn nở nhiệt, nhiệt dung và điểm nóng chảy của các hệ hợp kim thép (FeC, FeCrSi) và kim loại chịu nhiệt (TaSi, WSi, VWSi) hữu ích cho việc thiết kế vật liệu làm việc ở nhiệt độ và áp suất cao.
  • Các nhà khoa học địa vật lý và vật lý thiên thể: Các kết quả về đường cong nóng chảy của FeC và chuyển pha $\alpha - \gamma$ của sắt dưới áp suất cao cung cấp thông số tham chiếu để mô hình hóa trạng thái vật chất trong lõi Trái Đất.
  • Phần nội dung đáng tham khảo nhất: Hệ thống công thức giải tích năng lượng tự do và các hệ số nhiệt động tại Chương 2, Chương 3 và quy trình giải số phương trình trạng thái trên Maple.

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

  1. Mô hình hình học của ô cơ sở hợp kim xen kẽ nhị nguyên $AC$ cấu trúc FCC và BCC được phân bố như thế nào? Trong cấu trúc FCC, nguyên tử xen kẽ C nằm ở tâm khối của ô lập phương, nguyên tử kim loại chính $A_1$ nằm ở các tâm mặt và $A_2$ nằm ở các đỉnh. Trong cấu trúc BCC, nguyên tử xen kẽ C nằm ở các tâm mặt, nguyên tử kim loại chính $A_1$ nằm ở tâm khối và $A_2$ nằm ở các đỉnh.

  2. Phương pháp thống kê mômen (SMM) xử lý tính phi điều hòa của dao động mạng tinh thể bằng cơ chế nào? Phương pháp đưa trực tiếp các số hạng phi điều hòa bậc ba và bậc bốn của thế tương tác vào biểu thức năng lượng thông qua chuỗi các tham số tinh thể $\gamma_{1X}, \gamma_{2X}, \gamma_X$. Các đại lượng nhiệt động học được dẫn xuất từ năng lượng tự do Helmholtz thông qua chuỗi truy chứng mômen và các hệ số Bernoulli $B_{2m}$, không bị giới hạn bởi phép gần đúng chuẩn điều hòa.

  3. Nồng độ nút khuyết cân bằng trong hợp kim xen kẽ nhị nguyên $AC$ có dạng biểu thức như thế nào? Nồng độ nút khuyết cân bằng có dạng $n_v^{AC} = n_v^A \exp\left(-\frac{c_C g_v^f(C)}{k_{Bo}T}\right)$, trong đó $n_v^A$ là nồng độ nút khuyết của kim loại cơ sở A, $c_C$ là nồng độ nguyên tử xen kẽ, và $g_v^f(C)$ là độ biến thiên thế Gibbs tạo nút khuyết gắn với nguyên tử xen kẽ C, tính theo công thức $g_v^f(X) = n_1(\psi_X^{(1)} - \psi_X) + (B_X - 1)\psi_X$.

  4. Phạm vi áp suất và nồng độ nguyên tử được luận án giới hạn khảo sát là bao nhiêu? Luận án giới hạn phạm vi áp suất từ 0 đến dưới 100 GPa; khoảng nhiệt độ từ 0 K đến nhiệt độ nóng chảy của kim loại cơ sở; khoảng nồng độ nguyên tử thay thế từ 0 đến 10%; và khoảng nồng độ nguyên tử xen kẽ từ 0 đến 5%.

  5. Luận án đã tính số và đối chứng thực nghiệm cho những hệ hợp kim cụ thể nào? Luận án thực hiện tính số cho các hệ nhị nguyên FCC (AuSi, PtSi), nhị nguyên BCC (TaSi, WSi, FeC), tam nguyên FCC (AuCuSi, PtCuSi) và tam nguyên BCC (FeCrSi, VWSi). Kết quả nhiệt độ nóng chảy và chuyển pha được so sánh với thực nghiệm buồng kim cương DAC của Errandonea và cộng sự, giản đồ pha CALPHAD của Guo và cộng sự, cùng các tính toán ab initio và động lực học phân tử (MD).

Kết luận

Luận án của nghiên cứu sinh Lê Hồng Việt đã áp dụng phương pháp thống kê mômen (SMM) để xây dựng hệ thống biểu thức giải tích xác định các tính chất nhiệt động, nồng độ nút khuyết cân bằng, nhiệt độ nóng chảy và chuyển pha cấu trúc của hợp kim xen kẽ nhị nguyên và tam nguyên cấu trúc FCC và BCC. Các tính toán số trên các hệ hợp kim cụ thể (AuSi, PtSi, AuCuSi, PtCuSi, FeCrSi, VWSi, TaSi, WSi, FeC) trong dải áp suất đến 80–100 GPa cho thấy sự phù hợp với thực nghiệm và các mô hình lý thuyết khác, đồng thời chỉ rõ vai trò của nút khuyết nhiệt ở vùng nhiệt độ cao. Công trình cung cấp cơ sở dữ liệu giải tích và số liệu định hướng có giá trị cho nghiên cứu vật liệu chịu nhiệt, áp suất cao và vật lý chất rắn.