Tổng quan nghiên cứu

Vật liệu đa tính sắt multiferroic bismuth ferrite BiFeO3 (BFO) là một trong những hệ vật liệu điện từ nổi bật nhất trong vật lý chất rắn hiện đại nhờ khả năng đồng thời tồn tại trật tự sắt điện với nhiệt độ Curie chuyển pha rất cao khoảng 1103 K và trật tự phản sắt từ với nhiệt độ Neel khoảng 643 K ngay tại nhiệt độ phòng. Sự kết hợp độc đáo này mở ra triển vọng to lớn trong việc thiết kế các linh kiện spin điện tử, cảm biến từ trường và bộ nhớ truy cập ngẫu nhiên từ điện thế hệ mới. Tuy nhiên, việc đưa vật liệu BFO vào ứng dụng thực tế gặp phải hai rào cản kỹ thuật nghiêm trọng: dòng điện rò nội tại tương đối lớn do khuyết tật nút mạng oxy và từ tính vĩ mô rất yếu bị triệt tiêu bởi cấu trúc spin xoắn ốc cycloid dạng sóng không điều hòa với chu kỳ dài khoảng 620 Å.

Nhằm giải quyết triệt để bài toán này, đề tài nghiên cứu tập trung vào việc biến tính cấu trúc mạng nền BFO thông qua kỹ thuật pha tạp ion đất hiếm Europi (Eu3+) thay thế cho vị trí của Bismuth (Bi3+). Mục tiêu trọng tâm của nghiên cứu là tổng hợp thành công hệ vật liệu gốm Bi1-xEuxFeO3 với các nồng độ pha tạp biến thiên x = 0,00; 0,05; 0,10; 0,15 và 0,20 bằng phương pháp phản ứng pha rắn truyền thống. Trên cơ sở đó, nghiên cứu tiến hành phân tích sâu sắc mối liên hệ hữu cơ giữa sự biến dạng mạng tinh thể, quá trình chuyển pha cấu trúc, sự tiến hóa của các mốt dao động phonon và sự chuyển đổi trật tự từ tính ở thang vĩ mô.

Toàn bộ quá trình thực nghiệm, chế tạo và khảo sát đặc trưng vật lý được thực hiện hoàn chỉnh tại Khoa Vật lý, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội cùng sự phối hợp đo đạc từ tính tại Viện Khoa học Vật liệu vào năm 2014. Kết quả nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc: việc pha tạp Eu3+ với tỷ lệ từ x = 0,10 trở lên đã loại bỏ hoàn toàn các pha tạp chất thứ cấp, đồng thời nâng giá trị từ độ cực đại từ mức 0,07 emu/g ở mẫu BFO nguyên bản lên 0,4484 emu/g ở mẫu pha tạp x = 0,20 (tăng trưởng hơn 540%), mở đường cho việc phát triển các vật liệu chuyển đổi điện từ hiệu năng cao.

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 vững chắc trên các nguyên lý nhiệt động học và vật lý lượng tử chất rắn:

  • Cấu trúc mạng tinh thể Perovskite ABO3 lý tưởng và thừa số dung hạn Goldschmidt: Cấu trúc perovskite chuẩn thuộc hệ lập phương đối xứng cao với khối bát diện BO6 nằm ở tâm ô cơ sở. Độ bền vững và mức độ biến dạng mạng được định lượng qua thừa số dung hạn Goldschmidt t = (RA + RO) / [sqrt(2) * (RB + RO)]. Khi giá trị t sai lệch khỏi mức 1 lý tưởng do sự chênh lệch bán kính ion giữa Eu3+ (khoảng 1,07 Å) và Bi3+ (khoảng 1,17 Å), mạng tinh thể bị biến dạng méo, kéo theo sự quay và nghiêng góc của các bát diện FeO6 xung quanh các trục tinh thể.
  • Hiệu ứng méo mạng Jahn-Teller và trạng thái năng lượng: Sự hạ thấp tính đối xứng không gian dẫn tới việc giải tỏa trạng thái suy biến của các orbital điện tử d trong cation kim loại chuyển tiếp Fe3+, giúp hệ chuyển dịch về trạng thái có năng lượng tự do cực tiểu bền vững hơn.
  • Lý thuyết trật tự từ và tương tác siêu trao đổi gián tiếp (Superexchange Interaction): Tính chất từ của BFO bắt nguồn từ tương tác siêu trao đổi góc Fe3+-O2--Fe3+ qua trung gian anion oxy. Ở trạng thái chuẩn, tương tác này duy trì trật tự phản sắt từ loại G bị bao phủ bởi cấu trúc spin xoắn cycloid chu kỳ 620 Å. Sự méo mạng tinh thể do ion đất hiếm gây ra phá vỡ tính tuần hoàn xoắn ốc này, giải phóng các mômen từ spin chưa bù trừ và kích hoạt trật tự sắt từ vĩ mô.

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

  • Cỡ mẫu và quy trình chế tạo: Nghiên cứu khảo sát hệ 5 mẫu gốm đơn pha và đa pha Bi1-xEuxFeO3 với các tỷ lệ thành phần danh định x lần lượt là 0,00; 0,05; 0,10; 0,15 và 0,20. Mẫu được chế tạo theo phương pháp phản ứng pha rắn (gốm truyền thống) từ các tiền chất oxit ban đầu có độ sạch cao gồm Bi2O3 (độ tinh khiết 99%), Fe2O3 (độ tinh khiết 99%) và Eu2O3 (độ tinh khiết đạt 99,9%). Hỗn hợp bột oxit được nghiền cơ học liên tục trong 8 giờ bằng cối mã não để đảm bảo độ đồng nhất kích thước hạt. Sau đó, bột được phối trộn chất kết dính hữu cơ PVA 2% và ép tạo hình khối chữ nhật dưới áp suất nén 5 tấn/cm2. Mẫu ép được sấy khô ổn định ở 80 oC trong 10 giờ, tiếp tục nung thiêu kết tại nhiệt độ 825 oC trong thời gian 5 giờ với tốc độ gia nhiệt 5 oC/phút, cuối cùng làm nguội tự nhiên theo buồng lò về nhiệt độ phòng 25 oC.
  • Lý do lựa chọn hệ thống phương pháp phân tích:
    • Phép đo nhiễu xạ tia X (XRD): Thực hiện trên hệ Bruker D5005 với ống phát bức xạ Cu-Kalpha (bước sóng lambda = 1,54056 Å), kết hợp phần mềm tính toán Chekcell nhằm xác định chính xác kiểu ô mạng, hằng số mạng (a, b, c), thể tích ô cơ sở và nhận diện thành phần pha cấu trúc mà không phá hủy mẫu.
    • Kính hiển vi điện tử quét (SEM) và phổ tán sắc năng lượng tia X (EDS): Khảo sát trên thiết bị Nova NanoSEM 450 với độ phân giải siêu cao 1,4 nm trong chân không cao, cho phép đánh giá vi cấu trúc bề mặt, độ xốp liên kết hạt và định lượng tỷ lệ hợp thức nguyên tử thực nghiệm.
    • Phổ vi tán xạ Raman (Micro-Raman): Đo trên hệ quang phổ Labram HR800 sử dụng nguồn kích thích laser bước sóng 632,8 nm nhằm khảo sát các mốt dao động liên kết mạng cục bộ Bi-O và Fe-O.
    • Từ kế mẫu rung (VSM): Khảo sát sự phụ thuộc của từ độ vào từ trường ngoài trong dải từ trường quét rộng từ -11,5 kOe đến +11,5 kOe tại nhiệt độ phòng 300 K để thu thập các thông số từ độ dư Mr, từ độ cực đại Mmax và lực kháng từ Hc.

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 4 phát hiện cốt lõi:

  1. Sự triệt tiêu hoàn toàn pha tạp chất thứ cấp: Trên phổ XRD, mẫu BFO không pha tạp (x = 0,00) và mẫu pha tạp nồng độ thấp (x = 0,05) xuất hiện các đỉnh nhiễu xạ phụ đặc trưng cho pha thứ cấp Bi2Fe4O9 do sự bay hơi không đồng đều của oxit bismuth. Tuy nhiên, khi nâng nồng độ pha tạp Eu3+ lên mức x = 0,10; 0,15 và 0,20, toàn bộ các đỉnh tạp biến mất hoàn toàn, mẫu đạt trạng thái đơn pha 100%.
  2. Chuyển pha cấu trúc tinh thể sắc nét: Tại dải nồng độ thấp 0,00 <= x <= 0,10, vật liệu kết tinh ở cấu trúc mặt thoi (Rhombohedral) thuộc nhóm không gian đối xứng R3c. Thể tích ô cơ sở giảm dần từ 372,71 Å3 (ở mẫu x = 0,00 với a = 5,574 Å, c = 13,854 Å) xuống 369,28 Å3 (ở mẫu x = 0,10 với a = 5,560 Å, c = 13,792 Å). Khi nồng độ đạt x = 0,15 và 0,20, hệ vật liệu chuyển pha hoàn toàn sang cấu trúc trực thoi (Orthorhombic) thuộc nhóm không gian Pnma với sự thu hẹp thể tích ô cơ sở xuống 236,50 Å3 (x = 0,15) và 238,51 Å3 (x = 0,20). Kích thước tinh thể nano trung bình tính toán theo công thức Debye-Scherrer dao động trong khoảng 10 đến 30 nm.
  3. Sự biến đổi mốt dao động mạng phonon: Phổ tán xạ Raman của mẫu BiFeO3 thể hiện rõ 3 mốt dao động A1(LO) ở vùng tần số thấp (140, 172, 221 cm-1) gắn liền với liên kết Bi-O và các mốt E(TO) ở vùng tần số cao hơn (261 đến 520 cm-1). Khi nồng độ Eu3+ tăng, các đỉnh A1 dịch dần về vùng tần số cao và suy giảm mạnh cường độ; tại x = 0,20, hai mốt dao động A1-1 và A1-2 biến mất hoàn toàn, khẳng định ion Eu3+ đã thế chỗ thành công vào vị trí của cation Bi3+.
  4. Kích hoạt và gia tăng mạnh mẽ tính sắt từ ở nhiệt độ phòng: Mẫu không pha tạp thể hiện tính thuận từ yếu tuyến tính với từ độ cực đại Mmax chỉ đạt 0,07 emu/g và từ độ dư Mr xấp xỉ bằng 0. Khi pha tạp Eu3+, đường cong từ trễ M(H) mở rộng rõ rệt: mẫu x = 0,05 đạt Mr = 0,0669 emu/g, Hc = 865 Oe; mẫu x = 0,15 đạt Mr = 0,1262 emu/g, lực kháng từ đạt đỉnh 1422 Oe; mẫu x = 0,20 đạt giá trị từ độ cực đại Mmax = 0,4484 emu/g (tăng gấp 6,4 lần so với mẫu gốc), Mr = 0,1395 emu/g và Hc = 1212 Oe.

Thảo luận kết quả

Cơ chế vật lý sâu xa của sự biến đổi tính chất bắt nguồn từ việc cation Eu3+ có bán kính nhỏ hơn Bi3+, làm suy giảm thừa số Goldschmidt t. Để giải tỏa áp lực cơ học nội tại, khối bát diện FeO6 buộc phải xoay quanh trục phân cực, làm thay đổi góc liên kết Fe-O-Fe và độ dài liên kết Bi/Eu-O. Sự biến dạng mạng này phá hủy hoàn toàn cấu trúc spin xoắn ốc tuần hoàn 620 Å vốn có của BFO, khiến các mômen từ spin cục bộ không còn bị triệt tiêu bù trừ mà định hướng song song một phần dưới tác dụng của từ trường ngoài.

Để trực quan hóa các kết quả này, toàn bộ tham số cấu trúc tinh thể và dữ liệu từ tính có thể được tổng hợp trực quan thông qua bảng so sánh hằng số mạng và bảng thông số từ tính (Mr, Mmax, Hc). Trên đồ thị XRD, sự dịch chuyển góc 2theta của cặp đỉnh phản xạ (110) và (104) tại vùng góc 32 độ minh chứng rõ ràng cho sự suy giảm hằng số mạng, trong khi sự xuất hiện của đỉnh nhiễu xạ mới tại góc 2theta khoảng 26 độ đánh dấu ranh giới chuyển pha sang cấu trúc Orthorhombic. Đồng thời, biểu đồ biểu diễn tương quan giữa từ độ dư Mr theo nồng độ pha tạp x thể hiện quy luật tăng trưởng gần như tuyến tính.

So sánh với các công bố khoa học trước đây trên hệ mẫu pha tạp nồng độ thấp x <= 0,04 vốn chỉ quan sát thấy tính thuận từ, kết quả của luận văn đã chứng minh vùng nồng độ x từ 0,05 đến 0,20 là vùng chuyển biến then chốt để khai mở pha sắt từ. Khi đối chiếu với các nghiên cứu chế tạo màng mỏng Bi1-xEuxFeO3 bằng kỹ thuật lắng đọng xung laser (PLD) trên đế đơn tinh thể LaAlO3 (vốn chuyển sang pha tứ giác Tetragonal tại x = 0,20 do ứng suất epitaxy từ đế), mẫu gốm khối trong nghiên cứu này chuyển sang pha Orthorhombic, phản ánh đúng bản chất nhiệt động học của vật liệu đa tinh thể dạng khối.

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

Dựa trên các phát hiện thực nghiệm, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật và định hướng ứng dụng:

  1. Tối ưu hóa chế độ xử lý nhiệt và kiểm soát áp suất riêng phần khí oxy: Các đơn vị nghiên cứu và phòng thí nghiệm vật liệu bán dẫn cần tiến hành nung thiêu kết trong môi trường khí quyển oxy tinh khiết 99,99% tại dải nhiệt độ 820 đến 830 oC nhằm triệt tiêu tối đa các nút khuyết oxy, giảm thiểu dòng điện rò xuống dưới ngưỡng 10^-6 A/cm2 trong thời gian 6 đến 12 tháng tới.
  2. Mở rộng thử nghiệm cơ chế đồng pha tạp (Co-doping): Nhóm nghiên cứu tại các viện và trường đại học nên triển khai hướng đồng pha tạp kết hợp cation đất hiếm (Eu3+, La3+, Gd3+) tại vị trí A và kim loại chuyển tiếp (Mn3+, Ti4+, Sc3+) tại vị trí B với tỷ lệ 5% đến 15% trong giai đoạn 2026 - 2027 nhằm tăng cường đồng thời phân cực sắt điện và nâng từ độ bão hòa vượt mức 1,0 emu/g.
  3. Chuẩn hóa quy trình chế tạo màng mỏng chất lượng cao: Các cơ sở sản xuất vi điện tử cần ứng dụng công nghệ PLD hoặc phún xạ catot magnetron để chế tạo màng mỏng Bi0,80Eu0,20FeO3 với chiều dày từ 50 đến 200 nm trên đế Si/SiO2/Pt, tạo tiền đề tích hợp vào cấu trúc chip nhớ truy cập ngẫu nhiên MERAM thương mại trong lộ trình 2 đến 3 năm.
  4. Ứng dụng mô phỏng tính toán vi mô nguyên lý đầu (DFT): Các nhóm nghiên cứu vật lý lý thuyết cần đẩy mạnh tính toán cấu trúc dải năng lượng và mật độ trạng thái điện tử (DOS) có tính đến tham số tương quan Hubbard U đối với hệ Bi1-xEuxFeO3, hoàn thành bộ dữ liệu chuẩn hóa trong vòng 18 tháng để định hướng chính xác nồng độ pha tạp tối ưu.

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

Tài liệu luận văn mang lại giá trị học thuật và thực tiễn chuyên sâu cho 4 nhóm đối tượng chính:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Vật lý chất rắn, Khoa học vật liệu và Vật lý kỹ thuật: Nắm bắt phương pháp luận nghiên cứu chuẩn mực, quy trình tổng hợp mẫu gốm oxit phức hợp và kỹ năng xử lý dữ liệu nhiễu xạ tia X, phổ Raman và từ kế mẫu rung VSM.
  2. Kỹ sư nghiên cứu và phát triển (R&D) trong ngành công nghiệp bán dẫn và vi điện tử: Khai thác các dữ liệu thực nghiệm về lực kháng từ cao (Hc > 1200 Oe) và tính đơn pha của hệ gốm để thiết kế các phần tử nhớ từ điện FeRAM, cảm biến đo từ trường độ nhạy cao và linh kiện spin điện tử.
  3. Giảng viên và nhà khoa học tại các trường đại học chuyên ngành kỹ thuật: Sử dụng làm tài liệu tham khảo chất lượng cao phục vụ giảng dạy các chuyên đề cao học về Vật liệu multiferroic, Hiệu ứng Jahn-Teller, và Kỹ thuật phân tích quang phổ thực nghiệm.
  4. Chuyên gia phát triển công nghệ cảm biến và chuyển đổi năng lượng: Tìm kiếm giải pháp tối ưu hóa hiệu ứng ghép cặp từ - điện (ME effect) ứng dụng trong các thiết bị vi cơ điện tử (MEMS) và các bộ thu hồi năng lượng cơ từ trường.

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

  1. Vì sao vật liệu BiFeO3 nguyên bản lại có từ tính vĩ mô rất yếu ở nhiệt độ phòng? Vật liệu BiFeO3 nguyên bản sở hữu trật tự phản sắt từ loại G với các mômen từ của ion Fe3+ định hướng đối song song. Thêm vào đó, cấu trúc spin xoắn ốc cycloid không điều hòa với chu kỳ dài khoảng 620 Å bao phủ lên mạng tinh thể đã triệt tiêu hoàn toàn từ độ thực, khiến mẫu x = 0,00 chỉ thể hiện tính thuận từ với từ độ cực đại 0,07 emu/g.

  2. Pha tạp Eu3+ đóng vai trò gì trong việc triệt tiêu các pha tạp chất thứ cấp? Bismuth ferrite là pha không bền nhiệt và oxit bismuth rất dễ bay hơi trong quá trình nung ở nhiệt độ cao, tạo ra các pha tạp như Bi2Fe4O9. Việc đưa ion Eu3+ có liên kết mạng bền vững hơn vào thay thế vị trí của Bi3+ giúp ổn định hóa cấu trúc nhiệt động học, loại bỏ hoàn toàn pha tạp chất ở nồng độ pha tạp từ x = 0,10 đến 0,20.

  3. Dấu hiệu thực nghiệm nào chứng minh sự chuyển pha cấu trúc từ Rhombohedral sang Orthorhombic? Trên giản đồ nhiễu xạ tia X (XRD), khi nồng độ pha tạp tăng lên x = 0,15 và 0,20, các đỉnh nhiễu xạ mở rộng, xuất hiện đỉnh phản xạ mới tại góc 2theta khoảng 26 độ và có sự tách đỉnh rõ rệt tại góc 32 độ, đồng thời thể tích ô cơ sở tính toán co rút mạnh từ 369,28 Å3 xuống 236,50 Å3.

  4. Phổ tán xạ Raman phản ánh sự thế chỗ của ion Eu3+ vào mạng tinh thể như thế nào? Phổ Raman của BFO có 3 mốt dao động A1(LO) ở dải tần số thấp 140, 172 và 221 cm-1 gắn liền trực tiếp với liên kết Bi-O. Khi nồng độ Eu3+ tăng dần, các đỉnh A1 suy giảm mạnh cường độ và dịch về tần số cao; tại x = 0,20, hai đỉnh A1-1 và A1-2 biến mất hoàn toàn, chứng minh liên kết Bi-O đã bị thay thế bởi liên kết Eu-O.

  5. Vì sao có sự khác biệt về cấu trúc giữa mẫu gốm khối và mẫu màng mỏng tại nồng độ x = 0,20? Mẫu gốm khối kết tinh đa tinh thể tự do nên chuyển sang cấu trúc trực thoi Orthorhombic (Pnma) để giải tỏa ứng suất nội. Ngược lại, mẫu màng mỏng chế tạo bằng phương pháp lắng đọng xung laser (PLD) chịu lực căng mạng epitaxy đơn tinh thể mạnh mẽ từ đế LaAlO3, dẫn tới việc định hình cấu trúc tứ giác Tetragonal với tỷ số tham số mạng c/a khoảng 1,42.

Kết luận

  • Tổng hợp thành công hệ vật liệu đơn pha Bi1-xEuxFeO3 (x = 0,10; 0,15 và 0,20) bằng phương pháp phản ứng pha rắn truyền thống tại nhiệt độ nung thiêu kết 825 oC trong 5 giờ.
  • Khám phá ranh giới chuyển pha cấu trúc rõ nét từ hệ mặt thoi Rhombohedral (nhóm R3c) sang hệ trực thoi Orthorhombic (nhóm Pnma) khi nồng độ pha tạp Eu3+ đạt ngưỡng x = 0,15 đến 0,20.
  • Xác thực cơ chế thay thế của cation Eu3+ vào vị trí của Bi3+ thông qua sự biến mất của các mốt dao động mạng A1 trên phổ vi tán xạ Raman và sự biến thiên tỷ lệ nguyên tố trên phổ EDS.
  • Cải thiện vượt bậc tính chất sắt từ ở nhiệt độ phòng với từ độ cực đại đạt 0,4484 emu/g (tăng hơn 540%), từ độ dư đạt 0,1395 emu/g và lực kháng từ đạt 1212 Oe tại nồng độ pha tạp x = 0,20.
  • Cung cấp luận cứ khoa học quan trọng khẳng định cơ chế phá vỡ trật tự spin xoắn ốc bằng biến dạng mạng tinh thể, đặt nền tảng vững chắc cho lộ trình 12 đến 24 tháng tới trong việc chế tạo các linh kiện nhớ từ điện đa trạng thái hiệu năng cao. Các nhà nghiên cứu và kỹ sư vi điện tử hãy chủ động ứng dụng dữ liệu thực nghiệm này để mở rộng phát triển các thế hệ linh kiện chuyển đổi năng lượng và cảm biến thông minh.