Luận văn thạc sĩ HUS mô hình ising và ứng dụng cho các chất sắt từ 04

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Kim Oanh MƠ HÌNH ISING VÀ ỨNG DỤNG CHO CÁC CHẤT SẮT TỪ LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – Năm 2014 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Thị Kim Oanh MƠ HÌNH ISING VÀ ỨNG DỤNG CHO CÁC CHẤT SẮT TỪ Chuyên ngành: Vật lý chất rắn Mã số: 60440104 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: TS Bạch Hương Giang GS.TS Bạch Thành Công Hà Nội – Năm 2014 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    MỤC LỤC MỞ ĐẦU CHƯƠNG 1: SẮT TỪ VÀ MƠ HÌNH ISING 1.1: Đặc điểm chất sắt từ 1.2:Hiện tượng chuyển pha vật liệu sắt từ 1.2.1.Pha chuyển pha 1.2.2.Phân loại chuyển pha 1.2.3: Chuyển pha sắt từ-thuận từ 1.3 : Mơ hình Ising cho chất sắt từ 1.3.2: Lời giải xác cho mơ hình hai chiều 13 1.3.3 Mơ hình Ising ba chiều 15 1.3.4 Năng lượng tự , mô men từ , độ từ hóa mơ hình Ising .19 1.3.5 : Kết luận 21 CHƯƠNG 2: MƠ HÌNH ISING MẤT TRẬT TỰ VỚI TÍCH PHÂN TRAO ĐỔI THĂNG GIÁNG VÀ ỨNG DỤNG .22 2.1: Hệ thức Callen cho mơ hình Ising trật tự .22 2.1.1: Hệ thức Callen cho mơ hình Ising trật tự 22 2.1.2 Mô hình Ising trật tự với tích phân trao đổi thăng giáng hệ thức Callen 24 2.1.3: Phương trình đại số cho mơ men từ nút mạng nhận phương pháp biến đổi tích phân 26 2.2: Phương pháp Monte Carlo [5] 34 2.2.1: Thuật toán Metropolis 34 2.2.2: Áp dụng cho mơ hình Ising hai chiều .37 Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   CHƯƠNG : KẾT QUẢ VÀ THẢO LUẬN 40 3.1: Đường cong từ nhiệt m(t) có khơng có từ trường ngồi 41 3.1.1: Mạng hai chiều .41 3.1.2: Mạng ba chiều (z=6) 45 3.2 Đường biểu diễn phụ thuộc nhiệt độ chuyển pha Curie vào xác suất p .47 3.2.1 : Mạng hai chiều 47 3.2.2: Mạng ba chiều 48 3.3 : Sự phụ thuộc mơ men từ vào từ trường ngồi h nhiệt độ thấp 49 3.3.1: Mạng hai chiều .49 3.3.2: Mạng ba chiều 51 3.3.3:Áp dụng mơ hình Ising có tích phân trao đổi thăng giáng cho chuyển phameta từ 54 KẾT LUẬN 56 TÀI LIỆU THAM KHẢO 57 Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    LỜI CẢM ƠN Em xin bày tỏ lòng cảm ơn sâu sắc tới TS Bạch Hương Giang GS.TS Bạch Thành Cơng tận tình hướng dẫn động viên suốt trình thực luận văn để em hồn thành tốt đề tài “ Mơ hình Ising ứng dụng với chất sắt từ” Em xin gửi lời cảm ơn chân thành tới Th.S Nguyễn Văn Chinh –bộ môn Lý Sinh –Học viện Quân Y nhiệt tình giúp đỡ trình em thực tính tốn ngơn ngữ lập trình Matlab Em xin gửi lời cảm ơn PTN tính tốn KHVL, thầy mơn Vật lý chất rắn thầy cô Khoa Vật lý – Trường Đại học Khoa học Tự nhiên, thầy cô công tác trường Đại học Sư Phạm Hà Nội trang bị kiến thức chuyên môn tạo điều kiện thuận lợi giúp em hoàn thành luận văn Cám ơn đề tài NAFOSTED 103.02.2012.73 giúp đỡ tính tốn máy tính để thực thành công luận văn Cuối em xin gửi lời cảm ơn tới gia đình, bạn bè sát cánh, động viên suốt trình em học tập hồn thành luận văn tốt nghiệp Hà Nội ngày 01 tháng 08 năm 2014 Người thực Nguyễn Thị Kim Oanh Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    MỞ ĐẦU Vật liệu từ phát cách hàng nghìn năm ứng dụng tiêu biểu thời kì kim la bàn Chính la bàn tạo điều kiện cho ngành hàng hải phát triển, góp phần tìm lục địa Việc phát loại vật liệu với tính chất đặc biệt tạo bước ngoặt lớn tiến loài người Ngày nay, vật liệu từ ứng dụng rộng rãi thiết bị đại sống xung quanh điện thoại, la bàn, ổ cứng, ti vi… Song song với phát triển loại vật liệu từ phát triển ngành từ học nghiên cứu tính chất tượng vật liệu Một tượng quen thuộc nhận nhiều quan tâm nhà khoa học tượng chuyển pha vật liệu từ Các mơ hình lý thuyết giải thích tượng từ cách tượng luận đưa mơ hình lý thuyết trường phân tử Weiss (1907) giải thích tượng sắt từ, mơ hình Neel (1904-2000) giải thích tượng phản sắt từ feri từ ….Tuy nhiên việc phát triển mơ hình vi mơ để giải thích chất lượng tử tượng từ nhiệm vụ cần thiết Hiện nay, q trình từ hóa vật liệu có cạnh tranh tương tác trình cạnh tranh phản sắt từ sắt từ hợp chất, hợp kim perovskite pha tạp, sắt từ pha tạp …đang nhiều phịng thí nghiệm giới nghiên cứu Đặc biệt q trình từ hóa nhiệt độ thấp với ảnh hưởng từ trường vật liệu đa tinh thể Ví dụ cơng trình nghiên cứu R Mahendiran [15] khảo sát phụ thuộc mơ men từ vào trường ngồi vật liệucho kết lý thú tồn bước nhảy mô men từ nhiệt độ thấp gần độ Kelvin.Các nhảy bậc đường cong từ hóa nhiệt độ thấp mơi trường có tồn cạnh tranh tương tác khảo sát mơ hình Ising hai chiều [7] Mơ hình Ising (1920) mơ hình tốn học đơn giản cho tượng từ học thống kê Mơ hình bao gồm biến độc lập gọi spin nhận hai giá trị -1 Các biến spin xếp mạng tinh thể nút mạng tương tác với lân cận nhà khoa học Ersnt Ising (1900-1998) xây dựng với số lý thuyết nêu công Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   trình khoa học sở để giải thích cho q trình chuyển pha từ hệ từ pha tạp mạnh có cạnh tranh tương tác Trong luận văn này, tiếp tục phát triển lý thuyết khảo sát q trình từ hóa vật liệu sắt từ tác dụng trường khác cho hệ thực (hệ hai chiều, hệ ba chiều) trật tự so sánh kết lý thuyết với thực nghiệm Các tính tốn thực gần phương pháp trường trung bình dựa đẳng thức Callen khảo sát kết dựa phương pháp Monte Carlo phương pháp tính tốn lý thuyết kết hợp với mơ Phương pháp nghiên cứu - Dựa mơ hình Ising hệ thức Callen thực bước biến đổi giải tích theo học thống kê để xây dựng biểu thức mô men từ tỉ đối nút mạng phụ thuộc vào thông số nhiệt độ, từ trường ngồi đặt vào, xác suất thăng giáng… Từ sử dụng phần mềm hỗ trợ Mathlab tính tốn số thu kết phụ thuộc mô men từ tỉ đối vào nhiệt độ vào từ trường phù hợp với lý thuyết chuyển pha thực nghiệm đo - Ngoài sử dụng phương pháp Monte Carlo áp dụng cho số trường hợp cụ thể để thu kết tương tự so với phương pháp giải tích Cấu trúc luận văn Bên cạnh phần mục lục, mở dầu cấu trúc luận văn gồm ba phần sau: Chương 1: Tổng quan mơ hình Ising Chương 2: Ứng dụng mơ hình Ising cho trình chuyển pha Chương 3: Kết thảo luận Kết luận Danh mục tài liệu tham khảo phụ lục Luận văn tốt nghiệp  2   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mô hình Ising ứng dụng cho chất sắt từ    Danh mục hình có luận văn: Hình1.1: Đường cong từ trễ Hình 1.2: Sự thay đổi định hướng đám spin theo nhiệt độ Hình 1.3: Mơ hình Ising 1D Hình 1.4: Mơ hình Ising 2D Hình 1.5: Mơ hình Ising 3D Hình3.1: Đường cong từ nhiệt với tham số z=4, h=0, p=0.5 giá trị   0.6(1);   0.8(2);   0.98(3);   1(4);   1.001(5);   1.02(6);   1.1(7);   1.106(8) Hình 3.2 : Đường cong từ nhiệt với z=4, p=0.5 , h=0.002 giá trị    0.8(1);   1(2);   1.02(3);   1.106(4);   1.11(5);   1.19(6) Hình 3.3:Đường cong từ nhiệt với z=4, p=0.5, h=0.02 giá trị    0.8(1);   1(2);   1.02(3);   1.106(4);   1.19(5);   1.25(6) Hình 3.4 : Đường cong từ nhiệt với z=4 , h=0.2 , p=0.5 với giá trị    1.19(1);   1.3(2);   1.401(3);   1.43(4);   1.45(5);   1.5(6) Hình 3.5: Đường cong từ hóa với Z=4 , p=0.5 , delta=1.02 giá trị h Hình 3.6 : Đường biểu diễn phụ thuộc mô men từ mộtnút mạng m vào nhiệt độ t z=6,p=0.5, h=0 giá trị ∆   0.6(1);   0.8(2);   1(3);   1.2(4);   1.5(5);   1.56(6) Hình 3.7: Đường biểu diễn phụ thuộc mô men từ nút mạng mvào nhiệt độ t z=6, p=0.5, h=0.002 giá trị    0.6(1);   0.8(2);   1.2(3);   1.5(4);   1.56(5) Hình 3.8: Sự phụ thuộc nhiệt độ Curie vào xác suất thăng giáng p với z=4,   1.01   1.1 h=0 Hình 3.9 : Sự phụ thuộc nhiệt độphụ thuộc tc vào xác suất thăng giáng với z=4,  =1.15 h=0, h=1.2 h=1.5 Hình 3.10 :Đồ thị phụ thuộc (p-t) với z=6, h=0,  =1.005  =1.15 Hình 3.11:Sự phụ thuộc nhiệt độtc vào xác suất p với z=6,  =1.15, h=0, h=1.5, h=1.6 h=1.8 Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   Hình 3.12: Đồ thị biểu diễn phụ thuộc mômen từ vào từ trườngvới z=4, ∆ 1.03, t= 0.01 giá trị thăng giáng p=0.2 ; p=0.4; p=0.45 Hình 3.13: Đồ thị biểu diễn phụ thuộc mômen từ vào từ trường vớiz=4,p=0.2,t=0.01và giá trị ∆ tích phân trao đổi   1.02(1);   1.03(2);   1.04(3) Hình 3.14: Đường biểu diễn phụ thuộc mô men từ vào từ trường với z=4, p=0.2, ∆=1.02 giá trị nhiệt độ t=0.01, t=0.001, t=0.0001 Hình 3.15: Đường biểu diễn phụ thuộc mô men từ vào từ trường vớiz=6, t=0.01, p=0.2 giá trị ∆ tích phân trao đổi   1.02(1);   1.038(2);   1.04(3) Hình 3.16: Đồ thị (m-h) với z=6, Delta=1.04, t=0.01 p=0.1, p=0.3, p=0.5 Hình 3.17: Đường biểu diễn phụ thuộc mô men từ vào từ trường với z=6, p=0.1, ∆=1.04 giá trị nhiệt độ t=0.01, t=0.001, t=0.0001 Hình 3.18 : Đồ thị biểu diễn m theo h với z=4, t=0.01, delta=1.04,p=0.1 Hình 3.19:Đồ thị biểu diễn phụ thuộc mô ment từ tỷ đối m vào từ trường h z=4,t=0.01,delta=1.04,p=0.5 phương pháp Callen Monte Carlo Hình 3.20: Đồ thị so sánh lý thuyết thực nghiệm cho mơmen từ(trong đơn vị µB ) nút mạng theo từ trường Luận văn tốt nghiệp  4   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    CHƯƠNG 1: SẮT TỪ VÀ MÔ HÌNH ISING 1.1: Đặc điểm chất sắt từ - Sắt từ chất có từ tính mạnh có mơmen từ ngun tử lớn (ví dụ sắt 2,2 μB, Gd μB ) Nhờ tương tác trao đổi mômen từ nguyên tử định hướng song song với theo vùng (gọi đơmen từ tính) Mơmen từ đơn vị thể tích vùng gọi từ độ tự phát - có nghĩa chất sắt từ có từ tính nội khơng có từ trường Đây nguồn gốc tạo nên tính chất chất sắt từ -Hiện tượng từ trễ: Khi từ hóa khối chất sắt từ mơmen từ có xu hướng xếp trật tự theo hướng từ trường ngồi từ độ mẫu tăng dần đến độ bão hòa từ trường đủ lớn (khi mơmen từ hồn tồn song song với nhau) Khi ngắt từ trường khử từ theo chiều ngược, liên kết mômen từ đômen từ, mômen từ không bị quay trở lại trạng thái hỗn độn chất thuận từ mà giữ từ độ giá trị khác khơng Có nghĩa đường cong đảo từ khơng khớp với đường cong từ hóa ban đầu, từ hóa khử từ theo chu trình kín từ trường ngồi, ta có đường cong kín gọi đường cong từ trễ Và đường cong từ trễ, ta có đại lượng đặc trưng chất sắt từ sau: Hình1.1: Đường cong từ trễ Luận văn tốt nghiệp     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   1,0 1.19   1.3 0,8 1.401   1.43  1.45 0,6 m   1.5 0,4 0,2 0,0 t Hình 3.4: Đường cong từ nhiệt với z=4 , h=0.2 , p=0.5 với giá trị    1.19(1);   1.3(2);   1.401(3);   1.43(4);   1.45(5);   1.5(6) Các hình (3.2), (3.3), (3.4) khảo sát trình chuyển pha giá trị thăng giáng  tích phân trao đổi khoảng 0.8    1.5 từ trường tăng từ h=0.002 đến h=0.02 h=0.2 Ta thấy từ trường tăng lên, pha sắt từ ngày chiếm ưu so với pha phản sắt từ, vùng sắt từ ngày mở rộng Điểm xảy chuyển pha sắt từ sang thuận từ ngày “ nhòe” rộng giá trị nhiệt chuyển pha lúc xác định độ dốc đường cong mômen từ theo nhiệt độ Ngay nhiệt không tuyệt đối t=0, hệ tồn mômen từ dư tác dụng từ trường đặt vào làm cho số spin có xu hướng xếp theo hướng từ trường Từ trường lớn giá trị mơmen từ dư lớn Khi tiếp tục tăng nhiệt độ, lượng chưa đủ lớn để phá vỡ tác dụng từ trường ngồi có tác dụng làm cho liên kết spin lỏng lẻo có xu hướng quay dần theo hướng từ trường, giá trị mômen từ tăng đật đến giá trị cực đại Tiếp tục tăngnăng lượng nhiệt đủ lớn có tác dụng phá vỡ cấu trúc mơ hình, đa số spin xếp theo hướng hỗn độn khác làm cho mô men từ giảm Luận văn tốt nghiệp   44   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    Hình 3.5 mơ tả đường cong từ nhiệtkhi trường ngồi có giá trị khác 1.0 h=0 h=0.002 h=0.02 h=0.04 h=0.07 h=0.08 0.8 m 0.6 0.4 0.2 0.0 t Hình 3.5: Đường cong từ hóa với Z=4 , p=0.5 , delta=1.02 giá trị h Từ hình (3.5) ta thấy so sánh thấy rõkhi tăng giá trị h, vùng sắt từ mở rộng dần, giá trị mômen từ ban đầu nhiệt độ t=0 tăng gần điểm chuyển pha tồn mô men từ khác Tại h=0 ta tìm hai nhiệt độ chuyển pha Curie xác Với giá trị h , nhiệt độ chuyển pha xác định độ dốc đường cong mô men từ với trục nhiệt độ 3.1.2: Mạng ba chiều (z=6) Đường cong từ nhiệt khơng có trường ngồi h=0 có trường h  với p=0.5 độ thăng giáng tích phân trao đổi  có giá trị khác hình 3.6 hình 3.7 Các đường cong có dạng tương tự trường hợp hai chiều điểm chuyển pha cao (thấp) trường hợp chiều lớn (nhỏ hơn) so với trường hợp hai chiều lấy tham số p,  h Luận văn tốt nghiệp 45     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   1,0   0.6   0.8   0,8   1.2   1.5   1.56 m 0,6 0,4 0,2 0,0 t Hình 3.6 : Đường biểu diễn phụ thuộc mô men từ nút mạng m vào nhiệt độ t z=6 ,p=0.5, h=0 giá trị ∆   0.6(1);   0.8(2);   1(3);   1.2(4);   1.5(5);   1.56(6) 1,0  0.6  0.8   1.2 0,8   1.5   1.56 m 0,6 0,4 0,2 0,0 t Hình 3.7: Đường biểu diễn phụ thuộc mô men từ nút mạng m vào nhiệt độ t z=6, p=0.5, h=0.002 giá trị    0.6(1);   0.8(2);   1.2(3);   1.5(4);   1.56(5) Luận văn tốt nghiệp   46   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    3.2 Đường biểu diễn phụ thuộc nhiệt độ chuyển pha Curie vào xác suất p 3.2.1 : Mạng hai chiều Khảo sát phụ thuộc nhiệt độ chuyển pha vào xác suất thăng giáng p Ta có đường biểu diễn phụ thuộc sau (hình 3.8) trường hợp mạng hai chiều z=4: 2,5   0   1 2,0 tC 1,5 1,0 0,5 0,0 0,48 0,52 0,56 0,60 p Hình 3.8: Sự phụ thuộc nhiệt độ Curie vào xác suất thăng giáng p với z=4,   1.005   1.15 h=0 Trên đồ thị phụ thuộc nhiệt độ Curie vào xác suất thăng giáng (khi thăng giáng   vùng có cạnh tranh hai pha trạng thái sắt từ phản sắt từ) ta thấy có vùng giá trị p mà p  pc tồn hai giá trị của nhiệt độ chuyển pha TC1 TC2 với giá trị xác suất p Khi giá trị p nằm ngồi khoảng tìm nhiệt độ chuyển pha TC xảy trình chuyển pha Khi  tăng giá trị p tìm hai nhiệt độ chuyển tăng lên, với   1.005 xác suất tới hạn p 0.484,   1.15 xác suất p 0.545 Khi  tăng ta thấy nhiệt độ TC2 giảm đi, tức trình chuyển pha từ trạng thái sắt từ sang thuận từ xảy nhiệt dộ thấp tăng, vùng sắt từ bị thu hẹp dần ( phù hợp với kết đồ thị hình 3.1) thăng giáng phản sắt từ tăng cường Luận văn tốt nghiệp 47     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   Khi thay đổi từ trường ngồi đặt vào, ta có đồ thị biểu diễn mối liên hệ nhiệt độ chuyển pha xác suất thăng giáng 3.5 h=0 h=1.2 h=1.3 h=1.5 3.0 2.5 tC 2.0 1.5 1.0 0.5 0.0 0.5 0.6 0.7 0.8 0.9 1.0 p Hình 3.9 :Sự phụ thuộc nhiệt độ tc vào xác suất p                                            với z=4,  =1.15,h=0, h=1.2 h=1.5   3.2.2: Mạng ba chiều 3.5 3.0 1.005 1.15 2.5 tC 2.0 1.5 1.0 0.5 0.0 0.40 0.44 0.48 0.52 p Hình 3.10 : Đồ thị phụ thuộc (p-t )với z=6,h=0,   1.005   1.15 Luận văn tốt nghiệp   48   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    4.0 h=0 h=1.5 h=1.6 h=1.8 3.5 3.0 tC 2.5 2.0 1.5 1.0 0.5 0.0 0.4 0.5 0.6 0.7 0.8 0.9 1.0 p Hình 3.11: Sự phụ thuộc nhiệt độ tc vào xác suất p với z=6,  =1.15, h=0, h=1.5 ,h=1.6 h=1.8 3.3 : Sự phụ thuộc mơmen từ vào từ trường ngồi h nhiệt độ thấp 3.3.1: Mạng hai chiều Khảo sát thay đổi mơ men từ với trường ngồi trường hợp ∆ 1ở nhiệt độ thấp ta thu kết sau 1.0 0.8 m 0.6 0.4 p=0.2 p=0.4 p=0.45 0.2 0.0 0.00 0.04 0.08 0.12 0.16 h Hình 3.12: Đồ thị biểu diễn phụ thuộc mômen từ vào từ trường với z=4, ∆ 1.03, t= 0.01 giá trị thăng giáng p=0.2 ; p=0.4; p=0.45 Luận văn tốt nghiệp 49     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   0 m 0 0 1.02 1.03 1.04 0 0 h Hình 3.13: Đồ thị biểu diễn phụ thuộc mômen từ vào từ trường với z=4,p=0.2,t=0.01và giá trị ∆ tích phân trao đổi   1.02(1);   1.03(2);   1.04(3) Từ hình (3.12) (3.13) ta thấy xác định mômen từ tăng đột ngột vài giá trị hc từ trường, tượng nhảy bậc đặc trưng cho trình quay đột ngột đám spin phản sắt từ trở thành đám sắt từ Do hệ có tương tác sắt từ (FM) phản sắt từ(AF) ta giả thiết ban đầu mạng có đám spin phản sắt từ sắt từ Khi tăng từ trường đặt vào đám spin phản sắt từ dần quay theo hướng từ trường dọc theo hướng đám sắt từ (FM) làm mômen từ tăng dần, tiếp tục tăng từ trường đến giá trị tới hạn hc số lớn đám phản sắt từ đột ngột đảo hướng song song với từ trường làm mơmen từ tăng nhảy bậc Q trình q trình từ hóa loại I theo lý thuyết nhiệt động học vật lý thống kê Khi ta thay đổi độ thăng giáng  tích phân trao đổi từ trường tới hạn xảy tượng nhảy bậc tăng Hình 3.15 cho thấy mạng spin hai chiều ta thấy có hai bước nhảy mơmen từ hc1  0.032 hc  0.071 p=0.2,   1.02 bước nhảy sắc nét nhiệt độ thấp Luận văn tốt nghiệp   50   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    1.0 0.8 m 0.6 0.4 0.2 t=0.01 t=0.001 t=0.0001 0.0 0.00 0.02 0.04 0.06 0.08 0.10 0.12 h Hình 3.14: Đường biểu diễn phụ thuộc mơmen từ vào từ trường với z=4, p=0.2, ∆=1.02 giá trị nhiệt độ t=0.01, t=0.001, t=0.0001 3.3.2: Mạng ba chiều Đối với mạng ba chiều trình từ hóa loại I xảy nhiệt độ thấp với bước nhảy (xem hình 3.16) Khi độ thăng giáng tích phân trao đổi  tăng đường cong từ hóa từ trường tới hạn có xu hướng dịch phía từ trường cao 0 m   1.02   1.038   1.04 0 0 0 0 0 h Hình 3.15: Đường biểu diễn phụ thuộc mômen từ vào từ trường với z=6, t=0.01, p=0.2 giá trị ∆ tích phân trao đổi   1.02(1);   1.038(2);   1.04(3) Luận văn tốt nghiệp 51     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   1.0 0.8 m 0.6 0.4 p=0.1 p=0.3 p=0.5 0.2 0.0 0.00 0.05 0.10 0.15 0.20 0.25 0.30 h Hình 3.16: Đồ thị (m-h) với z=6, Delta=1.04, t=0.01 vàp=0.1, p=0.3, p=0.5 Hình 3.16 cho thấy   xác suất p tương tác sắt từ bé tăng độ lớn mơmen từ trường bé chưa bão hòa tăng giá trị từ tới hạn, số bước nhảy mômen từ cũ hay giá trị từ trường tới hạn không phụ thuộc xác suất thăng giáng 1.0 0.8 m 0.6 0.4 0.2 0.0 0.00 t=0.01 t=0.001 t=0.0001 0.05 0.10 0.15 0.20 0.25 0.30 h Hình 3.17: Đường biểu diễn phụ thuộc mômen từ vào từ trường với z=6, p=0.1, ∆=1.04 giá trị nhiệt độ t=0.01, t=0.001, t=0.0001 Luận văn tốt nghiệp   52   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    Sử dụng phương pháp Monte Carlo để tính thuộc thuộc mơmen từ vào từ trường ta thu kết tương tự tính tốn theo hệ thức Callen (xem hình 3.18) lập trình Mathlab 1.0 Callen Monte Carlo 0.8 m 0.6 0.4 0.2 0.0 0.00 0.05 0.10 0.15 0.20 h Hình 3.18 : Đồ thị biểu diễn m theo h với z=4, t=0.01,   1.04 ,p=0.1bằng phương pháp Monte Carlovà phương pháp Callen-Mathlab 1.0 0.8 m 0.6 0.4 0.2 Callen Monte Carlo 0.0 0.00 0.05 0.10 0.15 0.20 h Hình 3.19:Đồ thị biểu diễn phụ thuộc mômen từ tỷ đối m vào từ trường h z=4,t=0.01,delta=1.04,p=0.5bằng phương pháp Callen Monte Carlo Luận văn tốt nghiệp 53     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   3.3.3:Áp dụng mơ hình Ising có tích phân trao đổi thăng giáng cho chuyển pha meta từ Sử dụng lý thuyết trình từ hóa loại I cho mẫu Ising trật tự chúng tơi so sánh tính tốn lý thuyết kết thực nghiệm Y.Y.Wu et al [12] cho chuyển pha meta từ Pr 0.5 Ca0.5 Mn0.97 Ga0.03 O3 đo nhiệt độ T=2.5 (K) ly thuyet thuc nghiem T=2.5 K 2.0 m 1.5 1.0 0.5 0.0 H(T) Hình 3.20: Đồ thị so sánh lý thuyết thực nghiệm cho mômen từ (trong đơn vị µB ) nút mạng theo từ trường với z=6, t=0.01,   1.04 , p=0.34 Bảng số liệu so sánh bước nhảy mô men từ từ trường tới hạn m(B ) H c (T ) Lý thuyết Thực nghiệm Lý thuyết Thực nghiệm 0.80542 0.80617 1.93973 3.41884 0.49491 0.2389 4.84954 5.24766 0.41722 0.44822 7.43599 7.23436 Kết lý thuyết minh họa tốt ba bước nhảy quan sát thực nghiệm Nhìn chung độ lớn bước nhảy mômen từ tương đối phù hợp Hai từ trường tới hạn cao có giá trị phù hợp với cịn từ trường tới hạn thứ có Luận văn tốt nghiệp   54   LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ    khác biệt đáng kể lý thuyết thực nghiệm Điều nguyên nhân trạng thái ban đầu mẫu đo phụ thuộc mạnh vào điều kiện thực nghiệm Luận văn tốt nghiệp 55     LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình Ising ứng dụng cho chất sắt từ   KẾT LUẬN Luận văn khảo sát mơ hình Ising cho hệ spin với tích phân trao đổi thăng giáng có giá trị dương (trao đổi sắt từ), âm (trao đổi phản sắt từ) phụ thuộc vào độ thăng giáng ∆ Luận văn thu số kết sau: 1.Khi ∆ có cạnh tranh tương tác sắt từ phản sắt từ dẫn tới đường cong mô men từ phụ thuộc nhiệt độ khác vùng nhiệt độ hữu hạn (có hai điểm chuyển pha TC1 thấp TC2 cao) tượng gọi tượng chuyển pha từ trở lại (Reentrant phase transition) TC1 điểm ứng với trình chuyển từ pha sắt từ sang phản sắt từ giảm nhiệt độ TC2 điểm ứng với chuyển pha sắt từ sang thuận từ tăng nhiệt độ Hiện tượng Chuyển pha trở lại xảy hệ chiều Cho giá trị ∆ trường ngồi (h=0) tồn giá trị tới hạn pc p> pc xảy tượng chuyển pha từ trở lại, h tùy thuộc vào độ lớn h tồn hai giá trị tới hạn 2.Trong trường hợp có chuyển pha từ trở lại (trường hợp độ thăng giáng   ) nhiệt độ thấp gần hệ có q trình chuyển pha loại từ trường với xuất bước nhảy bậc mômen từ (độ lớn mô men từ tăng đột ngột) giá trị từ trường tới hạn hc khác Số bước nhảy độ lớn bước nhảy phụ thuộc vào xác suất thăng giáng số chiều mạng (thể qua số lân cậ gần nhất) Mạng (3) chiều có số bước nhảy xác xuất tương tác sắt từ p