Mẫu Ising và một số ứng dụng

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ĐỒNG MINH SƠN HUYỀN TRANG MẪU ISING VÀ MỘT SỐ ỨNG DỤNG TÓM TẮT LUẬN VĂN THẠC SĨ KHOA HỌC Chuyên ngành: Vật lý chất rắn Mã số: 60.44.01.04 HÀ NỘI - 2016 Cơng trình đƣợc hồn thành tại: Trƣờng Đại Học Khoa Học Tự Nhiên Ngƣời hƣớng dẫn khoa học: TS Bạch Hƣơng Giang GS TS Bạch Thành Công Phản biện 1: TS Phạm Ngọc Anh Huy Phản biện 2: TS Nguyễn Tiến Cường Luận văn bảo vệ trước Hội đồng chấm luận văn Thạc sĩ tại: Trường Đại Học Khoa Học Tự Nhiên Có thể tìm hiểu luận văn tại: - Thư viện Đại học Quốc gia Hà Nội TÓM TẮT LUẬN VĂN MẪU ISING VÀ MỘT SỐ ỨNG DỤNG Mẫu Ising mơ hình tốn học đơn giản mơ tả tượng học thống kê Mục đích ban đầu mẫu Ising, chủ đề luận án tiến sĩ Ising giải thích cấu trúc tính chất chất sắt từ Ở đây, Ising cố gắng giải thích số liệu thực nghiệm quan sát vật liệu sắt từ cách sử dụng mơ hình người thầy Lenz đề xuất năm 1920 Kể từ mẫu Ising cho phép đơn giản hóa tương tác phức tạp ứng dụng thành cơng lĩnh vực khoa học Thống kê cho thấy khoảng năm từ 1969 đến 1997 có 12.000 báo mẫu Ising công bố, số tăng khơng ngừng Có thể kể đến biến thể mẫu Ising giúp hiểu chất sâu xa nhiều tượng lý sinh đường cong bão hòa Hemoglobin, tốc độ phản ứng ban đầu enzyme allosteric hay mơ hình mạng thần kinh đặc trưng quan trọng màng lipid Ngồi mơ hình Ising ứng dụng lĩnh vực khác kinh tế học (nghiên cứu ảnh hưởng kinh tế - xã hội đến số kinh doanh, phân tích chuỗi thời gian tài thị trường kinh doanh), xã hội học (hành vi xã hội cá nhân thay đổi để phù hợp với hành vi cá nhân khác vùng lân cận họ) hay ngôn ngữ học (sự thay đổi ngôn ngữ), … Đối với ngành Vật lý, nhiều thập kỷ qua, mẫu Ising chủ yếu áp dụng để nghiên cứu vật liệu từ Gần đây, phát triển kỹ thuật màng mỏng mở nhiều hướng nghiên cứu mẫu Ising hai chiều, mẫu Ising cho màng mỏng sắt điện sắt từ trường Những hướng nghiên cứu thu hút nhà vật lý lý thuyết vật lý thực nghiệm Về mặt lý thuyết giúp xác định tính chất vĩ mô hệ vật chất Về mặt thực nghiệm, ứng dụng nhiều lĩnh vực cơng nghệ khác nhau, chẳng hạn lưu trữ liệu, xúc tác, điện tử, tạo bước ngoặt lớn tiến khoa học công nghệ Trong luận văn này, tiếp tục nghiên cứu phát triển mơ hình Ising mặt lý thuyết, ứng dụng việc khảo sát tham số nhiệt động mẫu Ising chiều, mẫu Ising trường dọc, cho màng mỏng có trật tự so sánh kết lý thuyết với thực nghiệm cho điểm Curie màng mỏng sắt điện Các tính tốn thực gần phương pháp trường trung bình lý thuyết Landau cho mẫu Ising đồng thời so sánh với kết dựa phương pháp Monte – Carlo cho mẫu Ising 2D (màng mỏng lớp) Mục tiêu nghiên cứu - Xây dưṇ g lý thuyết Landau cho màng mỏng có chuyển pha trật tự - khơng trật tự, tham số xuấ t phá t từ mô hiǹ h Ising vớ i spin tù y ý theo phương phá p trườ ng trung biǹ h - So sá nh phương phá p trườ ng trung binh ̀ phương phá p Monte - Carlo cho mà ng đơn lớ p - Áp dụng lý thuyết trường trung bìn h giả i thích sư ̣ phu ̣ thuôc̣ củ a nhiêṭ đ ộ Curie vào độ dày màng Phƣơng pháp nghiên cứu - Dựa mơ hình Ising lý thuyết trường trung bình, lý thuyết chuyển pha Landau 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 cho lượng tự do, độ từ hóa nhiệt dung hệ spin đặc trưng cho hệ có trật tự xa, khảo sát tượng chuyển pha xảy hệ khơng vàkhi có trường ngồi tác dụng Từ sử dụng phần mềm hỗ trợ tính tốn số thu kết áp dụng để phân tích kết thực nghiệm khác cho đại lượng nhiệt động tương ứng - Sử dụng phương pháp Monte – Carlo áp dụng cho số trường hợp màng lớp (mẫu 2D)có trật tự xa để so sánh với phương pháp giải tích gần trường trung bình Cấu trúc luận văn Bên cạnh phần mục lục mở đầu, cấu trúc luận văn gồm ba phần sau: Chương 1: Mẫu Ising lý thuyết chuyển pha Landau Chương 2: Áp dụng mẫu Ising lý thuyết chuyển pha Landau cho màng mỏng có trật tự xa Chương 3: Tính tốn Monte – Carlo cho mẫu Ising 2D (màng mỏng lớp) Kết luận CHƢƠNG 1: MẪU ISING VÀ LÝ THUYẾT CHUYỂN PHA LANDAU 1.1 Mẫu Ising hai trạng thái (S = 1/2) Mẫu Ising cách biểu diễn lý thuyết đơn giản cho tượng sắt từ, nhiên dùng để mơ tả hệ có trật tự khác trật tự sắt từ, sắt điện, hợp kim có trật tự,… Xuất phát tốn học mẫu Ising: Coi nút mạng có spin σivà có hai định hướng lên (spin up σi = +1) xuống (spin down σi = -1) Đối với vật liệu có trật tự khác nhau, spin đặc trưng cho độ phân cực từ (vật liệu từ) hay độ phân cực điện (vật liệu sắt điện) hay tỷ số nồng độ thành phần hợp kim đôi trật tự 1.2 Lý thuyết chuyển pha Landau 1.2.1 Lý thuyết chuyển pha Landau trường ngồi Lý thuyết chuyển pha Landau lý thuyết tượng luận cho chuyển pha trật tự - trật tự tham số trật tự η pha trật tự khác khơng khơng pha trật tự Biểu thức lượng tự Landau gần điểm chuyển pha có dạng: F = a(T – TC)� � +β + … (ở đây, η tham số trật tự độ từ hóa trung bình, trật tự xa, mật độ dòng siêu lỏng mật độ dòng siêu dẫn) Đặt: A = mω0 = a(T - TC) ; B = β 2− ; C = γ > ta biểu thức củad c' lượng tự hệ với trật tự đặc trưng tham số η chủ yếu: F =A + sgn B η (B)η η6 +C Đặt: α F0 = B ; (1.13) AC ; B C2 f = = F F0 ta biểu thức lượng tự không thứ nguyên: α f = ( B) ξ 2+ ξ 4+ sgn ξ (1.17) 6 Biểu thức (1.17) thuận lợi cho việc khảo sát phụ thuộc lượng tự vào tham số trật tự không thứ nguyên ξ cách thay đổi giá trị α Hình 1.1: Đường biểu diễn phụ thuộc lượng tự không thứ nguyên f theo tham số trật tự không thứ nguyên ξ với giá trị khác tham số α 1.2.2 Lý thuyết chuyển pha Landau có trường ngồi Biểu thức lượng tư F có dạng: F =A + sgn B η ( B )η η +C 6 (1.24) − h 'η B η ; 12 ⇒ Đặt ξ = B ⇒ η0 η =   η = η = ξ Với α = AC B  C  ξη0 ; h = h'C 32 B 52 ; B F0 = C C ta thu lượng tự không thứ nguyên: f= ( B) α ξ ξ 2+ + sgn ξ − hξ Trong đó: ξ độ từ hóa tỉ đối khơng thứ ngun η độ từ hóa η0 độ từ hóa bão hòa (1.27) Dựa vào biểu thức (1.27) ta khảo sát phụ thuộc lượng tự vào tham số trật tự khơng thứ ngun ξ Hình 1.2: Đường biểu diễn phụ thuộc tham số trật tự không thứ nguyên ξ theo từ trường h, tham số α = 0.1, sgn(B) = -1 Hình 1.2 đường biểu diễn phụ thuộc tham số trật tự không thứ nguyên ξ vào từ trường h tham số α = 0.1 vàsgn(B) = -1.Trạng thái giả bền (đường nét đứt) tồn trường ngồi bé 1.3 Lý thuyết trƣờng trung bình cho mẫu Ising spin -1/2 trƣờng dọc Xét hệ spin mạng tinh thể với Hamiltonian Ising sau: H=-J Nq/2 ∑ σi σj N - h ∑σ k =1 k (1.32) đó: σi = ±1/2 giá trị hình chiếu spin nguyên tử từ tính nút mạng thứ i Biến spin σi = ±1/2 xem thành phần z spin -1/2 chuẩn (ma trận z Pauli) với giá trị riêng σ i = ±1/2 heo bậc bé độ từ hóa K Curie dùng ) q JTkhi 4kC gần điểm Curie có dạng: F(h h a = 0, m→0) = N (A0+ A2m + i A4 m + … ) khai triển (1.47.a) r f i ể n n ă n g ( b ) Trong công thức trên: l g − β ự d o lượn g tự do: (1 48)  qJ 1  2  qJ  β qJ  A =24 (1.49)   Tại nhiệt độ chuyển pha loại 2, T = TC , β = βC = 1/kBTC ; A2(βC) = Từ (1.48) ta lại nhận t Trong lý thuyết chuyển pha Landau, chuyển pha liên tục bị ràng buộc điều kiện A2 = 0, A4 > Và theo lý thuyết điểm chuyển pha spin Ising -1/2 giá trị nhiệt độ chuyển pha phụ thuộc vào số phối vị q độ lớn tích phân trao đổi J spin (S = 1/2) t cho A0 = −kBT ln = q J ợ au = F/N = A0+ A2m + A4m + … A n Land Năng lượng tự tính hạt là: t B biểu thức (1.44) cho nhiệt độ CHƢƠNG 2:ÁP DỤNG MẪU ISING VÀ LÝ THUYẾT CHUYỂN PHA LANDAU CHO MÀNG MỎNG CÓ TRẬT TỰ XA 2.1 Mẫu Ising cho màng mỏng có trật tự xa, lời giải lý thuyết trƣờng trung bình Xét màng mỏng có trật tự xa gồm n lớp song song xếp chồng lên nhau, lớp có N spin Hamiltonian Ising cho hệ spin màng mỏng cấu trúc lập phương gồm n lớp viết sau: H=- ∑ Jν j ,ν ' j 'Sν j Sν j ' S - hμ ∑ z z ν j ,ν ' j ' Với: ν j (2.1) ν j Sνz toán tử spin nút rνj j μ magneton Bohr Jν j,ν ' j ' spin tích phân trao đổi S z ν S z hai nút mạng rνjvà rν’j’ ν' Trong gần trường trung bình, Hamiltonian biểu thức (2.1) viết dạng sau:  Hm f = − ∑ ∑ Jν j,ν ' j ' z  z z Sν ' S − νj j' hµ ∑ S ν j (2.3) ν j ν 'j'  ν j ký hiệu trung bình thống kê lấy theo Hamiltonian H m f Ta thu biểu thức tổng thống kê hệ spin gần trường trung bình sau: sh n ν Z= ∑ ν  α  S  +     =1 (2.9)   sh  αν    Năng lượng tự hệ spin màng mỏng là: N F = −k T ln Z + J (0) m m ∑ B νν ' νν ' Trong đó: β = (với kB số Boltzmann, T nhiệt độ tuyệt đối) r Jνν ' (0) ảnh Fourier tích phân trao đổi k = k BT Năng lượng tự màng mỏng spin n lớp phép gần trường trung bình là: N F= N J (0) m m − n ∑ νν ν ' ν n    sh  αν 2S + ln ∑ ν =1 ν ' β '  sh  αν   (2.10)   2 z Giá trị trung bình thống kê momen từ nút mạng νj với spin Sν j gần trường trung bình là: m = cth α ν  S+    ν S + − cth    αν     phép      Trong đó:   �� βµ h   + = µ  uur Jνν ' (0) = ∑ j Jνν '   ν ' (0)  ∑ ν Jνν ' (R j ) 2.2 Phƣơng trình xác định điểm Curie cho màng mỏng trật tự lý thuyết trƣờng trung bình Phương trình xác định điểm Curie cho màng mỏng có trật tự xa: kBTC JS 2S(S+1)  J 2+ P = π cos    n +1   JS Đặt η = JP JS (2.22.a) tham số đặc trưng cho dị hướng tương tác trao đổi spin lân cận lớp ( J S ) lớp ( J P ) ta được: kBTC JS = cos  2S(S+1)  π   +η (2.22.b)  n +1  Hình 2.2: Đường biểu diễn phụ thuộc nhiệt độ Curie ‫ح‬C = kBTC vào số lớp màng JS Hình 2.3: Đồ thị so sánh kết lý thuyết trường trung bình thực nghiệm phụ thuộc nhiệt độ Curie sắt điện perovskite PbTiO3 vào độ dày màng mỏng Hình 2.2 biểu diễn phụ thuộc nhiệt độ Curie vào độ dày màng mỏng (số lớp spin n) S = Đồ thị cho thấy trao đổi dị hướng tăng η lớn nhiệt độ Curie tăng theo xu hướng chung nhiệt độ Curie giảm độ dày màng giảm Hình 2.3 thể kết tính tốn cho điểm Curie theo lý thuyết trường trung bình so sánh với kết thực nghiệm cho điểm Curie sắt điện xác định từ thực nghiệm cho màng mỏng có độ dày khác Trong đó, đường liền nét đường biểu diễn tỉ số kBTC theo số lớp màng n Các điểm rời rạc màu đỏ điểm thực nghiệm JS Kết lý thuyết thực nghiệm tương đối phù hợp.Kết tốt cho mơ hình Ising trường dọc ngang 2.3 Khai triển Landau cho màng mỏng có trật tự mơ tả mơ hình Ising Khai triển Landau cho lượng tự có dạng: F = N Jν (0) ∑ ν ν' m ν ν mν ' N S +      '  + 14 + −+ + + )S  )1 3 l ν  n α  S ) (4 ( α 2( SS S + ++ 1 +( + …  )) αS − ( n β ν   ∑ ν =1    Năng lượng tự tính nút mạng (hay spin): f = F ∑ n β J ' (0) mν ' ν mν νν ν ∑ ' ν ν 90 − 1 ( S +1)α S − Đặt: A = S(S+1) (2S + 2S ν  ν 6! 63   vν ν ν ν f ν nν ν ν năn ν gν ν lượ ng t tính mộ spi n tron g m l m ng spi n f ó g: f = − ln 2S ( Nn nβ +1) + 2n (.0 b 23 ) ν + 1) α  S(S +1) +  (4S   + + ) S2 +  α +  (2.28)  −S ( S +1 )(2.29)  mν ν 2nβ  B  S(S+1) = (2S ν + 2S nβ  C = S ( S−+1) ν (4S 6! 63   nβ    α ν   f 1( + − )A lm α nm  + 1) ∑ 90  J ( 3  + 4S + ) +   B S α  Khi biểu thức lượng tự tính nút mạng có dạng: n + + 1+ )C +α + CĐặ Aα c ν β = h 04 ệ khơ cóβ J s tr m ố h A ν , B ν , C ν đ ợ c c h o b i ( ) ta nhận lại phương trình (2.15) xác định điểm Curie cho màng mỏng lớp kBT C J = 4S ( S +1) (2.31) 2.4 Nhiệt dung đẳng tích cho màng mỏng Cơng thức chung cho nhiệt dung đẳng tích vật liệu mà Hamiltonian H cho biểu thức sau:  H T  H £ V =  vớiθ = kB T = kB (2.32) Công thức chung cho nhiệt dung đẳng tích màng mỏng spin gồm n lớp CV θ = −ν kB νν ' n νν θ ν mν ' ∂ θ + ∑ J ( ) ∂αν ∂ ∂mν ' m  2 ∂ m ∂θ 2θ ∂ θ  +  ∑  n BS (αν ) ∂θ ν ' ∂  ∂ α )  Bα + ( 2'S ν   ∑ ∂θ n ν  α  + BS (αν ) ∂θ  ν (2.36) ν    2.5 Tính tốn số cho màng mỏng đơn lớp có trật tự 2.5.1 Phương trình xác định phụ thuộc vào nhiệt độ tham số trật tự Biểu thức độ từ hóa tính cho màng mỏng lớp sau: m  = S +  ct h α  S +  4m  − cth  2m     2  ‫ح‬      (2.39)  ‫ح‬     VkBT J =nhiệt độ khơng thứ ngun đo đơn vị tích phân trao đổi J lân cận gần Hình 2.4: Đường biểu diễn phụ thuộc momen từ nút mạng m vào nhiệt độ không thứ nguyên ‫ = ح‬kBT J TÀI LIỆU THAM KHẢO A- Tài liệutiếng Việt [1] Nguyễn Thị Kim Oanh, luận văn Thạc sĩ khoa học “Mơ hình Ising ứng dụng cho chất sắt từ”, ĐH Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, 2014 B- Tài liệu tiếng Anh [2] Clarendon Press, Statistical Mechanics of Phase Transitions, Oxford,1992 [3] C Kittel, Introduction to Solid State physics, chapter 16, eighth edition, John Wiley & Sons, Inc 2005 [4]Dillon D Fong, G Brien Stephenson, Stephen K Streiffer,Jeffrey A Eastman, Orlando Auciello, Paul H Fuoss, Carol Thompson ,Science304 (2004) 1650 [5] D K Khudier, Nabeil A Fawaz, Two dimensional Ising model application with Monte Carlo method7(2013) [6] J A Krumhansl, Solid State communication 84 (1992) 251 [7] J Borowska, L Lacinska, Jour of Appl Math Comput.Mech.14(2015) 11 [8]J Strecka, M Jascur, Acta physica slovaca65 (2015) 235 [9] K Binder and D W Heermann, Monte Carlo Simulation in Statistical Physics, Springer, Berlin, 1997 [10] M Hjorth Jensen, Computational physics, University of Oslo, 2003 [11] N Metropolis, A W Rosenbluth, M N Rosenbluth, A H Teller and E Teller,Journal of Chemical Physics, 1953 [12]Nguyen Tu Niem, Bach Huong Giang, Bach Thanh Cong,Journal of Science: Advanced Materials and Devices1 (2016) 531 [13] Onsager Lars, Physical Review, Series II, 65 (3–4): 117–149, (1944) [14] Sergio A Cannas, Pablo M Gleiser, Francisco A Tamarit, Two dimentional Ising model with long-range competing interactions, Transworld Research Network 2004 [15]S V Tyablikov, Method in the quantum theory of magnerism, Plenumpress NewYork 1967 [16] W Nolting, A Ramakanth, Quantum theory of Magnetism, Springer 2009 ... gia Hà Nội TÓM TẮT LUẬN VĂN MẪU ISING VÀ MỘT SỐ ỨNG DỤNG Mẫu Ising mô hình tốn học đơn giản mơ tả tượng học thống kê Mục đích ban đầu mẫu Ising, chủ đề luận án tiến sĩ Ising giải thích cấu trúc... đây, Ising cố gắng giải thích số liệu thực nghiệm quan sát vật liệu sắt từ cách sử dụng mơ hình người thầy Lenz đề xuất năm 1920 Kể từ mẫu Ising cho phép đơn giản hóa tương tác phức tạp ứng dụng. .. Vật lý, nhiều thập kỷ qua, mẫu Ising chủ yếu áp dụng để nghiên cứu vật liệu từ Gần đây, phát triển kỹ thuật màng mỏng mở nhiều hướng nghiên cứu mẫu Ising hai chiều, mẫu Ising cho màng mỏng sắt điện