Tổng quan nghiên cứu

Hiện tượng siêu dẫn được nhà vật lý học Kamerlingh Onnes phát hiện lần đầu tiên vào năm 1911 trên mẫu thủy ngân ở nhiệt độ 4,2 K. Đến tháng 1 năm 1986, bước ngoặt lịch sử diễn ra khi Bednorz và Müller công bố trạng thái siêu dẫn trong hệ gốm La-Ba-Cu-O dưới nhiệt độ 35 K, chính thức mở ra kỷ nguyên nghiên cứu siêu dẫn nhiệt độ cao. Khám phá này đã vượt qua giới hạn nhiệt độ hóa lỏng của Nitơ (77 K), mang đến tiềm năng ứng dụng khổng lồ trong công nghiệp năng lượng, giao thông vận tải và điện tử học.

Vấn đề then chốt trong vật lý chất rắn hiện đại là việc giải mã cơ chế dẫn điện vi mô và tìm hiểu sự suy giảm tính chất siêu dẫn khi cấu trúc tinh thể bị xáo trộn bởi các nguyên tử pha tạp. Nghiên cứu này tập trung phân tích sự biến đổi cấu trúc và tính chất điện, từ của hai hệ vật liệu siêu dẫn tiêu biểu là Y-123 và Bi-2223 khi thay thế một lượng nhỏ các kim loại chuyển tiếp từ tính gồm Sắt (Fe), Coban (Co) và Niken (Ni) vào vị trí của Đồng (Cu).

Mục tiêu cụ thể của đề tài là làm chủ quy trình công nghệ chế tạo gốm phản ứng pha rắn, tối ưu hóa các thông số xử lý nhiệt để đạt tỷ phần pha siêu dẫn cao trên 93%, đồng thời đánh giá định lượng sự biến thiên của nhiệt độ chuyển pha siêu dẫn (Tc) và các hằng số mạng khi nồng độ pha tạp thay đổi từ x = 0,00 đến 0,10. Nghiên cứu được thực hiện vào năm 2018 tại Phòng thí nghiệm 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. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc xác định rõ vai trò quyết định của các mặt phẳng dẫn CuO2 đối với trạng thái siêu dẫn có nhiệt độ chuyển pha đạt tới 112 K, cung cấp cơ sở dữ liệu thực nghiệm quan trọng để nâng cao mật độ dòng tới hạn (Jc) trong ứng dụng cáp truyền tải điện không tổn hao.

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 của các lý thuyết nhiệt động học và vật lý vi mô kinh điển về siêu dẫn.

Lý thuyết vi mô BCS (Bardeen-Cooper-Schrieffer) giải thích sự hình thành trạng thái siêu dẫn thông qua các cặp electron liên kết gián tiếp qua tương tác với dao động mạng tinh thể (phonon), gọi là các cặp Cooper. Ở dưới nhiệt độ tới hạn Tc, một khe năng lượng Delta xuất hiện ở mặt Fermi, ngăn các cặp Cooper va chạm và triệt tiêu hoàn toàn điện trở thuần. Khi đưa các ion kim loại chuyển tiếp 3d có từ tính như Fe, Co, Ni vào mạng tinh thể, tương tác trao đổi giữa mô-men từ cục bộ và spin của electron sẽ làm bẻ gãy các cặp Cooper, dẫn đến sự suy giảm nhanh chóng của nhiệt độ chuyển pha.

Lý thuyết hiện tượng luận Ginzburg-Landau mô tả quá trình chuyển pha thông qua thông số trật tự psi(r), đóng vai trò như một hàm sóng vĩ mô của các cặp Cooper. Lý thuyết này phân lập hoàn hảo các đặc trưng của siêu dẫn loại 2, nơi tồn tại trạng thái hỗn hợp với ba giá trị từ trường tới hạn là Hc1 (từ trường tới hạn dưới), Hc (từ trường nhiệt động) và Hc2 (từ trường tới hạn trên).

Phương trình London cùng hiệu ứng Meissner - Ochsenfeld giải thích tính chất nghịch từ lý tưởng của chất siêu dẫn khi toàn bộ từ trường ngoài bị đẩy ra khỏi khối vật liệu với độ thẩm thấu sâu đặc trưng lambda_L. Cùng với đó, mô hình lý thuyết Anderson về mối quan hệ giữa nhiệt độ chuyển pha Tc và số lượng lớp dẫn CuO2 trong ô cơ sở Perovskite phân lớp cung cấp khung đối chiếu quan trọng cho các hệ gốm siêu dẫn gốc Bismuth.

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

Hệ mẫu khảo sát gồm 12 mẫu gốm oxit thuộc hợp thức danh định Bi2Sr2Ca2(Cu1-xMx)3Oy, trong đó M lần lượt là Fe, Co, Ni với nồng độ thay thế x biến thiên từ 0,00 đến 0,10 (bước nhảy 0,02).

Phương pháp chế tạo mẫu dựa trên công nghệ phản ứng pha rắn (công nghệ gốm). Các bột oxit ban đầu gồm Bi2O3, SrCO3, CaCO3, CuO cùng các oxit pha tạp Fe2O3, CoO, Ni2O3 có độ tinh khiết cao từ 3N đến 4N (99,9% đến 99,99%) được cân chính xác đến 10^-4 gam theo tỷ lệ hợp thức tính toán tự động bằng phần mềm chuyên dụng. Hỗn hợp bột được nghiền mịn bằng cối mã não từ 60 đến 90 phút nhằm đồng đều kích thước hạt và tăng diện tích tiếp xúc. Bột sau đó được ép thủy lực thành các viên mẫu hình tròn đường kính 20 mm, dày 1 đến 2 mm dưới áp lực cao.

Quy trình xử lý nhiệt trải qua ba giai đoạn khắt khe: nung sơ bộ dưới 820 độ C, nung thiêu kết tạo pha trong khoảng 840 đến 870 độ C với thời gian duy trì từ 24 đến 240 giờ, và ủ nhiệt định hình trong dải 500 đến 600 độ C từ 24 đến 150 giờ trước khi để nguội chậm theo lò về nhiệt độ phòng trong môi trường kiểm soát.

Cấu trúc tinh thể và các hằng số mạng (a, b, c, thể tích ô cơ sở V) được phân tích bằng phương pháp nhiễu xạ tia X (XRD) trên máy quét góc 2-theta từ 20 đến 80 độ, bước quét 0,100 độ sử dụng bức xạ Cu-K_alpha. Các quá trình nhiệt động, điểm chuyển pha nhiệt và biến thiên khối lượng được xác định qua phép đo nhiệt lượng kế quét vi phân (DSC) kết hợp phân tích nhiệt trọng lượng (TGA/DTG). Tính chất dẫn điện và nhiệt độ chuyển pha siêu dẫn được đo bằng phương pháp bốn mũi dò mạ vàng tiêu chuẩn tiếp xúc trực tiếp trên bề mặt mẫu, kết nối hệ thống đo tự động Keithley qua giao tiếp IEEE488 trong buồng chân không 10^-2 đến 10^-4 mmHg, sử dụng cặp nhiệt điện Cu-Constantan chuẩn hóa tại nhiệt độ Nitơ lỏng 77 K.

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

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

Quá trình phân tích thực nghiệm đã mang lại các phát hiện cốt lõi về cấu trúc và tính chất chuyển pha của hệ siêu dẫn nhiệt độ cao:

Thứ nhất, phân tích nhiệt DSC và TGA chỉ ra rằng pha siêu dẫn Bi-2223 có khoảng nhiệt độ phân rã sắc nét từ 890,21 độ C đến 898,00 độ C với nhiệt lượng phân rã Q = 97,53 J/g và độ giảm trọng lượng do thoát oxy đạt 4,232%. Bằng quy trình công nghệ tối ưu, mẫu Bi-2223 không pha tạp đạt tỷ phần thể tích pha siêu dẫn lên tới 93%, thể hiện độ tinh khiết pha vượt trội so với các phương pháp gốm truyền thống.

Thứ hai, các phép đo điện trở bốn mũi dò xác định mẫu Bi-2223 chuẩn có nhiệt độ bắt đầu chuyển pha Tc1 = 112 K và bước chuyển pha thứ hai Tc2 = 90 K tương ứng với sự đồng tồn tại của pha Bi-2212. Khi tăng nồng độ pha tạp x từ 0,00 đến 0,10, nhiệt độ chuyển pha Tc suy giảm liên tục. Đối với tạp chất Ni, tính siêu dẫn duy trì đến nồng độ x = 0,08. Ngược lại, đối với tạp chất Fe và Co, trạng thái siêu dẫn bị triệt tiêu hoàn toàn khi nồng độ đạt x = 0,05, vật liệu chuyển hẳn sang tính chất bán dẫn ở dải nhiệt độ dưới 70 K.

Thứ ba, dữ liệu nhiễu xạ tia X khẳng định các ion Fe, Co, Ni khi pha tạp không làm thay đổi hệ tinh thể Perovskite trực thoi của pha Bi-2223. Tuy nhiên, các hằng số mạng a và b giảm nhẹ khoảng 0,8%, trục c duy trì ổn định quanh giá trị 37,1 Angstrom, trong khi thể tích ô cơ sở V giảm tuyến tính từ 540 Angstrom khối xuống còn khoảng 532 Angstrom khối khi x tăng từ 0,00 lên 0,10.

Thảo luận kết quả

Sự co rút của thể tích ô cơ sở V và các hằng số mạng a, b được giải thích trực tiếp bởi bán kính ion của các kim loại chuyển tiếp 3d thay thế (Fe3+, Co3+, Ni2+) nhỏ hơn đáng kể so với bán kính ion Cu2+. Trong cấu trúc tinh thể của hệ Bi-2223, nguyên tử Đồng chỉ chiếm duy nhất vị trí Cu(2) nằm bên trong các mặt phẳng dẫn CuO2. Do đó, các nguyên tử tạp chất Fe, Co, Ni có xu hướng chiếm trực tiếp các vị trí Cu(2) này mà không làm xáo trộn các lớp đệm BiO hay SrO.

Điều này tạo nên sự khác biệt căn bản khi so sánh với hệ siêu dẫn Y-123. Trong hệ Y-123, Đồng chiếm cả hai vị trí gồm Cu(1) trên chuỗi một chiều CuO và Cu(2) trên mặt phẳng hai chiều CuO2. Khi pha tạp Fe hoặc Co vào Y-123, tạp chất ưu tiên chiếm vị trí Cu(1), gây ra sự chuyển pha cấu trúc từ trực thoi sang lập phương và làm xáo trộn mạnh trật tự oxy. Ngược lại, ở hệ Bi-2223, sự thay thế chỉ diễn ra tại mặt phẳng Cu(2), bảo toàn cấu trúc mạng nhưng tác động trực tiếp vào cơ chế dẫn điện.

Sự suy giảm mạnh của nhiệt độ chuyển pha Tc trong Bi-2223 pha tạp Fe, Co, Ni bắt nguồn từ hiện tượng bẻ gãy cặp Cooper do tán xạ spin từ tính. Mô-men từ cục bộ của các ion 3d tương tác với spin của các electron trong cặp Cooper, phá vỡ tính đối xứng nghịch đảo thời gian và làm suy giảm mật độ trạng thái điện tử tại mức Fermi N(EF). Mẫu pha tạp Ni duy trì tính siêu dẫn tốt hơn Fe và Co là do cấu hình điện tử và trạng thái hóa trị của Ni2+ gần tương đồng với Cu2+, ít gây nhiễu loạn mật độ hạt tải lỗ trống trong mặt phẳng dẫn.

Về mặt trình bày dữ liệu, các kết quả này được minh họa tối ưu qua biểu đồ đường cong điện trở chuẩn hóa theo nhiệt độ R(T)/R(300K) để làm nổi bật hai bước chuyển pha Tc1 và Tc2, kết hợp bảng thống kê chi tiết các tham số mạng tinh thể (a, b, c, V) và đồ thị phân tích nhiệt trọng lượng DTA-TGA biểu diễn lượng oxy thoát ra theo nhiệt độ nung.

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

Nhằm phát triển và ứng dụng hiệu quả các kết quả nghiên cứu về vật liệu siêu dẫn nhiệt độ cao, bốn khuyến nghị chiến lược được đề xuất:

Thứ nhất, tối ưu hóa quy trình thiêu kết và ủ nhiệt đẳng nhiệt. Các phòng thí nghiệm công nghệ vật liệu cần kiểm soát tốc độ gia nhiệt ở mức 2 đến 3 độ C/phút và duy trì thời gian ủ nhiệt kéo dài từ 120 đến 150 giờ trong môi trường dòng khí Oxy tinh khiết nhằm mục tiêu nâng cao tỷ phần pha Bi-2223 đạt trên 96%, hoàn thành trước quý IV năm 2027 do các nhóm nghiên cứu công nghệ gốm chủ trì.

Thứ hai, triển khai nghiên cứu giải pháp đồng pha tạp (Co-doping). Khuyến nghị kết hợp pha tạp kim loại chuyển tiếp với các nguyên tố Chì (Pb) hoặc đất hiếm (RE) với hàm lượng kiểm soát từ 0,01 đến 0,03 nhằm tạo các tâm ghim từ thông (flux pinning centers), nâng cao mật độ dòng điện tới hạn Jc vượt mức 10^5 A/cm2 trong giai đoạn 2027-2028 thực hiện bởi các viện nghiên cứu chuyên ngành vật lý nhiệt độ thấp.

Thứ ba, ứng dụng các mô hình mô phỏng tính toán cấu trúc điện tử. Các chuyên gia vật lý tính toán cần áp dụng phương pháp phiếm hàm mật độ (DFT) kết hợp thuật toán làm khớp Rietveld trong vòng 6 tháng tới để xác định chính xác phân bố mật độ điện tích và cấu hình spin của các ion Fe, Co, Ni tại vị trí Cu(2) trong mạng tinh thể.

Thứ tư, hiện đại hóa trang thiết bị đo đạc điện từ trường cao. Đề xuất các trường đại học khối khoa học tự nhiên đầu tư nâng cấp hệ đo bốn mũi dò tích hợp từ trường ngoài đến 5 Tesla trong lộ trình 12 tháng, phục vụ việc khảo sát toàn diện các giá trị từ trường tới hạn Hc1, Hc2 và cơ chế chuyển pha hỗn hợp.

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

Nội dung và dữ liệu thực nghiệm của luận văn mang lại giá trị tham khảo thiết thực cho bốn nhóm đối tượng chuyên môn:

Nhóm học viên cao học và nghiên cứu sinh chuyên ngành Vật lý chất rắn, Vật lý nhiệt độ thấp: Luận văn cung cấp tài liệu chuẩn mực về quy trình tổng hợp gốm siêu dẫn họ BSCCO, kỹ thuật xử lý số liệu nhiễu xạ tia X và phương pháp đo điện trở bốn mũi dò chính xác ở nhiệt độ Nitơ lỏng.

Nhóm kỹ sư nghiên cứu và phát triển (R&D) trong ngành kỹ thuật điện và năng lượng: Nắm vững các giới hạn nồng độ pha tạp và nguyên nhân suy giảm dòng tới hạn giúp các kỹ sư lựa chọn vật liệu tối ưu khi chế tạo cáp siêu dẫn truyền tải điện cao thế và cuộn dây nam châm cho máy chụp cộng hưởng từ (MRI).

Nhóm giảng viên và nhà nghiên cứu tại các viện, trường đại học: Dữ liệu thực nghiệm về phân tích nhiệt DSC/TGA và các phép so sánh cấu trúc giữa hệ Y-123 và Bi-2223 đóng vai trò là học liệu mẫu phục vụ giảng dạy chuyên đề chuyển pha nhiệt động và cấu trúc Perovskite.

Nhóm chuyên gia kiểm chuẩn và phân tích vật liệu tiên tiến: Khai thác phương pháp tính toán hợp thức hóa học, kỹ thuật phân tích pha Rơn-ghen và quy trình kiểm soát độ mất oxy để ứng dụng vào công nghệ chế tạo màng mỏng và vật liệu gốm kỹ thuật cao.

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

Tại sao việc thay thế kim loại chuyển tiếp vào vị trí Cu lại làm suy giảm nhiệt độ chuyển pha siêu dẫn Tc? Nguyên nhân chủ yếu do các ion Fe, Co, Ni mang mô-men từ cục bộ gây ra hiện tượng tán xạ đảo spin điện tử, phá vỡ liên kết của các cặp Cooper theo lý thuyết vi mô BCS. Khi nồng độ Fe đạt x = 0,05, nhiệt độ Tc giảm từ 112 K về 0 K, khiến vật liệu chuyển sang trạng thái bán dẫn.

Sự khác biệt cơ bản về vị trí chiếm chỗ của tạp chất Cu giữa hệ Y-123 và Bi-2223 là gì? Trong hệ Y-123, tạp chất chiếm cả vị trí Cu(1) dạng chuỗi và Cu(2) dạng mặt phẳng, làm biến đổi đối xứng mạng từ trực thoi sang lập phương. Ở hệ Bi-2223, tạp chất chỉ chiếm vị trí Cu(2) trong mặt phẳng CuO2, làm giảm thể tích ô cơ sở mà không làm thay đổi kiểu cấu trúc tinh thể.

Phương pháp bốn mũi dò mạ vàng có ưu điểm gì khi đo điện trở vật liệu siêu dẫn? Phương pháp bốn mũi dò tách biệt hoàn toàn mạch cấp dòng điện và mạch đo hiệu điện thế, giúp triệt tiêu hoàn toàn điện trở tiếp xúc tại các đầu đo. Thiết bị đo Keithley kết nối tự động cho phép ghi nhận chính xác giá trị điện trở triệt tiêu về mức dưới 10^-26 Ohm.m tại điểm chuyển pha Tc.

Vì sao tạp chất Niken (Ni) duy trì tính siêu dẫn ở nồng độ cao hơn so với Sắt (Fe) và Coban (Co)? Ion Ni2+ có trạng thái hóa trị và cấu hình điện tử tương đồng với ion Cu2+ hơn so với các ion Fe3+ và Co3+. Do đó, mức độ xáo trộn mật độ hạt tải và cường độ tán xạ bẻ gãy cặp Cooper do Ni gây ra thấp hơn, cho phép trạng thái siêu dẫn tồn tại đến nồng độ x = 0,08.

Nhiệt độ thiêu kết tối ưu để hình thành pha siêu dẫn Bi-2223 là bao nhiêu? Pha siêu dẫn Bi-2223 đòi hỏi năng lượng hoạt hóa cao khoảng 485 kJ/mol. Quá trình tạo pha tối ưu yêu cầu nung thiêu kết trong khoảng 840 đến 870 độ C từ 24 đến 240 giờ, sau đó ủ nhiệt tại 500 đến 600 độ C để đạt tỷ phần pha trên 93%.

Kết luận

  • Làm chủ quy trình chế tạo gốm phản ứng pha rắn, tổng hợp thành công hệ siêu dẫn nhiệt độ cao Bi-2223 với tỷ phần pha đạt 93% và nhiệt độ bắt đầu chuyển pha Tc1 đạt 112 K.
  • Xác định quy luật biến thiên tinh thể học khi thể tích ô cơ sở V giảm tuyến tính từ 540 Angstrom khối về 532 Angstrom khối theo sự gia tăng nồng độ pha tạp Fe, Co, Ni từ x = 0,00 đến 0,10.
  • Làm sáng tỏ cơ chế vi mô của sự suy giảm tính siêu dẫn thông qua hiện tượng tán xạ bẻ gãy cặp Cooper do mô-men từ cục bộ của kim loại chuyển tiếp 3d tại mặt phẳng dẫn CuO2.
  • Phân biệt rõ nét bản chất chiếm chỗ ion Cu(2) ở hệ Bi-2223 so với việc chiếm cả hai vị trí Cu(1) và Cu(2) ở hệ Y-123 khi pha tạp vi lượng.
  • Chứng minh giới hạn tồn tại của trạng thái siêu dẫn đối với tạp chất Ni đạt tới x = 0,08, vượt trội so với giới hạn x = 0,05 của Fe và Co.

Đóng góp lớn nhất của công trình là thiết lập hệ thống dữ liệu thực nghiệm toàn diện về mối quan hệ giữa cấu trúc tinh thể, trật tự từ tính và tính chất dẫn điện của vật liệu gốm siêu dẫn nhiệt độ cao chứa Bismuth.

Lộ trình nghiên cứu tiếp theo sẽ tập trung vào việc thử nghiệm đồng pha tạp nguyên tố đất hiếm trong vòng 18 tháng tới nhằm gia tăng mật độ dòng tới hạn Jc trong môi trường từ trường cao.

Các viện nghiên cứu, trường đại học và doanh nghiệp công nghệ cao hãy đẩy mạnh hợp tác khai thác các kết quả này để sớm đưa vật liệu siêu dẫn nhiệt độ cao vào ứng dụng thực tiễn trong hệ thống lưới điện thông minh và thiết bị y tế hiện đại.