Tổng quan luận án

Luận án tiến sĩ với đề tài "Nghiên cứu đặc tính pha và chuyển pha dị thường trong mạng nano" được thực hiện bởi nghiên cứu sinh Nguyễn Thị Kim Oanh, dưới sự hướng dẫn khoa học của GS.TS Bạch Thành Công và TS. Đặng Đình Long, tại chuyên ngành Vật liệu và linh kiện nano (Mã số: 944012801.QTD), Trường Đại học Công nghệ phối hợp cùng Khoa Vật lý - Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội, với sự hỗ trợ kinh phí từ đề tài NAFOSTED mã số 103.324.

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

Trong vật lý chất rắn và vật lý vật chất ngưng tụ hiện đại, các hệ thống tương quan mạnh và các mô hình mạng spin đóng vai trò nền tảng để dự đoán, giải thích các hành vi vật lý phức tạp. Hai vấn đề nổi bật thu hút sự quan tâm lớn bao gồm:

  1. Hiện tượng chuyển pha từ loại I với các bước nhảy từ hóa dị thường dưới tác dụng của từ trường ngoài trong các hệ từ tính có cấu trúc nano phân tách pha (như perovskite Mangan $\text{Pr}_{1-x}\text{Ca}_x\text{MnO}_3$) hoặc các hệ có tính chất vấp spin hình học (như hợp chất $\text{SrCu}_2(\text{BO}_3)_2$ và tetraboride đất hiếm $\text{RB}_4$ trên mạng Shastry – Sutherland). Cơ chế vi mô dẫn đến các bước nhảy từ hóa phân số và sự biến đổi của điện trở suất trong các vật liệu này đòi hỏi phải được mô hình hóa định lượng thông qua sự cạnh tranh tương tác và thăng giáng của tích phân trao đổi.
  2. Sự hình thành và ổn định của pha siêu tinh thể (Supersolid - SS) trong hệ hạt boson tương tác mạnh. Pha siêu tinh thể là trạng thái lượng tử dị thường, nơi trật tự đường chéo tầm xa (DLRO - đặc trưng cho cấu trúc tinh thể định xứ) và trật tự ngoài đường chéo tầm xa (ODLRO - đặc trưng cho tính siêu chảy SF) cùng tồn tại đồng thời. Mặc dù các nghiên cứu trước đây trên mô hình Bose – Hubbard với tương tác tầm xa đã tìm thấy pha siêu tinh thể, nhưng trạng thái này chủ yếu chỉ xuất hiện ở phía trên mật độ tinh thể do khuyết tật giữa nút; ở phía dưới mật độ tinh thể, hệ thường rơi vào trạng thái phân tách pha. Khoảng trống đặt ra là việc xác định cơ chế ngoại vi, cụ thể là tác dụng của thế năng ghim tuần hoàn, nhằm ổn định pha siêu tinh thể ở cả hai phía mật độ tinh thể (khuyết tật lỗ trống và khuyết tật giữa nút).

Trong khuôn khổ luận án, thuật ngữ "mạng nano" không quy ước cho một kích thước cấu trúc cụ thể cố định, mà mô tả các màng mỏng hai chiều/giả hai chiều có độ dày cỡ nanomet hoặc các vật liệu khối tồn tại sự phân tách pha với các miền pha từ tính, pha trật tự điện tích có kích thước nanomet cạnh tranh nhau.

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

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

  1. Khảo sát ảnh hưởng của xác suất phân bố và độ thăng giáng trong tương tác lên đường cong từ hóa và đường cong từ trở trong mô hình spin Ising cho mạng vuông và mạng Shastry – Sutherland; từ đó mô tả các hiện tượng tương đồng trong giản đồ trạng thái giữa mô hình lý thuyết với kết quả thực nghiệm, góp phần giải thích nguồn gốc các bước nhảy từ trong lớp vật liệu đa tinh thể perovskite Mangan và tetraboride đất hiếm $\text{RB}_4$.
  2. Khảo sát giản đồ pha và các tham số trật tự đặc trưng trong mô hình Bose – Hubbard của các hạt boson dưới tác dụng của thế năng ghim tuần hoàn; từ đó tìm ra các điều kiện cần thiết để ổn định pha siêu tinh thể trong mô hình.

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

  • Đối tượng nghiên cứu: Các đặc tính pha, quá trình chuyển pha nhiệt động và chuyển pha lượng tử trong hai mô hình mạng hai chiều: mô hình spin Ising có cạnh tranh tương tác/mất trật tự và mô hình hạt boson Bose – Hubbard (lõi cứng và lõi mềm) dưới tác dụng của thế năng ghim tuần hoàn.
  • Phạm vi nghiên cứu:
    • Không gian: Các mô hình mạng phẳng hai chiều bao gồm mạng vuông hai chiều tương tác lân cận gần nhất (NN), mạng vấp spin Shastry – Sutherland với tương tác lân cận gần nhất (NN) và lân cận gần nhì (NNN), cùng mạng quang học hai chiều cho hạt boson.
    • Điều kiện vật lý và nhiệt động: Khảo sát trong vùng nhiệt độ từ nhiệt độ hữu hạn đến không độ tuyệt đối ($T \to 0\text{ K}$), dưới tác động biến thiên của các tham số điều khiển bao gồm từ trường ngoài không thứ nguyên ($h$), thế hóa học trung bình ($\mu$), thế năng ghim tuần hoàn ($\varepsilon$), biên độ thăng giáng tích phân trao đổi ($\Delta$), xác suất phân bố tương tác ($p$), và cường độ tương tác đẩy tại nút ($U$) cùng tương tác tầm xa ($V_{nn}, V_{nnn}$).

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

Luận án đã tổng thuật hệ thống các công trình nghiên cứu lý thuyết và thực nghiệm trong và ngoài nước liên quan đến nhiệt động lực học chuyển pha, mô hình spin và hệ boson tương quan mạnh:

Lĩnh vực nghiên cứu Các tác giả, công trình và mốc thời gian tiêu biểu Nội dung chính đã được ghi nhận trong y văn
Nhiệt động học và chuyển pha cổ điển H. Kamerlingh Onnes (1908, 1910); W. H. Keesom (1932); P. Ehrenfest (1933); L. Landau (1937); E. Stoner; F. Arrott Khám phá tính chất hóa lỏng và chuyển pha dị thường của Helium lỏng tại điểm $\lambda$ (2,17 K); phân loại chuyển pha loại I (có ẩn nhiệt, gián đoạn đạo hàm bậc nhất của thế nhiệt động) và loại II (liên tục, gián đoạn đạo hàm bậc hai); lý thuyết Landau về tham số trật tự $\eta$ và khai triển năng lượng tự do; phép dựng đồ thị Arrott xác định nhiệt độ Curie $T_C$.
Chuyển pha lượng tử và pha dị thường P. Kapitza (1938); L. Tisza và L. Landau; F. London; O. Penrose và L. Onsager (1956); A. F. Andreev, I. M. Lifshitz và G. V. Chester (1969); S. Sachdev Khám phá hiện tượng siêu chảy (SF) không ma sát ở $\text{He}^4$; mô hình hai chất lỏng ($\rho = \rho_s + \rho_n$); liên hệ giữa ngưng tụ Bose – Einstein với trật tự tầm xa ngoài đường chéo (ODLRO); định nghĩa pha siêu tinh thể (SS) đồng thời tồn tại DLRO và ODLRO; lý thuyết chuyển pha lượng tử tại điểm tới hạn lượng tử QCP ($T = 0\text{ K}, r = r_c$).
Mô hình Ising và perovskite Mangan W. Lenz (1920); E. Ising (1925); L. Onsager (1944); H. B. Callen (1963); R. Landauer; J. Inoue và S. Maekawa; E. Dagotto et al. Lời giải chính xác mô hình Ising 2D; hệ thức Callen cho giá trị trung bình spin; lý thuyết truyền đạn đạo và hiệu ứng từ trở khổng lồ (CMR) xuyên hầm phân cực spin; hiện tượng phân tách pha kích thước nano giữa các pha sắt từ (FM), phản sắt từ (AF) và trật tự điện tích (CDW) trong manganite $\text{Pr}_{1-x}\text{Ca}_x\text{MnO}_3$.
Mạng Shastry – Sutherland và vấp spin B. S. Shastry và B. Sutherland (1981); H. Kageyama et al. (1999); Z. Fisk et al. Mô hình mạng hình học vấp spin 2D; thực nghiệm trên $\text{SrCu}_2(\text{BO}_3)_2$ và các hợp chất tetraboride đất hiếm $\text{RB}_4$ ($\text{TbB}_4, \text{ErB}_4, \text{TmB}_4$) phát hiện chuỗi bước nhảy và cao nguyên từ hóa ở các phân số $1/3, 1/2...$ của mômen từ bão hòa.
Mô hình Bose – Hubbard và thuật toán Worm M. P. A. Fisher et al. (1989); D. Jaksch et al. (1998); M. Greiner et al. (2002); N. Prokof'ev, B. Svistunov và I. Tupitsyn (1998) Mô tả hạt boson trên mạng quang học với các pha SF, điện môi Mott (MI); phát triển thuật toán Worm trong Monte Carlo lượng tử (QMC) khắc phục triệt để hiện tượng nghẽn tới hạn và bẫy cấu hình kín; các nghiên cứu QMC chỉ ra pha SS bị phân tách pha ở vùng dưới mật độ tinh thể.

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

  1. Vận dụng phương pháp giải tích trường trung bình cải biên kết hợp biến đổi Fourier hàm delta Dirac trên hệ thức Callen để xây dựng biểu thức đóng cho mômen từ và điện trở suất của mạng spin Ising vuông có cạnh tranh ngẫu nhiên FM – AF, giải thích định lượng bước nhảy từ hóa loại I và từ trở trong vật liệu đa tinh thể perovskite Mangan $\text{Pr}{0.7}\text{Ca}{0.3}\text{MnO}_3$ tại nhiệt độ thấp.
  2. Thiết lập mô hình Ising mất trật tự trên mạng Shastry – Sutherland có xét đến thăng giáng tích phân trao đổi NN và NNN, sử dụng mô phỏng Monte Carlo Metropolis để khảo sát chi tiết sự suy biến và biến dạng của cao nguyên từ hóa phân số $1/3$ dưới các điều kiện nhiệt độ và xác suất phân bố khác nhau.
  3. Ứng dụng thuật toán Worm Monte Carlo lượng tử cho mô hình Bose – Hubbard mở rộng có mặt thế năng ghim tuần hoàn trên mạng vuông hai chiều, nhằm làm rõ điều kiện vi mô triệt tiêu hiện tượng phân tách pha và ổn định pha siêu tinh thể ở cả hai vùng mật độ khuyết tật lỗ trống ($\rho < 1/2$) và khuyết tật giữa nút ($\rho > 1/2$).

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

Cơ sở lý thuyết và các khái niệm sử dụng

  • Nhiệt động lực học và lý thuyết Landau: Mô tả hệ qua các hàm nhiệt động Helmholtz ($F = U - T S_e$), Gibbs ($G = F + PV$), entropy ($S_e = -(\partial G/\partial T)_P$). Trong chuyển pha từ loại II, năng lượng tự do được khai triển theo tham số trật tự mômen từ $M$: $$F_M = \frac{1}{2} A(T) M^2 + \frac{1}{4} B(T) M^4 - \mu_0 M H$$ với $A(T) = a(T - T_C)$. Lý thuyết trường trung bình Landau suy ra các số mũ tới hạn $\beta' = 1/2$ và $\delta = 3$. Đối với chuyển pha meta từ loại I, sự xuất hiện của hệ số $B(T) < 0$ dẫn đến cấu trúc hai cực tiểu của thế nhiệt động tại các giá trị mômen từ $M_1$ và $M_2$, gây ra bước nhảy gián đoạn khi từ trường đạt giá trị tới hạn $H_C$.
  • Chuyển pha lượng tử (QCP): Quá trình chuyển pha liên tục diễn ra tại $T = 0\text{ K}$ khi biến đổi tham số điều khiển $r$ qua giá trị tới hạn $r_c$, được thúc đẩy thuần túy bởi các dao động lượng tử theo nguyên lý bất định Heisenberg, với năng lượng dao động tỷ lệ theo $\hbar \omega_c \propto |r - r_c|^{\nu z}$ ($z$ là số mũ tới hạn nhiệt động).
  • Tiêu chuẩn nhận diện pha siêu chảy (SF) và siêu tinh thể (SS):
    • Trật tự ngoài đường chéo tầm xa (ODLRO): Xác định qua ma trận mật độ đơn hạt $n(\mathbf{r}, \mathbf{r}') = \langle \psi^\dagger(\mathbf{r}) \psi(\mathbf{r}') \rangle$, tiến tới giá trị hữu hạn $n_0 = N_0/V \neq 0$ khi $|\mathbf{r} - \mathbf{r}'| \to \infty$, tương ứng mật độ siêu chảy $\rho_s > 0$.
    • Trật tự đường chéo tầm xa (DLRO): Mật độ hạt địa phương tuần hoàn $\delta\rho(\mathbf{r}) = \delta\rho(\mathbf{r} + \mathbf{T})$, xác định qua hệ số cấu trúc tĩnh $S(\mathbf{Q}) = \frac{1}{N_s} |\sum_\mathbf{r} e^{-i\mathbf{Q}\cdot\mathbf{r}} \delta\rho(\mathbf{r})|^2 \neq 0$ tại vectơ sóng mạng đảo $\mathbf{Q}$.
    • Pha SF có $\rho_s > 0, S(\mathbf{Q}) = 0$; Pha rắn/CDW có $\rho_s = 0, S(\mathbf{Q}) > 0$; Pha SS có đồng thời $\rho_s > 0$ và $S(\mathbf{Q}) > 0$.

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

Luận án sử dụng kết hợp ba phương pháp tính toán và mô phỏng số:

  1. Phương pháp giải tích trường trung bình kết hợp tính toán số trên Matlab: Áp dụng cho mô hình Ising mạng vuông $S = 1/2$ có cạnh tranh tương tác. Hamiltonian mô hình: $$H = -\sum_{\langle ik \rangle} J_{ik} \sigma_i \sigma_k - h \sum_i \sigma_i \quad (\sigma_i = \pm 1)$$ với phân bố xác suất của tích phân trao đổi: $$P(J_{ik}) = p,\delta(J_{ik} - J_{FM}) + (1-p),\delta(J_{ik} - J_{AF})$$ trong đó $J_{FM} = J(1+\Delta), J_{AF} = J(1-\Delta)$, với $J=1$ và $\Delta > 1$. Sử dụng hệ thức Callen $\langle \sigma_i \rangle = \langle \tanh(\beta E_i) \rangle$ với trường cục bộ $E_i = \sum_k J_{ik}\sigma_k + h$, kết hợp biến đổi Fourier hàm delta Dirac: $$\tanh(\beta E_i) = \frac{1}{\pi} \int_{-\infty}^{+\infty} \text{Im}\left[\frac{\exp(it E_i)}{\sinh(t/2\beta)}\right] dt$$ Thiết lập phương trình đại số tự nhất quán xác định mômen từ tỷ đối trung bình trên một nút mạng $m = \langle \sigma_i \rangle_r$: $$m = \sum_{n=0}^z C_z^n A_n(\tau, p, \Delta, z, h) m^n$$ trong đó $z=4$ là số lân cận gần nhất, $\tau = k_B T/J$, và các hệ số $A_n$ được tích phân số chính xác trên phần mềm Matlab. Hệ số từ trở tỷ đối được tính toán qua mô hình xuyên hầm Landauer – Inoue – Maekawa: $$MR = \frac{\rho_h(\tau, h)}{\rho_0(\tau)} = \frac{1}{1 + P^2 m^2}$$ với độ phân cực spin của điện tử dẫn $P = \tanh\left[\frac{J_H}{k_B T}\left(S_{t2g} + \frac{1}{2}\right)\right]$ được xác định từ phương pháp hàm Green hai thời điểm trong giới hạn tương tác trao đổi Hund mạnh $J_H \gg t$ ($S_{t2g} = 3/2$).
  2. Phương pháp mô phỏng Monte Carlo cổ điển với thuật toán Metropolis: Áp dụng cho mô hình Ising trên mạng vấp spin Shastry – Sutherland hai chiều. Khảo sát năng lượng vi mô, sự lật spin theo xác suất $W = \min(1, \exp(-\beta \Delta E))$, kiểm tra điều kiện biên tuần hoàn và phân tích hiệu ứng kích thước hữu hạn trên các kích thước mạng khác nhau ($L \times L$).
  3. Phương pháp Monte Carlo lượng tử (QMC) với thuật toán Worm: Áp dụng cho mô hình Bose – Hubbard mở rộng: $$H = -t \sum_{\langle ij \rangle} (b_i^\dagger b_j + b_j^\dagger b_i) + \frac{U}{2} \sum_i n_i(n_i - 1) + V_{nn} \sum_{\langle ij \rangle} n_i n_j + V_{nnn} \sum_{\langle\langle ik \rangle\rangle} n_i n_k - \sum_i (\mu - \mu_i) n_i$$ với $\mu_i = \pm \varepsilon$ biểu diễn thế năng ghim tuần hoàn. Thuật toán Worm làm việc trên không gian cấu hình tích phân đường (đường thế giới - worldlines), mở rộng không gian trạng thái bằng việc đưa vào hai toán tử gián đoạn (Worm head và Worm tail) để lấy mẫu hàm Green đơn hạt, thực hiện các cập nhật sinh/hủy Worm, dịch chuyển thời gian ảo và dịch chuyển không gian, giải quyết triệt để sự đóng băng cấu hình tại nhiệt độ thấp.

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

CHƯƠNG 1: TỔNG QUAN VỀ PHA VÀ CHUYỂN PHA

Chương 1 hệ thống hóa cơ sở lý thuyết về trạng thái pha và các cơ chế chuyển pha nhiệt động học và lượng tử:

  • Trình bày định nghĩa pha nhiệt động dựa trên tính đồng nhất và tính đối xứng nội tại của hệ; thiết lập các hệ thức nhiệt động cho năng lượng tự do Helmholtz $F$, năng lượng tự do Gibbs $G$, nội năng $U$, enthalpy $H_E$, entropy $S_e$ và nhiệt dung riêng $C_p$.
  • Phân loại chuyển pha theo Ehrenfest: làm rõ sự khác biệt giữa chuyển pha loại I (có biến đổi gián đoạn của đạo hàm bậc nhất thế nhiệt động, kèm theo ẩn nhiệt và sự cùng tồn tại của hai pha tại $T_C$) và chuyển pha loại II (đạo hàm bậc nhất liên tục, gián đoạn ở đạo hàm bậc hai, không có ẩn nhiệt nhưng có sự biến đổi tính đối xứng và trật tự nội tại).
  • Phân tích lý thuyết Landau về tham số trật tự $\eta$ và ứng dụng cho chuyển pha từ; thiết lập phương trình trạng thái từ nhiệt động, giải thích hiện tượng chuyển pha meta từ cảm ứng bởi từ trường ngoài thông qua sự thay đổi dấu của hệ số khai triển bậc bốn $B(T)$, dẫn đến hai cực tiểu năng lượng tự do cạnh tranh tại các giá trị từ hóa $M_1$ và $M_2$.
  • Trình bày cơ chế chuyển pha lượng tử tại lân cận điểm tới hạn lượng tử QCP ($T = 0\text{ K}$), cấu trúc giản đồ pha gồm ba vùng đặc trưng: vùng mất trật tự nhiệt động, vùng mất trật tự lượng tử và vùng lượng tử tới hạn điều khiển bởi số mũ tới hạn động học $z$.
  • Phân tích bản chất vật lý của hai pha lượng tử dị thường: pha siêu chảy SF (đặc trưng bởi tham số trật tự ODLRO, mật độ ngưng tụ $n_0$ và mô hình hai thành phần chất lỏng của Tisza – Landau) và pha siêu tinh thể SS (đặc trưng bởi sự cùng tồn tại đồng thời của trật tự không gian tầm xa DLRO với đỉnh tán xạ cấu trúc tĩnh $S(\mathbf{Q}) \neq 0$ và trật tự pha tầm xa ODLRO với $\rho_s \neq 0$).

CHƯƠNG 2: QUÁ TRÌNH TỪ HÓA LOẠI I TRONG PEROVSKITE MANGAN

Chương 2 tập trung giải quyết bài toán chuyển pha từ loại I và hiệu ứng từ trở trong vật liệu perovskite Mangan pha tạp $\text{Pr}_{1-x}\text{Ca}_x\text{MnO}_3$:

  • Phân tích nguồn gốc vi mô của sự phân tách pha kích thước nano: sự cạnh tranh giữa tương tác siêu trao đổi phản sắt từ (AF) giữa các ion $\text{Mn}^{3+}-\text{Mn}^{3+}$ hoặc $\text{Mn}^{4+}-\text{Mn}^{4+}$ với tương tác trao đổi kép sắt từ (FM) giữa các ion $\text{Mn}^{3+}-\text{Mn}^{4+}$ qua cầu oxy.
  • Xây dựng mô hình giải tích trường trung bình cho mạng Ising vuông $S = 1/2$ với số lân cận $z = 4$, trong đó tích phân trao đổi mang tính ngẫu nhiên nhị phân với xác suất $p$ cho liên kết FM và $(1-p)$ cho liên kết AF.
  • Thiết lập hệ phương trình đại số phi tuyến bậc cao cho mômen từ tỷ đối $m$ và tính toán đường cong từ trở tỷ đối $\rho_h/\rho_0$ theo từ trường ngoài $h$ ở các giá trị nhiệt độ không thứ nguyên $\tau$, độ thăng giáng tích phân trao đổi $\Delta$, và xác suất phân bố $p$.
  • Phát hiện quy luật: khi độ thăng giáng $\Delta > 1$ và ở nhiệt độ thấp ($\tau$ nhỏ), đường cong từ hóa $m(h)$ xuất hiện hai bước nhảy gián đoạn tại các giá trị từ trường tới hạn xác định, tương ứng với sự sụt giảm bậc thang đột ngột của đường cong từ trở tỷ đối $\rho_h/\rho_0$. Khi nhiệt độ $\tau$ tăng, các bước nhảy từ dần biến mất và chuyển thành biến thiên liên tục.
  • So sánh kết quả tính toán lý thuyết với số liệu đo thực nghiệm trên mẫu đa tinh thể $\text{Pr}{0.7}\text{Ca}{0.3}\text{MnO}_3$ tại nhiệt độ $T = 3\text{ K}$:
Đại lượng vật lý so sánh Kết quả tính toán theo mô hình lý thuyết ($z=4, p=0.25, \Delta=1.35, \tau=0.01$) Số liệu đo thực nghiệm trên mẫu $\text{Pr}{0.7}\text{Ca}{0.3}\text{MnO}_3$ tại $T = 3\text{ K}$
Dạng đường cong từ hóa Xuất hiện 2 bước nhảy từ gián đoạn (chuyển pha loại I) Ghi nhận 2 bước nhảy từ hóa rõ nét dưới tác dụng của từ trường
Trường tới hạn thứ nhất ($H_{C1}$) $\approx 1{,}8\text{ T}$ $\approx 1{,}9\text{ T}$
Trường tới hạn thứ hai ($H_{C2}$) $\approx 3{,}8\text{ T}$ $\approx 4{,}0\text{ T}$
Hành vi đường cong từ trở Sụt giảm mạnh dạng nhảy bậc tại các điểm trường tới hạn Điện trở suất giảm đột ngột nhiều bậc độ lớn tại lân cận $H_{C1}, H_{C2}$

CHƯƠNG 3: BƯỚC NHẢY TỪ PHÂN SỐ TRONG MÔ HÌNH ISING CHO MẠNG SHASTRY – SUTHERLAND

Chương 3 trình bày các kết quả mô phỏng Monte Carlo cổ điển (thuật toán Metropolis) cho mô hình Ising mất trật tự trên mạng vấp spin Shastry – Sutherland:

  • Mô tả cấu trúc hình học mạng Shastry – Sutherland gồm các liên kết lân cận gần nhất $J$ (cạnh ô vuông) và liên kết đường chéo lân cận gần nhì $J'$ phân bố đối xứng. Nhận diện các trạng thái cơ bản cạnh tranh: pha Ising – dimer (hai spin đối song trên đường chéo), pha Néel, và pha phân số $1/3$ (cấu hình Up – Up – Down, ký hiệu UUD, hình thành trên các tam giác con chứa đường chéo Shastry – Sutherland).
  • Đánh giá hiệu ứng kích thước hữu hạn trên các mạng $L \times L$ và xác định điều kiện cân bằng nhiệt động cho hệ mô phỏng.
  • Phân tích ảnh hưởng của tỷ số tương tác $J/J'$:
    • Khi tỷ số nhỏ ($J/J' = 0{,}01$ đến $0{,}1$), đường cong từ hóa $m(h)$ biểu hiện bậc nhảy cao nguyên rõ rệt tại giá trị phân số $m/m_s = 1/3$.
    • Khi tỷ số tương tác tăng lên ($J/J' = 1$), cấu trúc bậc nhảy phân số bị biến dạng do sự cạnh tranh mạnh mẽ giữa các trật tự spin lân cận.
  • Khảo sát ảnh hưởng của sự mất trật tự và thăng giáng $\Delta$:
    • Thêm thăng giáng vào tương tác NN ($J$), tương tác NNN ($J'$), hoặc đồng thời cả hai loại tương tác với các xác suất phân bố $p = 0{,}2; 0{,}5; 0{,}8$.
    • Kết quả chỉ ra rằng độ thăng giáng $\Delta$ và phân bố ngẫu nhiên $p$ làm thay đổi độ rộng của cao nguyên $1/3$, làm dịch chuyển vị trí các điểm nhảy từ trường tới hạn $h_c$, và ở nhiệt độ hữu hạn ($\tau = 0{,}2; 0{,}5$) các bậc nhảy phân số bị làm mờ, chuyển dần sang dạng biến thiên trơn.
  • Tính toán năng lượng trung bình trên mỗi nút mạng $E/J'$ theo xác suất $p$ tại $J/J' = 0{,}01, h = 0$, chứng minh sự biến đổi tuyến tính của năng lượng trạng thái cơ bản khi thay đổi tỷ lệ thành phần tương tác ngẫu nhiên.

CHƯƠNG 4: CHUYỂN PHA TRONG MÔ HÌNH HẠT BOSON BOSE – HUBBARD DƯỚI TÁC DỤNG CỦA THẾ NĂNG GHIM TUẦN HOÀN

Chương 4 khảo sát chi tiết các pha lượng tử và điều kiện ổn định pha siêu tinh thể trong mô hình Bose – Hubbard bằng thuật toán Monte Carlo lượng tử Worm:

  • Thiết lập mô hình hạt boson tương tác mạnh trên mạng vuông hai chiều chịu tác dụng của thế năng ghim tuần hoàn $\mu_i = \pm \varepsilon$ phân bố luân phiên trên các nút mạng.
  • Khảo sát mô hình hạt boson lõi cứng (Hard-core bosons):
    • Khảo sát đường cong mật độ hạt $\rho$ theo thế hóa học $\mu$ với tương tác lân cận $V_{nn} = 6t$. Khi chưa có thế ghim ($\varepsilon = 0$), hệ chỉ xuất hiện pha siêu tinh thể ở vùng mật độ $\rho > 1/2$ (gây ra bởi khuyết tật giữa nút), trong khi vùng $\rho < 1/2$ xảy ra hiện tượng phân tách pha (nhảy bậc mật độ gián đoạn).
    • Khi áp đặt thế năng ghim tuần hoàn ($\varepsilon = 0{,}5t; 1{,}0t$), đường cong $(\rho, \mu)$ trở nên liên tục ở cả hai phía của mật độ tinh thể $\rho = 1/2$.
    • Tính toán đồng thời mật độ siêu chảy $\rho_s$ và hệ số cấu trúc tĩnh $S(\mathbf{Q})$ tại $\mathbf{Q} = (\pi, \pi)$, xác nhận sự tồn tại đồng thời của $\rho_s > 0$ và $S(\pi, \pi) > 0$ ở cả hai vùng $\rho < 1/2$ (khuyết tật lỗ trống) và $\rho > 1/2$ (khuyết tật giữa nút).
  • Khảo sát mô hình hạt boson lõi mềm (Soft-core bosons):
    • Xây dựng giản đồ pha trạng thái cơ bản $(U, \varepsilon)$ phân định ranh giới giữa pha siêu chảy SF và pha sóng mật độ điện tích (CDW) tại các mật độ hạt $\rho = 1/2, 1, 3/2$.
    • Phân tích các cấu hình trật tự cấu trúc tinh thể của hạt boson: dạng ô bàn cờ (checkerboard), dạng sọc (stripe), và dạng sao (star).
    • Khảo sát đường cong $(\rho, \mu)$ và các tham số trật tự $\rho_s, S(\pi, 0)$ dưới tác dụng của tương tác đẩy tại nút $U = 4t$ và $U = 8t$ kết hợp thế ghim $\varepsilon = 20t$.
    • Thực hiện phép ngoại suy kích thước vô hạn ($1/L \to 0$) cho tham số trật tự $\rho_s$ và $S(\pi, 0)$ tại mật độ hạt $\rho \approx 0{,}85$, chứng minh trật tự siêu tinh thể SS duy trì giá trị khác không trong giới hạn nhiệt động học vô hạ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 cụ thể về mặt lý luận và mô hình hóa tính toán:

Đóng góp về mặt lý luận khoa học

  1. Xây dựng và hoàn thiện phương pháp giải tích trường trung bình cải biên dựa trên hệ thức Callen kết hợp biến đổi tích phân Fourier hàm delta Dirac cho mô hình spin Ising có cạnh tranh tương tác ngẫu nhiên FM – AF; thiết lập biểu thức toán học xác định đồng thời đường cong từ hóa và hệ số từ trở tỷ đối phụ thuộc vào nhiệt độ, từ trường ngoài và độ thăng giáng tương tác.
  2. Làm rõ cơ chế hình thành các bước nhảy từ hóa loại I trong mạng nano perovskite Mangan phân tách pha, chỉ ra rằng sự tồn tại đồng thời của thăng giáng tích phân trao đổi $\Delta > 1$ và tỷ lệ tương tác cạnh tranh $p$ là điều kiện cần thiết để xuất hiện các bước nhảy từ hóa và sự sụt giảm điện trở suất gián đoạn tại các giá trị trường tới hạn $H_C$.
  3. Xác định quy luật biến dạng của cao nguyên từ hóa phân số $1/3$ (cấu hình spin UUD) trên mạng vấp spin Shastry – Sutherland dưới ảnh hưởng độc lập và đồng thời của sự mất trật tự trên các liên kết NN ($J$) và NNN ($J'$).
  4. Chứng minh bằng phương pháp Monte Carlo lượng tử (thuật toán Worm) rằng thế năng ghim tuần hoàn $\varepsilon$ có khả năng triệt tiêu trạng thái phân tách pha trong mô hình Bose – Hubbard, tạo ra cơ chế ổn định pha siêu tinh thể (SS) bền vững ở cả hai phía mật độ tinh thể tương ứng với khuyết tật lỗ trống và khuyết tật giữa nút trên cả hệ boson lõi cứng và boson lõi mềm.

Ý nghĩa thực tiễn và định hướng ứng dụng

  • Cung cấp mô hình tính toán lý thuyết định lượng phù hợp với các dữ liệu đo đạc thực nghiệm từ hóa và từ trở trên hệ vật liệu perovskite Mangan đa tinh thể $\text{Pr}{0.7}\text{Ca}{0.3}\text{MnO}_3$ và các vật liệu vấp spin họ tetraboride $\text{RB}_4$.
  • Đưa ra các giá trị ngưỡng tham số vi mô (tỷ số năng lượng $U/t, V_{nn}/t$, biên độ thế ghim $\varepsilon/t$) làm cơ sở định hướng cho các nghiên cứu thực nghiệm thiết lập thông số trên mạng quang học nguyên tử siêu lạnh nhằm quan sát và kiểm chứng pha siêu tinh thể.

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

Hạn chế của luận án

  • Các mô hình tính toán chủ yếu tập trung trên không gian mạng hai chiều phẳng ($2D$), chưa mở rộng đầy đủ sang cấu trúc không gian ba chiều ($3D$) của mạng tinh thể thực.
  • Trong mô hình Ising cho perovskite Mangan, việc sử dụng gần đúng trường trung bình với số phối trí lân cận cố định ($z=4$) và giả thiết spin định xứ cổ điển chưa xét đến đầy đủ các tương tác tầm xa hơn cũng như dao động mạng động học (hiệu ứng Jahn – Teller động).
  • Mô phỏng số Monte Carlo Metropolis và QMC Worm chịu giới hạn về kích thước mạng hữu hạn, đòi hỏi kỹ thuật ngoại suy kích thước để đánh giá giới hạn nhiệt động.

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

  • Mở rộng khảo sát mô hình spin mất trật tự có cạnh tranh tương tác trên các mạng hình học vấp spin ba chiều hoặc các mạng phẳng phức tạp khác như mạng Kagome, mạng Pyrochlore.
  • Khảo sát mô hình Bose – Hubbard với các dạng thế năng ngoài ngẫu nhiên (disorder potential) hoặc thế năng quasi-periodic (mô hình Aubry – André) để nghiên cứu sự cạnh tranh giữa pha siêu tinh thể và pha thủy tinh Bose (Bose glass).
  • Đưa thêm tương tác spin – quỹ đạo và các hiệu ứng lượng tử bậc cao vào mô hình giải tích cho màng mỏng từ tính nano.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị tham khảo cụ thể cho các đối tượng và nội dung sau:

  • Nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn, Vật lý tính toán, Vật liệu nano: Tham khảo phương pháp giải tích trường trung bình dựa trên hệ thức Callen cho hệ spin mất trật tự; quy trình xây dựng thuật toán mô phỏng Monte Carlo Metropolis cổ điển và kỹ thuật mô phỏng Monte Carlo lượng tử Worm giải bài toán tích phân đường cho hệ boson tương quan mạnh.
  • Giảng viên và nhà nghiên cứu lý thuyết: Sử dụng làm tài liệu giảng dạy và nghiên cứu chuyên sâu về lý thuyết chuyển pha Landau, chuyển pha lượng tử xung quanh điểm tới hạn lượng tử QCP, các tiêu chuẩn trật tự DLRO/ODLRO, và động lực học hạt boson trên mạng quang học.
  • Nhà nghiên cứu thực nghiệm vật liệu từ và mạng quang học: Tham khảo các phân tích về cơ chế bước nhảy từ hóa loại I, hành vi từ trở khổng lồ trong perovskite Mangan, cấu trúc trạng thái vấp spin Shastry – Sutherland, cùng các điều kiện tham số thế ghim tuần hoàn phục vụ thiết kế thí nghiệm nguyên tử siêu lạnh.

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

1. Cơ chế vi mô nào dẫn đến các bước nhảy từ hóa loại I trong hệ perovskite Mangan được mô hình hóa trong luận án?

Trong perovskite Mangan pha tạp $\text{Pr}_{1-x}\text{Ca}_x\text{MnO}_3$, sự cùng tồn tại của các ion $\text{Mn}^{3+}$ và $\text{Mn}^{4+}$ tạo ra sự cạnh tranh giữa tương tác siêu trao đổi phản sắt từ (AF) giữa các ion cùng hóa trị và tương tác trao đổi kép sắt từ (FM) giữa các ion khác hóa trị. Sự phân bố ngẫu nhiên của các đám pha nano FM và AF kết hợp với độ thăng giáng tích phân trao đổi $\Delta > 1$ tại nhiệt độ thấp làm xuất hiện chuyển pha meta từ loại I; khi từ trường ngoài đạt các giá trị tới hạn $H_C$, cấu trúc spin của hệ biến đổi đột ngột qua các bước nhảy gián đoạn của mômen từ tỷ đối $m$, kéo theo sự sụt giảm bậc thang của điện trở suất.

2. Cấu hình spin Up – Up – Down (UUD) trên mạng Shastry – Sutherland có ý nghĩa gì đối với bước nhảy từ phân số 1/3?

Mạng Shastry – Sutherland có cấu trúc vấp spin hình học tạo bởi các liên kết lân cận gần nhất $J$ và liên kết đường chéo $J'$. Ở tỷ số tương tác thích hợp và dưới tác dụng của từ trường ngoài, hệ không thể thỏa mãn đồng thời tất cả các liên kết phản sắt từ mà ưu tiên hình thành cấu hình gồm hai spin hướng lên và một spin hướng xuống (cấu hình UUD) trong mỗi tam giác con chứa đường chéo. Cấu hình này tạo nên trạng thái cơ bản có mômen từ bão hòa đạt chính xác giá trị phân số $m/m_s = 1/3$, biểu hiện thành cao nguyên từ hóa (plateau) trên đường cong từ hóa $m(h)$.

3. Hai điều kiện tham số trật tự bắt buộc để xác nhận sự xuất hiện của pha siêu tinh thể (SS) là gì?

Pha siêu tinh thể đòi hỏi sự cùng tồn tại đồng thời của hai loại trật tự vốn loại trừ lẫn nhau:

  • Trật tự đường chéo tầm xa (DLRO): Thể hiện tính chất cấu trúc tinh thể định xứ, được xác định khi hệ số cấu trúc tĩnh $S(\mathbf{Q}) > 0$ tại vectơ sóng mạng đảo $\mathbf{Q}$ đặc trưng.
  • Trật tự ngoài đường chéo tầm xa (ODLRO): Thể hiện tính chất siêu chảy không ma sát, được xác định khi ma trận mật độ đơn hạt ở khoảng cách vô hạn đạt giá trị hữu hạn, tương ứng mật độ siêu chảy $\rho_s > 0$.

4. Thế năng ghim tuần hoàn đóng vai trò gì trong việc ổn định pha siêu tinh thể trong mô hình Bose – Hubbard?

Trong mô hình Bose – Hubbard thông thường chỉ có tương tác tầm xa đồng nhất, pha siêu tinh thể chỉ ổn định ở vùng mật độ hạt $\rho > 1/2$ do khuyết tật giữa nút; vùng $\rho < 1/2$ bị hiện tượng phân tách pha phá hủy. Khi đặt thêm thế năng ghim tuần hoàn $\mu_i = \pm \varepsilon$, thế năng ngoài này đóng vai trò ghim cố định một phần cấu trúc hạt, triệt tiêu sự phân tách pha và cho phép các khuyết tật lỗ trống di chuyển tự do mà không phá vỡ trật tự tinh thể nền, giúp pha siêu tinh thể xuất hiện bền vững ở cả hai phía mật độ $\rho < 1/2$ và $\rho > 1/2$.

5. Thuật toán Worm trong phương pháp Monte Carlo lượng tử (QMC) khắc phục nhược điểm gì của các thuật toán lấy mẫu truyền thống?

Các thuật toán QMC truyền thống dựa trên không gian cấu hình vòng kín thường gặp hiện tượng "nghẽn tới hạn" (critical slowing down) và xác suất chấp nhận cực thấp khi mô phỏng tại nhiệt độ thấp hoặc gần điểm chuyển pha do các cấu hình đường thế giới bị khóa chặt. Thuật toán Worm mở rộng không gian trạng thái bằng cách cho phép các đường thế giới bị gián đoạn tại hai điểm (Worm head và Worm tail), tạo điều kiện lấy mẫu trực tiếp hàm Green đơn hạt và thực hiện các bước cập nhật cục bộ linh hoạt (dịch chuyển không gian, thời gian, sinh/hủy Worm), từ đó tính toán chính xác mật độ siêu chảy $\rho_s$ và cấu trúc pha mà không bị bẫy cấu hình.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Kim Oanh đã nghiên cứu toàn diện các đặc tính pha và chuyển pha dị thường trong các mô hình mạng spin Ising và mô hình hạt boson Bose – Hubbard hai chiều. Bằng việc kết hợp phương pháp giải tích trường trung bình cải biên, mô phỏng Monte Carlo cổ điển Metropolis và Monte Carlo lượng tử với thuật toán Worm, công trình đã giải thích thành công cơ chế vi mô của các bước nhảy từ hóa loại I trong perovskite Mangan $\text{Pr}{0.7}\text{Ca}{0.3}\text{MnO}_3$, làm rõ tính chất cao nguyên từ hóa phân số $1/3$ trên mạng vấp spin Shastry – Sutherland, đồng thời xác lập điều kiện thế năng ghim tuần hoàn để ổn định pha siêu tinh thể ở cả hai phía mật độ tinh thể. Các kết quả đạt được cung cấp hệ thống luận cứ lý thuyết và số liệu mô phỏng nhất quán, đóng góp dữ liệu khoa học phục vụ công tác nghiên cứu vật liệu từ nano và các hệ lượng tử trên mạng quang học.