Tổng quan luận án

Luận án tiến sĩ vật lý của tác giả Nguyễn Văn Hinh, hoàn thành năm 2019 tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, mang tiêu đề: "Phương pháp tiếp cận vi mô Ginzburg-Landau đối với sự đồng tồn tại pha trong hệ nhiều hạt". Công trình được thực hiện dưới sự hướng dẫn của TS. Nguyễn Trí Lân và GS. Nguyễn Toàn Thắng tại Bộ môn Vật lý lý thuyết – Vật lý toán, Viện Vật lý.

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

Kể từ khi lý thuyết vi mô Bardeen-Cooper-Schrieffer (BCS) ra đời năm 1957, quan niệm truyền thống trong vật lý chất rắn cho rằng tính siêu dẫn và tính sắt từ là hai trạng thái đối kháng, loại trừ lẫn nhau. Trong khuôn khổ lý thuyết BCS chuẩn, ngưng tụ siêu dẫn được thiết lập thông qua tương tác hút qua trung gian phonon liên kết các điện tử có spin đối song thành các cặp Cooper spin-singlet ($S = 0$). Khi đưa các nguyên tử tạp chất từ vào mạng tinh thể siêu dẫn, trường cục bộ của mô men từ tạp chất sẽ phá hủy các cặp singlet, dẫn đến sự sụt giảm nhanh chóng nhiệt độ chuyển pha siêu dẫn $T_{sc}$.

Mặc dù sự đồng tồn tại giữa siêu dẫn và trật tự phản sắt từ (AFM) đã được quan sát từ những năm 1980 (như trong các hợp chất fermion nặng chứa Cerium, nơi mô men từ trung bình bằng không), sự phát hiện tính siêu dẫn xuất hiện sâu bên trong pha sắt từ của hợp chất fermion nặng $UGe_2$ vào năm 2000 dưới áp suất cao (bởi Saxena và cộng sự), tiếp nối bởi các chất $UIr$, $URhGe$, $UCoGe$, đã đặt ra thách thức lý thuyết căn bản. Trong các hợp chất này, tính siêu dẫn xuất hiện ở nhiệt độ $T_{sc}$ thấp hơn hẳn nhiệt độ Curie $T_c$ mà không loại trừ trật tự sắt từ.

Khoảng trống nghiên cứu được tác giả xác định cụ thể như sau:

  • Cơ chế vi mô dẫn đến sự kết cặp siêu dẫn trong các hợp chất sắt từ vẫn chưa đạt được sự đồng thuận hoàn chỉnh: bản chất của cặp Cooper là spin-singlet hay spin-triplet; nguồn gốc tương tác hút do trao đổi phonon, magnon hay sóng mật độ; trật tự từ do các điện tử 5f định xứ hay điện tử linh động đảm nhận.
  • Các mô hình lý thuyết hiện tượng học trước đó (chẳng hạn mô hình Shopova, Dahl, Machida - Ohmi) chủ yếu dựa trên các lập luận đối xứng vĩ mô và dữ liệu nhiệt động lực học, thiếu sự suy dẫn hiển ngôn từ một Hamiltonian vi mô của hệ nhiều hạt.
  • Mối liên hệ giải tích giữa nhiệt độ chuyển pha, độ từ hóa tự phát và áp suất ngoài trong các giản đồ pha chưa được thiết lập một cách nhất quán từ lý thuyết trường lượng tử vi mô.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu bao trùm của luận án là hình thành một phiếm hàm năng lượng Ginzburg-Landau (GL) từ nền tảng vi mô để mô tả sự đồng tồn tại của nhiều pha trong hệ fermion tương quan mạnh. Các nhiệm vụ nghiên cứu được văn bản xác định gồm:

  1. Phát triển phương pháp biến đổi lý thuyết trường fermion thành lý thuyết trường hiệu dụng dựa trên các tham số trật tự thông qua phép biến đổi Hubbard-Stratonovich (HS) chia tách theo các kênh tương tác khả dĩ.
  2. Xây dựng biểu diễn giải tích chung của phiếm hàm Ginzburg-Landau cho hệ ba tham số trật tự bằng kỹ thuật tích phân phiếm hàm trên biến Grassmann và hàm Green.
  3. Rút ra biểu thức giải tích vi mô cho phiếm hàm Ginzburg-Landau hai tham số trật tự đại diện cho pha sắt từ và pha siêu dẫn spin-triplet trong các hợp chất fermion nặng gốc Uranium.
  4. Áp dụng phiếm hàm thiết lập được để nghiên cứu trạng thái cân bằng pha, thiết lập giản đồ pha nhiệt độ - áp suất ($T - P$) và khảo sát sự phụ thuộc tương hỗ giữa các tham số trật tự trong hợp chất $UGe_2$.

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

  • Đối tượng nghiên cứu: Trạng thái đồng tồn tại giữa trật tự sắt từ (FM) và trật tự siêu dẫn (SC), đặc biệt là trạng thái kết cặp spin-triplet ($S = 1$) trong hệ nhiều hạt fermion nặng; các mô hình Hamiltonian vi mô (mô hình Hubbard, mô hình Heisenberg, Hamiltonian tương tác kết cặp tổng quát).
  • Phạm vi nghiên cứu:
    • Về mặt lý thuyết: Tiếp cận bằng phương pháp lý thuyết trường giải tích vi mô, tích phân phiếm hàm trên đại số Grassmann, phương pháp hàm Green ở nhiệt độ hữu hạn, lý thuyết chuyển pha Ginzburg-Landau trong gần đúng trường trung bình.
    • Về hệ vật liệu thực tế: Tập trung khảo sát hệ fermion nặng $UGe_2$ trong vùng nhiệt độ thấp ($T < 60\text{ K}$) và dải áp suất thủy tĩnh từ $0$ đến khoảng $2\text{ GPa}$ tương ứng với các vùng chuyển pha lượng tử được công bố trong các nghiên cứu thực nghiệm.

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

Các hướng nghiên cứu thực nghiệm và lý thuyết tiền đề

Văn bản luận án điểm lại chi tiết lịch sử phát triển của lĩnh vực siêu dẫn và từ tính qua các mốc nghiên cứu then chốt:

  • Thực nghiệm siêu dẫn kinh điển và hiệu ứng đồng vị: Kamerlingh Onnes (1911) phát hiện tính dẫn điện hoàn hảo của thủy ngân ở $4,2\text{ K}$; Meissner và Ochsenfeld (1933) phát hiện tính nghịch từ hoàn hảo; Maxwell (1950) và Reynolds cùng cộng sự (1950) chứng minh thực nghiệm mối liên hệ hiệu ứng đồng vị $M^{1/2} T_c = \text{const}$, khẳng định vai trò của dao động mạng ion đối với cơ chế siêu dẫn.
  • Lý thuyết hiện tượng học và vi mô chuẩn: Heinz London và Fritz London (1937) thiết lập phương trình điện động lực học siêu dẫn và độ xuyên sâu $\lambda$; Ginzburg và Landau (1950) đề xuất lý thuyết hiện tượng học mô tả chuyển pha loại hai; Bardeen, Cooper và Schrieffer (1957) thiết lập lý thuyết vi mô BCS; Gor'kov (1959) chứng minh rằng lý thuyết Ginzburg-Landau có thể suy ra chính xác từ lý thuyết vi mô BCS với tham số trật tự $\Psi$ tỉ lệ với khe năng lượng $\Delta$; Abrikosov (1957) giải hệ phương trình GL và khám phá mạng xoáy trong siêu dẫn loại II mang lượng tử từ thông $\Phi_0 = hc/(2e)$.
  • Sự quan sát đồng tồn tại từ tính và siêu dẫn: Lynn cùng cộng sự (1981), Genicon cùng cộng sự, Rossat cùng cộng sự khảo sát tinh thể $HoMo_6S_8$, ghi nhận chuyển pha siêu dẫn tại $T = 1,82\text{ K}$ và chuyển pha từ tại $T = 0,67\text{ K}$, bão hòa ở $T = 0,4\text{ K}$. Frank Steglich và cộng sự (1979) phát hiện siêu dẫn trong kim loại fermion nặng $CeCu_2Si_2$ dưới $T_c \approx 0,5\text{ K}$. Jérome và cộng sự (1980) phát hiện siêu dẫn trong chất dẫn hữu cơ $(TMTSF)2PF_6$ tại $T_c = 1,1\text{ K}$. Kamihara và cộng sự (2008) khám phá siêu dẫn pnictide gốc sắt $LaFeAsO{1-x}F_x$ với $T_c \approx 26\text{ K}$.
  • Hệ siêu dẫn sắt từ gốc Uranium: Saxena cùng các cộng sự (2000) phát hiện tính siêu dẫn ngập sâu trong pha sắt từ của $UGe_2$ dưới áp suất; Akazawa và cộng sự (2004) phát hiện $UIr$; Aoki và cộng sự (2001) phát hiện $URhGe$; Huy và cộng sự (2007) phát hiện $UCoGe$.
  • Các nghiên cứu lý thuyết về siêu dẫn sắt từ:
    • Mô hình hiện tượng học: Shopova và cộng sự (2005) sử dụng phiếm hàm năng lượng tự do GL mở rộng có bổ sung số hạng ghép cặp tương tác dạng $f_I(\psi, M) = i\gamma_0 \vec{M} \cdot (\vec{\psi} \times \vec{\psi}^*)$ để nghiên cứu pha Meissner và các miền ổn định pha trong siêu dẫn spin-triplet.
    • Mô hình vi mô: Dahl và Sudbø đề xuất Hamiltonian ba thành phần gồm khí điện tử tự do, số hạng kết cặp BCS tổng quát và số hạng trao đổi sắt từ Heisenberg để rút ra phiếm hàm GL sóng p.
    • Hướng nghiên cứu tại Việt Nam: GS. Đỗ Trần Cát (2003, 2004) phân tích lý thuyết pha siêu dẫn và sắt từ trong $UGe_2$ dựa trên cấu trúc dải năng lượng kết ô (nesting); Nhóm Hệ tương quan mạnh tại Trung tâm Vật lý lý thuyết thuộc Viện Vật lý (từ năm 2006) khảo sát hai pha sắt từ FM1, FM2 và đề xuất cơ chế kết cặp siêu dẫn thông qua kích thích tách trường tinh thể của các mức 5f Uranium.

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

Mặc dù các mô hình hiện tượng học của Shopova và cộng sự mô tả định tính tương đối tốt một số tính chất nhiệt động, các hệ số trong phiếm hàm năng lượng hoàn toàn mang tính kinh nghiệm và không liên hệ được với các thông số vi mô của vật liệu. Mặt khác, các nghiên cứu vi mô trước đó thường xử lý riêng rẽ từng tương tác hoặc sử dụng các gần đúng làm triệt tiêu các kênh tương tác cạnh tranh. Luận án lựa chọn giải quyết bài toán: xây dựng một sơ đồ nhất quán sử dụng phép biến đổi tích phân phiếm hàm và biến đổi Hubbard-Stratonovich để khử các bậc tự do fermion, thu nhận phiếm hàm Ginzburg-Landau đa thành phần từ một Hamiltonian vi mô ban đầu, từ đó khảo sát một cách có hệ thống sự đồng tồn tại pha từ - siêu dẫn trong hợp chất $UGe_2$.

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

Các lý thuyết và mô hình nền tảng

  1. Lý thuyết Ginzburg-Landau và lý thuyết chuyển pha Landau: Tiếp cận hiện tượng chuyển pha thông qua việc phá vỡ đối xứng tự phát; xây dựng phiếm hàm năng lượng tự do $F = F_0(T) + F_L(T, \psi)$ dưới dạng khai triển đa thức giải tích theo tham số trật tự $\psi$. Tham số trật tự là đại lượng triệt tiêu ở pha mất trật tự đối xứng cao ($T > T_c$) và nhận giá trị khác không ở pha trật tự ($T < T_c$).
  2. Lý thuyết trật tự sắt từ:
    • Mô hình Heisenberg cho moment từ định xứ: Hamiltonian tương tác trao đổi $H = - \sum_{i \neq j} J_{ij} \hat{S}_i \hat{S}_j$, xuất phát từ tích phân trao đổi $J$ trong cơ học lượng tử (khảo sát hệ hai điện tử của nguyên tử Heli). Trạng thái spin-triplet đối xứng ($S = 1$) có mức năng lượng $E_2 = E^0 + K - J$, trạng thái spin-singlet phản đối xứng ($S = 0$) có năng lượng $E_1 = E^0 + K + J$. Điều kiện $J > 0$ ưu tiên sắp xếp spin song song (sắt từ).
    • Lý thuyết trường trung bình phân tử Weiss: Thay thế tương tác giữa các spin bằng từ trường hiệu dụng phân tử $H_W = 2 \sum_j J_{ij} \langle S_j \rangle$, dẫn đến liên hệ giữa nhiệt độ Curie và tích phân trao đổi: $k_B T_c = \frac{2S(S+1)}{3} \sum_j J_{ij} = \frac{2S(S+1)}{3} ZJ$.
    • Mô hình Hubbard cho hệ điện tử linh động: Xét Hamiltonian $H = -t \sum_{\langle j, l \rangle, \sigma} c_{j\sigma}^\dagger c_{l\sigma} + U \sum_j n_{j\uparrow} n_{j\downarrow} - \mu \sum_{j,\sigma} n_{j\sigma}$. Nghiên cứu các giới hạn giải được chính xác: giới hạn dải dẫn không tương tác ($U = 0$) với phổ tán sắc $\xi_k = -2t \sum_{\nu=1}^d \cos k_\nu - \mu$ và giới hạn nguyên tử định xứ ($t = 0$).
    • Tiêu chuẩn Stoner cho tính sắt từ dải dẫn: Xuất phát từ sự cạnh tranh giữa động năng gia tăng và thế năng tương tác Coulomb đẩy $U$, điều kiện làm xuất hiện độ từ hóa tự phát là $U N(E_F) > 1$, trong đó $N(E_F)$ là mật độ trạng thái tại mức Fermi.
  3. Lý thuyết BCS và BCS tổng quát:
    • Cơ chế hút qua trung gian phonon: Phân tích thời gian trễ của biến dạng mạng ion $\tau_D \sim 10^{-12}\text{ s}$ so với tốc độ chuyển động của điện tử $v_F \sim 10^8\text{ cm/s}$ (quãng đường dịch chuyển $\tau_D v_F \sim 10^{-4}\text{ cm}$ vượt xa tầm tác dụng chắn của lực đẩy Coulomb).
    • Trạng thái cơ bản BCS: $|\psi_{BCS}\rangle = \prod_k (u_k + v_k c_{k\uparrow}^\dagger c_{-k\downarrow}^\dagger)|0\rangle$ với chuẩn hóa $|u_k|^2 + |v_k|^2 = 1$.
    • Biến đổi Bogolyubov: Chéo hóa Hamiltonian thông qua phép biến đổi sang các toán tử chuẩn hạt mới $\gamma_{k0}, \gamma_{k1}$, xác định phổ năng lượng kích thích $E_k = \sqrt{\xi_k^2 + |\Delta_k|^2}$.
    • Lý thuyết BCS tổng quát ma trận $2 \times 2$: Biểu diễn ma trận hàm khe $\hat{\Delta}k = \begin{pmatrix} \Delta{k,\uparrow\uparrow} & \Delta_{k,\uparrow\downarrow} \ \Delta_{k,\downarrow\uparrow} & \Delta_{k,\downarrow\downarrow} \end{pmatrix}$. Tách thành phần quỹ đạo và spin: kênh spin-singlet đối xứng chẵn $\hat{\Delta}_k = i\sigma_y \Psi(k)$; kênh spin-triplet đối xứng lẻ biểu diễn qua véc tơ $\vec{d}(k)$ với $\hat{\Delta}_k = i (\vec{d}(k) \cdot \vec{\sigma}) \sigma_y$. Trạng thái unita thỏa mãn $\vec{q}(k) = i \vec{d}(k) \times \vec{d}^*(k) = 0$.

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

Luận án hoàn toàn sử dụng phương pháp tính toán giải tích vi mô của lý thuyết nhiều hạt:

  • Phương pháp hàm Green ở nhiệt độ hữu hạn: Thiết lập hệ phương trình chuyển động cho các hàm Green hạt đơn và hàm Green hai hạt để tính toán các giá trị trung bình nhiệt động và biểu diễn phương trình tự hợp cho các tham số trật tự.
  • Phương pháp tích phân phiếm hàm trên biến Grassmann: Biểu diễn hàm phân bố thống kê $Z = \text{Tr}, e^{-\beta(H - \mu N)}$ dưới dạng tích phân phiếm hàm đối với các trường Grassmann biến thiên theo thời gian ảo $\tau$.
  • Phép biến đổi Hubbard-Stratonovich (HS): Đưa các trường phụ liên tục (bosonic auxiliary fields) vào tích phân phiếm hàm để tách cặp số hạng tương tác bốn fermion thành dạng bậc hai theo các trường fermion. Luận án chia tách theo ba kênh vật lý:
    • Kênh mật độ: sinh ra tham số trật tự sóng mật độ điện tích (CDW).
    • Kênh kết cặp (Cooper): sinh ra tham số trật tự siêu dẫn (SC).
    • Kênh trao đổi: sinh ra tham số trật tự từ tính/sóng mật độ spin (SDW/FM).
  • Phương pháp lấy tích phân phiếm hàm Gauss và khai triển vết logarithm: Sau khi biến đổi HS, tích phân trên các biến Grassmann thuần nhất dạng Gauss được lấy chính xác, dẫn đến định thức fermion $\text{Det}[\hat{G}^{-1}]$. Khai triển $\text{Tr}\ln[\hat{G}^{-1}]$ theo các trường phụ gần điểm chuyển pha cho phép thu nhận dạng giải tích tường minh của phiếm hàm năng lượng Ginzburg-Landau.

Nguồn tư liệu và dữ liệu thực nghiệm đối chiếu

Luận án sử dụng các số liệu thực nghiệm đã công bố từ các nhóm nghiên cứu thực nghiệm quốc tế đối với hệ fermion nặng:

  • Giản đồ pha $T - P$ và các số đo độ từ hóa của hợp chất $UGe_2$ từ công trình của Saxena và cộng sự (2000), Huxley và cộng sự (2001), Pfleiderer và cộng sự (2002), Tateiwa và cộng sự (2001).
  • Giản đồ pha áp suất - nhiệt độ của hợp chất $CeRhIn_5$ từ các phép đo nhiệt dung riêng và tán xạ neutron của Hegger và cộng sự (2000), Knebel và cộng sự (2006), Park và cộng sự (2006).
  • Dữ liệu điện trở và nhiễu xạ neutron trên tinh thể $HoMo_6S_8$ của Lynn và cộng sự (1981), Rossat và cộng sự (1985).

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa toàn bộ bức tranh thực nghiệm và nền tảng lý thuyết liên quan đến từ tính, siêu dẫn và chuyển pha trong hệ tương quan mạnh.

Tác giả phân tích tính bất ổn định của chất lỏng Fermi trong kim loại dẫn tới hai xu hướng chính: sự hình thành cặp điện tử (siêu dẫn) và sóng mật độ (điện tích hoặc spin). Trong các kim loại fermion nặng chứa ion actinide và đất hiếm, lớp vỏ f lấp đầy không hoàn toàn (trạng thái $4f$ trong Ce hoặc $5f$ trong U) tương tác với dải dẫn tạo ra các chuẩn hạt nặng với khối lượng hiệu dụng $m^*$ lớn hơn khối lượng điện tử tự do $m$ hàng trăm đến hàng nghìn lần, dẫn đến độ rộng dải chuẩn hạt hẹp cỡ meV, đi kèm sự tăng vọt của nhiệt dung riêng và độ từ cảm Pauli.

Văn bản tổng hợp các tham số thực nghiệm đo đạc được trên các hợp chất tiêu biểu như sau:

Hợp chất Kiểu trật tự từ tại $P = 0$ Nhiệt độ chuyển pha từ ($P=0$) Đặc trưng pha siêu dẫn ($T_{sc}$) Vùng áp suất đồng tồn tại
$UGe_2$ Sắt từ (FM), trục dễ là trục $a$, $\mu_s = 1,4\mu_B$ $T_c = 52\text{ K}$ $T_{sc}^{\max} \approx 0,8\text{ K}$ (tại $P \approx 1,2\text{ GPa}$) $1,0\text{ GPa} \le P \le 1,6 - 1,7\text{ GPa}$ (nằm trọn trong pha FM)
$CeRhIn_5$ Phản sắt từ (AFM), mô men so le $\approx 0,5\mu_B$ $T_N = 3,8\text{ K}$ $T_c = T_N$ tại $P_c$; siêu dẫn không truyền thống đạt $T_{sc}$ cao khi $P > P_c^*$ $1,6\text{ GPa} \le P \le 1,95\text{ GPa}$ ($P_c^* = 1,95\text{ GPa}$)
$HoMo_6S_8$ Sắt từ (FM) $T_m = 0,67\text{ K}$ $T_{sc} = 1,82\text{ K}$ Tồn tại trạng thái tái nhập dẫn thường tại nhiệt độ thấp

Về mặt lý thuyết, chương này khảo sát chi tiết:

  • Nguồn gốc lượng tử của lực hút ghép cặp: Phân tích cơ chế Fröhlich trao đổi phonon và lý do lực hút gián tiếp chậm trễ vượt qua lực đẩy Coulomb tức thời.
  • Cơ chế thăng giáng spin trong ghép cặp triplet: Tại lân cận điểm tới hạn lượng tử sắt từ, các thăng giáng spin định hướng song song đóng vai trò môi giới ghép cặp các điện tử thành cặp Cooper spin-triplet với spin toàn phần $S = 1$ gồm ba trạng thái: $|\uparrow\uparrow\rangle$ ($L=1, S_z=1$), $|\downarrow\downarrow\rangle$ ($L=1, S_z=-1$), và $\frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle + |\downarrow\uparrow\rangle)$ ($L=1, S_z=0$).
  • Mô hình hiện tượng học Shopova cho siêu dẫn sắt từ với phiếm hàm: $$f_{GL}(\psi, M) = f_S(\psi) + f_F(M) + f_I(\psi, M) + \frac{1}{2}H^2$$ trong đó số hạng tương tác trực tiếp là $f_I(\psi, M) = i\gamma_0 \vec{M} \cdot (\vec{\psi} \times \vec{\psi}^*)$.

Chương 2: Một số phương pháp tiếp cận vi mô Ginzburg-Landau

Chương 2 trình bày các công cụ toán - lý vi mô để chuyển đổi từ Hamiltonian nhiều hạt sang phiếm hàm hiệu dụng.

  • Mục 2.1 - Phương pháp hàm Green: Thiết lập phương trình chuyển động cho hàm Green Matsubara hạt đơn. Tác giả áp dụng kỹ thuật hàm Green để giải bài toán trường trung bình cho hệ sắt từ Heisenberg và hệ siêu dẫn BCS, dẫn xuất hệ phương trình tự hợp xác định độ từ hóa và khe năng lượng $\Delta_k$.
  • Mục 2.2 - Phương pháp tích phân phiếm hàm và biến Grassmann:
    • Trình bày giải tích biến Grassmann: đại số phản giao hoán ${c_i, c_j} = 0$, tích phân Berezin trên biến Grassmann: $\int d\eta = 0$, $\int \eta d\eta = 1$.
    • Biểu diễn hàm phân bố thống kê của hệ fermion: $$Z = \int \mathcal{D}[\psi^\dagger, \psi] \exp\left{ - \int_0^\beta d\tau \left[ \sum_{k,\sigma} \psi_{k\sigma}^\dagger(\tau) \partial_\tau \psi_{k\sigma}(\tau) + H(\psi^\dagger, \psi) \right] \right}$$
    • Áp dụng phép biến đổi Hubbard-Stratonovich vào bài toán siêu dẫn BCS: Tương tác bốn fermion dạng $-V \psi_\uparrow^\dagger \psi_\downarrow^\dagger \psi_\downarrow \psi_\uparrow$ được tách cặp bằng cách tích hợp với trường phụ phức $\Delta(\tau), \Delta^*(\tau)$.
    • Tích phân Gauss trên các biến Grassmann cho ra định thức fermion. Khai triển chuỗi Taylor đối với số hạng logarithm $\text{Tr}\ln \hat{G}^{-1}[\Delta]$ theo bậc của tham số trật tự $\Delta$ xung quanh nhiệt độ chuyển pha $T_c$.
    • Thu nhận phiếm hàm năng lượng Ginzburg-Landau một thành phần chuẩn: $$F_{GL}[\Delta] = \int d^3r \left[ a(T) |\Delta|^2 + \frac{b}{2} |\Delta|^4 + c |\nabla \Delta|^2 \right]$$ với các hệ số $a(T) = \alpha (T - T_c)$, $b$, $c$ được biểu diễn giải tích hoàn toàn qua mật độ trạng thái tại mặt Fermi $N(0)$ và vận tốc Fermi $v_F$.

Chương 3: Thiết lập phiếm hàm năng lượng Ginzburg-Landau nhiều thành phần

Chương 3 là phần đóng góp lý thuyết mở rộng, tập trung vào việc mô hình hóa hệ nhiều tham số trật tự cạnh tranh.

  • Mô hình ba tham số trật tự: Xây dựng Hamiltonian tổng quát mô tả tương tác hai hạt trong hệ fermion nhiều hạt gồm động năng và các số hạng tán xạ tổng quát.
  • Biến đổi Hubbard-Stratonovich đa kênh: Tách tương tác hạt thành ba kênh độc lập:
    1. Kênh mật độ (Cooper-density/CDW): Đưa vào trường phụ vô hướng mô tả thăng giáng mật độ điện tích.
    2. Kênh kết cặp (Cooper/SC): Đưa vào trường phụ tensor/véc tơ mô tả sự ngưng tụ cặp Cooper.
    3. Kênh trao đổi (Exchange/FM/SDW): Đưa vào trường phụ véc tơ mô tả thăng giáng spin và độ từ hóa.
  • Khử tham số kết cặp và lấy tích phân trường fermion: Xây dựng ma trận nghịch đảo hàm Green trong không gian mở rộng (spin $\otimes$ Nambu). Thực hiện phép lấy vết trên các biến fermion để nhận được phiếm hàm tác dụng hiệu dụng $S_{eff}$ chỉ phụ thuộc vào ba trường tham số trật tự.
  • Phiếm hàm Ginzburg-Landau ba thành phần: Khai triển phiếm hàm hiệu dụng tới bậc bốn theo các trường và bậc hai theo gradient không gian. Tác giả thu được dạng tường minh của các số hạng tự tương tác của từng trường và các số hạng ghép cặp tương tác chéo giữa tham số trật tự siêu dẫn, từ tính và sóng mật độ điện tích.
  • Khảo sát các mô hình rút gọn: Phân tích các trường hợp giới hạn khi đóng băng một hoặc hai kênh tương tác, tập trung vào kênh trao đổi spin để phục vụ khảo sát tương tác từ tính trong chương tiếp theo.

Chương 4: Sự đồng tồn tại các trật tự sắt từ và siêu dẫn trong các hợp chất fermion nặng

Chương 4 áp dụng hình thức luận đã phát triển vào bài toán cụ thể của hợp chất fermion nặng gốc Uranium ($UGe_2$).

  • Nguồn gốc vi mô của phiếm hàm GL hai thành phần: Xuất phát từ Hamiltonian hạt kết hợp tương tác dải dẫn và tương tác trao đổi spin, thực hiện biến đổi HS giữ lại hai trường: trường véc tơ phức $\vec{\psi}$ (tham số trật tự siêu dẫn spin-triplet) và trường véc tơ thực $\vec{M}$ (tham số trật tự từ hóa sắt từ).
  • Thiết lập biểu thức giải tích phiếm hàm: Chứng minh từ cấp độ vi mô sự xuất hiện của số hạng ghép cặp tương tác bất đối xứng quay: $$f_{GL}(\vec{\psi}, \vec{M}) = a_s |\vec{\psi}|^2 + \frac{b_s}{2} |\vec{\psi}|^4 + a_f |\vec{M}|^2 + \frac{b_f}{2} |\vec{M}|^4 + i \gamma_0 \vec{M} \cdot (\vec{\psi} \times \vec{\psi}^*)$$ trong đó hệ số ghép cặp $\gamma_0$ bắt nguồn từ tích phân ba hàm Green chứa đỉnh tương tác spin - quỹ đạo và cấu trúc dải điện tử.
  • Khảo sát giản đồ pha và điều kiện đồng tồn tại trong $UGe_2$:
    • Cực tiểu hóa phiếm hàm $f_{GL}$ đối với các thành phần của $\vec{\psi}$ và $\vec{M}$.
    • Thiết lập giản đồ pha trong mặt phẳng quy dẫn $(t, r)$ và mặt phẳng thực nghiệm nhiệt độ - áp suất ($T - P$).
    • Kết quả chỉ ra rằng khi tham số ghép cặp $\gamma_0 \neq 0$, trạng thái siêu dẫn triplet có thể xuất hiện và ổn định nhiệt động lực học trong một miền nhiệt độ xác định bên dưới đường nhiệt độ Curie $T_c(P)$.
    • Phân tích sự chuyển đổi bản chất chuyển pha: Chuyển pha từ pha sắt từ FM1 sang pha thuận từ PM đổi từ bậc hai sang bậc một tại điểm ba $T_{cr}$ trên giản đồ $T(P)$. Vùng pha đồng tồn tại sắt từ - siêu dẫn (FS) được định vị chính xác trong miền lân cận áp suất $P \sim 1,2\text{ GPa}$, giải thích sự sụt giảm đồng thời của nhiệt độ Curie và nhiệt độ siêu dẫn quanh áp suất tới hạn $P_c \approx 1,6 - 1,7\text{ GPa}$.

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

Đóng góp mới về mặt lý luận

  1. Dẫn xuất vi mô phiếm hàm Ginzburg-Landau đa thành phần: Luận án đã thiết lập thành công một sơ đồ giải tích xuất phát từ Hamiltonian vi mô của hệ fermion tương tác, áp dụng tích phân phiếm hàm trên biến Grassmann và phép biến đổi Hubbard-Stratonovich đa kênh để thu nhận phiếm hàm Ginzburg-Landau nhiều thành phần. Phương pháp này khắc phục tính kinh nghiệm thuần túy của các lý thuyết hiện tượng học trước đó.
  2. Làm sáng tỏ nguồn gốc vi mô của số hạng tương tác từ - siêu dẫn: Luận án chứng minh một cách giải tích rằng số hạng ghép cặp tương tác $i\gamma_0 \vec{M} \cdot (\vec{\psi} \times \vec{\psi}^*)$ giữa trật tự sắt từ và siêu dẫn spin-triplet có thể được rút ra trực tiếp từ lý thuyết vi mô nhiều hạt khi tính đến các tương tác trao đổi spin, xác lập bản chất vi mô cho hệ số tương tác $\gamma_0$.

Đóng góp mới về mặt thực tiễn và tính toán giản đồ pha

  1. Giải thích giản đồ pha của hợp chất $UGe_2$: Áp dụng phiếm hàm GL hai thành phần đã thiết lập để tính toán và tái tạo định tính và định lượng giản đồ pha nhiệt độ - áp suất ($T - P$) của $UGe_2$. Luận án làm rõ điều kiện nhiệt động để pha siêu dẫn triplet tồn tại ổn định hoàn toàn bên trong miền sắt từ.
  2. Làm rõ sự biến thiên tham số trật tự theo áp suất: Mô hình giải thích được cơ chế chuyển đổi bậc chuyển pha (từ bậc hai sang bậc một tại điểm ba tới hạn) và sự biến mất đồng thời của trật tự từ tính và trật tự siêu dẫn tại lân cận áp suất tới hạn $P_c \sim 1,6 - 1,7\text{ GPa}$.

Kiến nghị và định hướng ứng dụng

  • Sử dụng khung lý thuyết vi mô Ginzburg-Landau đa thành phần đã xây dựng để tiếp tục khảo sát các hệ vật liệu fermion nặng sắt từ khác đã phát hiện như $UIr$, $URhGe$, $UCoGe$.
  • Vận dụng phương pháp phân tách kênh Hubbard-Stratonovich để nghiên cứu sự cạnh tranh giữa các sóng mật độ điện tích (CDW), sóng mật độ spin (SDW) và siêu dẫn trong các họ vật liệu tương quan mạnh mới như pnictide sắt và chalcogenide.

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

Văn bản luận án cho thấy một số giả thiết giới hạn trong quá trình tính toán giải tích:

  • Giới hạn cặp unita: Trong quá trình chéo hóa Hamiltonian Bogolyubov - de Gennes ở Chương 1 và Chương 4, tác giả giới hạn khảo sát ở trường hợp kết cặp unita, tức điều kiện vector $\vec{q}(k) = i \vec{d}(k) \times \vec{d}^*(k) = 0$ (đảm bảo năng lượng kích thích chuẩn hạt không bị phân cực thêm $\vec{E}{k\uparrow} = \vec{E}{k\downarrow}$). Trạng thái kết cặp không unita phức tạp hơn liên quan đến sự phá vỡ đối xứng đảo thời gian nội tại chưa được giải tích toàn diện.
  • Giả thiết đẳng hướng: Trong một số bước rút gọn phiếm hàm GL ở Chương 4, các ảnh hưởng của tính bất đẳng hướng tinh thể trực giao sâu sắc của mạng tinh thể $UGe_2$ được giản lược hóa để thu nhận các nghiệm giải tích đối xứng.
  • Xử lý thăng giáng trong gần đúng trường trung bình: Phiếm hàm Ginzburg-Landau được thiết lập bằng cách khai triển quanh nghiệm trường tĩnh trung bình; các hiệu ứng thăng giáng lượng tử phi Gauss ở sát điểm tới hạn lượng tử ($T \to 0$, $P \to P_c$) chưa được tích hợp bằng phương pháp nhóm tái chuẩn hóa (Renormalization Group).

Giá trị tham khảo

Bản luận án có giá trị tham khảo chuyên môn trong các lĩnh vực cụ thể:

  • Đối với nghiên cứu sinh và học viên cao học ngành Vật lý lý thuyết và Vật lý toán: Chương 2 và Chương 3 cùng ba Phụ lục (Phụ lục A: Tính hệ fermion hiệu dụng bậc hai kết cặp với các trường phụ; Phụ lục B: Khai triển hàm mũ chứa các đạo hàm; Phụ lục C: Các phép đạo hàm lên hàm mũ) là tài liệu tham khảo chi tiết về kỹ thuật tích phân phiếm hàm trên biến Grassmann và phép biến đổi Hubbard-Stratonovich trong lý thuyết trường vật chất ngưng tụ.
  • Đối với giảng viên và nhà nghiên cứu vật lý chất rắn / hệ tương quan mạnh: Cung cấp khung phương pháp luận vi mô để chuyển đổi từ một mô hình Hamiltonian nhiều hạt sang phiếm hàm nhiệt động Ginzburg-Landau nhiều tham số trật tự, có thể tái sử dụng cho các hệ có pha cạnh tranh phức tạp.
  • Đối với các nhóm nghiên cứu thực nghiệm: Cung cấp cơ sở giải thích định tính và định lượng cho các giản đồ pha đo đạc dưới điều kiện nhiệt độ cực thấp và áp suất cao trên các hợp chất liên kim loại gốc actinide và đất hiếm.

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

1. Tại sao theo lý thuyết BCS truyền thống tính siêu dẫn và sắt từ lại loại trừ lẫn nhau?

Theo lý thuyết vi mô BCS (1957), trạng thái siêu dẫn truyền thống được tạo thành từ các cặp Cooper spin-singlet gồm hai điện tử có động lượng và spin đối song ($k\uparrow, -k\downarrow$), với tổng spin $S = 0$. Trong chất sắt từ, tương tác trao đổi nội tại định hướng các spin điện tử song song cùng chiều để tạo ra độ từ hóa tự phát. Trường từ phân tử mạnh mẽ này tác động làm lệch mức năng lượng của các điện tử có spin ngược chiều (hiệu ứng Zeeman), bẻ gãy liên kết của cặp Cooper singlet và phá hủy trạng thái siêu dẫn. Do đó, siêu dẫn singlet và sắt từ không thể đồng thời tồn tại trong cùng một vùng không gian điện tử.

2. Trạng thái kết cặp siêu dẫn spin-triplet trong $UGe_2$ có cấu trúc như thế nào để cùng tồn tại được với sắt từ?

Trong hợp chất $UGe_2$, các điện tử 5f của Uranium chịu trách nhiệm cho cả từ tính và tính siêu dẫn. Do có mặt trật tự sắt từ mạnh, các điện tử kết cặp thành các cặp Cooper spin-triplet với spin toàn phần $S = 1$, trong đó các spin thành phần định hướng song song: $|\uparrow\uparrow\rangle$ ứng với hình chiếu spin $S_z = 1$, và $|\downarrow\downarrow\rangle$ ứng với $S_z = -1$. Vì các điện tử trong cặp Cooper triplet có spin song song, chúng không bị bẻ gãy bởi trường trao đổi sắt từ nội tại mà trái lại, các thăng giáng spin sắt từ đóng vai trò chính là lực hút hiệu dụng liên kết các điện tử thành cặp.

3. Phép biến đổi Hubbard-Stratonovich đóng vai trò gì trong việc thiết lập phiếm hàm Ginzburg-Landau của luận án?

Trong lý thuyết trường lượng tử nhiều hạt, số hạng tương tác điện tử - điện tử là tương tác bốn toán tử fermion, khiến việc tính toán tích phân phiếm hàm không thể thực hiện chính xác. Phép biến đổi Hubbard-Stratonovich là một đồng nhất thức tích phân Gauss cho phép chuyển hóa tương tác bốn fermion thành tương tác bậc hai của fermion liên kết với các trường phụ boson (auxiliary fields). Sau khi thực hiện biến đổi này, các bậc tự do fermion được tích phân chính xác (lấy vết định thức fermion), để lại một phiếm hàm tác dụng hiệu dụng chỉ phụ thuộc vào các trường phụ boson. Các trường phụ này ở giới hạn tĩnh chính là các tham số trật tự của lý thuyết Ginzburg-Landau.

4. Dạng toán học của số hạng tương tác giữa trật tự sắt từ và siêu dẫn triplet trong phiếm hàm Ginzburg-Landau được biểu diễn như thế nào?

Trong phiếm hàm năng lượng tự do Ginzburg-Landau của hệ siêu dẫn sắt từ triplet, tương tác trực tiếp giữa tham số trật tự véc tơ từ hóa thực $\vec{M}$ và tham số trật tự siêu dẫn véc tơ phức $\vec{\psi}$ được biểu diễn thông qua số hạng: $$f_I(\vec{\psi}, \vec{M}) = i\gamma_0 \vec{M} \cdot (\vec{\psi} \times \vec{\psi}^)$$ Số hạng này mô tả sự ghép cặp bất đối xứng quay giữa độ từ hóa tự phát với mô men từ riêng của cặp Cooper spin-triplet (tỉ lệ với biểu thức $i\vec{\psi} \times \vec{\psi}^$). Luận án chứng minh số hạng hiện tượng học này bắt nguồn trực tiếp từ tích phân vết của ba hàm Green fermion trong mô hình vi mô.

5. Giản đồ pha thực nghiệm của hợp chất $UGe_2$ dưới áp suất có những đặc điểm nổi bật nào được luận án mô hình hóa?

Tại áp suất thường ($P = 0$), $UGe_2$ là chất sắt từ linh động với nhiệt độ Curie $T_c = 52\text{ K}$, độ từ hóa tự phát $\mu_s = 1,4\mu_B$. Khi tăng áp suất:

  • Hệ trải qua chuyển pha giữa hai pha sắt từ có độ phân cực khác nhau: pha phân cực mạnh FM2 ($\mu \approx 1,5\mu_B$) chuyển sang pha phân cực yếu FM1 ($\mu \approx 0,9\mu_B$) tại đường ranh giới $T_x(P)$.
  • Pha siêu dẫn (FS) chỉ xuất hiện trong miền nhiệt độ rất thấp ($T_{sc} \le 0,8\text{ K}$) và nằm gọn hoàn toàn bên trong miền sắt từ ở dải áp suất trung gian $1,0\text{ GPa} \le P \le 1,6 - 1,7\text{ GPa}$.
  • Đỉnh nhiệt độ siêu dẫn cao nhất đạt được tại lân cận áp suất mà đường chuyển pha $T_x$ tiến tới triệt tiêu ($T_x \to T_{sc}$). Cả trật tự sắt từ và siêu dẫn cùng biến mất tại áp suất tới hạn $P_c \approx 1,6 - 1,7\text{ GPa}$.

Kết luận

Luận án của tác giả Nguyễn Văn Hinh đã giải quyết một vấn đề căn bản của vật lý lý thuyết chất ngưng tụ: xây dựng thành công phương pháp tiếp cận vi mô để thiết lập phiếm hàm Ginzburg-Landau đa thành phần mô tả sự đồng tồn tại pha trong hệ nhiều hạt tương quan mạnh. Thông qua việc kết hợp kỹ thuật tích phân phiếm hàm Grassmann, phép biến đổi Hubbard-Stratonovich đa kênh và phương pháp hàm Green, công trình đã xác lập cơ sở vi mô cho số hạng tương tác giữa trật tự sắt từ và siêu dẫn spin-triplet. Kết quả tính toán giải tích đã giải thích thành công cơ chế cân bằng pha và cấu trúc giản đồ pha nhiệt độ - áp suất của hợp chất fermion nặng $UGe_2$, đóng góp phương pháp luận có giá trị cho việc nghiên cứu các trạng thái vật chất lượng tử đồng tồn tại.