Tổng quan nghiên cứu

Kể từ cột mốc phát hiện hiệu ứng từ nhiệt khổng lồ trên hợp chất Gd5Ge2Si2 vào năm 1997, lĩnh vực nghiên cứu vật liệu làm lạnh từ trạng thái rắn đã ghi nhận tốc độ tăng trưởng học thuật vượt bậc với hơn 300% số lượng công bố khoa học trong hai thập kỷ qua. Công nghệ làm lạnh từ dựa trên hiệu ứng từ nhiệt đang trở thành giải pháp đột phá thay thế công nghệ nén khí gas truyền thống, giúp loại bỏ hoàn toàn các loại môi chất làm lạnh gây hiệu ứng nhà kính và tiết kiệm từ 20% đến 30% mức tiêu thụ điện năng. Tuy nhiên, thách thức lớn nhất hiện nay là phần lớn các vật liệu từ nhiệt có hiệu ứng khổng lồ thường gắn liền với quá trình chuyển pha từ loại một, gây ra hiện tượng trễ từ và trễ nhiệt lớn, làm tiêu hao đáng kể công suất làm lạnh tương đối.

Nhằm giải quyết bài toán trên, đề tài tập trung nghiên cứu sâu về cấu trúc tinh thể, tính chất từ nhiệt và cấu trúc từ vi mô của hệ vật liệu liên kim đất hiếm Ho2In và Tb2In trong vùng nhiệt độ thấp dưới 200 K. Mục tiêu cụ thể là xác định bản chất các chuyển pha từ, hằng số mạng, mômen từ nguyên tử và đánh giá tiềm năng ứng dụng của chúng trong các tầng làm lạnh nhiệt độ thấp. Nghiên cứu được thực hiện tại Bộ môn Vật lý Nhiệt độ thấp thuộc Trường Đại học Khoa học Tự nhiên Hà Nội phối hợp cùng Viện Khoa học và Công nghệ Tiên tiến Nhật Bản và Lò phản ứng Nghiên cứu Đại học Kyoto trong giai đoạn từ năm 2010 đến 2014. Kết quả của luận văn mang ý nghĩa khoa học lớn khi cung cấp bức tranh chi tiết về cấu trúc từ vi mô với mômen từ đạt 9,0 Bohr magneton tại 2 K và định hướng góc quay 41 độ so với trục c, đặt tiền đề tối ưu hóa công suất làm lạnh lên tới 660 J/kg trong dải nhiệt độ từ 1,8 K đến 170 K.

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 động học từ tính kết hợp mô hình trường tinh thể và tương tác trao đổi gián tiếp RKKY giữa các ion đất hiếm qua trung gian điện tử dẫn. Bốn khái niệm then chốt chi phối toàn bộ phân tích bao gồm:

Hiệu ứng từ nhiệt: Hiện tượng biến thiên nhiệt độ đoạn nhiệt hoặc biến thiên entropy từ của vật liệu từ tính dưới tác dụng của từ trường ngoài thay đổi từ 0 đến 70 kOe, được mô tả định lượng thông qua hệ thức nhiệt động Maxwell.

Chuyển pha từ loại một và loại hai: Phân loại dựa trên tính liên tục của đạo hàm hàm nhiệt động. Chuyển pha loại một thể hiện bước nhảy entropy đột ngột và dị thường dạng hàm số lambda trên đường nhiệt dung, trong khi chuyển pha loại hai có biến thiên trơn tru, giúp loại bỏ tổn hao trễ năng lượng.

Dị hướng từ trường tinh thể: Tương tác tĩnh điện giữa các đám mây điện tử lớp 4f của ion đất hiếm với điện trường không đồng nhất do các ion lân cận tạo ra, được mô tả toán học qua toán tử Stevens bậc 2, bậc 4 và bậc 6, quyết định trục từ hóa dễ của tinh thể.

Công suất làm lạnh tương đối: Đại lượng đo lường lượng nhiệt năng mà một chu trình làm lạnh từ có thể vận chuyển, được tính bằng tích phân của biến thiên entropy từ trên độ rộng nửa cực đại của đường cong nhiệt độ.

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

Nghiên cứu sử dụng quy trình thực nghiệm đa bước từ chế tạo mẫu đến khảo sát chuyên sâu:

Chuẩn bị mẫu: Mẫu hợp chất đa tinh thể Ho2In và Tb2In với cỡ mẫu 10 gam được chế tạo bằng phương pháp nóng chảy hồ quang điện cực vonfram không tiêu hao trên đế đồng làm mát bằng nước. Phương pháp chọn mẫu là nấu hợp kim từ các kim loại thành phần có độ sạch cao trên 99,9%, sắp xếp theo thứ tự nhiệt độ nóng chảy tăng dần. Quá trình nấu diễn ra trong môi trường khí argon bảo vệ ở áp suất 0,6 đến 0,8 bar sau khi buồng mẫu được hút chân không phân tử đạt áp suất 10^-5 Torr và nấu chảy viên titan mồi để hấp thụ triệt để oxy dư. Các thỏi mẫu được đảo mặt 3 lần để đảm bảo độ đồng nhất tuyệt đối, sau đó bọc trong ống thạch anh chân không và ủ nhiệt ở 1073 K trong 5 ngày (120 giờ).

Phương pháp phân tích cấu trúc: Phép đo nhiễu xạ bột tia X thực hiện trên thiết bị Rigaku RINT-2000 với bức xạ Cu-Kalpha, góc quét 2-theta từ 10 độ đến 100 độ, bước nhảy 0,01 độ trên khối lượng mẫu 0,1 gam nhằm kiểm tra độ đơn pha và xác định thông số mạng. Phép đo nhiễu xạ bột nơtron được tiến hành tại Lò phản ứng Đại học Kyoto với bước sóng 1,6430 Angstrom trên lượng mẫu 10 gam ở dải nhiệt độ 1,8 K đến 250 K. Lý do lựa chọn nhiễu xạ nơtron là vì nơtron mang mômen từ lưỡng cực, cho phép tương tác trực tiếp với mômen từ nguyên tử để phân tách chính xác trật tự từ tính vi mô mà tia X không thể quan sát được. Toàn bộ cấu trúc được làm mịn bằng phương pháp Rietveld qua chương trình Rietan-2000.

Phương pháp khảo sát từ tính và nhiệt dung: Sử dụng từ kế lượng tử SQUID với từ trường siêu dẫn lên đến 70 kOe và hệ thống đo tính chất vật lý PPMS của Quantum Design trên các mẫu kích thước chuẩn 2,0 x 2,0 x 1,0 mm3 với khối lượng khoảng 20 mg trong dải nhiệt độ từ 1,8 K đến 300 K.

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

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

Phân tích nhiễu xạ tia X tại nhiệt độ phòng khẳng định cả hai mẫu Ho2In và Tb2In đều đạt trạng thái đơn pha hoàn toàn, kết tinh ở cấu trúc lục giác kiểu Ni2In thuộc nhóm không gian P63/mmc với hai đơn vị công thức trên mỗi ô cơ sở. Các nguyên tử đất hiếm phân bố tại hai vị trí tinh thể học bất tương đương là 2a và 2d, trong khi nguyên tử In chiếm vị trí 2c. Hằng số mạng của Ho2In được xác định chính xác là a = 5,3216 Angstrom và c = 6,6667 Angstrom với hệ số tin cậy Rwp = 5,54% và hệ số phù hợp S = 1,45.

Khảo sát từ tính của Ho2In phát hiện sự tồn tại của hai quá trình chuyển pha từ liên tiếp: chuyển pha từ trạng thái thuận từ sang sắt từ tại nhiệt độ Curie TC = 87 K và một chuyển pha trật tự từ thứ hai tại Tt = 36 K khi đo ở từ trường 1,0 kOe. Trong pha thuận từ, nghịch đảo độ cảm từ tuân theo định luật Curie-Weiss với nhiệt độ Curie thuận từ đạt 90 K và mômen từ hiệu dụng thực nghiệm đạt 10,61 Bohr magneton, hoàn toàn phù hợp với giá trị lý thuyết của ion tự do Ho3+ là 10,60 Bohr magneton.

Phổ nhiễu xạ nơtron vi mô của Ho2In tại nhiệt độ siêu thấp 2 K đã xác định được mômen từ của ion Ho tại vị trí 2a là 9,0 Bohr magneton và tại vị trí 2d là 8,8 Bohr magneton, tạo một góc nghiêng 41 độ so với trục đối xứng c. Khi nhiệt độ tăng lên 45 K, cường độ đỉnh nhiễu xạ từ tính (102) suy giảm rõ rệt, mômen từ tương ứng tại hai vị trí giảm xuống còn 7,4 Bohr magneton và 7,0 Bohr magneton, đồng thời các mômen từ có xu hướng quay dần về phía trục c.

Đối với hợp chất Tb2In, kết quả ghi nhận quá trình chuyển pha thuận từ sang sắt từ tại TC = 160 K và chuyển pha từ sắt từ sang phản sắt từ tại nhiệt độ Neel TN = 50 K. Dữ liệu nhiệt dung cho thấy bước nhảy nhiệt dung tại 160 K đạt 18 J/mol.K với dạng đường cong đối xứng đặc trưng cho chuyển pha loại hai, mang lại tiềm năng công suất làm lạnh rất cao.

Thảo luận kết quả

Các kết quả thực nghiệm được biểu diễn trực quan thông qua đồ thị phụ thuộc nhiệt độ của từ độ M(T) và nghịch đảo độ cảm từ trong dải từ trường từ 0,1 kOe đến 1,0 kOe, kết hợp hệ đường cong từ hóa đẳng nhiệt M(H) trong dải nhiệt độ 1,8 K đến 150 K dưới từ trường quét lên đến 70 kOe. Đồ thị nhiệt dung C(T) minh họa rõ nét các dị thường nhiệt động học, trong đó mẫu Ho2In xuất hiện bước nhảy nhiệt dung lớn Delta C = 20 J/mol.K tại 85 K dạng lambda và bước nhảy nhỏ Delta C = 2 J/mol.K tại 36 K, minh chứng cho bản chất chuyển pha loại một. Bảng tổng hợp dữ liệu cấu trúc và từ tính chỉ ra sự tương thích cao giữa mô hình làm mịn Rietveld và thực nghiệm.

Từ độ bão hòa của Ho2In đo được tại 1,8 K dưới từ trường 70 kOe đạt 8,41 Bohr magneton, thấp hơn giá trị lý thuyết của ion tự do Ho3+ là 10,0 Bohr magneton khoảng 15,9%. Sự suy giảm này được giải thích là do tác động mạnh mẽ của trường tinh thể nội tại lên lớp vỏ điện tử 4f, làm suy biến và đóng băng một phần mômen quỹ đạo. Khi so sánh với các vật liệu cùng hệ như Dy2In (TC = 125 K, biến thiên entropy từ 7,1 J/kg.K, công suất làm lạnh 473 J/kg) và Er2In (TC = 46 K, biến thiên entropy từ 16 J/kg.K, công suất làm lạnh 490 J/kg), hợp chất Ho2In (TC = 85 K, biến thiên entropy từ 11,2 J/kg.K, công suất làm lạnh 360 J/kg) và Tb2In (TC = 165 K, biến thiên entropy từ 6,6 J/kg.K, công suất làm lạnh đạt mức kỷ lục 660 J/kg dưới từ trường 5 Tesla) chứng minh ưu thế vượt trội về khoảng nhiệt độ làm việc rộng và công suất tải nhiệt lớn.

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

Tối ưu hóa công nghệ luyện kim và chế độ ủ nhiệt: Các phòng thí nghiệm vật liệu từ và cơ sở luyện kim chuyên sâu cần áp dụng quy trình kiểm soát chân không cao dưới 10^-5 Torr kết hợp hệ thống tinh luyện hồ quang đa điện cực. Thời gian thực hiện trong vòng 6 đến 12 tháng, đặt mục tiêu nâng cao độ đồng nhất thành phần pha lên trên 99,8% và giảm hệ số tin cậy Rietveld Rwp xuống dưới mức 3,5% đối với các mẻ nấu quy mô lớn.

Thực hiện pha tạp nguyên tố để dịch chuyển vùng nhiệt độ làm việc: Các nhóm nghiên cứu khoa học vật liệu tại các viện nghiên cứu cần triển khai thử nghiệm pha tạp một phần các nguyên tố kim loại chuyển tiếp như sắt, coban hoặc đất hiếm nặng vào vị trí của Ho và Tb. Lộ trình thực hiện từ 12 đến 24 tháng, nhằm mục tiêu đưa nhiệt độ chuyển pha Curie TC tịnh tiến lên vùng 250 K đến 300 K, đồng thời duy trì công suất làm lạnh tương đối RCP trên mức 700 J/kg dưới từ trường biến thiên 5 Tesla.

Chế tạo cấu trúc vật liệu dạng màng mỏng và bột nano: Các trung tâm công nghệ cao cần nghiên cứu ứng dụng phương pháp nghiền năng lượng cao hoặc phún xạ catốt để tạo cấu trúc hạt kích thước nano. Kế hoạch triển khai từ 2 đến 3 năm, hướng tới mục tiêu tăng diện tích tiếp xúc truyền nhiệt lên 45% và triệt tiêu hoàn toàn hiện tượng trễ từ xuống dưới 2%, giúp cải thiện tốc độ chu trình làm lạnh.

Thiết kế mô hình thử nghiệm máy làm lạnh từ nhiệt độ thấp: Các đơn vị kỹ thuật nhiệt lạnh và năng lượng cần phối hợp thiết kế cụm buồng trao đổi nhiệt sử dụng vật liệu Ho2In và Tb2In cho các tầng làm lạnh khí hóa lỏng nitơ và heli. Dự án kéo dài từ 3 đến 5 năm, đặt chỉ tiêu cắt giảm 25% điện năng tiêu thụ so với hệ thống máy nén khí cổ điển tại các cơ sở nghiên cứu và y tế.

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

Nhà nghiên cứu và giảng viên chuyên ngành Vật lý chất rắn, Vật lý nhiệt độ thấp: Tài liệu cung cấp cơ sở dữ liệu thực nghiệm toàn diện về cấu trúc tinh thể Ni2In, phân tích phổ nhiễu xạ nơtron và lý thuyết trường tinh thể Stevens, rất hữu ích cho công tác giảng dạy chuyên đề sau đại học và phát triển các đề tài nghiên cứu cơ bản.

Kỹ sư thiết kế hệ thống nhiệt lạnh và công nghệ cryo: Luận văn là nguồn tham khảo quan trọng về các thông số nhiệt động học từ tính như nhiệt dung, bước nhảy nhiệt lượng và công suất làm lạnh tương đối, phục vụ trực tiếp cho việc tính toán thiết kế tầng tiền làm lạnh từ 20 K đến 170 K trong các hệ thống hóa lỏng khí công nghiệp.

Chuyên gia phát triển vật liệu tiên tiến tại các viện công nghệ: Hưởng lợi từ quy trình thực nghiệm chi tiết về kỹ thuật nấu hồ quang chân không, ủ nhiệt đẳng nhiệt 1073 K trong ống thạch anh và phương pháp xử lý dữ liệu Rietveld qua phần mềm Rietan-2000, giúp rút ngắn thời gian chuẩn hóa quy trình chế tạo hợp kim liên kim loại.

Học viên cao học và nghiên cứu sinh ngành Khoa học Vật liệu: Đây là hình mẫu tiêu chuẩn về phương pháp luận nghiên cứu khoa học thực nghiệm, kết hợp đồng bộ các kỹ thuật phân tích hiện đại từ XRD, NPD tại lò phản ứng hạt nhân, từ kế SQUID độ nhạy cao đến thiết bị đo nhiệt dung PPMS.

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

Cấu trúc tinh thể của hợp chất Ho2In và Tb2In có đặc trưng hình học như thế nào?

Cả hai hợp chất đều kết tinh ở hệ tinh thể lục giác loại Ni2In thuộc nhóm không gian P63/mmc với hai đơn vị công thức trên mỗi ô cơ sở. Điểm đặc biệt là các nguyên tử đất hiếm phân bố ở hai vị trí tinh thể bất tương đương là vị trí 2a và vị trí 2d với môi trường phối trí khác nhau, trong khi nguyên tử In chiếm vị trí 2c, tạo nên tính dị hướng từ trường tinh thể rất phức tạp.

Vì sao mômen từ bão hòa đo được của Ho2In tại 1,8 K lại nhỏ hơn giá trị lý thuyết ion tự do?

Tại nhiệt độ 1,8 K dưới từ trường 70 kOe, mômen từ bão hòa thực nghiệm của Ho2In đạt 8,41 Bohr magneton, thấp hơn mức lý thuyết 10,0 Bohr magneton của ion tự do Ho3+. Nguyên nhân chính là do điện trường nội tại của mạng tinh thể tác động lên đám mây điện tử lớp 4f, làm suy biến các mức năng lượng và dẫn đến hiện tượng đóng băng một phần mômen quỹ đạo.

Kỹ thuật nhiễu xạ nơtron đem lại ưu thế gì vượt trội so với nhiễu xạ tia X trong nghiên cứu này?

Tia X chỉ tương tác với đám mây điện tử nên chỉ xác định được cấu trúc tinh thể nguyên tử. Trong khi đó, nơtron mang mômen từ nội tại nên tán xạ mạnh với mômen từ của ion đất hiếm. Nhờ vậy, phép đo nơtron tại bước sóng 1,6430 Angstrom đã tách biệt được mômen từ cục bộ tại vị trí 2a (9,0 Bohr magneton) và vị trí 2d (8,8 Bohr magneton), đồng thời xác định chính xác góc nghiêng mômen từ 41 độ.

Bản chất chuyển pha từ của Ho2In và Tb2In khác nhau như thế nào qua phép đo nhiệt dung?

Hợp chất Ho2In thuộc loại chuyển pha từ loại một với đỉnh nhiệt dung sắc nhọn dạng hàm lambda tại nhiệt độ Curie 87 K có độ lớn bước nhảy 20 J/mol.K. Ngược lại, hợp chất Tb2In thể hiện chuyển pha từ loại hai tại nhiệt độ Curie 160 K với đường cong nhiệt dung đối xứng có bước nhảy 18 J/mol.K, sau đó tiếp tục chuyển pha sang trạng thái phản sắt từ tại 50 K.

Khả năng ứng dụng thực tế của hệ vật liệu R2In trong công nghệ làm lạnh từ ra sao?

Hệ hợp chất R2In sở hữu công suất làm lạnh tương đối rất lớn, điển hình như Tb2In đạt 660 J/kg và Ho2In đạt 360 J/kg dưới từ trường 5 Tesla. Các thông số này biến chúng thành vật liệu nhiệt từ lý tưởng để chế tạo các máy làm lạnh từ tính nhiều tầng hoạt động trong khoảng nhiệt độ từ 20 K đến 170 K, phục vụ hóa lỏng khí tự nhiên và làm mát thiết bị y tế siêu dẫn.

Kết luận

Tổng kết những đóng góp khoa học chính của luận văn:

  • Chế tạo thành công hệ hợp chất liên kim Ho2In và Tb2In đơn pha cấu trúc lục giác loại Ni2In thuộc nhóm không gian P63/mmc bằng phương pháp nóng chảy hồ quang chân không 10^-5 Torr kết hợp ủ nhiệt đẳng nhiệt 1073 K trong 120 giờ.
  • Phát hiện và làm sáng tỏ quy luật hai chuyển pha từ liên tiếp trong Ho2In (TC = 87 K, Tt = 36 K) và Tb2In (TC = 160 K, TN = 50 K) thông qua các phép đo từ độ SQUID từ 1,8 K đến 300 K.
  • Định lượng chính xác cấu trúc từ vi mô của Ho2In bằng nhiễu xạ bột nơtron tại 2 K với mômen từ vị trí 2a đạt 9,0 Bohr magneton, vị trí 2d đạt 8,8 Bohr magneton và góc nghiêng 41 độ so với trục c.
  • Xác định bản chất nhiệt động học của chuyển pha từ loại một trên Ho2In với bước nhảy nhiệt dung Delta C = 20 J/mol.K và chuyển pha loại hai trên Tb2In với bước nhảy Delta C = 18 J/mol.K.
  • Chứng minh tiềm năng ứng dụng vượt trội của hệ vật liệu R2In trong công nghệ làm lạnh từ tính nhiệt độ thấp với công suất làm lạnh đạt từ 360 J/kg đến 660 J/kg dưới từ trường 5 Tesla.

Đóng góp quan trọng nhất của công trình là cung cấp bức tranh toàn diện và chuẩn xác đầu tiên về mối liên hệ giữa cấu trúc tinh thể, cấu trúc từ vi mô và tính chất từ nhiệt của hệ R2In trong nước. Kế hoạch tiếp theo dự kiến triển khai mở rộng pha tạp đa nguyên tố trong 12 tháng tới và hoàn thiện thiết kế module thử nghiệm làm lạnh từ trong vòng 24 tháng. Các đơn vị nghiên cứu, viện công nghệ và doanh nghiệp thiết bị nhiệt lạnh hãy chủ động liên hệ hợp tác để cùng khai thác nguồn dữ liệu thực nghiệm giá trị này, chung tay phát triển các giải pháp làm lạnh từ tính xanh, hiệu quả và bền vững cho tương lai.