Tổng quan luận án

Luận án tiến sĩ 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" do nghiên cứu sinh Nguyễn Văn Hinh thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Toàn Thắng và TS. Nguyễn Trí Lân tại Học viện Khoa học và Công nghệ cùng bộ môn Vật lý lý thuyết – Vật lý toán thuộc Viện Vật lý (Viện Hàn lâm Khoa học và Công nghệ Việt Nam), hoàn thành vào tháng 9 năm 2019. Luận án thuộc chuyên ngành Vật lý lý thuyết và Vật lý toán.

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

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

Tuy nhiên, thực nghiệm bắt đầu ghi nhận sự đồng tồn tại của siêu dẫn với trật tự phản sắt từ (năm 1980) và đặc biệt là sự phát hiện ra tính siêu dẫn nằm ngập sâu trong pha sắt từ ($T_{sc} < T_c$) ở hợp chất fermion nặng $\text{UGe}_2$ vào năm 2000, tiếp theo là $\text{UIr}$, $\text{URhGe}$ và $\text{UCoGe}$. Các hiện tượng này cho thấy cơ chế kết cặp truyền thống không thể giải thích thỏa đáng sự xuất hiện của siêu dẫn trong môi trường sắt từ có các mô men từ linh động mạnh. Dù đã có các mô hình hiện tượng học dựa trên lý thuyết Ginzburg-Landau và đối xứng vĩ mô để mô tả sự cùng tồn tại của hai pha, bản chất vi mô của tương tác ghép cặp, cấu trúc của tham số trật tự siêu dẫn (spin-singlet hay spin-triplet), nguồn gốc của từ tính (từ điện tử định xứ hay điện tử linh động), cũng như sự phụ thuộc của nhiệt độ chuyển pha và độ từ hóa tự phát vào áp suất vẫn chưa có lời giải đáp lý thuyết trọn vẹn. Đây là khoảng trống nghiên cứu vật lý lý thuyết mà luận án đặt ra giải quyết.

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

Mục tiêu nghiên cứu của luận án được xác định gồm ba nội dung chính:

  1. Thực hiện phép biến đổi lý thuyết trường fermion thành lý thuyết hiệu dụng dựa trên các trường kết cặp biểu diễn dưới dạng các tham số trật tự theo các kênh khả dĩ; thiết lập phiếm hàm năng lượng tự do Ginzburg-Landau (GL) từ phương pháp tiếp cận vi mô thông qua phép biến đổi Hubbard-Stratonovich (HS) và hình thức luận tích phân phiếm hàm kết hợp kỹ thuật hàm Green.
  2. Thiết lập biểu diễn chung của phiếm hàm Ginzburg-Landau nhiều thành phần cho hệ ba tham số trật tự (kênh mật độ, kênh kết cặp Cooper và kênh trao đổi), từ đó rút gọn về các mô hình cho từng kênh riêng biệt.
  3. Rút ra biểu diễn chính thức cho phiếm hàm Ginzburg-Landau hai thành phần (mô tả trật tự sắt từ và trật tự siêu dẫn spin-triplet) áp dụng cho hệ fermion nặng dựa trên Uranium ($\text{UGe}_2$); từ đó khảo sát sự đồng tồn tại của hai pha, vẽ giản đồ pha và phân tích các miền pha cân bằng theo các thông số nhiệt độ và áp suất.

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

  • Đối tượng nghiên cứu: Sự đồng tồn tại của các trật tự từ (sắt từ, phản sắt từ) và trật tự siêu dẫn phi truyền thống (đặc biệt là siêu dẫn spin-triplet) trong các hệ điện tử tương quan mạnh, trọng tâm là các hợp chất fermion nặng.
  • Phạm vi nghiên cứu:
    • Về lý thuyết: Tập trung vào phương pháp tiếp cận vi mô dẫn xuất phiếm hàm năng lượng Ginzburg-Landau từ mô hình vi mô nhiều hạt (Hamiltonian dạng Hubbard và Heisenberg mở rộng), áp dụng gần đúng trường trung bình, hình thức luận tích phân phiếm hàm trên đại số Grassmann và biến đổi Hubbard-Stratonovich.
    • Về đối tượng vật liệu: Khảo sát trọng tâm hợp chất sắt từ siêu dẫn $\text{UGe}_2$, kết hợp đối chiếu với các dữ liệu thực nghiệm trên các hợp chất fermion nặng khác như $\text{CeRhIn}_5$ và hợp chất từ $\text{HoMo}_6\text{S}_8$.
    • Về điều kiện nhiệt động: Dải nhiệt độ thấp (từ 0 K đến 52 K) và dải áp suất ngoài từ 0 GPa đến trên 2 GPa, bao quanh các điểm tới hạn lượng tử của vật liệu.

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

Các hướng nghiên cứu trước đây được luận án điểm lại

Luận án hệ thống hóa các mốc nghiên cứu lý thuyết và thực nghiệm then chốt về sự tương tác giữa trật tự từ và siêu dẫn:

  • Các nghiên cứu thực nghiệm đầu tiên: Quan sát điện trở bằng không trong trạng thái sắt từ của $\text{HoMo}_6\text{S}_8$ được công bố bởi Lynn và cộng sự (1980), tiếp theo bởi Genicon và cộng sự (1984). Thí nghiệm nhiễu xạ neutron trên đơn tinh thể $\text{HoMo}_6\text{S}_8$ của Rossat và cộng sự khẳng định khi làm lạnh chậm dưới $T < 0{,}60\text{ K}$, pha sắt từ xuất hiện, cường độ từ trường bão hòa quanh $T = 0{,}4\text{ K}$; vật liệu có chuyển pha siêu dẫn ở $T = 1{,}82\text{ K}$ và chuyển pha từ ở $T = 0{,}67\text{ K}$.
  • Khám phá siêu dẫn sắt từ: Phát hiện của Saxena và cộng sự (2000) về tính siêu dẫn ngập sâu trong pha sắt từ ($T_{sc} < T_c$) ở hợp chất $\text{UGe}_2$ dưới áp suất; tiếp sau là các phát hiện trên $\text{UIr}$ (Akazawa và cộng sự, 2004), $\text{URhGe}$ (Aoki và cộng sự, 2001) và $\text{UCoGe}$ (Huy và cộng sự, 2007). Các vật liệu này có điểm chung là trật tự sắt từ do mô men từ của các điện tử lớp $5f$ của Uranium chi phối, có đặc tính linh động mạnh và tính siêu dẫn xuất hiện kề cận điểm không ổn định từ tính.
  • Cơ chế kết cặp spin-triplet và thăng giáng spin: Các công trình lý thuyết của Fay và Appel (1980), Kirkpatrick và cộng sự (2001) đề xuất mô hình thăng giáng spin: gần điểm tới hạn lượng tử sắt từ, các thăng giáng từ đóng vai trò trung gian ghép cặp điện tử thành các cặp Cooper spin-triplet với spin toàn phần $S = 1$, bao gồm ba trạng thái $|{\uparrow\uparrow}\rangle$ ($L=1, S_z=1$), $|{\downarrow\downarrow}\rangle$ ($L=1, S_z=-1$), và $(|{\uparrow\downarrow}\rangle + |{\downarrow\uparrow}\rangle)/\sqrt{2}$ ($L=1, S_z=0$). Suhl (2001) đề xuất vai trò kết hợp của thăng giáng spin ngang và dọc. Ngược lại, Abrikosov (2001) và Mineev cùng cộng sự (2001) cho rằng siêu dẫn sóng $s$ có thể xuất hiện từ tương tác điện tử qua trung gian các mô men định xứ sắt từ.
  • Mô hình hiện tượng học Ginzburg-Landau: Mô hình của Machida và Ohmi (2001), Shopova và cộng sự (2003, 2005) xây dựng phiếm hàm năng lượng tự do hiện tượng luận cho siêu dẫn sắt từ spin-triplet dựa trên lý thuyết đối xứng, kết hợp số hạng siêu dẫn $f_S(\boldsymbol{\psi})$, số hạng sắt từ $f_F(\mathbf{M})$ và số hạng tương tác $f_I(\boldsymbol{\psi}, \mathbf{M}) = i\gamma_0 \mathbf{M} \cdot (\boldsymbol{\psi} \times \boldsymbol{\psi}^*)$.
  • Mô hình vi mô: Dahl và Sudbø (2007) khảo sát Hamiltonian vi mô gồm động năng điện tử tự do, số hạng kết cặp BCS spin tổng quát và số hạng trao đổi sắt từ Heisenberg, từ đó rút ra phiếm hàm Ginzburg-Landau cho siêu dẫn sắt từ sóng $p$.
  • Các nghiên cứu tại Việt Nam: GS. Đỗ Trần Cát (2003, 2004) nghiên cứu lý thuyết pha siêu dẫn và pha sắt từ trong $\text{UGe}_2$ dựa trên cấu trúc vùng năng lượng; nhóm nghiên cứu Hệ tương quan mạnh tại Viện Vật lý từ năm 2006 khảo sát hai pha từ FM1, FM2 trong $\text{UGe}_2$ và đề xuất cơ chế siêu dẫn qua trao đổi kích thích từ do tách trường tinh thể của các mức Uranium định xứ.

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

Các mô hình hiện tượng học của Machida, Ohmi hay Shopova tuy mô tả thành công các miền ổn định của pha nhưng được xây dựng thuần túy trên lập luận đối xứng nhiệt động lực học, không liên kết trực tiếp với các thông số vi mô cơ bản. Trong khi đó, các mô hình vi mô hiện có thường tách biệt từng cơ chế (chỉ xét trao đổi magnon, sóng mật độ spin, hoặc sóng mật độ điện tích) và chưa giải thích được sự phụ thuộc đồng thời của nhiệt độ chuyển pha và độ từ hóa vào áp suất ngoài. Luận án chọn khoảng trống: xây dựng một phương pháp tiếp cận vi mô tổng quát, biến đổi lý thuyết trường fermion sang lý thuyết hiệu dụng cho các tham số trật tự thông qua tích phân phiếm hàm và biến đổi Hubbard-Stratonovich, từ đó suy dẫn tường minh các hệ số của phiếm hàm Ginzburg-Landau đa thành phần và áp dụng phân tích giản đồ pha của $\text{UGe}_2$.

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 chuyển pha Landau và Ginzburg-Landau: Luận án dựa trên tiên đề Landau về sự phá vỡ đối xứng tự phát khi hệ chuyển qua nhiệt độ tới hạn $T_c$, gắn liền với sự xuất hiện của tham số trật tự $\psi$. Năng lượng tự do được khai triển thành chuỗi giải tích của tham số trật tự tuân theo các tính chất đối xứng của hệ.
  2. Lý thuyết sắt từ trong kim loại: Sử dụng mô hình Heisenberg cho các mô men spin định xứ với Hamiltonian tương tác trao đổi $H = -\sum_{i \neq j} J_{ij} \mathbf{S}_i \cdot \mathbf{S}j$, lý thuyết trường trung bình phân tử Weiss dẫn tới biểu thức nhiệt độ Curie $k_B T_c = \frac{2S(S+1)}{3} \sum_j J{ij}$. Đối với hệ spin linh động, luận án sử dụng mô hình Hubbard và tiêu chuẩn Stoner cho sắt từ: $U N(E_F) > 1$, trong đó $U$ là thế năng đẩy đẩy tại chỗ và $N(E_F)$ là mật độ trạng thái tại mức Fermi.
  3. Lý thuyết BCS tổng quát: Mở rộng lý thuyết BCS từ kết cặp singlet truyền thống (hàm khe $\hat{\Delta}{\mathbf{k}} = i\sigma_y \Psi(\mathbf{k})$ với $\Psi(\mathbf{k}) = \Psi(-\mathbf{k})$) sang kết cặp triplet tính lẻ ($\hat{\Delta}{\mathbf{k}} = (\mathbf{d}(\mathbf{k}) \cdot \boldsymbol{\sigma}) i\sigma_y$ với $\mathbf{d}(-\mathbf{k}) = -\mathbf{d}(\mathbf{k})$). Cấu trúc hàm khe triplet liên hệ mật thiết với véc tơ $\mathbf{q}(\mathbf{k}) = i \mathbf{d}(\mathbf{k}) \times \mathbf{d}^*(\mathbf{k})$; khi $\mathbf{q} \neq 0$, trạng thái kết cặp là không unita và phá vỡ đối xứng đảo thời gian.

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

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

  • Phương pháp tích phân phiếm hàm trên đại số Grassmann: Biểu diễn hàm phân bố thống kê của hệ fermion $Z = \int \mathcal{D}(c^, c) \exp(-S[c^, c])$, trong đó các toán tử sinh và hủy fermion được thay bằng các biến số Grassmann phản giao hoán.
  • Phép biến đổi Hubbard-Stratonovich (HS): Kỹ thuật phân rã số hạng tương tác bốn fermion thành tương tác bậc hai của fermion kết cặp với các trường phụ boson. Các trường phụ này sau đó đóng vai trò là các tham số trật tự của hệ trong các kênh kết cặp Cooper, kênh mật độ, hoặc kênh trao đổi spin.
  • Phương pháp hàm Green: Thiết lập hệ phương trình chuyển động cho các hàm Green cấu tử, kết hợp kỹ thuật lấy vết trên không gian spin-Nambu để tích phân các bậc tự do fermion, thu được tác dụng hiệu dụng và phiếm hàm năng lượng tự do Ginzburg-Landau.

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

Dữ liệu đối chiếu trong luận án là các số đo thực nghiệm công bố về các hợp chất fermion nặng:

  • Số liệu đo độ từ hóa, nhiệt dung riêng, điện trở dưới áp suất cao của $\text{UGe}_2$ từ các nhóm Saxena (2000), Huxley (2001), Tateiwa (2001), Bauer (2001).
  • Số liệu giản đồ pha $P - T$ của $\text{CeRhIn}_5$ xác định từ nhiệt dung riêng khi không có từ trường ngoài từ Hegger (2000), Nicklas (2001), Llobet (2004), Park (2006).
  • Số liệu tán xạ neutron và độ từ cảm trên hợp chất $\text{HoMo}_6\text{S}_8$ của Lynn (1980), Rossat (1985).

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

Chương 1: Tổng quan về sự đồng tồn tại của trật tự từ và siêu dẫn trong hệ fermion nặng

Chương 1 cung cấp cơ sở tổng quan về thực nghiệm và lý thuyết chuyển pha trong các hệ fermion nặng.

  • Hiện tượng đồng tồn tại pha từ - siêu dẫn: Trình bày bản chất của trạng thái chất lỏng Fermi trong kim loại và hai tính bất ổn định cơ bản: siêu dẫn và sóng mật độ. Trong các hợp chất fermion nặng (chứa ion actinide hoặc đất hiếm có lớp vỏ $f$ lấp đầy không hoàn toàn), hiệu ứng làm hẹp vùng năng lượng do tương tác Coulomb mạnh tạo ra các chuẩn hạt có khối lượng hiệu dụng $m^*$ lớn hơn khối lượng điện tử tự do hàng trăm đến hàng nghìn lần. Độ rộng vùng năng lượng chuẩn hạt chỉ vào cỡ vài $\text{meV}$.
  • So sánh cơ chế khối lượng nặng giữa hợp chất Cerium và Uranium: Trong hợp chất Ce (trạng thái $4f^1$ định xứ hoàn toàn), cơ chế Kondo chiếm ưu thế; tương tác kết cặp xảy ra qua trao đổi thăng giáng spin của các chuẩn hạt linh động tại véc tơ sóng phản sắt từ. Trong hợp chất Uranium (trạng thái $5f$), các điện tử vừa có tính định xứ vừa có tính linh động (tính đối ngẫu); việc kết cặp siêu dẫn được điều phối chủ yếu bởi sự trao đổi các kích thích nội lan truyền (exciton từ) của hệ con $5f$ định xứ.
  • Dữ liệu thực nghiệm trên $\text{UGe}_2$ và $\text{CeRhIn}_5$: Luận án tổng hợp các tham số chuyển pha thực nghiệm, làm rõ sự khác biệt giữa hai hệ đại diện:
Hợp chất Trật tự từ tại $P = 0$ Nhiệt độ chuyển pha từ ($P = 0$) Mô men từ tự phát ($P = 0$) Áp suất tới hạn xuất hiện siêu dẫn Nhiệt độ siêu dẫn cực đại ($T_{sc}^{\max}$) Đặc trưng đồng tồn tại pha
$\text{UGe}_2$ Sắt từ (FM), trục dễ là trục $a$ $T_c = 52\text{ K}$ $\mu_s = 1{,}4\mu_B$ $P \approx 1{,}0\text{ GPa}$ (bắt đầu pha FM-SC) $T_{sc} \approx 0{,}7 - 0{,}8\text{ K}$ (tại $P \approx 1{,}2\text{ GPa}$) Siêu dẫn nằm ngập sâu trong pha sắt từ; tồn tại hai pha sắt từ $\text{FM1}$ ($\mu \approx 0{,}9\mu_B$) và $\text{FM2}$ ($\mu \approx 1{,}5\mu_B$); cả hai biến mất tại $P_c \approx 1{,}6 - 1{,}7\text{ GPa}$.
$\text{CeRhIn}_5$ Phản sắt từ (AFM) $T_N = 3{,}8\text{ K}$ Mô men so le $\approx 0{,}8\mu_B$ $P \approx 1{,}0\text{ GPa}$ (bắt đầu pha AFM+SC) $T_c \approx 2{,}3\text{ K}$ (tại điểm $T_c = T_N$) Đồng tồn tại AFM+SC trong dải hẹp ($1{,}0 - 1{,}95\text{ GPa}$); trên $P_c^* = 1{,}95\text{ GPa}$, trật tự AFM biến mất đột ngột qua chuyển pha loại một, chuyển sang siêu dẫn thuần tới $5\text{ GPa}$.
  • Lý thuyết sắt từ: Trình bày tương tác trao đổi lượng tử bắt nguồn từ tích phân trao đổi $J$ trong bài toán hai điện tử của nguyên tử Helium. Thiết lập mẫu Heisenberg, gần đúng trường trung bình Weiss và suy dẫn tiêu chuẩn Stoner $U N(E_F) > 1$ từ sự cạnh tranh giữa thế năng tĩnh điện đẩy $\delta P = -U(\delta n)^2$ và độ tăng động năng $\delta K = \frac{1}{N(E_F)}(\delta n)^2$.
  • Lý thuyết BCS và BCS tổng quát: Điểm lại lịch sử khám phá hiệu ứng Meissner (1933), hiệu ứng đồng vị ($M^{1/2} T_c = \text{const}$ năm 1950), cơ chế hút chậm trễ qua phonon ($\tau_D \sim 10^{-12}\text{ s}$ so với thời gian truyền của điện tử với vận tốc Fermi $v_F \sim 10^8\text{ cm/s}$ qua quãng đường $\sim 10^{-4}\text{ cm}$). Phân tích trạng thái cơ bản BCS, phép biến đổi Bogoliubov và phổ kích thích chuẩn hạt $E_{\mathbf{k}} = \sqrt{\xi_{\mathbf{k}}^2 + |\Delta_{\mathbf{k}}|^2}$. Xây dựng cấu trúc ma trận hàm khe $2 \times 2$ trong không gian spin cho trường hợp spin-singlet và spin-triplet.
  • Lý thuyết Ginzburg-Landau: Phân tích nguyên lý xây dựng phiếm hàm Landau, vai trò phá vỡ đối xứng của tham số trật tự, và giới thiệu mô hình hiện tượng học Machida-Shopova cho siêu dẫn sắt từ.

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 chi tiết các công cụ toán tử và tích phân hàm nhằm xây dựng cầu nối từ Hamiltonian vi mô sang lý thuyết hiệu dụng Ginzburg-Landau.

  • Phương pháp hàm Green: Thiết lập phương trình chuyển động cho hàm Green hai thời gian và ứng dụng tính toán các giá trị trung bình thống kê của toán tử mật độ spin và toán tử kết cặp trong trạng thái siêu dẫn và sắt từ.
  • Phương pháp tích phân phiếm hàm trên biến Grassmann: Định nghĩa tích phân trên đại số Grassmann đối với các toán tử trường fermion. Chuyển đổi hàm phân bố thống kê $Z = \text{Tr}, e^{-\beta(H - \mu N)}$ thành tích phân phiếm hàm liên tục theo thời gian ảo $\tau \in [0, \beta]$: $$Z = \int \mathcal{D}(c^, c) \exp\left( -\int_0^\beta d\tau \left[ \sum_{k\sigma} c_{k\sigma}^(\tau) \partial_\tau c_{k\sigma}(\tau) + H(c^*, c) \right] \right)$$
  • Phép biến đổi Hubbard-Stratonovich: Đưa vào trường phụ liên hợp để tuyến tính hóa số hạng tương tác bốn fermion $V c^\dagger c^\dagger c c$. Luận án chứng minh cách phân tách tương tác sang kênh kết cặp Cooper.
  • Dẫn xuất phiếm hàm Ginzburg-Landau một thành phần cho siêu dẫn BCS: Bằng cách tích phân chính xác các biến Grassmann của fermion, thu được tác dụng hiệu dụng $S_{\text{eff}}[\Delta, \Delta^*]$. Khai triển tác dụng hiệu dụng quanh điểm chuyển pha theo lũy thừa của tham số trật tự $\Delta$ đến bậc bốn và bậc hai của đạo hàm không gian, luận án thu được dạng chuẩn của phiếm hàm Ginzburg-Landau một thành phần: $$F_{GL}[\Delta] = \int d^3r \left[ a |\Delta|^2 + \frac{b}{2} |\Delta|^4 + c |\nabla \Delta|^2 \right]$$ Phương pháp vi mô cho phép xác định tường minh các hệ số $a, b, c$ thông qua mật độ trạng thái tại mức Fermi $N(0)$ và vận tốc Fermi $v_F$, khẳng định lại kết quả kinh điển của Gor'kov từ hình thức luận tích phân phiếm hàm hiện đại.

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 của tác giả, tổng quát hóa phương pháp vi mô cho các hệ có nhiều tham số trật tự cạnh tranh hoặc đồng tồn tại.

  • Mô hình ba tham số trật tự: Xuất phát từ Hamiltonian tổng quát của hệ fermion tương tác nhiều hạt, luận án đồng thời đưa vào ba kênh phân rã của phép biến đổi Hubbard-Stratonovich:
    1. Kênh "mật độ" (density channel): tương ứng với sóng mật độ điện tích (CDW).
    2. Kênh "kết cặp" (pairing channel): tương ứng với các cặp Cooper của trạng thái siêu dẫn (SC).
    3. Kênh "trao đổi" (exchange channel): tương ứng với độ từ hóa spin hoặc sóng mật độ spin (SDW/FM).
  • Khử tham số kết cặp và tích phân các bậc tự do fermion: Sau khi áp dụng phép biến đổi Hubbard-Stratonovich ba kênh, Hamiltonian chuyển thành dạng song tuyến tính đối với các trường fermion, ghép cặp với ba trường phụ. Sử dụng kỹ thuật ma trận Nambu mở rộng và hàm Green tương ứng, tác giả tích phân toàn bộ các trường fermion để thu được định thức phiếm hàm (functional determinant).
  • Thiết lập phiếm hàm Ginzburg-Landau ba thành phần: Khai triển logarit của định thức ma trận hàm Green theo các trường phụ, tác giả thiết lập biểu thức giải tích của phiếm hàm Ginzburg-Landau chứa đồng thời ba tham số trật tự, bao gồm các số hạng bậc hai, bậc bốn độc lập của từng tham số và các số hạng tương tác chéo giữa các tham số trật tự khác nhau.
  • Các mô hình đơn giản: Khảo sát 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. Trọng tâm là phân tích kênh trao đổi spin và sự kết hợp giữa kênh trao đổi với kênh kết cặp siêu dẫn, tạo tiền đề vi mô cho việc khảo sát hệ hai thành phần ở 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 vận dụng khung lý thuyết vi mô đã xây dựng để giải quyết bài toán đồng tồn tại pha sắt từ và siêu dẫn spin-triplet trong hợp chất $\text{UGe}_2$.

  • Nguồn gốc vi mô của phiếm hàm Ginzburg-Landau hai thành phần: Xuất phát từ Hamiltonian vi mô chứa tương tác trao đổi spin giữa các mô men $5f$ của Uranium và tương tác hút spin-triplet giữa các điện tử dẫn, tác giả sử dụng phép biến đổi HS hai kênh để suy dẫn phiếm hàm năng lượng GL hai thành phần phụ thuộc vào véc tơ từ hóa $\mathbf{M}$ và véc tơ trật tự siêu dẫn phức $\boldsymbol{\psi}$.
  • Biểu thức giải tích của phiếm hàm năng lượng tự do: Dạng phiếm hàm thu được bao gồm: $$f_{GL}(\boldsymbol{\psi}, \mathbf{M}) = a_s |\boldsymbol{\psi}|^2 + \frac{b_s}{2} |\boldsymbol{\psi}|^4 + a_f |\mathbf{M}|^2 + \frac{b_f}{2} |\mathbf{M}|^4 + f_I(\boldsymbol{\psi}, \mathbf{M})$$ với số hạng tương tác tương hỗ có dạng đối xứng: $$f_I(\boldsymbol{\psi}, \mathbf{M}) = i \gamma_0 \mathbf{M} \cdot (\boldsymbol{\psi} \times \boldsymbol{\psi}^*)$$ Sự xuất hiện của số hạng tương tác này phản ánh tính chất không unita của cặp Cooper spin-triplet trong môi trường phân cực từ: véc tơ từ hóa tự phát đóng vai trò định hướng trục lượng tử hóa của các cặp Cooper spin song song ($S=1$).
  • Tính toán giản đồ pha: Bằng cách cực tiểu hóa phiếm hàm tự do theo hai tham số trật tự $|\boldsymbol{\psi}|$ và $|\mathbf{M}|$, tác giả giải hệ phương trình trạng thái tự hợp để xác định các ranh giới pha.
    • Thiết lập giản đồ pha nhiệt độ - áp suất ($T - P$) cho $\text{UGe}2$ với các thông số đầu vào thực nghiệm: $T_s = 0$, $T{f0} = 52\text{ K}$, áp suất tới hạn chuẩn hóa $P_c = 1$. Kết quả tính toán thể hiện miền pha sắt từ - siêu dẫn đồng tồn tại (kí hiệu FS) được làm đậm, nằm hoàn toàn bên trong miền sắt từ.
    • Đường nét liền trên đồ thị biểu thị đường chuyển pha loại hai giữa pha sắt từ thông thường (FM) và pha sắt từ siêu dẫn (FS).
    • Phân tích sự phụ thuộc của cấu trúc pha trong mặt phẳng thông số rút gọn $(t, r)$ khi hệ số ghép cặp tương tác thay đổi từ $\gamma = 0$ đến $\gamma \neq 0$. Kết quả khẳng định tương tác chéo giữa tham số trật tự từ và siêu dẫn triplet là điều kiện cần thiết để ổn định hóa pha siêu dẫn dưới áp suất cao trước khi trật tự từ bị dập tắt hoàn toàn tại điểm tới hạn lượng tử.

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

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

  1. Luận án đã xây dựng thành công phương pháp tiếp cận vi mô dẫn xuất phiếm hàm năng lượng Ginzburg-Landau nhiều thành phần thông qua việc kết hợp phép biến đổi Hubbard-Stratonovich đa kênh và hình thức luận tích phân phiếm hàm Grassmann. Phương pháp này khắc phục hạn chế lớn của các mô hình hiện tượng học trước đó (vốn chỉ dựa vào đối xứng vĩ mô và gán ghép các hệ số tự do), thiết lập được mối liên hệ giải tích trực tiếp giữa các hệ số trong khai triển Ginzburg-Landau với các thông số vi mô của vật liệu (mật độ trạng thái tại mặt Fermi, vận tốc Fermi, cường độ tương tác trao đổi spin và tương tác hút kết cặp).
  2. Xây dựng biểu diễn hình thức luận cho hệ ba tham số trật tự (kênh mật độ, kết cặp Cooper và trao đổi spin) và suy dẫn phiếm hàm Ginzburg-Landau hai thành phần cho siêu dẫn sắt từ spin-triplet. Luận án làm rõ nguồn gốc vi mô của số hạng ghép cặp tương tác $f_I(\boldsymbol{\psi}, \mathbf{M}) = i \gamma_0 \mathbf{M} \cdot (\boldsymbol{\psi} \times \boldsymbol{\psi}^*)$, chứng minh rằng trật tự sắt từ nội tại có thể hỗ trợ và ổn định hóa trạng thái kết cặp spin-triplet không unita.

Đóng góp mới về mặt thực tiễn

  1. Áp dụng thành công mô hình vi mô để mô tả sự đồng tồn tại pha trong hợp chất fermion nặng $\text{UGe}_2$, tái hiện về mặt định tính và bán định lượng giản đồ pha nhiệt độ - áp suất ($T - P$).
  2. Xác định các miền pha cân bằng và bản chất các đường ranh giới chuyển pha (loại một và loại hai) giữa các pha thuận từ (PM), sắt từ (FM), và sắt từ - siêu dẫn (FS), cung cấp cơ sở lý thuyết giải thích tại sao pha siêu dẫn lại xuất hiện ngập sâu trong miền trật tự từ ở nhiệt độ thấp dưới áp suất cao.

Kiến nghị và đề xuất

  • Mở rộng hình thức luận vi mô bằng cách đưa thêm các số hạng bất đẳng hướng tinh thể và tương tác spin-quỹ đạo nội tại của các ion actinide nhằm phản ánh chính xác cấu trúc tinh thể trực giao của $\text{UGe}_2$.
  • Áp dụng phương pháp phân tách kênh Hubbard-Stratonovich đa thành phần để khảo sát các hệ vật liệu sắt từ siêu dẫn liên quan như $\text{URhGe}$ và $\text{UCoGe}$, nơi tính siêu dẫn xuất hiện ngay tại áp suất thường hoặc được kích thích bởi từ trường ngoài mạnh.

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

Hạn chế được chỉ ra từ văn bản

  • Các phép tính toán dẫn xuất phiếm hàm và giải giản đồ pha trong luận án chủ yếu thực hiện trong gần đúng trường trung bình (mean-field approximation), do đó đã bỏ qua đóng góp của các thăng giáng lượng tử mạnh ở sát điểm tới hạn lượng tử sắt từ ($P \approx P_c$).
  • Mô hình tương tác vi mô được đơn giản hóa bằng cách gộp các bậc tự do quỹ đạo nội của điện tử lớp $5f$ (vốn bị phân tách phức tạp bởi trường tinh thể và liên kết spin-quỹ đạo) thành các tham số hiệu dụng.
  • Luận án giới hạn khảo sát ở trạng thái kết cặp unita trong một số bước tính gần đúng và chưa giải bài toán vi mô đầy đủ cho các hàm khe có cấu trúc nút bất đẳng hướng phức tạp trên mặt Fermi.

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

  • Phát triển phương pháp tích phân phiếm hàm vượt ra ngoài gần đúng trường trung bình bằng cách tính đến các số hạng thăng giáng bậc cao của các trường phụ xung quanh điểm yên ngựa (saddle-point).
  • Khảo sát chi tiết sự cạnh tranh giữa hai pha sắt từ $\text{FM1}$ và $\text{FM2}$ trong $\text{UGe}_2$ dưới tác dụng của từ trường ngoài theo các trục tinh thể khác nhau.
  • Tính toán cụ thể các đại lượng nhiệt động có thể đo đạc trực tiếp như độ xuyên sâu từ trường $\lambda(T)$, nhiệt dung riêng $C(T)$ và độ từ hóa $M(T, P)$ từ phiếm hàm vi mô để so sánh định lượng với thực nghiệm.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị học thuật trong lĩnh vực vật lý lý thuyết và vật lý chất rắn, cụ thể:

  • Đối với nghiên cứu sinh và giảng viên chuyên ngành Vật lý lý thuyết và Vật lý toán: Chương 2 và Chương 3 cung cấp một quy trình mẫu mực về việc ứng dụng tích phân phiếm hàm trên đại số Grassmann và biến đổi Hubbard-Stratonovich để chuyển đổi bài toán nhiều hạt tương tác mạnh thành lý thuyết trường hiệu dụng. Các phụ lục A, B, C (chi tiết về tính toán hệ fermion hiệu dụng bậc hai, khai triển hàm mũ chứa đạo hàm và các phép vi phân toán tử) là tài liệu tham khảo trực tiếp cho các phép tính giải tích phức tạp trong lý thuyết chất ngưng tụ.
  • Đối với các nhà nghiên cứu thực nghiệm và lý thuyết về vật liệu từ - siêu dẫn: Chương 1 và Chương 4 cung cấp bức tranh tổng thể về các hợp chất fermion nặng ($\text{UGe}_2$, $\text{CeRhIn}_5$, $\text{HoMo}_6\text{S}_8$), cách phân tích ranh giới pha trên giản đồ $T - P$ và cơ chế vi mô của siêu dẫn spin-triplet dưới áp suất cao.

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

1. Tại sao trạng thái siêu dẫn truyền thống theo lý thuyết BCS lại không thể đồng tồn tại với trật tự sắt từ?

Trong lý thuyết BCS truyền thống, các điện tử kết cặp thành cặp Cooper spin-singlet với tổng spin $S = 0$ (hai điện tử có spin đối song). Khi đặt trong một 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 một trường phân cực từ rất mạnh. Trường này phá hủy cấu hình spin đối song của cặp Cooper singlet, dẫn đến sự suy giảm nhanh chóng và triệt tiêu trạng thái siêu dẫn.

2. Trạng thái kết cặp Cooper trong hợp chất $\text{UGe}_2$ thuộc loại nào và tại sao nó tồn tại được trong pha sắt từ?

Trong $\text{UGe}_2$, các cặp Cooper thuộc loại spin-triplet có spin toàn phần $S = 1$ với hàm sóng không gian có tính đối xứng lẻ (sóng $p$). Các trạng thái spin bao gồm $|!\uparrow\uparrow\rangle$ ($S_z = 1$), $|!\downarrow\downarrow\rangle$ ($S_z = -1$) và $(|!\uparrow\downarrow\rangle + |!\downarrow\uparrow\rangle)/\sqrt{2}$ ($S_z = 0$). Do các điện tử trong cặp có spin song song cùng chiều, chúng không bị phá hủy bởi trường phân cực từ nội tại mà ngược lại, tương tác trao đổi và thăng giáng spin sắt từ đóng vai trò như một môi trường trung gian kích thích sự hình thành cặp Cooper.

3. Phép biến đổi Hubbard-Stratonovich đóng vai trò gì trong phương pháp luận của luận án?

Phép biến đổi Hubbard-Stratonovich là một công cụ toán học trong tích phân phiếm hàm, cho phép chuyển đổi số hạng tương tác bốn toán tử fermion (tương tác hai hạt phức tạp) thành các số hạng bậc hai của fermion liên kết với các trường phụ boson (auxiliary fields). Sau khi chuyển đổi, các trường fermion có thể được tích phân chính xác bằng đại số Grassmann, để lại một tác dụng hiệu dụng chỉ còn phụ thuộc vào các trường phụ. Các trường phụ này chính là các tham số trật tự của hệ (như độ từ hóa và khe năng lượng siêu dẫn).

4. Dưới tác dụng của áp suất ngoài, trật tự từ và siêu dẫn trong $\text{UGe}_2$ biến đổi như thế nào theo thực nghiệm?

Ở áp suất $P = 0$, $\text{UGe}2$ là chất sắt từ có nhiệt độ Curie $T_c = 52\text{ K}$ và độ từ hóa bão hòa $\mu_s = 1{,}4\mu_B$. Khi tăng áp suất lên $P \approx 1{,}0\text{ GPa}$, hệ xuất hiện pha siêu dẫn đồng tồn tại với sắt từ (FM-SC). Nhiệt độ siêu dẫn đạt cực đại $T{sc} \approx 0{,}7 - 0{,}8\text{ K}$ quanh áp suất $1{,}2\text{ GPa}$. Trong miền sắt từ, hệ trải qua chuyển pha giữa hai pha từ $\text{FM2}$ ($\mu \approx 1{,}5\mu_B$) và $\text{FM1}$ ($\mu \approx 0{,}9\mu_B$) tại $P_x \approx 1{,}2\text{ GPa}$. Khi áp suất tiếp tục tăng đến áp suất tới hạn $P_c \approx 1{,}6 - 1{,}7\text{ GPa}$, cả tính sắt từ và tính siêu dẫn đều biến mất đồng thời, hệ chuyển sang trạng thái dẫn thường thuận từ (PM).

5. Tiêu chuẩn Stoner cho tính sắt từ trong hệ spin linh động được phát biểu như thế nào trong luận án?

Theo mô hình Hubbard trong luận án, tiêu chuẩn Stoner xuất hiện từ sự cạnh tranh năng lượng khi làm lệch mật độ điện tử giữa hai trạng thái spin hướng lên và hướng xuống một lượng $\delta n$. Thế năng tương tác giảm một lượng $\delta P = -U(\delta n)^2$, trong khi động năng tăng do phải chiếm các mức năng lượng cao hơn mức Fermi một lượng $\delta K = \frac{1}{N(E_F)}(\delta n)^2$. Hệ sẽ ưu tiên chuyển sang trạng thái sắt từ khi tổng biến thiên năng lượng $\delta E = \delta P + \delta K < 0$, dẫn tới điều kiện Stoner: $U N(E_F) > 1$.

Kết luận

Luận án của tác giả Nguyễn Văn Hinh là một công trình nghiên cứu hệ thống thuộc chuyên ngành Vật lý lý thuyết và Vật lý toán, tập trung giải quyết bài toán phức tạp về sự đồng tồn tại giữa trật tự từ và siêu dẫn trong các hệ tương quan mạnh. Bằng việc kết hợp chặt chẽ hình thức luận tích phân phiếm hàm Grassmann, phép biến đổi Hubbard-Stratonovich đa kênh và kỹ thuật hàm Green, tác giả đã suy dẫn thành công phiếm hàm Ginzburg-Landau nhiều thành phần trực tiếp từ các mô hình vi mô. Kết quả áp dụng cho hợp chất fermion nặng $\text{UGe}_2$ không chỉ giải thích định tính và bán định lượng giản đồ pha nhiệt độ - áp suất mà còn làm sáng tỏ cơ chế ghép cặp spin-triplet trong môi trường sắt từ, tạo tiền đề phương pháp luận quan trọng cho các nghiên cứu tiếp theo về các trạng thái ngưng tụ phi truyền thống.