Tổng quan nghiên cứu

Hiệu ứng từ nhiệt là hiện tượng biến đổi nhiệt độ của vật liệu từ dưới tác dụng của từ trường ngoài trong điều kiện đoạn nhiệt, hoặc biến đổi entropy trong điều kiện đẳng nhiệt. Kể từ khi được Warburg phát hiện vào năm 1881 và bước đột phá công nghệ vào năm 1997 với máy làm lạnh từ kim loại Gadolinium đạt công suất 600 W hoạt động bền bỉ suốt 14 năm, công nghệ làm lạnh từ tính đã khẳng định vị thế vượt trội. Công nghệ này có thể đạt hiệu suất tới 60% giới hạn lý thuyết Carnot, cao hơn hẳn mức 40% của hệ thống nén khí truyền thống, đồng thời loại bỏ hoàn toàn môi chất làm lạnh gây hiệu ứng nhà kính.

Vấn đề nghiên cứu trọng tâm hiện nay là tìm kiếm các hệ vật liệu từ nhiệt có hiệu ứng từ nhiệt khổng lồ, giá thành thấp, không độc hại và có nhiệt độ chuyển pha từ gần nhiệt độ phòng trong điều kiện từ trường thấp. Mục tiêu cụ thể của luận văn là nghiên cứu quy trình công nghệ chế tạo, khảo sát cấu trúc tinh thể loại NaZn13 và phân tích chi tiết các tính chất từ của hệ hợp chất liên kim loại La(Fe1-xSix)13 với các tỷ lệ nồng độ pha tạp Silic x biến thiên từ 0,12 đến 0,21.

Phạm vi nghiên cứu được thực hiện tại Bộ môn Vật lý Nhiệt độ thấp, Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội, với các phép đo thực nghiệm từ tính và cấu trúc trong dải nhiệt độ từ 1,8 K đến 300 K và từ trường cực đại lên đến 7 T. Ý nghĩa khoa học và thực tiễn của công trình thể hiện qua việc xác lập mẫu hợp chất tối ưu đạt biến thiên entropy từ cực đại lên tới 31,5 J/kg.K và hiệu suất làm lạnh tương đối 229,95 J/kg ở từ trường 7 T, đóng góp trực tiếp vào mục tiêu cắt giảm 20% đến 30% năng lượng tiêu thụ cho các thiết bị nhiệt lạnh thế hệ mới.

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

Khung lý thuyết áp dụng

Khung lý thuyết của đề tài được xây dựng trên nền tảng nhiệt động lực học hiệu ứng từ nhiệt và lý thuyết chuyển pha từ tính Landau kết hợp mô hình điện tử linh động. Mối liên hệ nhiệt động lực học giữa các thông số trạng thái được xác định thông qua hàm năng lượng tự do Gibbs và hệ thức Maxwell:

$$\Delta S_{mag} = \int_{0}^{H} \left(\frac{\partial M}{\partial T}\right)_{H, P} dH$$

Lý thuyết chuyển pha Landau và mô hình Wohlfarth - Rhodes được áp dụng để giải thích hiện tượng chuyển pha từ giả bền điện tử linh động (IEM) bắt nguồn từ cấu trúc dải năng lượng của các điện tử 3d trong Fe. Khi hệ số bậc bốn trong khai triển Landau nhận giá trị âm, hệ vật liệu xảy ra chuyển pha loại một dưới tác dụng của từ trường ngoài.

Các khái niệm chính được vận dụng bao gồm: cấu trúc lập phương kiểu NaZn13 thuộc nhóm không gian Fm3c chứa 8 đơn vị công thức trong một ô nguyên tố với tỷ lệ vị trí phối trí CoI và CoII là 1:12; nhiệt độ chuyển pha Curie (TC) xác định điểm chuyển tiếp giữa trạng thái sắt từ và thuận từ; biến thiên entropy từ (ΔSmag) và biến thiên nhiệt độ đoạn nhiệt (ΔTad); hiệu suất làm lạnh tương đối (RCP) được tính bằng tích số giữa biến thiên entropy cực đại và nửa độ rộng cực đại của đường cong nhiệt (TFWHM).

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

Nghiên cứu khảo sát cỡ mẫu gồm 5 thành phần hợp chất danh định La(Fe1-xSix)13 với các giá trị x = 0,12; 0,14; 0,15; 0,18 và 0,21. Phương pháp chọn mẫu dựa trên tính toán tỷ lượng hợp kim chính xác từ các đơn chất ban đầu có độ sạch cao 99,99% đối với La, Fe và 99,999% đối với Si. Mỗi thỏi mẫu có khối lượng 5,0 gam, trong đó thành phần đất hiếm La được tính dư thêm 2% khối lượng nhằm bù trừ lượng hao hụt do bay hơi ở nhiệt độ cao.

Quy trình thực nghiệm được triển khai theo timeline chuẩn hóa gồm ba giai đoạn chính:

  • Chế tạo hợp kim bằng phương pháp nóng chảy hồ quang trong môi trường khí trơ Argon áp suất 10 Torr, buồng nấu được hút chân không sâu tới 10^-5 Torr và khử oxy dư bằng việc nấu chảy viên Titan mồi. Mẫu được lật đảo 3 lần để đảm bảo độ đồng nhất pha.
  • Xử lý nhiệt luyện bằng cách hàn kín mẫu trong ống thạch anh chân không 10^-5 Torr, ủ nhiệt đẳng nhiệt ở 1100°C trong 7 ngày (1 tuần) hoặc 14 ngày (2 tuần), sau đó tôi nhanh trong nước đá để đóng băng pha cấu trúc 1:13.
  • Lý do lựa chọn phương pháp phân tích: Sử dụng phương pháp nhiễu xạ bột tia X (XRD) trên máy RINT-2000 với bức xạ Cu-Kalpha bước sóng 1,540598 Å để xác định độ đơn pha và hằng số mạng của các hạt bột kích thước 50-100 μm. Đo đạc từ tính bằng Giao thoa kế lượng tử siêu dẫn (SQUID) do thiết bị này sở hữu độ nhạy từ cực cao tốt hơn 10^-9 emu, cho phép xác định chính xác từ độ bão hòa, nhiệt độ chuyển pha và biến thiên entropy từ trong dải từ trường từ 0 đến 7 T.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thực nghiệm đã đem lại bốn phát hiện khoa học quan trọng có giá trị định lượng cao:

  • Thứ nhất, điều kiện ủ nhiệt 1100°C trong 7 ngày là yếu tố quyết định để hình thành pha 1:13 đơn pha đối với các hợp chất có nồng độ Si thấp (x = 0,12; 0,14; 0,15). Ở trạng thái đúc ban đầu chưa xử lý nhiệt, mẫu tồn tại chủ yếu ở pha alpha-Fe, nhưng sau khi ủ nhiệt, pha alpha-Fe phân rã hoàn toàn và chuyển thành pha 1:13 với lượng pha tạp còn lại dưới 3%. Với nồng độ Si cao (x >= 0,18), pha 1:13 chiếm ưu thế ngay ở trạng thái đúc ban đầu.
  • Thứ hai, hằng số mạng tinh thể giảm tuyến tính khi nồng độ Si tăng trong dải lập phương: từ 11,5902 Å (ở x = 0,12) giảm xuống 11,5346 Å (ở x = 0,18), tương đương mức giảm 0,48%. Khi nồng độ Si đạt x = 0,21, cấu trúc tinh thể chuyển đổi từ hệ lập phương NaZn13 sang hệ tứ diện với hằng số mạng a = b = 7,9316 Å và c = 11,7783 Å.
  • Thứ ba, nhiệt độ chuyển pha Curie tăng tuyến tính 28,7% theo nồng độ Si, từ 202 K (với x = 0,12) lên 260 K (với x = 0,21). Ngược lại, mômen từ bão hòa tại 1,8 K giảm tuyến tính 24,5%, từ 2,1907 μB/Fe nguyên tử (với x = 0,12) xuống 1,6542 μB/Fe nguyên tử (với x = 0,21).
  • Thứ tư, hợp chất La(Fe0,88Si0,12)13 thể hiện hiệu ứng từ nhiệt khổng lồ với biến thiên entropy từ cực đại đạt 31,5 J/kg.K dưới biến thiên từ trường 7 T, đem lại hiệu suất làm lạnh tương đối RCP đạt 229,95 J/kg ở 7 T và 103,5 J/kg ở 3 T.

Thảo luận kết quả

Cơ chế vật lý của sự suy giảm hằng số mạng khi nồng độ Si tăng được giải thích rõ ràng bởi bán kính ion của Si (0,11 μm) nhỏ hơn bán kính ion của Fe (0,13 μm), dẫn đến sự co rút thể tích ô mạng tinh thể. Sự co mạng này làm biến đổi khoảng cách giữa các nguyên tử trung gian và ion Fe ở các vị trí phối trí 8b và 96i, làm thay đổi cấu trúc dải điện tử và trực tiếp làm tăng nhiệt độ chuyển pha Curie TC. Ngược lại, sự suy giảm mômen từ bão hòa xuất phát từ việc mật độ nguyên tử sắt từ Fe trên mỗi đơn vị công thức bị suy giảm khi Si thế chỗ.

Hiện tượng chuyển pha từ giả bền điện tử linh động (IEM) được khẳng định rõ nét ở mẫu La(Fe0,88Si0,12)13 ngay phía trên nhiệt độ Curie 202 K. Bản chất của chuyển pha loại một này là do sự tách dải năng lượng của các điện tử 3d dưới tác dụng của từ trường ngoài, tạo ra sự biến thiên từ độ đột ngột.

Toàn bộ dữ liệu thực nghiệm được trình bày nhất quán thông qua các bảng số liệu hằng số mạng, biểu đồ phổ nhiễu xạ tia X, đồ thị từ hóa phụ thuộc nhiệt độ M(T) và họ đường cong Arrott plots. Đồ thị Arrott của mẫu x = 0,12 thể hiện rõ độ dốc âm và điểm uốn chữ S đặc trưng cho chuyển pha loại một, trong khi ở các mẫu nồng độ Si cao (x = 0,15 và 0,18), đường cong Arrott có độ dốc dương đồng đều của chuyển pha loại hai. Đồ thị biến thiên entropy từ ΔSmag theo nhiệt độ thể hiện dạng đỉnh nhọn đối xứng quanh nhiệt độ TC, khẳng định tiềm năng ứng dụng vượt bậc của vật liệu.

Đề xuất và khuyến nghị

Dựa trên các kết quả đạt được, luận văn đưa ra bốn nhóm giải pháp hành động cụ thể:

  • Chuẩn hóa chế độ nhiệt luyện công nghiệp: Các cơ sở nghiên cứu và xưởng luyện kim đặc biệt cần áp dụng quy trình ủ nhiệt mẫu La(Fe,Si)13 ở nhiệt độ chuẩn 1100°C trong thời gian đúng 7 ngày (168 giờ) dưới môi trường chân không 10^-5 Torr đối với các tỷ lệ Si từ 0,12 đến 0,15. Giải pháp này giúp loại bỏ trên 97% pha tạp alpha-Fe, ổn định pha 1:13 đơn pha với mục tiêu hoàn thiện quy trình trong 6 tháng đầu năm 2025 do đội ngũ kỹ thuật vật liệu thực hiện.
  • Mở rộng nghiên cứu pha tạp nguyên tố thứ tư để nâng nhiệt độ làm việc: Các nhóm nghiên cứu vật lý chất rắn cần tiến hành pha tạp bổ sung các nguyên tố như Co, Ni hoặc chèn nguyên tử xen kẽ H (hydrô hóa) vào mạng La(Fe1-xSix)13 nhằm đẩy nhiệt độ chuyển pha Curie từ dải 202-260 K lên dải nhiệt độ phòng 295-310 K, đồng thời duy trì độ biến thiên entropy từ trên 15 J/kg.K tại từ trường 2 T trong giai đoạn 2025-2026.
  • Hoàn thiện hệ thống đo kiểm và chuẩn hóa đánh giá hiệu ứng từ nhiệt: Các phòng thí nghiệm trọng điểm quốc gia phối hợp với viện nghiên cứu vật lý xây dựng quy trình đo từ nhiệt gián tiếp qua từ kế SQUID kết hợp phép đo nhiệt dung trực tiếp, kiểm soát sai số đo đạc dưới 2% trước quý 4 năm 2025.
  • Thiết kế và chế tạo mô hình thử nghiệm máy làm lạnh từ: Các kỹ sư cơ điện lạnh và doanh nghiệp công nghệ năng lượng triển khai thiết kế khối trao đổi nhiệt từ tính sử dụng vật liệu La(Fe,Si)13 dạng tấm mỏng vi kênh, hướng tới mục tiêu chế tạo máy làm lạnh từ công suất 500 W với hệ số hiệu quả COP tăng 25% trong giai đoạn 2026-2028.

Đối tượng nên tham khảo luận văn

  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Vật lý chất rắn, Vật lý nhiệt: Tiếp cận tài liệu chuyên sâu về lý thuyết chuyển pha Landau, cơ chế chuyển pha từ giả bền điện tử linh động (IEM) và phương pháp xác định hiệu ứng từ nhiệt qua từ kế SQUID trong dải đo 1,8 K đến 300 K.
  • Kỹ sư nghiên cứu và phát triển (R&D) ngành nhiệt lạnh và thiết bị gia dụng: Vận dụng dữ liệu thực nghiệm về vật liệu La(Fe,Si)13 để thiết kế các chu trình làm lạnh từ tính không phát thải khí nhà kính, khai thác hiệu suất làm việc 60% của chu trình từ nhiệt thay thế máy nén gas.
  • Kỹ sư công nghệ luyện kim và vật liệu từ tính: Ứng dụng quy trình công nghệ nấu nóng chảy hồ quang và chế độ xử lý nhiệt 1100°C trong 7 ngày để tối ưu hóa việc chế tạo hợp kim đơn pha cấu trúc NaZn13 với độ sạch cao 99,99%.
  • Chuyên gia tư vấn chính sách năng lượng và môi trường: Sử dụng các số liệu khoa học của luận văn để xây dựng báo cáo đánh giá tiềm năng tiết kiệm 20% đến 30% điện năng tiêu thụ quốc gia thông qua việc chuyển đổi từ công nghệ nén khí sang công nghệ làm lạnh từ tính xanh.

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

Cấu trúc tinh thể NaZn13 đóng vai trò gì đối với hiệu ứng từ nhiệt của hợp chất La(Fe,Si)13?

Cấu trúc lập phương NaZn13 thuộc nhóm không gian Fm3c chứa 8 đơn vị công thức với các vị trí phối trí 8b và 96i tạo nên mật độ nguyên tử sắt từ rất cao. Mạng tinh thể này là điều kiện nền tảng để tạo ra mômen từ lớn và làm xuất hiện hiện tượng chuyển pha từ giả bền, mang lại hiệu ứng từ nhiệt khổng lồ.

Vì sao cần phải ủ nhiệt ở 1100°C trong 7 ngày đối với các hợp chất có nồng độ Si thấp?

Với các mẫu x từ 0,12 đến 0,15, trạng thái sau khi đúc chứa chủ yếu là pha alpha-Fe không mong muốn. Quá trình ủ nhiệt ở 1100°C trong 7 ngày cung cấp đủ động học khuếch tán để hòa tan hoàn toàn pha alpha-Fe, chuyển hóa trên 97% thành phần sang pha đơn tinh thể 1:13 đồng nhất.

Nồng độ Silic ảnh hưởng như thế nào đến nhiệt độ Curie và mômen từ bão hòa?

Khi nồng độ Si tăng từ 0,12 lên 0,21, nhiệt độ chuyển pha Curie tăng tuyến tính từ 202 K lên 260 K do mạng tinh thể bị co rút bởi bán kính ion Si nhỏ (0,11 μm). Ngược lại, mômen từ bão hòa giảm tuyến tính từ 2,1907 xuống 1,6542 μB/Fe nguyên tử do mật độ Fe giảm.

Chuyển pha từ giả bền điện tử linh động mang lại lợi ích gì cho vật liệu từ nhiệt?

Chuyển pha từ giả bền điện tử linh động là chuyển pha loại một kèm theo sự biến đổi đột ngột của từ độ dưới tác dụng từ trường. Hiện tượng này tạo ra bước nhảy entropy từ cực lớn, giúp mẫu La(Fe0,88Si0,12)13 đạt biến thiên entropy từ lên tới 31,5 J/kg.K ở từ trường 7 T.

Biến thiên entropy từ trong luận văn được xác định bằng phương pháp thực nghiệm nào?

Biến thiên entropy từ được xác định gián tiếp bằng cách đo hệ đường cong từ hóa đẳng nhiệt M(H) trên từ kế SQUID trong dải nhiệt độ lân cận điểm chuyển pha. Diện tích giới hạn giữa hai đường từ hóa liên tiếp được tích phân theo hệ thức Maxwell để tính ra giá trị biến thiên entropy chính xác.

Kết luận

  • Luận văn đã xác định thành công quy trình công nghệ tối ưu chế tạo hệ vật liệu từ nhiệt La(Fe1-xSix)13 bằng phương pháp nóng chảy hồ quang kết hợp ủ nhiệt tại 1100°C trong 7 ngày, đảm bảo cấu trúc đơn pha 1:13 đạt độ tinh khiết cao.
  • Làm sáng tỏ quy luật chuyển dịch cấu trúc tinh thể từ dạng lập phương NaZn13 với hằng số mạng 11,5902 Å sang dạng tứ diện có thông số a = 7,9316 Å và c = 11,7783 Å khi nồng độ Si tăng đến x = 0,21.
  • Chứng minh quy luật nhiệt độ Curie tăng tuyến tính từ 202 K lên 260 K đi kèm sự suy giảm mômen từ bão hòa từ 2,1907 xuống 1,6542 μB/Fe nguyên tử do hiệu ứng co mạng tinh thể.
  • Phát hiện chuyển pha từ giả bền điện tử linh động loại một trên hợp chất La(Fe0,88Si0,12)13, đem lại biến thiên entropy từ kỷ lục 31,5 J/kg.K dưới từ trường 7 T.
  • Xác lập hiệu suất làm lạnh tương đối RCP đạt 229,95 J/kg, khẳng định tiềm năng ứng dụng to lớn trong kỹ thuật làm lạnh từ tính không phát thải.
  • Kế hoạch nghiên cứu tiếp theo trong 6 tháng tới sẽ tập trung vào kỹ thuật pha tạp Co và xử lý hydrô hóa để đưa nhiệt độ Curie lên vùng nhiệt độ phòng 300 K.
  • Các viện nghiên cứu và doanh nghiệp công nghệ năng lượng quan tâm xin vui lòng liên hệ hợp tác nghiên cứu phát triển mẫu thử nghiệm máy làm lạnh từ tính ứng dụng thực tế.