Tổng quan luận án

Luận án tiến sĩ vật lý với 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 thực hiện bởi nghiên cứu sinh Nguyễn Văn Hinh, dưới sự hướng dẫn khoa học của TS. Nguyễn Trí Lân và GS. Nguyễn Toàn Thắng, hoàn thành 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).

Tính cấp thiết của luận án xuất phát từ sự biến đổi trong nhận thức vật lý chất rắn về mối quan hệ giữa tính siêu dẫn và từ tính. Theo 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 thông qua lực hút gián tiếp giữa các điện tử nhờ dao động mạng tinh thể (phonon), tạo nên các cặp Cooper spin-singlet với spin đối song song ($S = 0$). Dưới tác dụng của trường từ cục bộ sinh ra bởi các mômen từ hoặc tạp chất từ, các cặp Cooper singlet bị bẻ gãy, dẫn đến sự suy giảm nhanh chóng của nhiệt độ chuyển pha siêu dẫn $T_{sc}$. Quan niệm hai pha này loại trừ nhau tồn tại trong thời gian dài cho đến thập niên 1980 khi hiện tượng siêu dẫn đồng tồn tại với trật tự phản sắt từ được phát hiện trong các hợp chất fermion nặng. Sự kiện phát hiện tính siêu dẫn xuất hiện sâu trong pha sắt từ dưới áp suất cao ở hợp chất $UGe_2$ vào năm 2000 (Saxena và cộng sự), tiếp nối bởi các vật liệu $UIr$, $URhGe$ và $UCoGe$, đã đặt ra yêu cầu phải xây dựng cơ sở lý thuyết mới.

Khoảng trống nghiên cứu được tác giả xác định là:

  1. Bản chất ghép cặp của pha siêu dẫn đồng tồn tại với sắt từ (thuộc loại spin-singlet hay spin-triplet) và cơ chế hình thành trật tự từ (do mômen định xứ của điện tử $5f$ hay do hệ điện tử linh động) vẫn còn nhiều bất đồng và chưa có lời giải lý thuyết vi mô trọn vẹn.
  2. Nguyên nhân thực sự dẫn dắt tính siêu dẫn trong các hệ này (do tương tác trao đổi phonon, magnon, kích thích exciton từ trường tinh thể hay sóng mật độ điện tích) chưa được định lượng rõ ràng.
  3. Sự phụ thuộc của nhiệt độ chuyển pha và độ từ hóa tự phát vào áp suất chưa được tính toán một cách hệ thống từ một mô hình vi mô nhất quán; các mô hình hiện tượng học Ginzburg-Landau hiện có chủ yếu dựa trên các lập luận đối xứng vĩ mô mà thiếu mối liên hệ tường minh với các tham số vi mô của vật liệu.

Mục tiêu nghiên cứu của luận án được xác định gồm ba nhiệm vụ:

  1. Thực hiện phép biến đổi lý thuyết trường fermion sang một lý thuyết trường 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ĩ khác nhau.
  2. Xây dựng một phiếm hàm năng lượng Ginzburg-Landau đa thành phần mô tả sự đồng tồn tại của nhiều pha từ Hamiltonian vi mô thông qua phép biến đổi Hubbard-Stratonovich và kỹ thuật hàm Green.
  3. Thiết lập biểu diễn toán học chính thức cho phiếm hàm Ginzburg-Landau hai tham số trật tự (sắt từ và siêu dẫn spin-triplet), áp dụng khảo sát cụ thể giản đồ pha cân bằng nhiệt độ – áp suất của hợp chất fermion nặng $UGe_2$.

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

  • Đối tượng: Hệ nhiều hạt fermion tương tác mạnh, tập trung vào các hợp chất fermion nặng gốc actinide và đất hiếm, trọng tâm là hợp chất siêu dẫn sắt từ $UGe_2$ và hợp chất phản sắt từ $CeRhIn_5$.
  • Phạm vi không gian: Cấu trúc tinh thể trực giao của $UGe_2$ với trục từ hóa dễ dọc theo trục $a$, xét ở quy mô mạng tinh thể vĩ mô thông qua các trường vi mô của điện tử linh động và điện tử định xứ $5f$.
  • Phạm vi thời gian và điều kiện vật lý: Trạng thái cân bằng nhiệt động tại vùng nhiệt độ thấp (từ 0 K đến nhiệt độ Curie $T_c = 52\text{ K}$, nhiệt độ siêu dẫn $T_{sc} \approx 0,7\text{ K}$) và dải áp suất thủy tĩnh mở rộng từ 0 đến trên $1,7\text{ GPa}$.

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

Luận án điểm lại các mốc thực nghiệm và lý thuyết then chốt định hình hướng nghiên cứu:

Về mặt thực nghiệm:

  • Quan sát đầu tiên về trạng thái điện trở bằng không trong pha sắt từ ở hợp chất $HoMo_6S_8$ được công bố bởi Lynn và cộng sự (1981) [57], tiếp nối bởi Genicon và cộng sự [31]. Các phép đo nhiễu xạ neutron trên đơn tinh thể $HoMo_6S_8$ của Rossat và cộng sự (1985) [71] xác nhận trật tự sắt từ xuất hiện ở $T < 0,60\text{ K}$, độ từ hóa bão hòa quanh $T = 0,4\text{ K}$, với 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}$.
  • Phát hiện siêu dẫn trong kim loại fermion nặng $CeCu_2Si_2$ ở $T_c \approx 0,5\text{ K}$ của Steglich và cộng sự (1979) [81], chứng minh các chuẩn hạt có khối lượng hiệu dụng $m^* \gg m_e$ tham gia vào trạng thái ngưng tụ.
  • Phát hiện siêu dẫn sắt từ trong $UGe_2$ dưới áp suất của Saxena và cộng sự (2000) [74], kèm theo các khảo sát chi tiết của Huxley và cộng sự (2001) [43], Tateiwa và cộng sự (2001) [85], Kobayashi và cộng sự (2002) [38].
  • Các quan sát thực nghiệm trên các chất siêu dẫn sắt từ khác: $UIr$ (Akazawa và cộng sự, 2004 [7]), $URhGe$ (Aoki và cộng sự, 2001 [10]), và $UCoGe$ (Huy và cộng sự, 2007 [44]).
  • Khảo sát thực nghiệm trên hợp chất phản sắt từ $CeRhIn_5$ của Hegger và cộng sự (2000) [40], Knebel và cộng sự (2004) [49], Park và cộng sự (2006) [67], Mito và cộng sự (2001) [92], Christianson và cộng sự (2004) [11], Llobet và cộng sự (2004) [32].

Về mặt lý thuyết:

  • Lý thuyết vi mô BCS của Bardeen, Cooper, Schrieffer (1957) [12, 13] và cơ chế hút gián tiếp qua phonon của Fröhlich (1950) [30].
  • Gor'kov (1959) [34] chứng minh lý thuyết hiện tượng học Ginzburg-Landau (1950) [33] có thể được suy dẫn vi mô từ lý thuyết BCS thông qua kỹ thuật hàm Green.
  • Abrikosov (1957) [6] phát triển lý thuyết siêu dẫn loại II và mạng xoáy lượng tử từ thông $\Phi_0 = hc/2e$.
  • Các mô hình hiện tượng học Ginzburg-Landau cho siêu dẫn sắt từ của Machida và Ohmi (2001) [58], Shopova và cộng sự (2003, 2005) [76, 77, 78, 79], khảo sát số hạng ghép cặp dạng $i\gamma_0 \mathbf{M} \cdot (\boldsymbol{\psi} \times \boldsymbol{\psi}^*)$.
  • Mô hình vi mô kết hợp tương tác BCS tổng quát và tương tác trao đổi Heisenberg của Dahl và Sudbø (2006) [23] cho siêu dẫn sóng $p$.
  • Các nghiên cứu trong nước: GS. Đỗ Trần Cát và cộng sự (2004, 2006) [19, 20] nghiên cứu pha siêu dẫn và sắt từ trong $UGe_2$ dựa trên cấu trúc dải kết ô; nhóm Hệ tương quan mạnh tại Viện Vật lý nghiên cứu các pha $FM_1$, $FM_2$ [87] và cơ chế exciton từ do tách mức trường tinh thể [88].

Khoảng trống luận án lựa chọn giải quyết: Luận án không đi theo hướng tiếp cận hiện tượng luận đơn thuần mà xây dựng cầu nối vi mô: sử dụng hình thức luận tích phân phiếm hàm trên biến Grassmann và biến đổi Hubbard-Stratonovich để khử các bậc tự do fermion, từ đó thiết lập biểu thức giải tích tường minh cho các hệ số của phiếm hàm Ginzburg-Landau đa thành phần, giải thích trạng thái đồng tồn tại sắt từ và siêu dẫn triplet dưới ảnh hưởng của nhiệt độ và áp suất.

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

Cơ sở lý thuyết và mô hình toán lý

Luận án vận dụng hệ thống các lý thuyết vật lý chất rắn và lý thuyết trường nhiều hạt:

  1. Lý thuyết tương tác trao đổi và mẫu Heisenberg: Xuất phát từ phương trình Schrödinger cho nguyên tử Helium hai điện tử với hàm sóng không gian đối xứng $\Phi_s$ và phản đối xứng $\Phi_a$. Năng lượng trạng thái tương ứng là $E_1 = E^0 + K + J$ (spin-singlet, $S = 0$) và $E_2 = E^0 + K - J$ (spin-triplet, $S = 1$), trong đó $K$ là tích phân tương tác tĩnh điện Coulomb và $J$ là tích phân trao đổi sinh ra do tính không phân biệt của hạt đồng nhất. Toán tử Hamiltonian Heisenberg trong không gian mạng tinh thể được viết dưới dạng: $$\hat{H} = - \sum_{i \neq j} J_{ij} \hat{\mathbf{S}}i \cdot \hat{\mathbf{S}}j$$ với $J{ij} > 0$ đặc trưng cho trật tự sắt từ và $J{ij} < 0$ cho phản sắt từ. Trong lý thuyết trường phân tử Weiss, từ trường hiệu dụng được xác định bởi $\mathbf{H}W = 2 \sum_j J{ij} \langle \mathbf{S}j \rangle$, dẫn đến độ từ hóa tính theo hàm Brillouin $B_S(x)$ và hệ thức Curie: $$k_B T_c = \frac{2S(S+1)}{3} \sum_j J{ij}$$
  2. Mô hình Hubbard cho hệ spin linh động: Miêu tả động năng nhảy nút của điện tử giữa các vị trí lân cận và lực đẩy tĩnh điện tại chỗ: $$H = -t \sum_{\langle jl \rangle, \sigma} c_{j\sigma}^\dagger c_{l\sigma} + U \sum_j n_{j\uparrow} n_{j\downarrow} - \mu \sum_j (n_{j\uparrow} + n_{j\downarrow})$$ Khảo sát ở hai giới hạn giải được chính xác: giới hạn không tương tác $U = 0$ với tán sắc $\xi_\mathbf{k} = -2t \sum_{\nu=1}^d \cos k_\nu - \mu$ và giới hạn nguyên tử $t = 0$. Tiêu chuẩn Stoner cho sắt từ linh động xuất hiện khi tương tác đẩy vượt qua mật độ trạng thái tại mức Fermi: $$U N(E_F) > 1$$
  3. Lý thuyết BCS và BCS tổng quát: Mô tả sự hình thành cặp Cooper thông qua phân cực mạng tinh thể chậm trễ (thời gian hồi đáp ion $\tau_D \sim 10^{-12}\text{ s}$, vận tốc đ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 bán kính chắn Coulomb $1\text{ - }2\text{ \AA}$). Trong lý thuyết tổng quát, ma trận khe năng lượng $2 \times 2$ trong không gian spin được tham số hóa: $$\hat{\Delta}\mathbf{k} = \begin{pmatrix} \Delta{\mathbf{k},\uparrow\uparrow} & \Delta_{\mathbf{k},\uparrow\downarrow} \ \Delta_{\mathbf{k},\downarrow\uparrow} & \Delta_{\mathbf{k},\downarrow\downarrow} \end{pmatrix}$$ Đối với trạng thái chẵn spin-singlet, $\hat{\Delta}_\mathbf{k} = i\hat{\sigma}y \Psi(\mathbf{k})$. Đối với trạng thái lẻ spin-triplet, ma trận được biểu diễn qua véctơ ba thành phần $\mathbf{d}(\mathbf{k})$: $$\hat{\Delta}\mathbf{k} = i (\mathbf{d}(\mathbf{k}) \cdot \hat{\boldsymbol{\sigma}}) \hat{\sigma}_y$$ với điều kiện trạng thái kết cặp unita là $\mathbf{q}(\mathbf{k}) = i \mathbf{d}(\mathbf{k}) \times \mathbf{d}^*(\mathbf{k}) = 0$.
  4. Lý thuyết chuyển pha Landau và Ginzburg-Landau: Khai triển năng lượng tự do $F = F_0(T) + F_L(T, \psi)$ theo các bất biến đối xứng của tham số trật tự $\psi$, cho phép xác định bản chất chuyển pha loại một hoặc loại hai mà không phụ thuộc vào chi tiết vi mô tức thời.

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

  • 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ố nhiệt động $Z$ dưới dạng tích phân lộ trình trên không gian các biến phản giao hoán Grassmann $\bar{\psi}, \psi$.
  • Phép biến đổi Hubbard-Stratonovich (HS): Tách số hạng tương tác bốn hạt fermion thành các trường phụ hiệu dụng liên kết qua ba kênh: kênh mật độ (charge-density), kênh kết cặp Cooper (superconductivity), và kênh trao đổi (magnetism).
  • Phương pháp hàm Green hai thời gian: Thiết lập phương trình chuyển động của hàm Green, kết hợp phép biến đổi Bogoliubov-Valatin chéo hóa Hamiltonian hiệu dụng.
  • Kỹ thuật khai triển đạo hàm (Gradient expansion): Khai triển vết logarit ($\mathrm{Tr}\ln$) của ma trận truyền Green và các toán tử hàm mũ chứa vi phân không gian để thu nhận các số hạng vi phân bậc hai trong phiếm hàm năng lượng Ginzburg-Landau.

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

Dữ liệu sử dụng trong luận án là dữ liệu thứ cấp từ các kết quả đo thực nghiệm đã được công bố quốc tế về hợp chất $UGe_2$ và $CeRhIn_5$, bao gồm các số đo độ từ hóa dưới áp suất cao, nhiệt dung riêng, điện trở suất, giản đồ pha nhiệt độ – áp suất, và các thông số mạng tinh thể.

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 này hệ thống hóa các kiến thức cơ sở thực nghiệm và lý thuyết về hiện tượng chuyển pha trong kim loại tương quan mạnh:

  • Phân tích hiện tượng đồng tồn tại pha: Giới thiệu trạng thái chất lỏng Fermi nặng với 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 lần, xuất phát từ sự lai hóa giữa điện tử dẫn rộng và các lớp vỏ $4f$ ($Ce$) hoặc $5f$ ($U$). Dẫn xuất các cơ chế kích thích năng lượng thấp, exciton từ trường tinh thể, và các bất ổn định sóng mật độ spin (SDW), sóng mật độ điện tích (CDW).
  • Số liệu thực nghiệm trên $UGe_2$ và $CeRhIn_5$: Tổng hợp số liệu thực nghiệm từ các giản đồ pha của hai hợp chất đại diện:
Hợp chất Thông số áp suất ($P$) Nhiệt độ chuyển pha Trật tự từ / Mômen từ Trạng thái pha
$UGe_2$ $P = 0\text{ GPa}$ $T_c = 52\text{ K}$ Sắt từ trục $a$, $\mu_s = 1,4\mu_B$ Pha sắt từ ($FM$)
$P = 1,0\text{ GPa}$ Chuyển pha lượng tử Phân cực spin Chuyển sang $FM\text{-}SC$
$P \approx 1,2\text{ GPa}$ $T_{sc}^{\max} = 0,7\text{ K}$ $FM_2 \to FM_1$ tại $T_x = T_{sc}$ Siêu dẫn cực đại trong $FM$
$P \approx 1,2\text{ GPa}$ $T_x \approx 30\text{ K}$ (ở $P=0$) $FM_2$ ($\mu = 1,5\mu_B$), $FM_1$ ($\mu = 0,9\mu_B$) Chuyển pha giữa 2 pha từ
$P_c \approx 1,7\text{ GPa}$ $T_c \to 0\text{ K}$, $T_{sc} \to 0\text{ K}$ Mất từ tính tự phát Chuyển sang thuận từ ($PM$)
$CeRhIn_5$ $P = 0\text{ GPa}$ $T_N = 3,8\text{ K}$ Phản sắt từ so le Pha $AFM$ thuần túy
$P \sim 0,5\text{ GPa}$ $T_N$ đạt cực đại trơn Trật tự $AFM$ $AFM$ ổn định
$P < P_c^*$ $T_c < T_N$ Đồng tồn tại Pha $AFM + SC$
$P_c^* = 1,95\text{ GPa}$ $T_c = T_N$ Chuyển pha bậc nhất Mất trật tự $AFM$
$1,95\text{ - }5,0\text{ GPa}$ $T_{sc}$ duy trì Không có từ tính Pha $SC$ phi truyền thống
  • Cơ chế vi mô của từ tính và siêu dẫn: Trình bày chi tiết phép giải bài toán Helium dẫn tới tích phân trao đổi $J$, mô hình Heisenberg, lý thuyết trường phân tử Weiss, mô hình Hubbard với tiêu chuẩn Stoner $UN(E_F) > 1$. Phân tích chi tiết cơ chế tương tác hút điện tử do trao đổi phonon trễ thời gian theo Fröhlich.
  • Lý thuyết BCS tổng quát và ma trận khe: Dẫn xuất hệ phương trình Bogoliubov-Valatin và phương trình khe tự hợp cho kênh singlet và triplet. Thiết lập phiếm hàm Ginzburg-Landau hiện tượng học cho hệ siêu dẫn sắt từ với số hạng ghép đôi $f_I(\psi, \mathbf{M}) = i\gamma_0 \mathbf{M} \cdot (\boldsymbol{\psi} \times \boldsymbol{\psi}^*)$.

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

Chương 2 thiết lập công cụ toán học và kỹ thuật lý thuyết trường lượng tử:

  • Phương pháp hàm Green: Xây dựng hệ phương trình chuyển động cho hàm Green hai thời gian đối với các toán tử fermion. Áp dụng giải bài toán phổ chuẩn hạt trong chất sắt từ và chất siêu dẫn BCS, minh họa khe năng lượng đóng vai trò tham số trật tự.
  • Phương pháp tích phân phiếm hàm: Thiết lập giải tích tích phân trên đại số Grassmann. Xây dựng hàm phân bố thống kê: $$Z = \int \mathcal{D}[\bar{\psi}, \psi] \exp\left( - \int_0^\beta d\tau \int d\mathbf{r} , \mathcal{L}[\bar{\psi}, \psi] \right)$$
  • Áp dụng cho hệ siêu dẫn BCS: Thực hiện biến đổi Hubbard-Stratonovich đưa vào trường phụ $\Delta(\mathbf{r}, \tau)$ tương ứng với khe siêu dẫn. Tích phân trực tiếp các biến fermion để thu được tác dụng hiệu dụng phụ thuộc vào trường $\Delta$. Thực hiện khai triển chuỗi perturbation và khai triển gradient đối với số hạng vết logarit ($\mathrm{Tr}\ln$), từ đó dẫn xuất vi mô phiếm hàm Ginzburg-Landau một thành phần chuẩn bao gồm số hạng năng lượng thế pha và số hạng gradient không gian.

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 mở rộng hình thức luận tích phân phiếm hàm cho các hệ nhiều tham số trật tự:

  • Mô hình ba tham số trật tự: Xây dựng Hamiltonian tổng quát chứa tương tác hai hạt có thể phân rã đồng thời theo ba kênh vật lý:
    1. Kênh mật độ điện tích (CDW).
    2. Kênh kết cặp Cooper (siêu dẫn).
    3. Kênh trao đổi spin (từ tính).
  • Biến đổi Hubbard-Stratonovich đa kênh: Tách tương tác bốn fermion bằng ba trường phụ độc lập.
  • Khử tham số kết cặp: Sử dụng kỹ thuật tích phân phiếm hàm fermion trong không gian Nambu mở rộng. Bằng cách lấy tích phân trên các trường Grassmann, hệ fermion được tích phân ra ngoài hoàn toàn, để lại một tác dụng hiệu dụng chỉ chứa các trường tham số trật tự.
  • Dẫn xuất phiếm hàm GL ba thành phần: Khai triển tác dụng hiệu dụng theo các bậc của các tham số trật tự quanh điểm chuyển pha. Xác định tường minh các hệ số của các số hạng bậc hai, bậc bốn và các số hạng tương tác chéo giữa các tham số trật tự.
  • Khảo sát các mô hình đơn giản: Phân tích trường hợp giới hạn khi chỉ có kênh trao đổi, làm rõ cách thức các tham số vi mô trong Hamiltonian gốc quy định dấu và độ lớn của các hệ số trong phiếm hàm Landau.

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 các kết quả giải tích vào bài toán vật liệu thực tế $UGe_2$:

  • Nguồn gốc vi mô của phiếm hàm Ginzburg-Landau hai thành phần: Thiết lập mô hình vi mô cho các điện tử $5f$ của Uranium vừa có đặc tính linh động vừa có đặc tính định xứ. Tách kênh tương tác thành hai tham số trật tự: độ từ hóa $\mathbf{M}$ (đặc trưng cho sắt từ) và véctơ $\boldsymbol{\psi}$ (đặc trưng cho siêu dẫn spin-triplet).
  • Thiết lập phiếm hàm năng lượng cho siêu dẫn sắt từ: Dẫn xuất biểu thức giải tích của phiếm hàm Ginzburg-Landau hai thành phần từ tác dụng hiệu dụng vi 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})$$ trong đó các hệ số $a_s, b_s, a_f, b_f$ và hằng số ghép đôi $\gamma_0$ trong số hạng tương tác $f_I$ được tính toán trực tiếp từ các tích phân trên mặt Fermi và mật độ trạng thái của dải năng lượng điện tử $5f$.
  • Xây dựng giản đồ pha: Cực tiểu hóa phiếm hàm tự do theo hai tham số trật tự. Tính toán giản đồ pha trong mặt phẳng nhiệt độ – áp suất $(T, P)$ và mặt phẳng thông số rút gọn $(t, r)$.
  • Phân tích kết quả: Kết quả tính toán giải tích chỉ ra sự tồn tại của nghiệm đồng thời $\mathbf{M} \neq 0$ và $\boldsymbol{\psi} \neq 0$ trong một vùng áp suất xác định, giải thích sự xuất hiện của pha siêu dẫn nằm trọn vẹn bên trong pha sắt từ như thực nghiệm trên $UGe_2$ đã ghi nhận.

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

Luận án tự xác định các đóng góp mới về mặt khoa học lý thuyết và thực tiễn:

  1. Về mặt lý thuyết vi mô:

    • Xây dựng thành công quy trình giải tích vi mô nhất quán dựa trên tích phân phiếm hàm Grassmann và biến đổi Hubbard-Stratonovich để chuyển từ một Hamiltonian fermion nhiều hạt tổng quát sang phiếm hàm năng lượng tự do Ginzburg-Landau đa thành phần.
    • Thiết lập công thức liên hệ tường minh giữa các hệ số hiện tượng học của phiếm hàm Ginzburg-Landau (hệ số thế bậc hai, bậc bốn, hệ số gradient và hệ số tương tác chéo giữa các tham số trật tự) với các tham số cấu trúc vi mô như mật độ trạng thái ở mức Fermi $N(E_F)$, thế đẩy Coulomb tại chỗ $U$, tích phân nhảy nút $t$, và hằng số tương tác kết cặp.

  2. Về mặt thực tiễn mô hình hóa hệ vật liệu fermion nặng:

    • Cung cấp cơ sở vi mô giải thích sự đồng tồn tại của hai trật tự sắt từ và siêu dẫn spin-triplet trong hợp chất fermion nặng $UGe_2$.
    • Tái lập bằng con đường giải tích giản đồ pha định tính của $UGe_2$, làm rõ mối liên hệ giữa sự chuyển pha từ trạng thái phân cực mạnh ($FM_2$) sang phân cực yếu ($FM_1$) và sự xuất hiện cực đại của nhiệt độ chuyển pha siêu dẫn $T_{sc}$ dưới áp suất cao.

  3. Kiến nghị và giải pháp đề xuất:

    • Đề xuất tiếp tục mở rộng hình thức luận này để nghiên cứu các hiệu ứng chui hầm lượng tử trong cả vùng điện tích và vùng spin đối với các tiếp xúc dị thể chứa chất siêu dẫn sắt từ.
    • Ứng dụng khung giải tích đã thiết lập để khảo sát các vật liệu siêu dẫn sắt từ có cùng cơ chế ghép cặp như $URhGe$, $UCoGe$, và các hệ có sóng mật độ điện tích hoặc sóng mật độ spin đồng tồn tại với siêu dẫn.

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

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

  • Các tính toán trong luận án được thực hiện trong gần đúng trường trung bình vi mô dựa trên khai triển tiệm cận Ginzburg-Landau, do đó độ chính xác cao nhất đạt được ở vùng lân cận nhiệt độ chuyển pha; các hiệu ứng thăng giáng lượng tử mạnh tại điểm không độ tuyệt đối ($T = 0\text{ K}$) chưa được tính toán đầy đủ.
  • Mô hình dải năng lượng sử dụng trong các bước tính toán giải tích vẫn là mô hình đơn giản hóa, chưa tính toán trọn vẹn toàn bộ các ma trận tách mức trường tinh thể phức tạp và tương tác spin-quỹ đạo bất đẳng hướng của các mức năng lượng $5f$ trong mạng tinh thể trực giao thực tế.

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

  • Khảo sát các trạng thái kết cặp spin-triplet không unita ($\mathbf{q}(\mathbf{k}) \neq 0$), nơi tính đối xứng nghịch đảo thời gian bị phá vỡ hoàn toàn do từ trường nội.
  • Mở rộng phiếm hàm để tính toán cấu trúc không gian của mạng xoáy Abrikosov và hiệu ứng Meissner bất đẳng hướng khi có từ trường ngoài đặt dọc theo các trục tinh thể khác nhau.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo dành cho:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Vật lý lý thuyết và Vật lý toán: Có thể sử dụng Chương 2 cùng ba Phụ lục (A, B, C) như tài liệu hướng dẫn kỹ thuật tính toán về tích phân phiếm hàm trên đại số Grassmann, phép biến đổi Hubbard-Stratonovich và khai triển gradient toán tử hàm mũ.
  • Giảng viên và các nhà nghiên cứu vật lý chất rắn, vật lý vật chất ngưng tụ: Tham khảo Chương 3 và Chương 4 để nắm bắt phương pháp xây dựng phiếm hàm năng lượng Ginzburg-Landau đa tham số trật tự từ mô hình vi mô, phục vụ mô hình hóa các pha cạnh tranh hoặc đồng tồn tại trong các hệ điện tử tương quan mạnh.
  • Các nhóm nghiên cứu thực nghiệm vật liệu mới: Cung cấp bức tranh lý thuyết về vùng tham số áp suất và nhiệt độ thuận lợi cho việc tìm kiếm trạng thái siêu dẫn trong các hợp chất từ tính mới.

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

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

Trong lý thuyết BCS chuẩn, ngưng tụ siêu dẫn được hình thành bởi các cặp Cooper spin-singlet gồm hai điện tử có spin ngược chiều nhau ($S = 0$). Trong vật liệu sắt từ dưới nhiệt độ Curie, tương tác trao đổi nội định hướng các spin điện tử song song cùng chiều, tạo ra một trường cục bộ mạnh. Trường này phá vỡ cấu hình đối song của cặp Cooper singlet, làm giảm nhanh và triệt tiêu nhiệt độ chuyển pha siêu dẫn $T_{sc}$.

2. Trạng thái kết cặp siêu dẫn nào cho phép đồng tồn tại với trật tự sắt từ trong $UGe_2$?

Trạng thái kết cặp phù hợp là siêu dẫn spin-triplet ($S = 1$). Ở trạng thái này, hai điện tử trong cặp Cooper có spin song song cùng chiều (các trạng thái $|{\uparrow\uparrow}\rangle$ với $S_z = 1$, $|{\downarrow\downarrow}\rangle$ với $S_z = -1$, hoặc $(|{\uparrow\downarrow}\rangle + |{\downarrow\uparrow}\rangle)/\sqrt{2}$ với $S_z = 0$). Do spin của cặp không triệt tiêu, cặp Cooper triplet có thể tồn tại bền vững trong môi trường có sự phân cực từ hóa tự phát của vật liệu sắt từ.

3. Biến đổi Hubbard-Stratonovich giữ vai trò gì trong phương pháp tiếp cận của luận án?

Biến đổi Hubbard-Stratonovich là công cụ toán học cho phép chuyển đổi số hạng tương tác bốn hạt fermion phức tạp trong Hamiltonian vi mô thành các số hạng bậc hai ghép cặp với các trường phụ liên hợp. Nhờ đó, người nghiên cứu có thể tích phân chính xác các biến trường fermion (biến Grassmann), đưa bài toán nhiều hạt về một lý thuyết trường hiệu dụng chỉ phụ thuộc vào các trường tham số trật tự vĩ mô tương ứng với các kênh mật độ, kết cặp siêu dẫn và trao đổi từ.

4. Giản đồ pha thực nghiệm của $UGe_2$ dưới áp suất thể hiện những đặc trưng cấu trúc nào?

Giản đồ pha thực nghiệm của $UGe_2$ dưới áp suất chỉ ra rằng ở $P = 0$, vật liệu là chất sắt từ với $T_c = 52\text{ K}$, $\mu_s = 1,4\mu_B$. Dưới áp suất tăng dần, hệ trải qua hai pha sắt từ phân biệt: $FM_2$ (phân cực mạnh, $\mu = 1,5\mu_B$) và $FM_1$ (phân cực yếu, $\mu = 0,9\mu_B$), phân tách bởi ranh giới $T_x$. Pha siêu dẫn xuất hiện ngập sâu trong vùng sắt từ quanh $P \approx 1,2\text{ GPa}$ với $T_{sc}^{\max} \approx 0,7\text{ K}$ ngay tại điểm $T_x = T_{sc}$. Khi áp suất đạt giá trị tới hạn $P_c \approx 1,7\text{ GPa}$, cả trật tự sắt từ và siêu dẫn đều biến mất, hệ chuyển sang pha thuận từ ($PM$).

5. Tiêu chuẩn Stoner cho trật tự sắt từ trong mô hình Hubbard được phát biểu như thế nào?

Dựa trên sự cạnh tranh giữa động năng tiêu tốn để làm lệch mật độ điện tử giữa hai trạng thái spin và thế năng liên kết Coulomb được giảm thiểu khi các spin cùng chiều, tiêu chuẩn Stoner được phát biểu bằng bất đẳng thức: $$U N(E_F) > 1$$ trong đó $U$ là năng lượng đẩy tĩnh điện Coulomb tại chỗ và $N(E_F)$ là mật độ trạng thái điện tử tại mức năng lượng Fermi. Khi điều kiện này thỏa mãn, độ thay đổi năng lượng toàn phần $\delta E < 0$, hệ ưu tiên trạng thái phân cực spin tự phát (sắt từ).

Kết luận

Luận án của nghiên cứu sinh Nguyễn Văn Hinh đã giải quyết vấn đề xây dựng cơ sở vi mô cho lý thuyết Ginzburg-Landau mô tả sự đồng tồn tại giữa trật tự sắt từ và trật tự siêu dẫn trong hệ fermion nặng. 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 và kỹ thuật hàm Green, tác giả đã thiết lập thành công phiếm hàm năng lượng tự do đa thành phần với các hệ số được tính trực tiếp từ các tham số vi mô của vật liệu. Kết quả giải tích áp dụng cho hợp chất $UGe_2$ đã tái hiện thành công giản đồ pha nhiệt độ – áp suất, giải thích sự tồn tại của trạng thái siêu dẫn spin-triplet nằm sâu trong pha sắt từ. Công trình đóng góp một khung phương pháp luận hoàn chỉnh cho lĩnh vực nghiên cứu lý thuyết các hệ điện tử tương quan mạnh.