Nghiên Cứu Mô Phỏng Sự Chuyển Pha Trong Một Số Hệ Spin Gián Đoạn

Nghiên cứu mô phỏng sự chuyển pha trong các hệ spin gián đoạn, khám phá các hiện tượng vật lý và ứng dụng trong công nghệ hiện đại.

• Tác giả: Lương Minh Tuấn
• Người hướng dẫn: PGS. Đào Xuân Việt
• Trường học: Trường Đại học Bách khoa Hà Nội
• Chuyên ngành: Khoa học Vật liệu
• Đề tài: Nghiên Cứu Mô Phỏng Sự Chuyển Pha Trong Hệ Spin Gián Đoạn
• Loại tài liệu: luận án tiến sĩ
• Năm xuất bản: 2021
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• MỤC LỤC
• DANH MỤC CÁC KÍ TỰ VIẾT TẮT
• DANH MỤC BẢNG BIỂU
• DANH MỤC HÌNH VẼ, ĐỒ THỊ
• MỞ ĐẦU
• 0.1. Giới thiệu và lý do chọn đề tài
• 0.2. Mục tiêu của Luận án
• 0.3. Nội dung nghiên cứu
• 0.4. Đối tượng nghiên cứu
• 0.5. Phương pháp nghiên cứu
• 0.6. Ý nghĩa khoa học và thực tiễn của đề tài. Những đóng góp mới của Luận án
• 0.7. Cấu trúc của Luận án
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Vật liệu từ
• 1.2. Giới thiệu và lịch sử của từ học
• 1.3. Bản chất từ tính của vật liệu
• 1.4. Phân loại vật liệu từ
• 1.5. Pha và chuyển pha trong vật liệu từ
• 1.6. Pha và chuyển pha
• 1.7. Phân loại chuyển pha
• 1.8. Hiện tượng chuyển pha trong các mô hình
• 1.8.1. Mô hình 2D Ising
• 1.8.2. Mô hình 2D XY
• 1.8.3. Mô hình q-state clock
• 1.8.4. Mô hình 2D XYhq
• 1.9. Tổng quan nghiên cứu cho các mô hình
• 1.9.1. Mô hình 2D q-state clock
• 1.9.2. Mô hình 2D XYhq
• 1.10. Bán dẫn từ pha loãng - DMS
• 1.11. Lý do và nội dung nghiên cứu
• 2. CHƯƠNG 2: PHƯƠNG PHÁP MÔ PHỎNG VÀ CÁC ĐẠI LƯỢNG VẬT LÍ
• 2.1. Phương pháp Monte Carlo
• 2.2. Phương pháp lấy mẫu
• 2.2.1. Lấy mẫu đơn giản
• 2.2.2. Lấy mẫu quan trọng
• 2.3. ‘Ergodic’ và ‘detailed balance’
• 2.4. Thuật toán Metropolis
• 2.5. Thuật toán Wolff
• 2.6. Thuật toán Over-relaxation
• 2.6.1. Thuật toán general Over-relaxation
• 2.7. Thuật toán Parallel-tempering
• 2.8. Chương trình chính
• 2.9. Hệ vật liệu hai chiều
• 2.10. Hệ bán dẫn từ pha loãng
• 2.11. Phân tích số liệu. Kiểm tra cân bằng
• 2.12. Các đại lượng vật lí
• 2.12.1. Nhiệt dung riêng
• 2.12.2. Tham số Binder
• 2.12.3. Đạo hàm của tham số Binder
• 2.12.4. Mô đun xoắn
• 2.12.5. Đạo hàm của mô đun xoắn
• 2.12.6. Tỉ số chiều dài tương quan
• 2.13. Phương pháp phân tích sai số
• 2.13.1. Cơ sở của phương pháp
• 2.13.2. Độ lệch xuất hiện như thế nào
• 2.13.3. Ước tính sai số khi có độ lệch bằng sử dụng Jackknife
• 2.14. Lí thuyết tỉ lệ kích thước hữu hạn
• 2.15. Tính nhiệt độ chuyển pha bậc 2
• 2.16. Tính nhiệt độ chuyển pha KT
• 3. CHƯƠNG 3: NGHIÊN CỨU CHUYỂN PHA TRONG MÔ HÌNH q-STATE CLOCK BẰNG PHƯƠNG PHÁP MÔ PHỎNG MONTE CARLO
• 3.1. Mô hình 2D q-state clock
• 3.2. Các tham số mô phỏng cho mô hình 2D q-state clock
• 3.3. Kết quả mô phỏng
• 3.4. Mô đun xoắn và đạo hàm mô đun xoắn
• 3.5. Tỉ số chiều dài tương quan và tham số Binder
• 3.6. Đạo hàm của tham số Binder
• 3.7. Tổng hợp và so sánh
• 4. CHƯƠNG 4: NGHIÊN CỨU PHÂN TÍCH CÁC CHUYỂN PHA VÀ XÂY DỰNG GIẢN ĐỒ PHA CHO MÔ HÌNH XYhq
• 4.1. Mô hình 2D XYhq
• 4.2. Các tham số mô phỏng
• 4.3. Kết quả mô phỏng
• 4.4. Nghiên cứu ảnh hưởng của h3 đến sự chuyển pha của mô hình XY hai chiều
• 4.5. Các tham số mô phỏng
• 4.6. Giản đồ pha của mô hình 2D XY với trường phá vỡ đối xứng 5-hướng
• 4.7. Các tham số mô phỏng
• 4.8. Nhiệt độ chuyển pha như hàm logarit của trường tinh thể bất đẳng hướng trong mô hình XYh6
• 4.9. Các tham số mô phỏng
• 4.10. Giản đồ pha T – h6 đầy đủ của mô hình XYh6
• 5. CHƯƠNG 5: NGHIÊN CỨU MÔ PHỎNG MONTE CARLO VỀ TRẬT TỰ SẮT TỪ Ở NHIỆT ĐỘ PHÒNG TRONG CHẤT BÁN DẪN TỪ PHA LOÃNG
• 5.1. Bán dẫn từ pha loãng
• 5.2. Các tham số mô phỏng
• 5.3. Mô hình spin liên lục
• 5.4. Mô hình spin rời rạc
• 5.5. Cân bằng nhiệt
• 5.6. Kết quả mô phỏng
• 5.7. Mô hình Heisenberg
• 5.8. Mô hình spin rời rạc
• KẾT LUẬN CHUNG VÀ KIẾN NGHỊ
• DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ CỦA LUẬN ÁN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Nghiên Cứu Chuyển Pha Tổng Quan Hệ Spin Gián Đoạn

Nghiên cứu về vật liệu từ luôn thu hút sự quan tâm lớn từ cả lý thuyết và thực nghiệm. Đặc biệt, các vật liệu từ màng mỏng và vật liệu từ pha loãng đang là tâm điểm do sự khác biệt về từ tính so với vật liệu từ khối. Câu hỏi đặt ra là liệu màng mỏng từ có duy trì trật tự từ, hay tồn tại các trật tự hoặc chuyển pha khác? Đối với vật liệu từ pha loãng, cơ chế tương tác nào tạo nên từ tính khi chỉ có một lượng nhỏ spin? Các nhà khoa học đã xác định màng mỏng từ có pha trật tự (sắt từ), pha mất trật tự (thuận từ) và pha giả trật tự, cùng với bằng chứng về chuyển pha bậc 2, bậc 1 và Kosterlitz-Thouless. Chuyển pha sắt từ ở nhiệt độ phòng cũng đã được ghi nhận ở vật liệu từ loãng. Về mặt mô phỏng lý thuyết, mô hình q-state clock và mô hình XY được sử dụng để mô tả các vật liệu từ có spin phân cực mạnh và yếu tương ứng. Nghiên cứu này tập trung vào hiện tượng chuyển pha trong các mô hình này, đặc biệt là các tranh luận về chuyển pha KT và sự phù hợp của mô hình cho vật liệu từ pha loãng.

1.1. Vật Liệu Từ Tính Giới Thiệu Chung và Ứng Dụng

Vật liệu từ tính đóng vai trò quan trọng trong nhiều ứng dụng công nghệ hiện đại, từ lưu trữ dữ liệu đến cảm biến và thiết bị y tế. Nghiên cứu về tính chất từ của vật liệu giúp chúng ta hiểu rõ hơn về cấu trúc và tương tác bên trong, từ đó mở ra cơ hội phát triển các vật liệu mới với tính chất từ ưu việt hơn. Các tính chất nhiệt động của vật liệu từ cũng là một yếu tố quan trọng cần xem xét trong quá trình thiết kế và ứng dụng.

1.2. Chuyển Pha Trong Vật Liệu Từ Các Loại Chuyển Pha

Chuyển pha là hiện tượng vật lý quan trọng, thể hiện sự thay đổi đột ngột trong cấu trúc và tính chất từ của vật liệu khi nhiệt độ hoặc áp suất thay đổi. Có nhiều loại chuyển pha khác nhau, bao gồm chuyển pha bậc một, chuyển pha bậc hai và chuyển pha Kosterlitz-Thouless (KT). Mỗi loại chuyển pha có những đặc điểm riêng và ảnh hưởng khác nhau đến tính chất từ của vật liệu.

II. Thách Thức Nghiên Cứu Chuyển Pha Trong Hệ Spin Gián Đoạn

Mặc dù đã có nhiều tiến bộ, nghiên cứu về chuyển pha trong hệ spin gián đoạn vẫn còn nhiều thách thức. Các tranh luận về loại chuyển pha KT và đại lượng vật lý phù hợp để mô tả chính xác các chuyển pha vẫn tiếp diễn. Mô hình q-state clock, mặc dù được chấp nhận rộng rãi, vẫn còn những điểm chưa rõ ràng về chuyển pha của nó. Mô hình XYhq, với trường tinh thể bất đẳng hướng, cũng đặt ra câu hỏi về tính phổ biến của chuyển pha và ảnh hưởng của các giá trị nhiễu loạn. Việc xác định mô hình phù hợp để mô tả từ tính cho vật liệu từ pha loãng cũng là một thách thức lớn. Các nghiên cứu cần tập trung vào việc giải quyết những tranh luận này và cung cấp cái nhìn sâu sắc hơn về chuyển pha trong hệ spin gián đoạn.

2.1. Mô Hình Q State Clock Tranh Luận Về Chuyển Pha KT

Mô hình q-state clock, với đối xứng rời rạc Zq, là một tổng quát hóa của mô hình Ising. Tuy nhiên, các nghiên cứu vẫn tranh luận về việc liệu chuyển pha của mô hình này có phải là loại KT hay không, tương ứng với mỗi giá trị q. Việc xác định đại lượng vật lý nào phù hợp để mô tả và xác định chính xác các chuyển pha của mô hình này vẫn là một vấn đề mở.

2.2. Mô Hình XYhq Ảnh Hưởng Của Trường Tinh Thể Bất Đẳng Hướng

Mô hình XYhq, với trường tinh thể bất đẳng hướng q hướng, đặt ra câu hỏi về tính phổ biến của chuyển pha và ảnh hưởng của các giá trị nhiễu loạn. Việc xác định rõ ràng các giá trị nhiễu loạn và ảnh hưởng của chúng đến chuyển pha là một thách thức quan trọng. Nghiên cứu cần tập trung vào việc phân biệt rõ ràng các giá trị nhiễu loạn và ảnh hưởng của chúng đến tính chất từ của vật liệu.

III. Phương Pháp Mô Phỏng Monte Carlo Nghiên Cứu Chuyển Pha

Để giải quyết các thách thức trên, nghiên cứu này sử dụng phương pháp mô phỏng Monte Carlo. Phương pháp mô phỏng Monte Carlo là một công cụ mạnh mẽ để nghiên cứu các hệ vật lý phức tạp, đặc biệt là các hệ có nhiều tương tác và chuyển pha. Thuật toán Metropolis và thuật toán Wolff là hai trong số các thuật toán mô phỏng Monte Carlo phổ biến nhất. Nghiên cứu này sẽ sử dụng các thuật toán này để mô phỏng các hệ spin gián đoạn và nghiên cứu chuyển pha của chúng. Phân tích thống kê kết quả mô phỏng sẽ giúp xác định các tính chất nhiệt động và tính chất từ của hệ.

3.1. Thuật Toán Metropolis Cơ Chế Hoạt Động và Ưu Điểm

Thuật toán Metropolis là một thuật toán mô phỏng Monte Carlo phổ biến, được sử dụng để lấy mẫu từ một phân bố xác suất phức tạp. Thuật toán này dựa trên việc tạo ra các cấu hình mới bằng cách thay đổi ngẫu nhiên cấu hình hiện tại, và sau đó chấp nhận hoặc từ chối cấu hình mới dựa trên một tiêu chí xác suất. Ưu điểm của thuật toán Metropolis là đơn giản và dễ thực hiện.

3.2. Thuật Toán Wolff Hiệu Quả Trong Mô Phỏng Hệ Spin

Thuật toán Wolff là một thuật toán mô phỏng Monte Carlo hiệu quả, đặc biệt phù hợp cho việc mô phỏng các hệ spin. Thuật toán này dựa trên việc tạo ra các cụm spin liên kết với nhau, và sau đó lật toàn bộ cụm spin. Ưu điểm của thuật toán Wolff là giảm thiểu sự tương quan giữa các bước mô phỏng, giúp tăng tốc quá trình hội tụ.

IV. Ứng Dụng Mô Hình Q State Clock Phân Tích Chuyển Pha Chi Tiết

Nghiên cứu này tập trung vào việc ứng dụng mô hình q-state clock để phân tích chuyển pha chi tiết. Mô hình q-state clock sở hữu đối xứng rời rạc Zq và là tổng quát hóa của mô hình Ising. Mô hình này được dự đoán bằng lí thuyết và sau đó được xác nhận bằng các nghiên cứu mô phỏng số và bằng học máy cùng đồng thuận rằng mô hình q-state clock chỉ có một chuyển pha bậc hai cho q < 5 và hai chuyển pha riêng biệt ở nhiệt độ xác định T1 và T2 (T2 > T1) cho q ≥ 5. Pha ở giữa các nhiệt độ này là pha giả trật tự giống như mô hình XY dưới nhiệt độ chuyển pha KT, pha trên T2 là pha mất trật tự, và pha dưới T1 là pha trật tự.

4.1. Xác Định Nhiệt Độ Chuyển Pha Phương Pháp và Kết Quả

Việc xác định chính xác nhiệt độ tới hạn là một bước quan trọng trong nghiên cứu chuyển pha. Nghiên cứu này sử dụng các phương pháp phân tích thống kê để xác định nhiệt độ tới hạn từ kết quả mô phỏng Monte Carlo. Các kết quả mô phỏng cho thấy sự phụ thuộc của nhiệt độ tới hạn vào kích thước hệ và các tham số của mô hình.

4.2. Phân Tích Các Đại Lượng Vật Lý Mô Đun Xoắn và Tham Số Binder

Các đại lượng vật lý như mô đun xoắn và tham số Binder cung cấp thông tin quan trọng về chuyển pha. Nghiên cứu này phân tích sự phụ thuộc của mô đun xoắn và tham số Binder vào nhiệt độ và kích thước hệ. Các kết quả phân tích cho thấy sự thay đổi đột ngột của các đại lượng này tại nhiệt độ tới hạn, cho thấy sự xuất hiện của chuyển pha.

V. Nghiên Cứu Mô Hình XYhq Giản Đồ Pha và Tính Phổ Biến

Mô hình XY với trường tinh thể bất đẳng hướng q hướng (XYhq) các giá trị hq xuất hiện dưới dạng một trường ngoài phá vỡ đối xứng cơ bản U(1). Giới hạn hq → ∞ tái tạo mô hình q-state clock với các góc quay hữu hạn, nhưng đối với bất kì giá trị hữu hạn nào của hq, góc quay đều liên tục và thành phần thứ hai của biểu thức có thể được coi là nhiễu loạn. Một hệ quả của điều này là trong khi tính bất đẳng hướng của mô hình q-state clock có thể liên quan hoặc không liên quan đến các lớp phổ biến (như các lớp Ising, Potts ba trạng thái, Ising-like, và XY), người ta phải cho trường tinh thể bất đẳng hướng hq phân biệt rõ ràng ứng với các giá trị nhiễu loạn khác nhau.

5.1. Xây Dựng Giản Đồ Pha Ảnh Hưởng Của Trường Bất Đẳng Hướng

Việc xây dựng giản đồ pha là một bước quan trọng để hiểu rõ hơn về chuyển pha trong mô hình XYhq. Nghiên cứu này xây dựng giản đồ pha của mô hình XYhq bằng cách mô phỏng Monte Carlo và phân tích thống kê. Giản đồ pha cho thấy sự phụ thuộc của nhiệt độ tới hạn vào giá trị của trường bất đẳng hướng.

5.2. Tính Phổ Biến Của Chuyển Pha So Sánh Với Các Mô Hình Khác

Nghiên cứu này so sánh tính chất từ của mô hình XYhq với các mô hình khác, như mô hình Ising và mô hình XY thông thường. So sánh này giúp xác định tính phổ biến của chuyển pha trong mô hình XYhq. Các kết quả so sánh cho thấy sự tương đồng và khác biệt giữa các mô hình, giúp hiểu rõ hơn về vai trò của trường bất đẳng hướng.

VI. Kết Luận và Hướng Nghiên Cứu Tương Lai Về Hệ Spin

Nghiên cứu này đã cung cấp cái nhìn sâu sắc hơn về chuyển pha trong hệ spin gián đoạn, đặc biệt là trong mô hình q-state clock và mô hình XYhq. Phương pháp mô phỏng Monte Carlo đã được chứng minh là một công cụ hiệu quả để nghiên cứu các hệ vật lý phức tạp. Các kết quả nghiên cứu này có thể được sử dụng để thiết kế và phát triển các vật liệu từ mới với tính chất từ ưu việt hơn. Hướng nghiên cứu tương lai có thể tập trung vào việc nghiên cứu các hệ spin phức tạp hơn, như các hệ spin có tương tác tầm xa hoặc các hệ spin không cân bằng.

6.1. Tổng Kết Các Kết Quả Nghiên Cứu Chính

Nghiên cứu này đã đạt được một số kết quả quan trọng, bao gồm việc xác định nhiệt độ tới hạn của mô hình q-state clock và mô hình XYhq, xây dựng giản đồ pha của mô hình XYhq, và so sánh tính chất từ của các mô hình khác nhau. Các kết quả này cung cấp thông tin quan trọng về chuyển pha trong hệ spin gián đoạn.

6.2. Đề Xuất Các Hướng Nghiên Cứu Tiếp Theo

Nghiên cứu này đề xuất một số hướng nghiên cứu tiếp theo, bao gồm việc nghiên cứu các hệ spin phức tạp hơn, như các hệ spin có tương tác tầm xa hoặc các hệ spin không cân bằng. Nghiên cứu cũng có thể tập trung vào việc phát triển các thuật toán mô phỏng Monte Carlo hiệu quả hơn, giúp nghiên cứu các hệ vật lý lớn hơn và phức tạp hơn.

Nội dung chính

Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu mô phỏng sự chuyển pha trong một số hệ spin gián đoạn" của nghiên cứu sinh Lương Minh Tuấn, chuyên ngành Khoa học Vật liệu (mã số: 9440122), được thực hiện tại Trường Đại học Bách khoa Hà Nội dưới sự hướng dẫn của PGS. Đào Xuân Việt và TS. Nguyễn Đức Trung Kiên. Nghiên cứu thuộc hướng vật lý tính toán và vật lý vật liệu từ, tập trung vào hiện tượng chuyển pha từ tính trong các hệ thấp chiều và vật liệu bán dẫn từ pha loãng.

Tính cấp thiết và khoảng trống nghiên cứu

Trong lĩnh vực từ học hiện đại, các hệ vật liệu từ thấp chiều (màng mỏng hai chiều - 2D) và vật liệu bán dẫn từ pha loãng (Diluted Magnetic Semiconductors - DMS) thu hút sự quan tâm lớn do bộc lộ các tính chất vật lý khác biệt so với vật liệu từ khối. Ở màng mỏng 2D, các câu hỏi cốt lõi đặt ra là sự tồn tại của trật tự từ, các dạng trật tự giả tầm xa và cơ chế chuyển pha Kosterlitz-Thouless (KT), chuyển pha bậc một hay bậc hai.

Để mô tả các hệ này, giới nghiên cứu sử dụng hai cách tiếp cận: mô hình spin liên tục (mô hình 2D XY) khi spin định hướng tự do trong mặt phẳng và mô hình spin gián đoạn (mô hình q-state clock) khi spin bị phân cực mạnh theo một số hướng đối xứng rời rạc. Tuy nhiên, nhiều vấn đề lý thuyết và mô phỏng số vẫn tồn tại bất đồng:

1. Đối với mô hình 2D q-state clock sở hữu đối xứng rời rạc $Z_q$, các tranh luận tập trung vào bản chất chuyển pha ở $q = 5$: liệu hai chuyển pha xảy ra ở các nhiệt độ $T_1$ và $T_2$ ($T_2 > T_1$) có cùng thuộc loại chuyển pha KT hay không, và cần sử dụng đại lượng vật lý nào để xác định chính xác các điểm chuyển pha này khi pha trung gian rất hẹp.
2. Đối với mô hình 2D XY chịu tác động của trường tinh thể bất đẳng hướng $q$-hướng ($XYh_q$) - mô hình trung gian chuyển tiếp giữa mô hình liên tục ($h_q = 0$) và mô hình gián đoạn ($h_q \rightarrow \infty$): quy luật phụ thuộc của nhiệt độ chuyển pha và bản chất chuyển pha vào cường độ trường tinh thể $h_q$ chưa được làm sáng tỏ, đặc biệt là sự thiếu vắng các giản đồ pha $T - h_q$ hoàn chỉnh cho các trường hợp $q = 3, 5, 6, 8$.
3. Đối với hệ bán dẫn từ pha loãng (DMS): mục tiêu thu nhận trật tự sắt từ ở nhiệt độ phòng với nồng độ tạp chất thấp có ý nghĩa quan trọng cho kỹ thuật spintronics. Các mô phỏng Monte Carlo trước đây chủ yếu dùng mô hình spin liên tục và thường phải đưa vào một từ trường ngoài cưỡng bức để tái lập nhiệt độ Curie ($T_C$) gần nhiệt độ phòng; khi không có từ trường ngoài, nhiệt độ $T_C$ tính toán được lại thấp hơn nhiều. Câu hỏi đặt ra là trạng thái spin thực tế của tạp chất nhẹ (như Carbon) trong mạng bán dẫn là liên tục hay rời rạc.

Mục tiêu nghiên cứu của luận án

Văn bản luận án xác định ba mục tiêu cụ thể:

1. Mô hình 2D q-state clock: Làm rõ hai chuyển pha trong mô hình 5-state clock là chuyển pha loại KT hay không phải loại KT; đề xuất đại lượng vật lý thích hợp để xác định chính xác các nhiệt độ chuyển pha.
2. Mô hình 2D $XYh_q$: Nghiên cứu ảnh hưởng của trường tinh thể đến chuyển pha trong mô hình 2D $XYh_q$ với $q = 3, 5, 6, 8$ nhằm làm rõ chuyển pha là loại KT hay không phải loại KT; xác định quy luật phụ thuộc của nhiệt độ chuyển pha vào cường độ trường tinh thể và xây dựng các giản đồ pha tương ứng.
3. Áp dụng mô hình spin vào chất bán dẫn từ pha loãng: Thực hiện mô phỏng Monte Carlo về trật tự sắt từ của hệ vật liệu CdS pha tạp Carbon ($CdS:C$ nồng độ 5.55%) bằng cả hai mô hình spin liên tục và spin rời rạc; từ đó đánh giá mô hình nào phù hợp để mô tả trạng thái spin trong hệ.

Đối tượng và phạm vi nghiên cứu

• Đối tượng nghiên cứu: Các mô hình spin lý thuyết bao gồm mô hình 2D q-state clock, mô hình 2D $XYh_q$, mô hình Heisenberg và mô hình spin rời rạc; các đại lượng vật lý thống kê mô tả hiện tượng chuyển pha từ tính; hệ vật liệu bán dẫn từ pha loãng CdS pha tạp Carbon.
• Phạm vi không gian và cấu trúc mạng:
  • Các hệ 2D được mô phỏng trên mạng tinh thể hình vuông hai chiều kích thước $L \times L$ với điều kiện biên tuần hoàn.
  • Hệ bán dẫn từ pha loãng $CdS:C$ được khảo sát trên cấu trúc mạng tinh thể wurtzite với nồng độ pha tạp thay thế vị trí nguyên tử lưu huỳnh (S) bằng carbon (C) ở mức 5.55%, xét tương tác trao đổi trong phạm vi 8 lớp lân cận gần nhất.
• Phạm vi phương pháp: Phương pháp mô phỏng số Monte Carlo cân bằng nhiệt kết hợp lý thuyết tỉ lệ kích thước hữu hạn (Finite-size Scaling - FSS).

---

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước được điểm lại

Trong Chương 1, tác giả hệ thống hóa các mốc nghiên cứu lý thuyết, thực nghiệm và mô phỏng:

• Nền tảng mô hình spin và chuyển pha: Khởi đầu từ công trình của Ernst Ising (1925) giải chính xác mô hình một chiều chỉ ra không có chuyển pha, Lars Onsager (1944) giải chính xác mô hình Ising 2D tìm ra chuyển pha bậc hai tại nhiệt độ Curie $T_C$. J. M. Kosterlitz và D. J. Thouless (1973, 1974) phát hiện cơ chế chuyển pha topo KT trong mô hình 2D XY thông qua sự phân ly của các cặp xoáy - phản xoáy (vortex - antivortex) tại nhiệt độ $T_{KT} \approx 0.89$. Tiêu chuẩn Nelson - Kosterlitz xác định bước nhảy phổ quát của mô đun xoắn tại điểm chuyển pha KT theo hệ thức $\Upsilon = 2T/\pi$.
• Mô hình 2D q-state clock: J. V. José và các cộng sự (1977) cùng Elitzur và cộng sự (1979) dự đoán sự tồn tại của hai chuyển pha cho $q \ge 5$. Nhiều nghiên cứu mô phỏng Monte Carlo sau đó (như Borisenko 2011, Kumano 2013, Tomita 2002) và nghiên cứu bằng học máy xác nhận sự tồn tại của pha giả trật tự tầm xa (Quasi-long-range-order - QLRO) nằm giữa hai nhiệt độ chuyển pha $T_1$ (chuyển pha từ trật tự sang QLRO) và $T_2$ (từ QLRO sang mất trật tự). Tuy nhiên, Lapilli và cộng sự (2006) lại cho rằng chuyển pha KT chỉ xuất hiện khi $q \ge 8$. Baek, Minnhagen và Kim (2009, 2011) cho rằng với $q = 5$, chuyển pha tại nhiệt độ cao $T_2$ không phải là chuyển pha loại KT dựa trên phân tích mô đun xoắn góc xoắn liên tục. Ngược lại, các công trình của Chatterjee (2018) và các nhóm khác sử dụng mô đun xoắn góc xoắn rời rạc $2\pi/q$ hoặc chiều dài tương quan từ lại ủng hộ bản chất KT nhưng thừa nhận tín hiệu bị mờ nhạt do dải nhiệt độ pha QLRO quá hẹp.
• Mô hình 2D $XYh_q$: José và cộng sự (1977) đề xuất mô hình; các nghiên cứu giải tích và mô phỏng của Taroni (2008), Lou (2007) chỉ ra tính chất chuyển pha phụ thuộc vào số hướng $q$ và trường tinh thể $h_q$. Mặc dù vậy, vùng trường yếu ($h_q \ll 1$) và quy luật tiệm cận của nhiệt độ chuyển pha theo độ lớn của trường chưa được khảo sát định lượng đầy đủ.
• Hệ bán dẫn từ pha loãng DMS: Sau phát hiện về tính sắt từ ở nhiệt độ phòng trong DMS (Dietl 2000, Matsumoto 2001), các tính toán nguyên lý đầu tiên kết hợp mô phỏng Monte Carlo được tiến hành cho các hệ pha tạp kim loại chuyển tiếp (Mn, Fe, Co) và tạp chất phi từ tính nhẹ (C, N) bởi Sato, Katayama-Yoshida (2002, 2010), Pan (2007), Zhou (2008). Trong hệ CdS pha tạp C, các tác giả đi trước chủ yếu dùng mô hình spin liên tục Heisenberg kết hợp từ trường ngoài nhỏ để thu được trật tự sắt từ gần 300 K; khi thiếu từ trường ngoài, nhiệt độ trật tự suy giảm nghiêm trọng.

Khoảng trống luận án lựa chọn giải quyết

Từ các bất đồng trên, luận án tập trung giải quyết:

1. Giải quyết dứt điểm tranh cãi về bản chất hai chuyển pha ở $q = 5$ trong mô hình 2D q-state clock bằng cách bổ sung đại lượng phân tích mới có độ nhạy cao.
2. Thiết lập bức tranh hoàn chỉnh về giản đồ pha $T - h_q$ và hàm giải tích mô tả sự biến thiên của các nhiệt độ chuyển pha theo trường $h_q$ cho mô hình 2D $XYh_q$ với $q = 3, 5, 6, 8$.
3. Kiểm tra tính xác đáng của giả định spin liên tục so với spin rời rạc trong việc giải thích trật tự sắt từ nhiệt độ phòng ở hệ bán dẫn từ pha loãng $CdS:C$ (5.55%) khi không sử dụng từ trường ngoài cưỡng bức.

---

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và các mô hình toán học

Luận án vận dụng lý thuyết vật lý thống kê cổ điển, lý thuyết chuyển pha và hiện tượng tới hạn, cùng các biểu thức Hamiltonian mô tả tương tác từ:

1. Mô hình 2D Ising: $$H = - \sum_{\langle i, j \rangle} J_{ij} S_i S_j$$ Trong đó spin chỉ nhận hai giá trị $S_i = \pm 1$, tương tác lân cận gần nhất $J_{ij} = J > 0$.

2. Mô hình 2D XY: $$H = - J \sum_{\langle i, j \rangle} \cos(\theta_i - \theta_j)$$ Trong đó vector spin phẳng hai chiều $\vec{S}_i = (\cos\theta_i, \sin\theta_i)$ với góc quay liên tục $\theta_i \in [0, 2\pi)$.

3. Mô hình 2D q-state clock: $$H = - J \sum_{\langle i, j \rangle} \cos(\theta_i - \theta_j)$$ Với các góc định hướng gián đoạn: $\theta_i = \frac{2\pi n}{q}$, trong đó $n \in {0, 1, \dots, q - 1}$. Khi $q = 2$, mô hình tương đương Ising; khi $q \rightarrow \infty$, mô hình tiệm cận mô hình XY.

4. Mô hình 2D $XYh_q$: $$H = - J \sum_{\langle i, j \rangle} \cos(\theta_i - \theta_j) - h_q \sum_i \cos(q\theta_i)$$ Số hạng thứ hai đóng vai trò trường tinh thể phá vỡ đối xứng liên tục $U(1)$, ưu tiên $q$ hướng đối xứng trong mặt phẳng.

5. Mô hình Heisenberg và mô hình spin rời rạc trong DMS: Tương tác trao đổi giữa các vị trí pha tạp C trong mạng CdS: $$H = - \sum_{\langle i, j \rangle} J_{ij} \vec{S}_i \cdot \vec{S}j$$ Hệ số tương tác trao đổi $J{ij}$ được lấy từ các tính toán nguyên lý đầu tiên cho 8 lớp lân cận gần nhất trong mạng wurtzite. Trong mô hình spin liên tục, $\vec{S}_i$ là vector đơn vị ba chiều tự do; trong mô hình spin rời rạc, $\vec{S}_i$ bị ràng buộc định hướng theo các trục tinh thể cụ thể.

Phương pháp mô phỏng Monte Carlo và các thuật toán

Nghiên cứu sử dụng kỹ thuật lấy mẫu quan trọng (Importance Sampling) dựa trên phân bố xác suất Boltzmann $p_i \propto e^{-E_i / k_B T}$, bảo đảm hai điều kiện tiên quyết là tính Ergodic và cân bằng chi tiết (Detailed Balance):

Thuật toán	Phạm vi ứng dụng trong luận án	Đặc điểm kỹ thuật và vai trò
Metropolis	Áp dụng cho mọi hệ (clock, $XYh_q$, Heisenberg, DMS)	Cập nhật đơn spin với xác suất chấp nhận $W_{i \rightarrow f} = \min(1, e^{-\Delta E / k_B T})$. Đóng vai trò thuật toán cơ sở đảm bảo tính ergodic.
Đám Wolff (Cluster)	Mô hình 2D clock và 2D $XYh_q$	Lật một đám spin cùng lúc qua mặt phẳng phản xạ ngẫu nhiên với xác suất kết nạp liên kết $p = 1 - e^{-2\beta (\vec{n}\cdot\vec{S}_i)(\vec{n}\cdot\vec{S}_j)}$; triệt tiêu hiện tượng chậm tới hạn (critical slowing down).
Over-relaxation (OR)	Mô hình Heisenberg cho DMS	Quay vector spin đối xứng qua hướng từ trường địa phương $\vec{h}i = \sum_j J{ij} \vec{S}_j$; giữ nguyên năng lượng ($\Delta E = 0$), giúp hệ vượt rào thế năng nhanh ở nhiệt độ thấp.
General Over-relaxation	Mô hình spin rời rạc và $XYh_q$	Quay spin trên mặt nón quanh từ trường địa phương với góc ngẫu nhiên thỏa mãn $\Delta E = 0$. Phải kết hợp cùng Metropolis để đảm bảo tính ergodic.
Parallel-tempering (Replica Exchange)	Mô hình spin rời rạc cho CdS:C	Chạy đồng thời $N$ bản sao ở tập nhiệt độ $[T_{\min}, T_{\max}]$; trao đổi cấu hình giữa hai nhiệt độ lân cận với xác suất $W = \min(1, \exp(\Delta\beta \Delta E))$. Sử dụng thuật toán lặp điều chỉnh tập nhiệt độ để tỉ lệ chấp nhận trao đổi đạt giá trị tối ưu đồng đều.

Quy trình tính toán cho hệ 2D kết hợp 5 bước Metropolis với 1 bước Wolff thành một bước Monte Carlo (MCs). Chương trình được viết bằng ngôn ngữ Fortran, biên dịch bằng Intel Fortran trên hệ điều hành Linux. Dãy số ngẫu nhiên được sinh bằng chương trình con tự xây dựng theo hàm phân bố đều.

Các đại lượng vật lý thống kê và kỹ thuật phân tích

Để nhận diện loại chuyển pha và định vị điểm tới hạn, luận án tính toán:

• Nhiệt dung riêng ($C$): Tính qua hai phương pháp độc lập để kiểm tra cân bằng nhiệt: $$C = \frac{\Delta E}{\Delta T} \quad \text{và} \quad C = \frac{\langle E^2 \rangle - \langle E \rangle^2}{k_B T^2}$$
• Độ từ hóa ($m$) và độ cảm từ ($\chi$): $$m = \frac{1}{N} \left| \sum_i \vec{S}_i \right|, \quad \chi = \frac{N}{k_B T} \left( \langle m^2 \rangle - \langle m \rangle^2 \right)$$
• Tham số Binder ($g$ hoặc $U_2$): $$g = \frac{3}{2} \left( 1 - \frac{\langle m^4 \rangle}{3 \langle m^2 \rangle^2} \right)$$ Các đường $g(T)$ với các kích thước $L$ khác nhau cắt nhau tại một điểm duy nhất biểu thị chuyển pha bậc hai; hành vi chồng chập trên một khoảng nhiệt độ biểu thị pha KT.
• Đạo hàm của tham số Binder theo nhiệt độ ($dg/dT$): Được tính thông qua thăng giáng năng lượng và độ từ hóa: $$\frac{dg}{dT} = \frac{1}{T^2} \left[ \frac{\langle m^4 H \rangle - \langle m^4 \rangle \langle H \rangle}{\langle m^2 \rangle^2} - \frac{2 \langle m^4 \rangle (\langle m^2 H \rangle - \langle m^2 \rangle \langle H \rangle)}{\langle m^2 \rangle^3} \right]$$ Đại lượng này xuất hiện các cực trị sắc nét tại các điểm chuyển pha, khắc phục nhược điểm của sự chồng chập tham số Binder.
• Mô đun xoắn ($\Upsilon$) và đạo hàm ($d\Upsilon/dT$): Đo khả năng kháng lại biến dạng xoắn spin nhỏ dọc theo một trục; rơi gián đoạn về 0 tại nhiệt độ chuyển pha KT khi $L \rightarrow \infty$.
• Tỉ số chiều dài tương quan từ ($\xi/L$): Xác định thông qua biến đổi Fourier của hàm tương quan spin.
• Xử lý số liệu và sai số: Sử dụng kỹ thuật Jackknife để khử độ lệch và ước lượng sai số thống kê; áp dụng lý thuyết tỉ lệ kích thước hữu hạn (FSS) để ngoại suy nhiệt độ chuyển pha về giới hạn nhiệt động học ($L \rightarrow \infty$). Quá trình lấy mẫu chỉ thực hiện trên nửa sau của chuỗi MCs sau khi hệ đã đạt cân bằng nhiệt hoàn toàn.

---

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 khái quát lịch sử phát triển của từ học từ khoáng vật Magnetite ($Fe_3O_4$), các định luật điện từ cổ điển đến thuyết lượng tử về mô-men từ nguyên tử (mô-men từ quỹ đạo, mô-men từ spin, mô-men từ hạt nhân). Tác giả phân loại 5 nhóm vật liệu từ cơ bản: nghịch từ, thuận từ, sắt từ, phản sắt từ và ferit từ.

Tiếp đó, chương luận giải các khái niệm nhiệt động học về pha, tham số trật tự, phân loại chuyển pha theo Ehrenfest (chuyển pha bậc một với bước nhảy gián đoạn của đạo hàm bậc nhất năng lượng tự do; chuyển pha bậc hai với điểm kỳ dị ở đạo hàm bậc hai) và chuyển pha topo Kosterlitz-Thouless với điểm kỳ dị cốt yếu. Chương 1 giới thiệu chi tiết 4 mô hình spin lý thuyết (2D Ising, 2D XY, 2D q-state clock, 2D $XYh_q$) và vật liệu bán dẫn từ pha loãng DMS, đồng thời tổng kết toàn diện các công trình nghiên cứu đi trước, chỉ rõ các mâu thuẫn học thuật để làm căn cứ hình thành nhiệm vụ nghiên cứu của luận án.

Chương 2: Phương pháp mô phỏng và các đại lượng vật lí

Chương 2 mô tả chi tiết phương pháp Monte Carlo và cơ sở thống kê cổ điển. Tác giả phân tích sự hạn chế của phương pháp lấy mẫu đơn giản và tính tất yếu của phương pháp lấy mẫu quan trọng thông qua xích Markov. Các thuật toán Metropolis, thuật toán đám Wolff, thuật toán Over-relaxation, thuật toán general Over-relaxation và thuật toán Parallel-tempering được thiết lập với đầy đủ biểu thức toán học xác suất chuyển đổi.

Chương này trình bày chi tiết thuật toán tối ưu hóa tỷ lệ chấp nhận trong Parallel-tempering qua ba bước lặp, kiểm tra thời gian tự tương quan $\tau$ (ước tính $\tau \approx 60$ tại $T_{\max} = 522.225\text{ K}$ cho hệ $CdS:C$). Tác giả thiết lập hệ thống biểu thức giải tích vi phân cho các đại lượng vật lý thống kê ($C, m, \chi, g, dg/dT, \Upsilon, d\Upsilon/dT, \xi/L$), xây dựng tiêu chuẩn kiểm tra cân bằng nhiệt qua sự trùng khớp giữa hai cách tính nhiệt dung riêng, và trình bày cơ sở toán học của kỹ thuật Jackknife cùng lý thuyết FSS cho chuyển pha bậc hai và chuyển pha KT.

Chương 3: Nghiên cứu chuyển pha trong mô hình q-state clock bằng phương pháp mô phỏng Monte Carlo

Chương 3 tập trung giải quyết các mô hình 2D q-state clock trên mạng vuông với kích thước $L$ thay đổi, trọng tâm là các trường hợp $q = 4, 5, 6$:

• Với $q = 4$: Tỉ số chiều dài tương quan $\xi/L$ và tham số Binder $g(T)$ với các kích thước $L$ khác nhau cắt nhau tại một điểm duy nhất, khẳng định mô hình chỉ có một chuyển pha bậc hai duy nhất tại $T_C \approx 1.13$.
• Với $q = 5$: Tác giả chỉ ra rằng mô đun xoắn liên tục không thể hiện rõ chuyển pha KT, nhưng đại lượng đạo hàm tham số Binder $dg/dT$ xuất hiện cực tiểu âm rõ rệt tại nhiệt độ chuyển pha thấp $T_1$, trong khi tỉ số $\xi/L$ và tham số Binder thể hiện vùng chồng chập đặc trưng giữa $T_1$ và $T_2$. Áp dụng FSS đối với vị trí các đỉnh, luận án xác định hai chuyển pha riêng biệt đều thuộc loại KT với các nhiệt độ chuyển pha $T_1 \approx 0.90$ và $T_2 \approx 0.95$. Kết quả này bác bỏ nhận định của Baek và cộng sự về việc $T_2$ không phải là chuyển pha KT.
• Với $q = 6$: Dải pha trung gian QLRO mở rộng rõ rệt. Cả đạo hàm tham số Binder $dg/dT$ và đạo hàm mô đun xoắn $d\Upsilon/dT$ đều phân kỳ tại hai nhiệt độ chuyển pha KT riêng biệt $T_1 \approx 0.68$ và $T_2 \approx 0.92$.

Chương 4: Nghiên cứu phân tích các chuyển pha và xây dựng giản đồ pha cho mô hình $XYh_q$

Chương 4 khảo sát mô hình 2D $XYh_q$ với các giá trị trường tinh thể $h_q$ biến thiên từ rất nhỏ đến rất lớn:

• Trường hợp $q = 3$: Nhiệt dung riêng và tỉ số chiều dài tương quan xác nhận mô hình chỉ có duy nhất một chuyển pha bậc hai. Giản đồ pha chia làm hai vùng:
  • Vùng trường mạnh ($h_3 > 0.01$): khi $h_3$ giảm, nhiệt độ chuyển pha $T_C$ giảm dần và tiệm cận về giá trị $T_C \rightarrow 1.01$.
  • Vùng trường yếu ($h_3 \le 0.01$): khi $h_3$ tiếp tục giảm, nhiệt độ chuyển pha không phụ thuộc vào $h_3$ và giữ giá trị không đổi xấp xỉ nhiệt độ chuyển pha KT của mô hình XY thuần túy ($T_C \approx T_{KT} \approx 0.89$).
• Trường hợp $q = 5$: Chứng minh chuyển pha tại nhiệt độ cao $T_2$ là chuyển pha KT. Giản đồ pha $T - h_5$ được thiết lập đầy đủ và phân chia thành ba vùng rõ rệt:
  • Khi $h_5 > 2$: Cả hai nhiệt độ chuyển pha $T_1$ và $T_2$ bão hòa, không phụ thuộc vào $h_5$.
  • Khi $1 \le h_5 \le 2$: $T_1$ suy giảm theo quy luật hàm logarit của $h_5$, trong khi $T_2$ giảm theo quy luật hàm lũy thừa.
  • Khi $h_5 < 1$: $T_1$ suy giảm dần về 0 theo quy luật logarit ($\sim \ln h_5$), trong khi $T_2$ không đổi và giữ nguyên giá trị chuyển pha KT thông thường ($T_2 \approx 0.89$).
• Trường hợp $q = 6$ và $q = 8$: Xây dựng giản đồ pha $T - h_6$ hoàn chỉnh từ trường lớn đến trường nhỏ. Với $h_6 \ge 1$, hai nhiệt độ chuyển pha không đổi; với $h_6 < 1$, $T_2$ giữ nguyên là chuyển pha KT thông thường, còn $T_1$ lần đầu tiên được xác định chính xác quy luật suy giảm tỷ lệ với logarit của cường độ trường tinh thể ($T_1 \propto \ln h_6$).

Chương 5: Nghiên cứu mô phỏng Monte Carlo về trật tự sắt từ ở nhiệt độ phòng trong chất bán dẫn từ pha loãng

Chương 5 áp dụng phương pháp Monte Carlo quy mô lớn cho hệ bán dẫn từ pha loãng CdS pha tạp Carbon ($CdS:C$) với nồng độ 5.55% trên mạng wurtzite, sử dụng bộ thông số tương tác trao đổi $J_{ij}$ trên 8 lớp lân cận:

• Mô phỏng với mô hình spin liên tục (Heisenberg): Không áp dụng từ trường ngoài cưỡng bức, kết quả cho thấy độ từ hóa tự phát $m$ giảm nhanh về 0 khi nhiệt độ còn rất thấp; tham số Binder không hội tụ về điểm chuyển pha ở nhiệt độ cao; nhiệt độ Curie ước tính thấp hơn nhiều so với nhiệt độ phòng. Điều này chứng minh mô hình Heisenberg liên tục cổ điển không phản ánh đúng thực nghiệm sắt từ nhiệt độ phòng khi thiếu từ trường ngoài.
• Mô phỏng với mô hình spin rời rạc: Nhờ kết hợp thuật toán Metropolis và Parallel-tempering, hệ vượt qua được các cực tiểu năng lượng địa phương. Độ từ hóa, độ cảm từ và tham số Binder thể hiện rõ nét chuyển pha sắt từ - thuận từ với nhiệt độ Curie $T_C$ nằm sát nhiệt độ phòng mà không cần đưa vào bất kỳ từ trường ngoài nhân tạo nào. Kết quả khẳng định mô hình spin rời rạc mô tả chính xác bản chất tương tác từ trong hệ $CdS:C$.

---

Kết quả và những đóng góp mới

Luận án đã đạt được các kết quả khoa học cụ thể, được tác giả tổng hợp trong phần đóng góp mới:

Đóng góp về lý luận và mô hình hóa

1. Đối với mô hình 2D q-state clock:
   • Chứng minh mô hình 5-state clock tồn tại hai chuyển pha KT riêng biệt ($T_1 \approx 0.90$ và $T_2 \approx 0.95$), giải quyết dứt điểm các quan điểm trái chiều trước đây.
   • Đề xuất và tính toán đại lượng đạo hàm tham số Binder theo nhiệt độ ($dg/dT$). Đây là đại lượng cho tín hiệu cực trị rõ ràng, cho phép xác định chính xác các nhiệt độ chuyển pha $T_1, T_2$ cho các trường hợp $q = 4, 5, 6$.
2. Đối với mô hình 2D $XYh_q$:
   • Với $q = 3$: Khẳng định mô hình chỉ có một chuyển pha bậc hai; phát hiện hiện tượng phân vùng ảnh hưởng của trường tinh thể tại ngưỡng $h_3 = 0.01$ (với $h_3 > 0.01$, $T_C$ giảm dần về 1.01; với $h_3 \le 0.01$, $T_C$ ổn định ở giá trị $T_{KT} \approx 0.89$).
   • Với $q = 5$: Chứng minh chuyển pha $T_2$ là chuyển pha loại KT; xây dựng hoàn chỉnh giản đồ pha $T - h_5$ chia thành ba vùng ứng xử (vùng bão hòa $h_5 > 2$; vùng trung gian $1 \le h_5 \le 2$ với $T_1 \sim \ln h_5$ và $T_2$ biến thiên theo hàm lũy thừa; vùng trường yếu $h_5 < 1$ với $T_1 \rightarrow 0$ theo quy luật logarit còn $T_2$ không đổi).
   • Với $q = 6$: Xây dựng thành công giản đồ pha $T - h_6$ đầy đủ từ trường lớn đến trường rất nhỏ; xác định quy luật nhiệt độ chuyển pha thấp suy giảm theo hàm logarit của cường độ trường ($T_1 \sim \ln h_6$) khi $h_6 < 1$.
3. Đối với chất bán dẫn từ pha loãng (DMS):
   • Đưa ra bằng chứng mô phỏng số chỉ rõ mô hình spin liên tục (Heisenberg) không phù hợp để mô tả từ tính của spin Carbon trong mạng tinh thể CdS khi không có từ trường ngoài.
   • Đề xuất mô hình spin rời rạc loại bỏ sự thiếu hiệu quả của mô hình spin liên tục, giải thích thỏa đáng sự xuất hiện của trật tự sắt từ ở nhiệt độ phòng trong vật liệu $CdS:C$ pha tạp thấp (5.55%).

Bảng tổng hợp các kết quả chuyển pha chính

Hệ mô hình	Số trạng thái / Cường độ trường	Bản chất chuyển pha	Nhiệt độ chuyển pha ($T_C$ hoặc $T_1, T_2$)	Quy luật phụ thuộc vào trường tinh thể
2D q-state clock	$q = 4$	Một chuyển pha bậc hai	$T_C \approx 1.13$	Không có trường ngoài
2D q-state clock	$q = 5$	Hai chuyển pha KT	$T_1 \approx 0.90; T_2 \approx 0.95$	Không có trường ngoài
2D q-state clock	$q = 6$	Hai chuyển pha KT	$T_1 \approx 0.68; T_2 \approx 0.92$	Không có trường ngoài
2D $XYh_3$	$h_3 > 0.01$	Một chuyển pha bậc hai	$T_C$ giảm dần về $1.01$	Phụ thuộc vào $h_3$
2D $XYh_3$	$h_3 \le 0.01$	Một chuyển pha bậc hai	$T_C \approx 0.89$ ($T_{KT}$)	Không đổi theo $h_3$
2D $XYh_5$	$h_5 > 2$	Hai chuyển pha KT	$T_1, T_2$ bão hòa	Không đổi theo $h_5$
2D $XYh_5$	$1 \le h_5 \le 2$	Hai chuyển pha KT	$T_1$ giảm, $T_2$ giảm	$T_1 \sim \ln h_5$, $T_2 \sim h_5^\alpha$
2D $XYh_5$	$h_5 < 1$	Hai chuyển pha KT	$T_1 \rightarrow 0; T_2 \approx 0.89$	$T_1 \sim \ln h_5$, $T_2$ không đổi
2D $XYh_6$	$h_6 < 1$	Hai chuyển pha KT	$T_1$ giảm; $T_2 \approx 0.89$	$T_1 \sim \ln h_6$, $T_2$ không đổi
DMS $CdS:C$ (5.55%)	Mô hình Heisenberg	Không trật tự nhiệt độ phòng	$T_C \ll 300\text{ K}$	Không áp dụng
DMS $CdS:C$ (5.55%)	Mô hình spin rời rạc	Chuyển pha sắt từ	$T_C \approx 300\text{ K}$ (nhiệt độ phòng)	Không áp dụng

Kiến nghị và giải pháp đề xuất

• Về mặt tính toán mô phỏng: Khuyến nghị sử dụng đại lượng đạo hàm tham số Binder ($dg/dT$) kết hợp lý thuyết tỉ lệ kích thước hữu hạn khi phân tích các chuyển pha loại Kosterlitz-Thouless có vùng giả trật tự hẹp; áp dụng thuật toán Parallel-tempering có tối ưu hóa bước nhiệt để giải quyết sự đóng băng cấu hình ở nhiệt độ thấp trong các hệ spin pha loãng.
• Về mặt mô hình hóa vật liệu từ: Đối với các hệ bán dẫn từ pha loãng có pha tạp phi kim loại nhẹ như C, cần chuyển dịch từ giả định spin liên tục sang mô hình spin gián đoạn/rời rạc phản ánh bất đẳng hướng từ tinh thể thực tế.

---

Hạn chế và hướng nghiên cứu tiếp

Văn bản luận án giới hạn phạm vi khảo sát trong các điều kiện tính toán cụ thể:

• Các mô phỏng 2D clock và 2D $XYh_q$ mới thực hiện trên mạng tinh thể hình vuông hai chiều đơn giản với điều kiện biên tuần hoàn; chưa mở rộng sang mạng tam giác, lục giác hay các mạng có tương tác phản sắt từ gây ức chế hình học (frustration).
• Đối với chất bán dẫn từ pha loãng $CdS:C$, tính toán cố định ở một mức nồng độ pha tạp duy nhất là 5.55% và xét tương tác trao đổi trong phạm vi 8 lớp lân cận của cấu trúc wurtzite; chưa khảo sát sự phụ thuộc liên tục của nhiệt độ Curie theo dải nồng độ pha tạp biến thiên dưới 10% hoặc ảnh hưởng của các khuyết tật mạng kèm theo.
• Hướng nghiên cứu tiếp tự nhiên được rút ra từ phạm vi luận án là mở rộng khảo sát các mô hình $XYh_q$ cho các giá trị $q$ chẵn lớn hơn ($q = 8$), khảo sát các hệ spin có tương tác tầm xa và áp dụng mô hình spin rời rạc cho các hệ DMS nền oxit hoặc hợp chất $A^{III}B^V$ khác.

---

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cho:

• Nghiên cứu sinh và học viên cao học ngành Vật lý lý thuyết, Vật lý tính toán và Khoa học Vật liệu: Cung cấp tài liệu hoàn chỉnh về quy trình xây dựng mã mô phỏng Monte Carlo (từ thuật toán Metropolis, thuật toán đám Wolff đến Parallel-tempering) viết bằng Fortran chạy trên môi trường Linux.
• Các nhà nghiên cứu về hiện tượng chuyển pha và từ học: Cung cấp phương pháp luận sử dụng đại lượng đạo hàm tham số Binder ($dg/dT$) để nhận diện và định vị chính xác chuyển pha KT - một kỹ thuật khắc phục được sự mờ nhạt của tham số Binder truyền thống và sự thiếu nhạy bén của mô đun xoắn liên tục.
• Các nhóm nghiên cứu về vật liệu Spintronics và bán dẫn từ: Phần phân tích so sánh giữa mô hình spin Heisenberg và spin rời rạc trong Chương 5 là cơ sở tham khảo quan trọng để lựa chọn mô hình tính toán từ tính cho các vật liệu pha tạp nguyên tố nhẹ mà không cần đưa vào các tham số trường ngoài thiếu căn cứ thực nghiệm.

---

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

1. Tại sao chuyển pha ở mô hình 5-state clock từng gây tranh cãi và luận án đã giải quyết như thế nào?

Tranh cãi xuất phát từ việc pha giả trật tự (QLRO) ở $q = 5$ rất hẹp, khiến các đại lượng vật lý truyền thống phản ánh mâu thuẫn: mô đun xoắn với góc xoắn liên tục cho dấu hiệu không phải chuyển pha KT (nghiên cứu của Baek), trong khi mô đun xoắn rời rạc cho dấu hiệu KT nhưng tín hiệu mờ nhạt. Luận án giải quyết vấn đề này bằng cách lần đầu tiên tính toán đại lượng đạo hàm của tham số Binder ($dg/dT$), làm xuất hiện cực tiểu âm rõ nét tại $T_1$, kết hợp với sự chồng chập của tỉ số chiều dài tương quan $\xi/L$ giữa $T_1$ và $T_2$, qua đó chứng minh mô hình có hai chuyển pha đều thuộc loại KT.

2. Đại lượng đạo hàm tham số Binder ($dg/dT$) có ưu điểm gì so với tham số Binder thông thường ($g$)?

Ở các chuyển pha loại KT, tham số Binder thông thường $g(T)$ biểu hiện hành vi chồng chập giữa các kích thước mạng khác nhau trên toàn bộ dải pha giả trật tự, gây khó khăn cho việc định vị chính xác nhiệt độ chuyển pha. Đạo hàm $dg/dT$ theo nhiệt độ khuếch đại sự biến thiên của các mô-men từ hóa bậc hai và bậc bốn theo năng lượng, tạo thành các đỉnh/cực trị phụ thuộc kích thước sắc nét gần nhiệt độ chuyển pha, cho phép áp dụng lý thuyết tỉ lệ kích thước hữu hạn (FSS) để xác định chính xác nhiệt độ chuyển pha.

3. Quy luật biến thiên của nhiệt độ chuyển pha theo trường tinh thể trong mô hình $XYh_5$ diễn ra như thế nào?

Luận án chỉ ra giản đồ pha $T - h_5$ chia thành ba vùng: khi $h_5 > 2$, nhiệt độ chuyển pha $T_1$ và $T_2$ bão hòa, không phụ thuộc vào trường; khi $1 \le h_5 \le 2$, $T_1$ giảm theo quy luật hàm logarit ($\sim \ln h_5$) và $T_2$ giảm theo quy luật hàm lũy thừa; khi $h_5 < 1$, $T_1$ giảm dần về 0 theo hàm logarit ($\sim \ln h_5$) trong khi $T_2$ không thay đổi và giữ nguyên giá trị nhiệt độ chuyển pha KT của mô hình XY thuần túy ($T_2 \approx 0.89$).

4. Vì sao mô hình Heisenberg cổ điển không mô tả được trật tự sắt từ nhiệt độ phòng của hệ $CdS:C$ trong mô phỏng?

Khi không đưa vào từ trường ngoài cưỡng bức, mô hình spin liên tục Heisenberg cho phép vector spin dao động tự do theo mọi hướng trong không gian 3 chiều. Ở nồng độ pha tạp thấp (5.55%), các thăng giáng nhiệt ở nhiệt độ thấp dễ dàng phá vỡ trật tự định hướng của spin, khiến độ từ hóa suy giảm nhanh về 0 và $T_C$ tính toán được thấp hơn nhiều so với nhiệt độ phòng. Khi chuyển sang mô hình spin rời rạc, bất đẳng hướng tinh thể được phản ánh, hạn chế các bậc tự do thăng giáng và duy trì được trật tự sắt từ ổn định ở nhiệt độ phòng.

---

Kết luận

Luận án tiến sĩ của tác giả Lương Minh Tuấn là công trình nghiên cứu chuyên sâu về mô phỏng số Monte Carlo trong vật lý thống kê và khoa học vật liệu từ tính. Tác giả đã giải quyết dứt điểm các tranh cãi học thuật về bản chất chuyển pha Kosterlitz-Thouless trong mô hình 2D 5-state clock thông qua việc đề xuất sử dụng đại lượng đạo hàm tham số Binder; xây dựng hoàn chỉnh giản đồ pha và quy luật phụ thuộc trường tinh thể cho mô hình 2D $XYh_q$ ($q = 3, 5, 6, 8$); đồng thời chứng minh tính ưu việt của mô hình spin rời rạc so với mô hình spin liên tục trong việc giải thích trật tự sắt từ nhiệt độ phòng ở chất bán dẫn từ pha loãng $CdS:C$. Các kết quả tính toán có tính hệ thống cao, đóng góp dữ liệu định lượng quan trọng cho lý thuyết chuyển pha trong các hệ spin thấp chiều và định hướng tính toán vật liệu spintronics.

Trích đoạn nội dung tài liệu

CHƯƠNG 1: TỔNG QUAN 1. Vật liệu từ 1. Giới thiệu và lịch sử của từ học Câu chuyện về từ tính bắt đầu từ một loại khoáng vật gọi là Magnetite (Fe3O4) (Hình 1. 1), đây là vật liệu từ tính đầu tiên được con người biết đến.

Lịch sử bắt đầu của nó rất mơ hồ, nhưng khả năng hút sắt của nó chắc chắn đã được biết đến 2500 năm trước công nguyên. Trong thế giới cổ đại, khoáng vật này xuất hiện nhiều ở khu vực Magnesia, của Thổ Nhĩ Kì ngày nay, và từ ”sắt từ” có nguồn gốc từ tiếng Hy Lạp, được cho là tên của khu vực này. Những đối thoại đầu tiên về từ học và sử dụng nam châm vĩnh cửu (đá thiên nhiên có từ tính) cho các mục đích khác nhau xuất hiện ở Hy Lạp [39], và sử dụng đá nam châm làm la bàn ở Trung Quốc [40]. Vật liệu từ tính Magnetite [41].

Nghiên cứu khoa học thực thụ đầu tiên về từ tính mãi đến năm 1600 mới được thực hiện bởi William Girlbert trong cuốn sách ông viết về điện và từ xuất bản ở Anh. Tiếp đó là phát hiện vĩ đại của Oersted (1775-1851) năm 1820 về tương tác giữa giữa nam châm với dòng điện và dòng điện với nhau đã dẫn đến kết luận dòng điện tạo ra từ trường. Ngay sau đó vào năm 1825 nam châm điện đầu tiên được chế tạo. Sau thí nghiệm Oersted, hàng loạt các nhà khoa học đã tiến hành thí nghiệm và những công trình nghiên cứu của Ampere, Gauss, Faraday,.

đã đưa ra những kiến thức nền tảng về từ học. Maxwell đã thống nhất lí thuyết về từ trường, điện trường, và quang học để phát triển thành lí thuyết tổng quát về trường điện từ. Để đến thế kỉ 20, từ học và các vật liệu từ được phát triển mạnh mẽ cả về mặt nghiên cứu lí thuyết và công nghệ chế tạo vật liệu, từ đó xây dựng các lí thuyết vi mô về hiện tượng từ dựa trên các lí thuyết 6 của cơ học lượng tử và vật lí chất rắn. Đi kèm với nó là sự phát triển của nhiều kĩ thuật chụp ảnh cấu trúc từ và đo đạc các tính chất từ của vật liệu.

Cuối thế kỉ 20, đầu thế kỉ 21, các vật liệu từ hệ thấp chiều được tập trung nghiên cứu và ngành mới spintronics ra đời dựa trên những thành tựu của từ học và điện tử học [42, 41]. Bản chất từ tính của vật liệu Ngay từ những năm 1820 Ampere đã giả thiết rằng từ tính của vật liệu liên quan đến sự tồn tại các dòng điện tròn không tắt dần trong nó. Tới đầu thế kỉ 20 Rutherford xây dựng mô hình hành tinh nguyên tử có các điện tử quay xung quanh một hạt nhân nặng, mang điện tích dương. Theo quan niệm này thì các dòng điện tròn của Ampere sinh ra do các điện tử quay trên các quỹ đạo quanh hạt nhận.

Sau đó, Planck, Bohr, Broglie, Schrodinger, Heisenberg, và nhiều người khác đã đưa ra thuyết lượng tử hoàn thiện thêm về cấu tạo vật chất, trên cơ sở đó làm sáng tỏ hơn bản chất từ tính của vật liệu. Sự hình thành của từ tính: a) mô-men từ quỹ đạo và b) spin của electron [42]. Nếu coi nguyên tử là phần tử nhỏ bé nhất cấu tạo nên các vật thể thì sự hình thành từ tính của nguyên tử chính là nguồn gốc tính chất từ của vật liệu. Vậy chúng ta hãy khảo sát từ tính của nguyên tử, xuất phát từ tính chất từ của điện tử và hạt nhân.

Tính chất từ của nguyên tử được đóng góp bởi ba thành phần. Thứ nhất, mô-men từ quỹ đạo của electron, thành phần từ tính này xuất hiện khi electron quay quanh hạt nhân trên quỹ đạo là một đường tròn bán kính r (Hình 1.