Tổng quan luận án

Luận án tiến sĩ vật lý với đề tài "Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS" được tác giả Trịnh Phi Hiệp thực hiện tại Trường Đại học Sư phạm Hà Nội 2 dưới sự hướng dẫn khoa học của PGS. Nguyễn Văn Thụ và PGS. Nguyễn Bá Đức, thuộc chuyên ngành Vật lý lý thuyết và Vật lý toán (mã số: 9 44 01 03), hoàn thành năm 2023.

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

Phổ cấu trúc tinh tế hấp thụ tia X mở rộng (Extended X-ray Absorption Fine Structure - EXAFS) cung cấp các thông tin cấu trúc vi mô quan trọng như số nguyên tử phối trí trong từng lớp nguyên tử và bán kính của các lớp nguyên tử thông qua ảnh biến đổi Fourier của phổ. Phương pháp EXAFS là công cụ lý thuyết và thực nghiệm then chốt để nghiên cứu vật liệu cấu trúc tinh thể, chất vô định hình và vật liệu nano. Ở vùng nhiệt độ thấp, biên độ dao động nhiệt của các nguyên tử là nhỏ, dao động mạng được coi là dao động điều hòa quanh vị trí cân bằng và mô hình điều hòa cho kết quả phù hợp với thực nghiệm. Tuy nhiên, khi nhiệt độ tăng cao, các hiệu ứng phi điều hòa xuất hiện do tương tác phonon-phonon làm cho thế năng tương tác giữa các nguyên tử trở nên bất đối xứng; nếu bỏ qua các hiệu ứng này thì các thông tin cấu trúc trích xuất từ phổ EXAFS sẽ phát sinh sai số lớn. Phép khai triển gần đúng cumulant đã được thiết lập để hiệu chỉnh các sai số phi điều hòa đó.

Mặt khác, áp suất cao tác động trực tiếp lên khoảng cách liên kết và độ dời nguyên tử, khiến phổ EXAFS biến đổi rõ rệt. Mặc dù một số mô hình lý thuyết như mô hình Debye tương quan phi điều hòa đã được sử dụng để khảo sát ảnh hưởng của áp suất lên các cumulant, các nghiên cứu trước đó chưa làm rõ sự phụ thuộc đồng thời của các cumulant và các tham số nhiệt động vào nhiệt độ, tỷ lệ chất pha tạp đối với các hệ kim loại nguyên chất và hợp kim liên kim loại có cấu trúc tinh thể lập phương dưới tác dụng của áp suất cao. Đây chính là khoảng trống khoa học mà luận án tập trung giải quyết.

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

Mục tiêu tổng quát của luận án là áp dụng lý thuyết EXAFS hiện đại để nghiên cứu ảnh hưởng của hiệu ứng phi điều hòa và áp suất lên các đại lượng nhiệt động của các tinh thể và hợp kim có cấu trúc lập phương. Luận án đặt ra 4 mục tiêu cụ thể được đánh số theo văn bản:

  1. Xây dựng biểu thức của các tham số nhiệt động cho vật liệu có cấu trúc tinh thể lập phương. Áp dụng mô hình tương quan phi điều hòa để khai triển và tính gần đúng các cumulant của tinh thể có cấu trúc lập phương. Áp dụng các biểu thức nhiệt động và các cumulant cho vật liệu có cấu trúc tinh thể lập phương.
  2. Áp dụng các mô hình Einstein tương quan phi điều hòa (ACEM) và mô hình Debye tương quan phi điều hòa (ACDM) để nghiên cứu sự phụ thuộc vào nhiệt độ, áp suất và tỉ lệ pha tạp của các cumulant, hệ số giãn nở nhiệt, hệ số đàn hồi, pha dao động trong phổ EXAFS phi điều hòa. Áp dụng đối với các kim loại nguyên chất và hợp kim liên kim loại có cấu trúc tinh thể lập phương.
  3. Nghiên cứu sự phụ thuộc của các tham số nhiệt động và các cumulant vào nhiệt độ, tỷ lệ pha tạp dưới ảnh hưởng của áp suất cao. Giải thích sự khác biệt đáng kể giữa mô hình dao động tương quan phi điều hòa và mô hình dao động đơn hạt điều hòa. Dựa vào các hệ thức được xây dựng, tiến hành tính số và so sánh với thực nghiệm và các lý thuyết khác.
  4. Xây dựng phương pháp xác định các tham số thế Morse bằng lý thuyết qua năng lượng phát xạ, khả năng nén và hằng số mạng tinh thể.

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

  • Đối tượng nghiên cứu: Các đại lượng nhiệt động (hệ số giãn nở nhiệt khối và tuyến tính, hệ số Debye-Waller DWF, hằng số lực đàn hồi hiệu dụng $k_{eff}$, các hệ số phi điều hòa $k_3, k_4$, độ dịch pha, độ dịch chuyển trung bình bình phương MSD, độ dịch chuyển tương đối trung bình bình phương MSRD, hàm dịch chuyển tương quan DCF) và các cumulant từ bậc một đến bậc bốn ($\sigma^{(1)}, \sigma^{(2)}, \sigma^{(3)}, \sigma^{(4)}$) trong phổ EXAFS phi điều hòa.
  • Đối tượng vật liệu: Các tinh thể và hợp kim có cấu trúc mạng lập phương, bao gồm lập phương tâm khối (bcc), lập phương tâm diện (fcc), lập phương đơn giản (sc), tinh thể cấu trúc kiểu kim cương (Si, Ge, hợp kim SiGe) và các kim loại, hợp kim liên kim loại lập phương (Cu, Ag, $\text{Cu}x\text{Ag}{1-x}$).
  • Phạm vi nghiên cứu: Khảo sát lý thuyết trong phạm vi phương pháp thống kê lượng tử và lý thuyết nhiễu loạn; phạm vi áp suất từ $0\text{ GPa}$ đến $14\text{ GPa}$; dải nhiệt độ từ vùng nhiệt độ thấp ($T \to 0\text{ K}$) đến vùng nhiệt độ cao ($T \to \infty$); tỷ lệ pha tạp $x$ của hợp kim nhị phân.

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 được điểm lại trong luận án

Luận án đã hệ thống hóa lịch sử phát triển của lý thuyết quang phổ hấp thụ tia X và các mô hình tính toán cumulant:

  • Hai trường phái lý luận về phổ XAFS: Trường phái trật tự tầm xa (Long Range Order - LRO) đặc trưng qua mật độ trạng thái vùng năng lượng với quãng đường tự do vô hạn [2, 31, 71], và trường phái trật tự tầm gần (Short Range Order - SRO) coi phổ XAFS là hiệu ứng giao thoa trạng thái cuối của sóng quang electron tán xạ ngược từ các nguyên tử lân cận [2, 62, 63, 68]. Lý thuyết hiện đại phát triển theo hướng SRO do khả năng sử dụng biến đổi Fourier để trích xuất cấu trúc nguyên tử cục bộ.
  • Các mô hình tính cumulant giải tích:
    • Mô hình tương quan đơn cặp (Single-Bond Model - SBM) do Fornasini và cộng sự phát triển [25], chỉ xét dao động của cặp nguyên tử riêng lẻ.
    • Phương pháp gần đúng nhiệt động toàn mạng (Full Lattice Dynamical Approach - FLDA) [35], tính toán động lực học toàn mạng tinh thể nhưng đòi hỏi khối lượng tính toán lớn.
    • Phương pháp thế phi điều hòa đơn hạt (Anharmonic Single-Particle - ASP) [73], xem xét hạt chuyển động trong trường thế trung bình độc lập.
    • Mô hình Einstein tương quan phi điều hòa (Anharmonic-Correlated Einstein Model - ACEM) [44] và Mô hình Debye tương quan phi điều hòa (Anharmonic-Correlated Debye Model - ACDM) [49], do nhóm nghiên cứu của Nguyễn Văn Hùng và cộng sự xây dựng. Các mô hình này khắc phục hạn chế của mô hình đơn hạt bằng cách đưa vào tương tác tương quan của cụm nguyên tử lân cận gần nhất và hiệu ứng phi điều hòa.
  • Nghiên cứu ảnh hưởng của áp suất trong EXAFS: Các công trình của Hung, Duc, Thu, Frahm và cộng sự [45, 50, 52] đã sử dụng mô hình ACDM để phân tích sự thay đổi của cumulant bậc hai dưới tác dụng của áp suất.
  • Hàm phân bố hiệu dụng và kỹ thuật khai triển cumulant: Dựa trên các công trình nền tảng của Bunker [12], Tranquada và Ingalls [72], các cumulant được định nghĩa thông qua các mômen của hàm phân bố cặp bất đối xứng $P(r, \gamma)$.
  • Các dạng thế tương tác cặp: Thế Lennard-Jones (phù hợp cho khí trơ), thế Madelung (tinh thể ion), và thế Morse phi điều hòa [17, 48] được sử dụng cho liên kết cộng hóa trị và kim loại.

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

Các công trình trước đây chủ yếu khảo sát riêng rẽ ảnh hưởng của nhiệt độ hoặc áp suất lên kim loại nguyên chất đơn lẻ. Luận án lựa chọn giải quyết bài toán tổng quát hơn: xây dựng hệ thức giải tích tích hợp đồng thời ba yếu tố (nhiệt độ, áp suất cao lên đến $14\text{ GPa}$, và tỷ lệ pha tạp $x$) lên các đại lượng nhiệt động và các cumulant trong cấu trúc lập phương; đồng thời thiết lập phương pháp lý thuyết thuần túy để tính toán các tham số của thế tương tác Morse mà không phụ thuộc vào các tham số thực nghiệm đo đạc trước.

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

Khung lý thuyết và khái niệm sử dụng

  1. Lý thuyết tán xạ đơn trong phổ EXAFS phi điều hòa: Phổ EXAFS cận K phụ thuộc nhiệt độ được biểu diễn qua hàm dao động tổng quát: $$\chi(k, T) = \sum_j \frac{S_0^2 N_j}{k R_j^2} F_j(k) \exp\left[-2k^2\left(\sigma_H^2(T) + \sigma_A^2(T)\right) - \frac{2 R_j}{\lambda(k)}\right] \sin\left[2k R_j + \phi_j(k) + \phi_A(k, T)\right]$$ Trong đó, $S_0^2$ là thừa số suy giảm nhiều hạt, $N_j$ là số nguyên tử phối trí lớp thứ $j$, $F_j(k)$ là biên độ tán xạ ngược, $\lambda(k)$ là quãng đường tự do trung bình của electron quang, $\sigma_H^2(T)$ và $\sigma_A^2(T)$ lần lượt là đóng góp điều hòa và phi điều hòa vào cumulant bậc hai (MSRD), $\phi_A(k, T)$ là độ lệch pha do hiệu ứng phi điều hòa.

  2. Thế tương tác Morse mở rộng: Thế cặp phi điều hòa Morse giữa hai nguyên tử có dạng: $$U(r) = D \left[ e^{-2\alpha(r-r_0)} - 2e^{-\alpha(r-r_0)} \right]$$ Với $D$ là năng lượng phân ly liên kết, $\alpha$ là tham số độ rộng thế, $r_0$ là khoảng cách cân bằng. Khai triển Taylor theo độ lệch tức thời $x = r - r_0$ đến số hạng bậc bốn: $$U(x) \approx D \left( -1 + \alpha^2 x^2 - \alpha^3 x^3 + \frac{7}{12} \alpha^4 x^4 \right)$$

  3. Mô hình Einstein tương quan phi điều hòa (ACEM): Bỏ qua tán sắc phonon, xét tương tác của một nhóm nhỏ nguyên tử với khối lượng rút gọn $\mu = \frac{M_1 M_2}{M_1 + M_2}$. Thế tương tác hiệu dụng được biểu diễn: $$U_E(x) = U(x) + \sum_{i=0,1; i \neq j} U\left( \frac{\mu}{M_i} x \hat{\Re}{01} \cdot \hat{\Re}{ij} \right) \approx \frac{1}{2} k_{eff} x^2 + k_3 x^3 + k_4 x^4$$ Chuyển sang biến dao động quanh vị trí giãn nở nhiệt $y = x - a$ (với $a = \langle x \rangle = \sigma^{(1)}$ và $\langle y \rangle = 0$), thế hiệu dụng trở thành: $$U_E(y) = U_E(a) + \frac{1}{2} k_{eff} y^2 + \delta U_E(y)$$

  4. Mô hình Debye tương quan phi điều hòa (ACDM): Có xét đến tán sắc phonon trong hệ một chiều với quan hệ tán sắc $\omega(q) = 2 \sqrt{\frac{k_{eff}}{M}} \left| \sin\left(\frac{qa}{2}\right) \right|$, tần số Debye cực đại $\omega_D = 2 \sqrt{\frac{k_{eff}}{M}}$, nhiệt độ Debye $\theta_D = \frac{\hbar \omega_D}{k_B}$. Độ dịch chuyển $x_n = u_{n+1} - u_n$ được lượng tử hóa qua toán tử phonon $A_q = \hat{a}q + \hat{a}{-q}^+$.

Phương pháp nghiên cứu thực tế

  • Phương pháp thống kê lượng tử và lý thuyết nhiễu loạn: Toán tử Hamiltonian toàn phần $H = H_0 + \delta U_E(y)$ gồm phần điều hòa $H_0 = \frac{P^2}{2\mu} + \frac{1}{2} k y^2$ và phần nhiễu loạn phi điều hòa $\delta U_E(y)$. Các trạng thái dao động tử điều hòa $|n\rangle$ có mức năng lượng $E_n = n \hbar \omega_E$.
  • Kỹ thuật ma trận mật độ: Trị trung bình nhiệt của toán tử $y^m$ được tính qua ma trận mật độ $\rho = \rho_0 + \delta\rho$ với $\beta = \frac{1}{k_B T}$, tổng thống kê không nhiễu loạn $Z_0 = \text{Tr}(\rho_0) = \frac{1}{1 - z}$, trong đó $z = \exp(-\beta \hbar \omega_E) = \exp(-\theta_E / T)$ là biến số nhiệt độ không thứ nguyên.
  • Tính toán số trị và so sánh: Lập trình giải tích và tính toán các đại lượng nhiệt động dưới áp suất biến thiên $0 - 14\text{ GPa}$, nhiệt độ $0 - 1000\text{ K}$, tỷ lệ pha tạp $x$; so sánh đối chứng kết quả với số liệu thực nghiệm EXAFS và các mô hình lý thuyết khác từ các công trình công bố.

Nguồn tư liệu và dữ liệu

Luận án sử dụng các bộ dữ liệu thực nghiệm EXAFS và hằng số mạng, môđun nén, hằng số đàn hồi đã được công bố của tinh thể Si, Ge, hợp kim SiGe, kim loại Cu, Ag và hợp kim Cu-Ag (từ các công trình trích dẫn [28, 57, 60, 68]).

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

Chương 1: Tổng quan về phổ cấu trúc tinh tế hấp thụ tia X

Chương 1 trình bày cơ sở vật lý thực nghiệm và lý thuyết của quang phổ hấp thụ tia X:

  • Cơ chế tạo bức xạ tia X: Trình bày nguyên lý ống tia X tạo bức xạ hãm liên tục (Bremsstrahlung) tuân theo định luật Duane-Hunt ($\lambda_{min} = \frac{hc}{eV}$) và bức xạ đặc trưng ($K_\alpha, K_\beta$) do dịch chuyển electron giữa các phân lớp $L \to K, M \to K$. Trình bày 6 đặc tính ưu việt của bức xạ Synchrotron (cường độ lớn, phổ liên tục rộng, độ ổn định cao, tính chuẩn trực lớn, phân cực phẳng, cấu trúc xung micro-giây, kích thước nguồn nhỏ).
  • Hiện tượng hấp thụ và cấu trúc tinh tế: Phân tích tương tác photon-vật rắn gồm tán xạ đàn hồi Rayleigh, tán xạ không đàn hồi Compton và hiệu ứng quang điện. Hệ số hấp thụ toàn phần $\epsilon(E) = \epsilon_a(E)[1 + \chi(E)]$ được phân tách thành phần hấp thụ nguyên tử cô lập $\epsilon_a(E)$ và phần cấu trúc tinh tế $\chi(E) = \frac{\epsilon(E) - \epsilon_a(E)}{\epsilon_a(E)}$. Phân loại các miền quang phổ: EXAFS ($E > 50\text{ eV}$ sau cận hấp thụ), XANES/NEXAFS ($E < 50\text{ eV}$), cùng các biến thể bề mặt SXANES, SEXAFS.
  • Biến đổi Fourier và hàm phân bố hiệu dụng: Thiết lập phép biến đổi Fourier $F(r) = \frac{1}{\sqrt{2\pi}} \int \chi(k) k^n e^{-2ikr} dk$ ($n = 1, 2, 3$) chuyển phổ từ không gian số sóng $k$ sang không gian tọa độ thực $r$, chứng minh vị trí các đỉnh của $|F(r)|$ phản ánh khoảng cách đến các lớp nguyên tử lân cận $R_j$. Trình bày hàm phân bố hiệu dụng $P(r, \gamma) = \rho_1(r) e^{-2\gamma r} / r^2$ (với $\gamma = \lambda^{-1}$) và phép khai triển mômen để thu được các cumulant $\sigma^{(n)}$.
  • Hệ số Debye-Waller (DWF) và thế tương tác: Phân tích vai trò của thừa số suy giảm Debye-Waller $\exp(-2k^2\sigma_j^2)$, trong đó $\sigma_j^2$ là độ dịch chuyển tương đối trung bình bình phương (MSRD), liên hệ trực tiếp với độ dịch chuyển trung bình bình phương (MSD, $u_j^2$) và hàm dịch chuyển tương quan (DCF, $C_R$) qua hệ thức $\sigma_j^2 = \langle (\vec{u}_j - \vec{u}_0)^2 \rangle = u_j^2 + u_0^2 - 2 C_R$. Trình bày biểu thức và sự khai triển của các thế tương tác Lennard-Jones, Madelung và thế Morse phi điều hòa.
+-----------------------------------------------------------------------------+
| Bảng 1.1: Tổng hợp các khái niệm cơ bản trong lý thuyết phổ EXAFS (Chương 1) |
+------------------------------------+----------------------------------------+
| Ký hiệu / Khái niệm                | Ý nghĩa vật lý                         |
+------------------------------------+----------------------------------------+
| EXAFS                              | Phổ cấu trúc tinh tế hấp thụ tia X     |
|                                    | mở rộng (động năng electron > 50 eV)   |
| DWF (Debye-Waller Factor)          | Hệ số suy giảm biên độ do dao động     |
|                                    | nhiệt và nhiễu loạn cấu trúc           |
| MSRD (sigma^2)                     | Độ dịch chuyển tương đối trung bình    |
|                                    | bình phương giữa 2 nguyên tử           |
| MSD (u^2)                          | Độ dịch chuyển trung bình bình phương  |
|                                    | của một nguyên tử                      |
| DCF (C_R)                          | Hàm dịch chuyển tương quan dao động    |
| P(r, gamma)                        | Hàm phân bố hiệu dụng của cặp nguyên tử|
+------------------------------------+----------------------------------------+

Chương 2: Lý thuyết về phổ cấu trúc tinh tế hấp thụ tia X mở rộng phi điều hòa

Chương 2 tập trung xây dựng mô hình giải tích dựa trên ACEM và ACDM kết hợp thống kê lượng tử:

  • Biểu diễn phổ EXAFS phi điều hòa qua cumulant: Hàm số mũ Debye-Waller phi điều hòa được khai triển đến bậc bốn: $$\omega(k) = 2ik\sigma^{(1)}(T) - 2k^2\sigma^{(2)}(T) - \frac{4ik\sigma^2(T)}{R}\left(1 - \frac{R}{\lambda(k)}\right) - \frac{4}{3}ik^3\sigma^{(3)}(T) + \frac{2}{3}k^4\sigma^{(4)}(T)$$ Các cumulant bậc chẵn ($\sigma^{(2)}, \sigma^{(4)}$) đóng góp vào sự suy giảm biên độ, các cumulant bậc lẻ ($\sigma^{(1)}, \sigma^{(3)}$) đóng góp vào độ dịch pha dao động.
  • Thiết lập công thức giải tích cho cumulant theo mô hình ACEM: Sử dụng lý thuyết nhiễu loạn ma trận mật độ bậc một cho hệ lập phương tâm diện (fcc), tác giả tính toán chi tiết các yếu tố ma trận chuyển mức $\langle n | y | n' \rangle, \langle n | y^3 | n' \rangle$ qua toán tử sinh hủy phonon $\hat{a}, \hat{a}^+$:
    • Cumulant bậc một (độ giãn nở nhiệt mạng): $$a(T) = \sigma^{(1)}(T) = \frac{3\alpha \hbar \omega_E}{40 D \alpha^2} \frac{1 + z}{1 - z} = \frac{3\alpha}{40 D \alpha^2} \sigma^{(2)}(T)$$
    • Cumulant bậc hai (MSRD / hệ số Debye-Waller): $$\sigma^{(2)}(T) = \sigma_0^2 \frac{1 + z}{1 - z}, \quad \text{với } \sigma_0^2 = \frac{\hbar}{2\mu\omega_E}$$
    • Cumulant bậc ba: $$\sigma^{(3)}(T) = \frac{5\alpha (\sigma_0)^6}{\hbar\omega_E} \frac{1 + 4z + z^2}{(1 - z)^3} = \frac{5\alpha (\hbar\omega_E)^2}{400 D^2 \alpha^3} \frac{1 + 4z + z^2}{(1 - z)^2}$$
  • Hệ số giãn nở nhiệt và giới hạn lượng tử: Thiết lập biểu thức hệ số giãn nở nhiệt tuyến tính $\alpha_T = \frac{1}{r_0}\frac{da}{dT}$ và hệ số giãn nở nhiệt khối $\alpha_V = 3\alpha_T$. Phân tích hành vi tiệm cận tại các vùng nhiệt độ giới hạn:
    • Khi $T \to 0\text{ K}$ ($z \to 0$): $\sigma^{(1)} \to \frac{3\alpha\hbar\omega_E}{40D\alpha^2}$, $\sigma^{(2)} \to \sigma_0^2$, $\sigma^{(3)} \to \frac{5\alpha\sigma_0^6}{\hbar\omega_E}$ (đặc trưng cho dao động không điểm lượng tử).
    • Khi $T \to \infty$ ($z \approx 1 - \frac{\hbar\omega_E}{k_B T}$): các cumulant phụ thuộc tuyến tính hoặc bậc cao vào nhiệt độ cổ điển: $\sigma^{(1)} \sim T$, $\sigma^{(2)} \sim T$, $\sigma^{(3)} \sim T^2$.
+-----------------------------------------------------------------------------+
| Bảng 2.1: Biểu thức giải tích các đại lượng nhiệt động theo ACEM (Chương 2) |
+----------------------------+------------------------------------------------+
| Đại lượng                  | Biểu thức giải tích                            |
+----------------------------+------------------------------------------------+
| Cumulant bậc 1 (a)         | a = [3*alpha*hbar*omega_E / (40*D*alpha^2)] *  |
|                            |     [(1 + z) / (1 - z)]                        |
| Cumulant bậc 2 (sigma^2)   | sigma^2 = sigma_0^2 * [(1 + z) / (1 - z)]      |
| Cumulant bậc 3 (sigma^(3)) | sigma^(3) = [5*alpha*(sigma_0)^6 /             |
|                            |     (hbar*omega_E)] * [(1 + 4z + z^2)/(1 - z)^3|
| Hệ thức tương quan         | a = [3*alpha / (40*D*alpha^2)] * sigma^2       |
+----------------------------+------------------------------------------------+

Chương 3: Sự phụ thuộc vào áp suất, tỷ lệ pha tạp và nhiệt độ của các đại lượng nhiệt động trong phổ EXAFS phi điều hòa

Chương 3 mở rộng lý thuyết để khảo sát ảnh hưởng đồng thời của áp suất và tỷ lệ pha tạp:

  • Sự phụ thuộc vào áp suất của các cumulant: Dưới tác dụng của áp suất ngoài $P$, thể tích tinh thể co lại $V(P)$, hằng số mạng $a(P)$ và khoảng cách cân bằng $r_0(P)$ giảm, dẫn đến sự tăng của hằng số lực đàn hồi hiệu dụng $k_{eff}(P)$, tần số dao động Einstein $\omega_E(P)$ và Debye $\omega_D(P)$. Kết quả là độ dịch chuyển tương đối $\sigma^2(T, P)$ giảm khi áp suất tăng tại cùng một nhiệt độ.
  • Sự phụ thuộc vào tỷ lệ pha tạp $x$ trong hợp kim nhị phân: Đối với hợp kim thay thế $\text{A}x\text{B}{1-x}$, các đại lượng hiệu dụng (khối lượng rút gọn $\mu(x)$, hằng số lực $k_{eff}(x)$, tham số thế phi điều hòa $k_3(x), k_4(x)$) được xác định thông qua tổ hợp có trọng số của các liên kết đơn cặp A-A, B-B và A-B. Tác giả xây dựng các hệ thức cho hệ số giãn nở nhiệt $\alpha_T(x, T, P)$, hệ số phi điều hòa và pha dao động $\phi(k, x, T, P)$.
  • Biểu thức giải tích cho các cấu trúc mạng lập phương: Thiết lập tường minh các tham số hiệu dụng cho 3 kiểu mạng tinh thể:
    • Mạng lập phương tâm khối (bcc): Phối trí 8 nguyên tử lân cận gần nhất, hình chiếu liên kết dọc trục.
    • Mạng lập phương tâm diện (fcc): Phối trí 12 nguyên tử lân cận gần nhất.
    • Mạng lập phương đơn giản (sc): Phối trí 6 nguyên tử lân cận.

Chương 4: Áp dụng tính số đối với một số tinh thể và hợp kim liên kim loại

Chương 4 kiểm chứng mô hình lý thuyết thông qua tính số trị cụ thể trên các hệ tinh thể Si, Ge, hợp kim SiGe và hợp kim $\text{Cu}x\text{Ag}{1-x}$ dưới áp suất lên đến $14\text{ GPa}$:

  • Phương pháp lý thuyết xác định tham số thế Morse: Tính toán các tham số $D, \alpha, r_0$ cho tinh thể Si, Ge và SiGe trực tiếp từ năng lượng kết khối, suất nén và hằng số mạng mà không cần số liệu khớp thực nghiệm trước. Kết quả so sánh trong Bảng 4.1 và Bảng 4.2 cho thấy sự phù hợp cao giữa lý thuyết và thực nghiệm.
+-----------------------------------------------------------------------------+
| Bảng 4.1: Tham số thế Morse hiệu dụng lý thuyết (LT) và thực nghiệm (TN)    |
|           của tinh thể Si, Ge và hợp kim SiGe [28, 60]                      |
+-------------------+--------------------+------------------------------------+
| Tinh thể          | D (eV) (LT / TN)   | alpha (A^-1) (LT / TN)             |
+-------------------+--------------------+------------------------------------+
| Si                | 1,8400 / 1,8400    | 1,4600 / 1,4600                    |
| Ge                | 1,6300 / 1,6300    | 1,3900 / 1,3900                    |
| SiGe              | 1,7350 / -         | 1,4250 / -                         |
+-------------------+--------------------+------------------------------------+

+-----------------------------------------------------------------------------+
| Bảng 4.2: Hằng số lực đàn hồi hiệu dụng (x 10^-11 N/m) của Si và Ge [68]   |
+-------------------+-------------------------+-------------------------------+
| Tinh thể          | Lý thuyết hiện tại (LT) | Thực nghiệm (TN)              |
+-------------------+-------------------------+-------------------------------+
| Si                | 4,7200                  | 4,7000                        |
| Ge                | 3,8500                  | 3,8300                        |
+-------------------+-------------------------+-------------------------------+
  • Khảo sát hệ hợp kim $\text{Cu}x\text{Ag}{1-x}$ dưới áp suất $0 - 14\text{ GPa}$:
    • Tính số các tham số thế Morse và hằng số lực đàn hồi hiệu dụng $k_{eff}$ dưới các mức áp suất $0, 2, 4, 6, 8, 10, 12, 14\text{ GPa}$ (Bảng 4.3, 4.4, 4.5, 4.6).
    • Khảo sát sự phụ thuộc nhiệt độ của độ dịch chuyển trung bình bình phương MSD $u^2(T, P)$, độ dịch chuyển tương đối trung bình bình phương MSRD $\sigma^2(T, P)$ và hàm dịch chuyển tương quan DCF $C_R(T, P)$ cho các thành phần $x = 0,72$ và $x = 0,50$ ($\text{Cu}{50}\text{Ag}{50}$).
    • Tỷ lệ tương quan $\sigma^2 / (2 u^2 C_R)$ (hoặc $\sigma^2 / u^2$) phản ánh mức độ đóng góp tương quan giữa các nguyên tử lân cận. Áp suất cao làm tăng độ cứng mạng tinh thể (tăng $k_{eff}$), làm giảm biên độ dao động nhiệt $u^2$ và $\sigma^2$, đồng thời làm dịch chuyển nhiệt độ Einstein và Debye về phía giá trị cao hơn.
    • Phân tích cơ chế mất ổn định và phá vỡ cấu trúc mạng của hợp kim liên kim loại khi chịu đồng thời nhiệt độ cao và áp suất lớn.

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

Luận án đã đạt được các kết quả khoa học chính và những đóng góp mới sau đây (được nêu trực tiếp trong phần đóng góp mới của văn bản):

  1. Xây dựng hệ thức giải tích hoàn chỉnh: Thiết lập thành công biểu thức giải tích mô tả sự phụ thuộc đồng thời vào nhiệt độ, tỷ lệ pha tạp và áp suất ngoài đối với các cumulant ($\sigma^{(1)}, \sigma^{(2)}, \sigma^{(3)}, \sigma^{(4)}$) và các tham số nhiệt động (hệ số giãn nở nhiệt, hệ số Debye-Waller, hằng số lực đàn hồi) trong phổ EXAFS phi điều hòa. Các biểu thức này áp dụng cho kim loại nguyên chất và hợp kim liên kim loại có các cấu trúc lập phương khác nhau (bcc, fcc, sc).
  2. Làm rõ bản chất tương quan phi điều hòa: Xác định và phân tích rõ sự khác biệt bản chất giữa mô hình dao động tương quan phi điều hòa (ACEM, ACDM) và mô hình dao động đơn hạt điều hòa/phi điều hòa; giải thích cơ chế biến dạng và phá vỡ cấu trúc của hợp kim liên kim loại dưới tác động kết hợp của áp suất cao (đến $14\text{ GPa}$) và sự gia tăng nhiệt độ.
  3. Phương pháp lý thuyết xác định tham số thế Morse: Xây dựng thành công phương pháp tính các tham số thế tương tác Morse ($D, \alpha, r_0$) thuần túy bằng lý thuyết thông qua các đại lượng vật lý vĩ mô (năng lượng phát xạ/kết khối, hệ số nén và hằng số mạng tinh thể), cho phép tính toán định lượng cho các hệ vật liệu mới hoặc hệ hợp kim chưa có sẵn số liệu đo đạc thực nghiệm.

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

Hạn chế của nghiên cứu

Căn cứ trên phạm vi và nội dung trình bày trong văn bản:

  • Mô hình lý thuyết chủ yếu tập trung khảo sát các vật liệu có cấu trúc mạng đối xứng cao thuộc hệ tinh thể lập phương (bcc, fcc, sc, kiểu kim cương); chưa mở rộng khảo sát chi tiết cho các hệ tinh thể có tính bất đối xứng cao hơn (như lục giác, trực thoi, đơn tà).
  • Khai triển thế năng tương tác hiệu dụng dừng lại ở phép gần đúng đến số hạng bậc bốn; các đóng góp phi điều hòa bậc cao hơn bị lược bỏ.
  • Đối tượng tính số tập trung vào hợp kim nhị phân tinh thể; chưa xét đến các hiệu ứng cấu trúc phức tạp trong vật liệu vô định hình hoặc hệ cấu trúc nano đa tầng.

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

  • Phát triển và mở rộng mô hình ACEM và ACDM dưới áp suất cho các tinh thể phi lập phương, các cấu trúc dị thể và vật liệu màng mỏng nano.
  • Tiếp tục hoàn thiện lý thuyết khai triển cumulant bậc cao cho các điều kiện môi trường siêu áp suất và dải nhiệt độ cực hạn.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao đối với các đối tượng và lĩnh vực sau:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý lý thuyết, Vật lý toán, Vật lý chất rắn: Cung cấp khung lý thuyết mẫu mực về ứng dụng phương pháp thống kê lượng tử, lý thuyết nhiễu loạn ma trận mật độ và phép lượng tử hóa thứ cấp để giải quyết bài toán dao động nhiệt phi điều hòa trong vật rắn (đặc biệt là nội dung Chương 2 và Chương 3).
  • Các nhà nghiên cứu thực nghiệm Quang phổ tia X (XAFS/EXAFS/XANES): Cung cấp các công thức giải tích chuẩn xác để làm khớp phổ (fitting), hiệu chỉnh sai số phi điều hòa ở nhiệt độ cao và trích xuất chính xác các thông số cấu trúc vi mô cục bộ (khoảng cách liên kết, số phối trí, hệ số Debye-Waller) dưới điều kiện nhiệt độ và áp suất thay đổi.
  • Các nhà khoa học tính toán vật liệu: Cung cấp quy trình xác định tham số thế tương tác Morse bằng lý thuyết ở Chương 4 và hệ thống số liệu tham chiếu về hằng số lực hiệu dụng, độ dịch chuyển tương đối MSRD của các hệ Si, Ge, SiGe, Cu, Ag và Cu-Ag dưới áp suất cao.

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

1. Tại sao ở nhiệt độ cao bắt buộc phải xét đến hiệu ứng phi điều hòa trong phân tích phổ EXAFS?

Ở nhiệt độ thấp, biên độ dao động nhiệt của nguyên tử nhỏ, thế năng tương tác có dạng đối xứng quanh vị trí cân bằng (gần đúng điều hòa). Khi nhiệt độ tăng cao, tương tác phonon-phonon trở nên mạnh mẽ, làm thế tương tác trở nên bất đối xứng. Nếu tiếp tục dùng mô hình điều hòa, việc phân tích phổ EXAFS sẽ dẫn đến sai số nghiêm trọng: khoảng cách liên kết biểu kiến bị rút ngắn giả tạo và số phối trí thực nghiệm bị đánh giá thấp hơn thực tế. Phép khai triển cumulant phi điều hòa cho phép hiệu chỉnh chính xác các sai số này.

2. Sự khác biệt căn bản giữa mô hình ACEM và ACDM trong luận án là gì?

Cả hai mô hình đều tính đến tương tác tương quan của cụm nguyên tử lân cận gần nhất và hiệu ứng phi điều hòa. Tuy nhiên:

  • Mô hình ACEM (Einstein tương quan phi điều hòa) giả định các dao động có cùng một tần số Einstein $\omega_E$ và bỏ qua sự tán sắc phonon, giúp đơn giản hóa việc tính toán giải tích các cumulant bậc cao ($\sigma^{(1)}, \sigma^{(2)}, \sigma^{(3)}, \sigma^{(4)}$).
  • Mô hình ACDM (Debye tương quan phi điều hòa) xét đến sự tán sắc của phonon trong dải tần số liên tục từ $0$ đến tần số Debye cực đại $\omega_D$, phản ánh sát hơn phổ dao động mạng thực tế, thường áp dụng hiệu quả cho cumulant bậc hai $\sigma^{(2)}$ dưới áp suất.

3. Phương pháp xác định tham số thế Morse bằng lý thuyết trong luận án dựa trên nguyên lý nào?

Luận án thiết lập mối liên hệ giải tích giữa 3 tham số của thế Morse ($D, \alpha, r_0$) với các đặc trưng vĩ mô của tinh thể: năng lượng kết khối (tương ứng năng lượng phân ly $D = U(r_0)$), suất nén tĩnh học $K$ (liên quan đến đạo hàm bậc hai của thế và tham số độ rộng $\alpha$), và hằng số mạng tinh thể đo được ở trạng thái cân bằng ($r_0$). Phương pháp này cho phép tính toán các tham số thế mà không cần dựa vào dữ liệu thực nghiệm đo phổ EXAFS từ trước.

4. Áp suất cao lên đến 14 GPa tác động như thế nào đến các đại lượng nhiệt động trong phổ EXAFS?

Khi áp suất tăng từ $0\text{ GPa}$ đến $14\text{ GPa}$, thể tích tinh thể bị nén, khoảng cách giữa các nguyên tử co ngắn lại làm tăng độ dốc của thế tương tác hiệu dụng, dẫn đến hằng số lực đàn hồi hiệu dụng $k_{eff}$ tăng lên. Hệ quả là tần số dao động riêng và nhiệt độ Debye/Einstein tăng, làm biên độ dao động nhiệt, độ dịch chuyển trung bình bình phương (MSD, $u^2$) và độ dịch chuyển tương đối trung bình bình phương (MSRD, $\sigma^2$) đều giảm đáng kể tại cùng một nhiệt độ xác định.

5. Ý nghĩa vật lý của các cumulant từ bậc một đến bậc bốn trong phổ EXAFS là gì?

  • Cumulant bậc một $\sigma^{(1)} = a(T)$: Độ giãn nở nhiệt của mạng tinh thể (độ dịch chuyển vị trí cân bằng trung bình).
  • Cumulant bậc hai $\sigma^{(2)} = \sigma^2$: Độ dịch chuyển tương đối trung bình bình phương (MSRD), đóng vai trò chính trong hệ số Debye-Waller điều hòa và phi điều hòa, làm suy giảm biên độ phổ.
  • Cumulant bậc ba $\sigma^{(3)}$: Thể hiện độ bất đối xứng của hàm phân bố khoảng cách nguyên tử do tính phi điều hòa, đóng góp trực tiếp vào độ dịch pha của phổ EXAFS.
  • Cumulant bậc bốn $\sigma^{(4)}$: Thể hiện độ nhọn/dẹt của hàm phân bố, đóng góp vào sự hiệu chỉnh biên độ phổ ở nhiệt độ cao.

Kết luận

Luận án tiến sĩ của tác giả Trịnh Phi Hiệp đã giải quyết trọn vẹn bài toán lý thuyết về ảnh hưởng đồng thời của hiệu ứng phi điều hòa, nhiệt độ, tỷ lệ pha tạp và áp suất cao lên các đại lượng nhiệt động và các cumulant trong phổ EXAFS của các tinh thể có cấu trúc lập phương. Bằng việc kết hợp phương pháp thống kê lượng tử, lý thuyết nhiễu loạn với các mô hình tương quan ACEM và ACDM, luận án đã xây dựng thành công hệ thống biểu thức giải tích tường minh và phương pháp xác định tham số thế Morse thuần túy lý thuyết. Các kết quả tính số cho Si, Ge, SiGe và hợp kim $\text{Cu}x\text{Ag}{1-x}$ dưới áp suất lên đến $14\text{ GPa}$ có độ tương thích cao với thực nghiệm, đóng góp hoàn thiện cơ sở lý thuyết của quang phổ EXAFS hiện đại.