Tổng quan nghiên cứu

Khủng hoảng năng lượng và biến đổi khí hậu toàn cầu đang đặt ra yêu cầu cấp thiết về việc cải tiến các thiết bị nhiệt lạnh dân dụng và công nghiệp. Hiện nay, các hệ thống làm lạnh nén khí truyền thống chỉ đạt hiệu suất nhiệt cực đại khoảng 40%, đồng thời phát thải lượng lớn môi chất gây suy giảm tầng ozone và gia tăng hiệu ứng nhà kính. Trong bối cảnh đó, công nghệ làm lạnh từ tính dựa trên hiệu ứng từ nhiệt khổng lồ của các hợp chất liên kim loại nổi lên như một giải pháp thay thế đột phá, mang lại hiệu suất năng lượng vượt trội và thân thiện với môi trường.

Trong số các hệ vật liệu tiềm năng, hợp chất gốc La(Fe,Si)13 nhận được sự quan tâm đặc biệt nhờ giá thành vật liệu nền hợp lý và hiệu ứng từ nhiệt lớn ở vùng nhiệt độ phòng. Tuy nhiên, Lantan là kim loại đất hiếm đắt đỏ và ngày càng khan hiếm, đặt ra bài toán kinh tế then chốt là phải tối ưu hóa, cắt giảm tỷ lệ Lantan mà vẫn bảo toàn cấu trúc và các đặc tính truyền dẫn nhiệt điện. Luận văn thạc sĩ chuyên ngành Vật lý Nhiệt tại Đại học Quốc gia Hà Nội đã tập trung giải quyết bài toán này thông qua đề tài nghiên cứu cấu trúc tinh thể, các thông số nhiệt điện và ảnh hưởng của áp suất thủy tĩnh lên điện trở suất trong hệ hợp chất thiếu Lantan La1-δ(Fe0,845Si0,155)13 với nồng độ khuyết thiếu δ = 0,03; 0,06; 0,09.

Nghiên cứu được thực hiện trong giai đoạn 2014–2015 tại Phòng thí nghiệm Bộ môn Vật lý Nhiệt độ thấp. Đề tài đã làm sáng tỏ mối liên hệ mật thiết giữa vi cấu trúc lập phương loại NaZn13 với các thông số truyền tải nhiệt và điện, tiêu biểu là độ dẫn nhiệt đạt 9,6 W/m.K tại 300 K và độ giảm điện trở suất 2,2% dưới áp suất 16,5 kbar. Kết quả này cung cấp cơ sở dữ liệu thực nghiệm quan trọng cho việc thiết kế các chu trình làm lạnh từ tính tái sinh chủ động tiết kiệm chi phí.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết nhiệt điện cổ điển và hiện đại, bao gồm ba hiệu ứng tương hỗ cốt lõi: hiệu ứng Seebeck (chuyển đổi gradient nhiệt thành suất điện động, phát hiện năm 1821), hiệu ứng Peltier (tỏa hoặc thu nhiệt tại mối nối khi có dòng điện chạy qua, phát hiện năm 1834) và hiệu ứng Thomson (sự thay đổi nhiệt lượng dọc theo vật dẫn có dòng điện và chênh lệch nhiệt độ, năm 1851). Mối liên hệ nhiệt động lực học giữa các hệ số này được mô tả chặt chẽ qua hai hệ thức Kelvin: dS/dT = μ/T và Π = S.T.

Hiệu năng chuyển đổi năng lượng của vật liệu nhiệt điện được lượng hóa thông qua hệ số phẩm chất ZT theo công thức ZT = (S^2.T) / (ρ.κ), trong đó S là hệ số Seebeck, ρ là điện trở suất, κ là hệ số dẫn nhiệt toàn phần và T là nhiệt độ tuyệt đối. Hệ số dẫn nhiệt được cấu thành từ hai thành phần chính: độ dẫn nhiệt điện tử κ_el và độ dẫn nhiệt mạng tinh thể κ_ph (κ = κ_el + κ_ph). Trong đó, κ_el tuân theo định luật Wiedemann–Franz với hệ số Lorenz tiêu chuẩn L0 = 2,45 x 10^-8 W/K^2 ở nhiệt độ cao hơn nhiệt độ Debye.

Bên cạnh đó, mô hình trường phân tử và lý thuyết chuyển pha từ loại một dạng giả bền của electron lưu động (IEM) được áp dụng để phân tích hành vi từ tính. Hệ số trường phân tử n_Fe-Fe quyết định nhiệt độ chuyển pha Curie Tc thông qua tương tác giữa các moment từ nguyên tử Sắt (có mômen từ hiệu dụng 3,5 μB ở trạng thái thuận từ). Cấu trúc lập phương loại NaZn13 thuộc nhóm không gian Fm-3c với 8 đơn vị công thức trên một ô mạng cơ sở là nền tảng hình học giải thích cơ chế phân bố ion và các biến dạng mạng tinh thể.

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

Quy trình thực nghiệm được thiết kế đồng bộ từ khâu chế tạo phôi mẫu đến các phép đo đặc tính chuyên sâu:

  • Nguồn dữ liệu và chế tạo mẫu: Sử dụng các đơn chất có độ sạch cao gồm La (99,9%), Fe (99,9%) và Si (99,99%). Cỡ mẫu chế tạo danh định là 5 gam cho mỗi mẻ hợp phần La1-δ(Fe0,845Si0,155)13 (δ = 0,03; 0,06; 0,09). Lượng Lantan được cân dư 2% khối lượng danh định để bù trừ lượng hao hụt do bay hơi trong quá trình nấu hồ quang.
  • Quy trình xử lý nhiệt: Các hợp kim được nấu nóng chảy hồ quang trong môi trường khí Argon tinh khiết (hút chân không sơ bộ đạt 3 x 10^-2 Torr và chân không cao đạt 10^-5 Torr, sử dụng bẫy khí Titan hoạt hóa). Mẫu được lật đảo 3 lần để đảm bảo độ đồng nhất pha. Sau đó, phôi mẫu được bọc kín trong ống thạch anh hút chân không 10^-5 Torr, ủ nhiệt ở 1100°C liên tục trong 7 ngày (168 giờ) và tôi nhanh trong nước đá. Ngoài ra, kỹ thuật phun băng nguội nhanh đơn trục với tốc độ làm nguội đạt 10^5 đến 10^6 K/s cũng được khảo sát.
  • Phương pháp phân tích cấu trúc: Cấu trúc tinh thể được xác định bằng phương pháp nhiễu xạ bột tia X (XRD) trên thiết bị sử dụng bức xạ Cu-Kα có bước sóng λ = 1,540598 Å, kết hợp phần mềm tinh chỉnh Rietveld để xác định chính xác các chỉ số mặt mạng (hkl) và hằng số mạng a.
  • Phương pháp đo đặc tính điện và nhiệt: Điện trở suất, hệ số Seebeck và hệ số dẫn nhiệt trong dải nhiệt độ từ 4 K đến 300 K được đo đồng thời trên hệ thiết bị đo tính chất vật lý đa năng PPMS Evercool II (sử dụng module TTO và ACT).
  • Phép đo điện trở suất dưới áp suất cao: Thực hiện bằng phương pháp bốn mũi dò một chiều (DC) trên mẫu cắt hình hộp chữ nhật kích thước 1,1 mm x 1,0 mm x 3,0 mm, nối bốn dây dẫn vàng đường kính 0,050 mm bằng máy hàn siêu âm. Mẫu đặt trong buồng Teflon chứa môi trường truyền áp suất thủy tĩnh n-pentane và isoamyl alcohol tỷ lệ 1:1, tạo áp suất lên tới 1,2 GPa (12 kbar) qua hệ piston-cylinder từ 80 K đến 300 K.
  • Đo từ tính: Tiến hành trên từ kế SQUID siêu dẫn với độ nhạy đạt mức 10^-9 emu trong từ trường lên đến 9 T ở dải nhiệt độ 1,8 K đến 190 K.

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

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

  1. Tính đơn pha bền vững của cấu trúc mạng NaZn13: Giản đồ nhiễu xạ tia X tại nhiệt độ phòng khẳng định tất cả các mẫu La1-δ(Fe0,845Si0,155)13 với δ = 0,03; 0,06; 0,09 đều kết tinh đơn pha trong cấu trúc lập phương tâm mặt loại NaZn13 (nhóm không gian Fm-3c). Hằng số mạng thực tế dao động rất nhỏ từ 11,449 Å đến 11,458 Å, thể tích ô mạng cơ sở biến thiên trong khoảng 1501,033 Å^3 đến 1503,252 Å^3 (mức chênh lệch chỉ khoảng 0,5% giữa các mẫu).
  2. Nhiệt độ Curie tăng theo độ khuyết thiếu Lantan: Đường cong từ độ theo nhiệt độ M(T) dưới từ trường 1 kOe cho thấy nhiệt độ chuyển pha từ trạng thái sắt từ sang thuận từ (Tc) của hệ hợp chất thiếu Lantan tăng đều từ 237 K (ở δ = 0,03) lên 245 K (ở δ = 0,09). Mức nhiệt độ này cao hơn rõ rệt so với hợp chất đủ Lantan tương ứng (Tc = 232 K, tăng khoảng 5,6%).
  3. Đặc trưng dẫn điện và truyền nhiệt của mẫu La0,91(Fe0,845Si0,155)13: Phép đo trên hệ PPMS xác định điện trở suất mang đặc tính kim loại điển hình, đạt giá trị 134 μΩ.cm ở nhiệt độ phòng 300 K và điện trở suất dư 105,2 μΩ.cm ở vùng nhiệt độ cận 0 K. Hệ số dẫn nhiệt toàn phần tăng dần theo nhiệt độ và đạt cực đại 9,6 W/m.K tại 300 K. Trong đó, độ dẫn nhiệt điện tử κ_el chiếm tới 6,9 W/m.K (tương đương 72% tổng độ dẫn nhiệt), trong khi độ dẫn nhiệt mạng dao động ổn định quanh mức 2,5 W/m.K.
  4. Hệ số Seebeck âm và sự phụ thuộc tuyến tính của điện trở suất vào áp suất: Hệ số Seebeck mang giá trị âm trong toàn dải đo, khẳng định hạt tải điện đa số là điện tử (bán dẫn loại n), đạt giá trị cực tiểu -7,5 μV/K tại 200 K và đạt -5,5 μV/K ở 300 K, mang lại hệ số ZT cực đại 1,55 x 10^-3 tại 200 K. Khi gia tải áp suất thủy tĩnh từ 0 kbar lên 16,5 kbar ở nhiệt độ phòng, điện trở suất giảm tuyến tính khoảng 2,2%, chứng minh khả năng tinh chỉnh mật độ trạng thái điện tử bằng áp suất bên ngoài.

Thảo luận kết quả

Sự gia tăng nhiệt độ chuyển pha Curie Tc khi tăng hàm lượng thiếu Lantan (từ 237 K lên 245 K) được giải thích thỏa đáng thông qua mô hình trường phân tử. Khi khuyết thiếu ion La tại vị trí 8a, trật tự mạng tinh thể bị xáo trộn nhẹ làm tăng khoảng cách tương tác Fe-Fe giữa các vị trí 8b và 96i, kéo theo sự gia tăng của hệ số trường phân tử n_Fe-Fe và làm bền vững trạng thái sắt từ ở nhiệt độ cao hơn.

Đặc tính truyền nhiệt vượt trội với κ = 9,6 W/m.K tại nhiệt độ phòng giúp La0,91(Fe0,845Si0,155)13 thể hiện ưu thế lớn so với các vật liệu từ nhiệt truyền thống như Gd5Si2Ge2 hay hợp chất MnAs (vốn có độ dẫn nhiệt thấp hơn nhiều). Tỷ trọng đóng góp 72% của thành phần điện tử κ_el phù hợp hoàn hảo với định luật Wiedemann–Franz khi nhiệt độ Debye của mẫu được xác định là 377 K.

Dưới góc độ thực nghiệm, mối quan hệ giữa điện trở suất và áp suất thủy tĩnh thể hiện rõ hiệu ứng nén mạng tinh thể. Khi áp suất tăng đến 16,5 kbar, thể tích ô mạng V co lại làm mật độ electron tự do n = N/V tăng lên, đồng thời thời gian hồi phục τ tăng tỷ lệ thuận với nhiệt độ Debye. Sự kết hợp này giải thích cho xu hướng giảm tuyến tính 2,2% của điện trở suất.

Các dữ liệu thực nghiệm trong nghiên cứu được mô phỏng trực quan qua biểu đồ phụ thuộc nhiệt độ của điện trở suất ρ(T) và từ độ M(T), làm nổi bật bước nhảy pha từ học tại vùng lân cận nhiệt độ Curie và sự dịch chuyển đều đặn của các đường đẳng áp trong bảng số liệu đo đạc bốn mũi dò.

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

  1. Tối ưu hóa thành phần pha tạp phụ gia: Khuyến nghị các nhóm nghiên cứu vật liệu từ nhiệt thực hiện pha tạp bổ sung các nguyên tố kim loại chuyển tiếp như Coban (Co) hoặc Nhôm (Al) với tỷ lệ từ 2% đến 6% nguyên tử. Mục tiêu là nâng cao hệ số Seebeck S từ -5,5 μV/K lên trên mức -50 μV/K trong vòng 12 tháng tới.
  2. Áp dụng công nghệ cấu trúc nano để giảm dẫn nhiệt mạng: Viện Khoa học Vật liệu và các trung tâm nghiên cứu chuyên sâu cần triển khai kỹ thuật tán xạ phonon ở biên hạt thông qua phương pháp nghiền năng lượng cao kết hợp ép nóng đẳng tĩnh. Giải pháp này hướng tới việc giảm hệ số dẫn nhiệt mạng κ_ph xuống dưới 1,5 W/m.K trong lộ trình 18 tháng.
  3. Mở rộng quy mô công nghệ phun băng nguội nhanh: Các doanh nghiệp sản xuất vật liệu tiên tiến cần phối hợp với phòng thí nghiệm đại học để hoàn thiện dây chuyền tạo băng nguội nhanh đơn trục với tốc độ làm nguội 10^5 đến 10^6 K/s. Quy trình này cho phép sản xuất dải băng mỏng 20 đến 60 μm đồng nhất, giúp cắt giảm 30% chi phí nguyên liệu Lantan đắt đỏ trong vòng 24 tháng.
  4. Thử nghiệm tích hợp vào máy làm lạnh từ tái sinh chủ động (AMRR): Đề xuất các kỹ sư nhiệt lạnh và viện nghiên cứu ứng dụng chế tạo khối trao đổi nhiệt dạng xốp từ hợp chất La0,91(Fe0,845Si0,155)13, phấn đấu nâng hiệu suất làm lạnh thực tế của hệ thống AMRR lên trên 55% chu trình Carnot trong giai đoạn 3 năm tới.

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

  • Nhà khoa học và giảng viên ngành Vật lý chất rắn, Vật lý nhiệt độ thấp: Luận văn cung cấp hệ thống số liệu thực nghiệm chuẩn xác về cấu trúc mạng NaZn13, cơ chế chuyển pha từ giả bền và các công thức tính toán hằng số mạng aRwZ-Ds hỗ trợ công tác giảng dạy và nghiên cứu chuyên sâu.
  • Kỹ sư R&D trong ngành công nghệ nhiệt lạnh và năng lượng sạch: Cung cấp bức tranh toàn diện về các đặc tính dẫn nhiệt (9,6 W/m.K) và độ biến thiên điện trở suất theo áp suất, phục vụ trực tiếp cho việc tính toán thiết kế buồng trao đổi nhiệt từ tính không phát thải khí nhà kính.
  • Học viên cao học và nghiên cứu sinh ngành Khoa học Vật liệu: Tài liệu là hình mẫu tham khảo chuẩn mực về quy trình chế tạo mẫu bằng phương pháp nấu hồ quang kết hợp ủ nhiệt chân không 1100°C và kỹ thuật đo đạc nhiệt điện bốn mũi dò chính xác.
  • Chuyên gia quản lý công nghệ và hoạch định chiến lược vật liệu mới: Hỗ trợ việc thẩm định các dự án phát triển vật liệu thay thế đất hiếm, giúp định hình lộ trình chuyển đổi công nghệ làm lạnh hiệu suất cao nhằm tiết kiệm từ 15% đến 25% điện năng tiêu thụ quốc gia.

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

Tại sao việc thiếu hụt Lantan lại được tập trung nghiên cứu trong hợp chất La(Fe,Si)13?

Lantan là kim loại đất hiếm đắt đỏ và có trữ lượng khai thác hạn chế. Việc tạo ra độ khuyết thiếu Lantan từ 3% đến 9% nguyên tử (δ = 0,03 đến 0,09) giúp giảm đáng kể giá thành chế tạo vật liệu mà vẫn duy trì đơn pha lập phương NaZn13, thậm chí còn giúp nâng nhiệt độ Curie từ 232 K lên 245 K.

Cấu trúc tinh thể của hợp chất biến đổi như thế nào khi hàm lượng thiếu Lantan thay đổi?

Kết quả phân tích nhiễu xạ tia X cho thấy cấu trúc mạng tinh thể lập phương NaZn13 hoàn toàn ổn định trong dải khuyết thiếu δ từ 0,03 đến 0,09. Hằng số mạng a chỉ biến thiên cực nhỏ từ 11,449 Å đến 11,458 Å (chênh lệch dưới 0,5%), chứng minh độ bền pha vững chắc của hợp chất.

Ảnh hưởng của áp suất thủy tĩnh lên điện trở suất của vật liệu diễn ra theo cơ chế nào?

Khi gia tăng áp suất thủy tĩnh từ 0 kbar lên 16,5 kbar, điện trở suất tại nhiệt độ phòng giảm tuyến tính 2,2%. Áp suất nén thể tích ô mạng làm tăng mật độ electron tự do n và tăng thời gian hồi phục τ, từ đó cải thiện độ dẫn điện tổng thể của kim loại.

Đâu là yếu tố chính quyết định khả năng dẫn nhiệt của hợp chất La0,91(Fe0,845Si0,155)13?

Tại nhiệt độ phòng 300 K, độ dẫn nhiệt toàn phần đạt 9,6 W/m.K, trong đó độ dẫn nhiệt do điện tử đóng góp tới 6,9 W/m.K (chiếm 72%). Sự truyền nhiệt chủ yếu qua cơ chế electron tự do tuân theo định luật Wiedemann–Franz với nhiệt độ Debye đạt 377 K.

Quy trình nhiệt luyện nào đóng vai trò quyết định để thu được mẫu đơn pha?

Phôi mẫu sau khi nấu chảy hồ quang cần được đưa vào ampul thạch anh hút chân không cao 10^-5 Torr, sau đó tiến hành ủ nhiệt đồng pha hóa ở nhiệt độ cao 1100°C liên tục trong 7 ngày (168 giờ) rồi tôi nhanh trong nước đá để cố định cấu trúc lập phương mong muốn.

Kết luận

  • Luận văn đã tổng hợp và chế tạo thành công hệ hợp chất thiếu Lantan La1-δ(Fe0,845Si0,155)13 (δ = 0,03; 0,06; 0,09) bằng phương pháp nóng chảy hồ quang và ủ nhiệt chân không 1100°C trong 7 ngày.
  • Khẳng định tất cả các mẫu nghiên cứu đều kết tinh đơn pha trong cấu trúc lập phương kiểu NaZn13 (nhóm không gian Fm-3c) với hằng số mạng a ổn định quanh mức 11,449 Å – 11,458 Å.
  • Xác định nhiệt độ chuyển pha Curie Tc tăng từ 237 K lên 245 K khi nồng độ thiếu Lantan tăng, cao hơn mẫu đủ hợp thức (232 K) nhờ sự gia tăng tương tác trao đổi Fe-Fe.
  • Đo đạc toàn diện bốn thông số nhiệt điện của La0,91(Fe0,845Si0,155)13 tại 300 K: điện trở suất 134 μΩ.cm, độ dẫn nhiệt 9,6 W/m.K (dẫn nhiệt điện tử chiếm 72%), hệ số Seebeck -5,5 μV/K và hệ số phẩm chất ZT đạt cực đại 1,55 x 10^-3 ở 200 K.
  • Phát hiện hiệu ứng áp suất thủy tĩnh làm giảm tuyến tính 2,2% điện trở suất ở 16,5 kbar, mở ra hướng tinh chỉnh tính chất truyền dẫn hạt tải điện.

Đóng góp khoa học cốt lõi của công trình là chứng minh tính khả thi của việc cắt giảm hàm lượng đất hiếm Lantan mà vẫn bảo toàn các đặc tính vật lý then chốt của vật liệu từ nhiệt. Để đưa kết quả nghiên cứu vào ứng dụng thực tiễn trong 3–5 năm tới, các nhà khoa học và doanh nghiệp công nghệ cần đẩy mạnh hợp tác pha tạp vi lượng và thử nghiệm chế tạo linh kiện trao đổi nhiệt trong máy lạnh từ tính thế hệ mới. Hãy liên hệ và kết nối ngay với nhóm nghiên cứu để thúc đẩy các dự án chuyển giao công nghệ vật liệu xanh tiên tiến.