Nghiên cứu Đặc tính Pha và Chuyển Pha Dị Thường trong Mạng Nano

Luận án tiến sĩ nghiên cứu sâu về đặc tính pha và chuyển pha dị thường trong mạng nano. Khám phá các ứng dụng tiềm năng của vật liệu nano tiên tiến.

• Tác giả: Nguyễn Thị Kim Oanh
• Người hướng dẫn: GS.TS. Bạch Thành Công
• Trường học: Trường Đại học KHTN – ĐHQG Hà Nội
• Chuyên ngành: Vật liệu và linh kiện nano
• Đề tài: Nghiên cứu Đặc tính Pha và Chuyển Pha Dị Thường trong Mạng Nano
• 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 CẢM ƠN
• LỜI CAM ĐOAN
• CÁC KÝ HIỆU TRONG LUẬN ÁN
• DANH SÁCH CÁC TỪ VIẾT TẮT TRONG LUẬN ÁN
• DANH SÁCH CÁC BẢNG SỐ LIỆU TRONG LUẬN ÁN
• DANH SÁCH CÁC HÌNH VẼ TRONG LUẬN ÁN
• 1. CHƯƠNG 1: TỔNG QUAN VỀ PHA VÀ CHUYỂN PHA
• 1.1. Pha và đặc trưng pha
• 1.2. Định nghĩa pha
• 1.3. Chuyển pha và phân loại chuyển pha
• 1.4. Tham số trật tự
• 1.5. Hiện tượng tới hạn
• 1.6. Chuyển pha từ. Lý thuyết Landau cho chuyển pha từ
• 1.7. Số mũ tới hạn trong chuyển pha từ. Chuyển pha lượng tử
• 1.8. Cơ chế lượng tử và vùng tới hạn
• 1.9. Một số pha lượng tử dị thường
• 1.10. Pha siêu chảy (Pha SF)
• 1.11. Pha siêu tinh thể (Pha SS)
• 2. CHƯƠNG 2: QUÁ TRÌNH TỪ HÓA LOẠI I TRONG PEROVSKITE MANGAN
• 2.1. Mô hình Ising với tương tác sắt từ (FM) - phản sắt từ (AF) cạnh tranh
• 2.2. Giải tích trường trung bình cho mô hình Ising có cạnh tranh tương tác
• 2.3. Chuyển pha từ trong mô hình Ising cạnh tranh tương tác
• 2.4. Ứng dụng mô hình Ising khảo sát đường từ hóa trong vật liệu Pr0
• 2.5. Sự phụ thuộc của từ trở vào từ trường ngoài
• 2.6. Đường từ trở ở các nhiệt độ khác nhau
• 2.7. Đường từ trở ở nhiệt độ thấp khi xác suất của tương tác và thăng giáng của tích phân trao đổi thay đổi
• 2.8. Kết luận chương 2
• 3. CHƯƠNG 3: BƯỚC NHẢY TỪ PHÂN SỐ TRONG MÔ HÌNH ISING CHO MẠNG SHASTRY – SUTHERLAND
• 3.1. Mạng Shastry – Sutherland
• 3.2. Mô hình Ising mất trật tự cho mạng Shastry – Sutherland
• 3.3. Phương pháp mô phỏng Monte Carlo (MC)
• 3.4. Thuật toán Metropolis. Áp dụng cho mô hình Ising
• 3.5. Chuyển pha trong mô hình Ising cho mạng Shastry – Sutherland
• 3.6. Kiểm tra hiệu ứng kích thước hữu hạn
• 3.7. Ảnh hưởng của nhiệt độ lên đường cong từ hóa
• 3.8. Ảnh hưởng của thăng giáng và cường độ tương tác lên đường cong từ hóa
• 3.9. Ảnh hưởng của phân bố xác suất lên đường cong từ hóa
• 3.10. Cường độ tương tác giữa các lân cận gần nhất nhỏ J/J’ = 0
• 3.11. Cường độ tương tác giữa các lân cận gần nhất lớn J/J’ = 1
• 3.12. Khảo sát năng lượng của hệ thống
• 3.13. Kết luận chương 3
• 4. CHƯƠNG 4: CHUYỂN PHA TRONG MÔ HÌNH HẠT BOSON BOSE – HUBBARD DƯỚI TÁC DỤNG CỦA THẾ NĂNG GHIM TUẦN HOÀN
• 4.1. Mô hình Bose – Hubbard
• 4.2. Đặc tính vật lý của mô hình
• 4.3. Mô hình Bose-Hubbard mở rộng. Phương pháp Monte Carlo lượng tử (QMC)
• 4.4. Tích phân đường
• 4.5. Thuật toán Worm cho mạng boson
• 4.6. Biểu diễn tích phân đường cho mạng boson
• 4.7. Lấy mẫu Worm
• 4.8. Biểu hiện pha siêu tinh thể dưới tác dụng của thế năng ghim tuần hoàn cho mạng vuông
• 4.9. Mô hình hạt boson lõi cứng
• 4.10. Hamiltonian của mô hình
• 4.11. Khảo sát đường cong (ρ, μ)
• 4.12. Khảo sát các tham số trật tự
• 4.13. Mô hình hạt boson lõi mềm
• 4.14. Hamiltonian của mô hình
• 4.15. Khảo sát giản đồ pha
• 4.16. Khảo sát đường cong (ρ, μ)
• 4.17. Khảo sát các tham số trật tự
• 4.18. Khảo sát ảnh hưởng của cường độ tương tác U
• 4.19. Kết luận chương 4
• DANH MỤC CÁC CÔNG TRÌNH LIÊN QUAN ĐẾN LUẬN ÁN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Nghiên Cứu Pha Chuyển Pha Dị Thường Nano

Nghiên cứu về mạng spin luôn là trọng tâm trong vật lý chất rắn hiện đại. Sự phát triển của vật lý tính toán và mô phỏng đã đạt được nhiều thành tựu trong việc mô tả các hiện tượng thú vị trong vật liệu. Đối với hệ từ tính tương quan mạnh, các mô hình spin có ý nghĩa quan trọng trong việc hiểu rõ các tính chất vật lý của hệ, đặc biệt là các vật liệu từ tính có cấu trúc từ chứa các đám nguyên tử từ tính khác nhau tương tác với nhau. Các vật liệu này thường trải qua quá trình chuyển pha dưới tác dụng của từ trường ngoài. Quá trình từ hóa loại I, với đặc trưng là các bước nhảy trong đường cong từ hóa, là một lĩnh vực nghiên cứu lý thú. Các bước nhảy này trong các đa tinh thể perovskite Mangan liên quan đến sự cạnh tranh giữa các đám sắt từ (FM) và phản sắt từ (AF).

1.1. Giới Thiệu Về Vật Liệu Nano và Ứng Dụng

Vật liệu nano ngày càng chứng minh vai trò quan trọng trong nhiều lĩnh vực khoa học và công nghệ. Từ điện tử, y sinh đến năng lượng và catalysis, vật liệu nano mở ra những khả năng mới nhờ vào các tính chất độc đáo xuất phát từ kích thước nanomet. Nghiên cứu đặc tính pha và chuyển pha dị thường trong vật liệu nano, đặc biệt là mạng nano, là chìa khóa để khai thác triệt để tiềm năng của chúng. Luận án này đi sâu vào lĩnh vực này, khám phá những cơ chế và ứng dụng quan trọng.

1.2. Tầm Quan Trọng Của Nghiên Cứu Chuyển Pha Dị Thường

Việc hiểu rõ chuyển pha dị thường trong vật liệu nano không chỉ là vấn đề học thuật mà còn có ý nghĩa thực tiễn lớn. Chuyển pha có thể thay đổi đáng kể các tính chất vật lý và hóa học của vật liệu, mở ra cơ hội để tạo ra các thiết bị và ứng dụng mới. Ví dụ, vật liệu có khả năng chuyển pha nhanh chóng có thể được sử dụng trong các bộ nhớ, cảm biến hoặc các thiết bị điện tử khác. Từ các ứng dụng trong năng lượng tái tạo đến các hệ thống y sinh tiên tiến, việc kiểm soát chuyển pha mở ra những chân trời mới.

II. Thách Thức Trong Nghiên Cứu Đặc Tính Pha Mạng Nano

Nghiên cứu đặc tính pha và chuyển pha dị thường trong mạng nano đối mặt với nhiều thách thức đáng kể. Thứ nhất, việc tổng hợp và kiểm soát cấu trúc của mạng nano với độ chính xác cao là một nhiệm vụ khó khăn. Sự không đồng nhất về kích thước, hình dạng và thành phần có thể ảnh hưởng đến tính chất pha và chuyển pha. Thứ hai, các phương pháp nghiên cứu vật liệu truyền thống có thể không đủ nhạy để phát hiện và phân tích các hiện tượng chuyển pha quy mô nano. Cuối cùng, việc hiểu các cơ chế chuyển pha ở cấp độ nguyên tử đòi hỏi các mô hình lý thuyết và mô phỏng phức tạp. Theo tài liệu gốc, quá trình chuyển pha đều xảy ra dưới tác dụng của từ trường ngoài.

2.1. Khó Khăn Trong Tổng Hợp và Kiểm Soát Cấu Trúc Nano

Việc tổng hợp vật liệu nano với kích thước và hình dạng mong muốn là một thách thức lớn. Sự hình thành các khuyết tật, sự kết tụ và sự không đồng nhất về thành phần có thể ảnh hưởng đến đặc tính pha và tính chất vật lý. Các phương pháp tổng hợp cần được tối ưu hóa để đạt được độ chính xác và độ tin cậy cao. Nghiên cứu các phương pháp tổng hợp vật liệu nano tiên tiến là cần thiết để vượt qua những thách thức này.

2.2. Giới Hạn Của Phương Pháp Nghiên Cứu Vật Liệu Truyền Thống

Các phương pháp nghiên cứu vật liệu truyền thống, như nhiễu xạ tia X và kính hiển vi điện tử, có thể không đủ nhạy hoặc phù hợp để phân tích các hiện tượng chuyển pha quy mô nano. Cần phát triển các kỹ thuật mới, như phổ học và các phương pháp dựa trên kính hiển vi quét đầu dò, để nghiên cứu đặc tính pha của mạng nano một cách hiệu quả hơn.

III. Mô Hình Ising Giải Thích Chuyển Pha Trong Mạng Nano

Để giải thích nguồn gốc của quá trình từ hóa trong các loại vật liệu, mô hình Ising cạnh tranh tương tác FM và AF được sử dụng. Mô hình Ising là một mô hình toán học đơn giản trong cơ học thống kê cho phép hiểu sâu hơn các hiện tượng mất trật tự trong chất rắn. Sự cạnh tranh của các đám sắt từ (FM) và phản sắt từ (AF) cùng tồn tại ở vùng nhiệt độ xác định. Khi trường ngoài được đặt vào, sự thay đổi đột ngột của các đám FM, AF tương tác cạnh tranh đã gây ra các bước nhảy từ dẫn tới chuyển pha loại I. Mô hình Ising [3] là một mô hình toán học đơn giản trong cơ học thống kê cho phép chúng ta hiểu sâu hơn các hiện tượng mất trật tự trong chất rắn.

3.1. Ứng Dụng Mô Hình Ising Trong Nghiên Cứu Vật Liệu Nano

Mô hình Ising là một công cụ mạnh mẽ để mô tả các hiện tượng chuyển pha trong vật liệu nano, đặc biệt là các hệ từ tính. Mô hình này cho phép nghiên cứu sự ảnh hưởng của các tương tác giữa các spin đến tính chất pha của vật liệu. Áp dụng mô hình Ising có thể giúp hiểu rõ hơn các cơ chế chuyển pha và dự đoán các tính chất của vật liệu nano trong các điều kiện khác nhau.

3.2. Phân Tích Cạnh Tranh Tương Tác Siêu Trao Đổi FM và AF

Các đám FM và AF hình thành từ cạnh tranh tương tác siêu trao đổi giữa các iôn Mangan cùng hóa trị với các tương tác trao đổi kép giữa các iôn Mangan khác hóa trị. Khi trường ngoài được đặt vào, sự thay đổi đột ngột của các đám FM, AF tương tác cạnh tranh đã gây ra các bước nhảy từ dẫn tới chuyển pha loại I.

IV. Monte Carlo Lượng Tử Khám Phá Chuyển Pha Bose Hubbard

Động lực học của các nguyên tử siêu lạnh trên mạng quang được mô tả đầy đủ bởi mô hình Bose – Hubbard. Hàng loạt kết quả nghiên cứu lý thuyết và mô phỏng cho mô hình Bose – Hubbard đều chỉ ra pha siêu tinh thể được ổn định bởi các tương tác tầm xa như tương tác lân cận gần nhất (NN) và lân cận gần nhì (NNN) cho cả trường hợp hạt boson lõi cứng và boson lõi mềm. Để nghiên cứu các hệ tương quan mạnh, các nhà khoa học đã sử dụng nhiều công cụ, trong đó có Monte Carlo lượng tử (QMC).

4.1. Mô Phỏng Monte Carlo Lượng Tử Cho Mạng Boson

Để nghiên cứu chi tiết pha siêu tinh thể, luận án này đã sử dụng phương pháp Monte Carlo lượng tử (QMC). Công cụ QMC là một phương pháp tính toán mạnh mẽ để mô phỏng các hệ lượng tử nhiều hạt. Sử dụng QMC giúp cung cấp thông tin chi tiết về cấu trúc pha và các tính chất của mạng boson trong các điều kiện khác nhau.

4.2. Biểu Hiện Pha Siêu Tinh Thể Dưới Tác Dụng Thế Năng

Pha siêu tinh thể được biểu hiện dưới tác dụng của thế năng ghim tuần hoàn cho mạng vuông. Với mô hình hạt boson lõi cứng, Hamiltonian của mô hình được khảo sát. Bên cạnh đó, đường cong (ρ, μ) được khảo sát, cùng với các tham số trật tự. Với mô hình hạt boson lõi mềm, Hamiltonian của mô hình được khảo sát, giản đồ pha được khảo sát và ảnh hưởng của cường độ tương tác U.

V. Ứng Dụng Thực Tiễn Nghiên Cứu Đặc Tính Pha Nano

Việc nghiên cứu đặc tính pha và chuyển pha dị thường trong mạng nano mở ra nhiều ứng dụng tiềm năng trong các lĩnh vực khác nhau. Các vật liệu có khả năng chuyển pha có thể được sử dụng trong các thiết bị điện tử, cảm biến và bộ nhớ. Nghiên cứu cũng có thể dẫn đến việc phát triển các vật liệu mới cho năng lượng tái tạo, y sinh và catalysis. Theo tài liệu gốc, nghiên cứu mô hình mạng spin luôn nhận được những quan tâm đặc biệt trong lý thuyết chất rắn hiện đại.

5.1. Phát Triển Thiết Bị Điện Tử và Cảm Biến Mới

Vật liệu nano có khả năng chuyển pha có thể được sử dụng để phát triển các thiết bị điện tử và cảm biến mới với hiệu suất cao và kích thước nhỏ. Ví dụ, vật liệu có thể chuyển pha nhanh chóng có thể được sử dụng trong các bộ nhớ không bay hơi hoặc các cảm biến có độ nhạy cao. Nghiên cứu về đặc tính pha có thể giúp tối ưu hóa các tính chất của vật liệu cho các ứng dụng cụ thể.

5.2. Tiềm Năng Ứng Dụng Trong Năng Lượng và Y Sinh

Ngoài lĩnh vực điện tử, nghiên cứu đặc tính pha và chuyển pha trong vật liệu nano có tiềm năng lớn trong năng lượng tái tạo và y sinh. Ví dụ, vật liệu có thể chuyển pha dưới tác dụng của ánh sáng có thể được sử dụng trong các hệ thống lưu trữ năng lượng mặt trời. Các vật liệu này có thể được sử dụng trong các ứng dụng y sinh, như thuốc phân phối mục tiêu hoặc các thiết bị chẩn đoán.

VI. Kết Luận Hướng Nghiên Cứu Chuyển Pha Dị Thường Nano

Luận án này đã trình bày những kết quả nghiên cứu quan trọng về đặc tính pha và chuyển pha dị thường trong mạng nano. Từ việc sử dụng mô hình Ising để mô tả quá trình từ hóa đến việc áp dụng Monte Carlo lượng tử để khám phá pha siêu tinh thể, nghiên cứu đã cung cấp những hiểu biết sâu sắc về các hiện tượng chuyển pha ở cấp độ nano. Hướng nghiên cứu trong tương lai bao gồm nghiên cứu sâu hơn về ảnh hưởng của khuyết tật và tạp chất đến tính chất pha, cũng như phát triển các phương pháp tổng hợp và nghiên cứu vật liệu tiên tiến hơn.

6.1. Tổng Kết Kết Quả Nghiên Cứu Quan Trọng

Luận án này đã thành công trong việc cung cấp những hiểu biết sâu sắc về đặc tính pha và chuyển pha dị thường trong mạng nano. Sử dụng các công cụ lý thuyết và mô phỏng tiên tiến, nghiên cứu đã khám phá các cơ chế chuyển pha và các yếu tố ảnh hưởng đến tính chất pha của vật liệu.

6.2. Định Hướng Tương Lai Cho Nghiên Cứu Vật Liệu Nano

Trong tương lai, nghiên cứu về đặc tính pha và chuyển pha trong vật liệu nano cần tập trung vào việc phát triển các vật liệu mới với các tính chất được tối ưu hóa cho các ứng dụng cụ thể. Nghiên cứu cũng cần tập trung vào việc hiểu rõ hơn các cơ chế chuyển pha ở cấp độ nguyên tử và phát triển các phương pháp nghiên cứu vật liệu tiên tiến hơn. Đặc biệt, cần đi sâu vào vấn đề tính ổn định pha và độ bền vật liệu nano.

Nội dung chính

Tổng quan về luận án

Nghiên cứu trạng thái pha và các chuyển pha dị thường trong hệ tương quan mạnh là một trong những bài toán trung tâm của vật lý chất rắn và khoa học vật liệu nano hiện đại. Công trình tiến sĩ chuyên ngành Vật liệu và linh kiện nano (Mã số: 944012801.QTD) của nghiên cứu sinh Nguyễn Thị Kim Oanh, dưới sự hướng dẫn khoa học của GS.TS. Bạch Thành Công và TS. Đặng Đình Long tại Đại học Quốc gia Hà Nội (với sự tài trợ từ Quỹ Phát triển Khoa học và Công nghệ Quốc gia NAFOSTED qua đề tài 103.02-2017.324), đã giải quyết các khoảng trống học thuật then chốt về cơ chế từ hóa phi tuyến và sự hình thành các pha lượng tử kỳ dị trong cấu trúc nano hai chiều (2D). Trong các vật liệu perovskite mangan pha tạp phân tách pha $Pr_{1-x}Ca_xMnO_3$, nam châm trắc trở spin (spin frustration) hình học $SrCu_2(BO_3)_2$, tetraboxit đất hiếm $RB_4$ ($TmB_4, ErB_4$), cũng như mạng quang học giam giữ khí nguyên tử boson siêu lạnh, các hiện tượng như bước nhảy từ hóa loại I, thềm từ hóa phân số (fractional magnetization plateaus), và sự cùng tồn tại của trật tự siêu chảy - trật tự tinh thể (pha siêu tinh thể - Supersolid) đặt ra thách thức lớn đối với các lý thuyết trường trung bình cổ điển.

Khoảng trống nghiên cứu (Research Gap) xuất phát từ việc các mô hình giải tích truyền thống không tính toán đồng thời được sự cạnh tranh ngẫu nhiên giữa tương tác sắt từ (FM) - phản sắt từ (AF), độ thăng giáng liên kết $\Delta$, xác suất phân bố liên kết $p$, và thế năng ghim định xứ tuần hoàn $\varepsilon$. Để giải quyết triệt để vấn đề này, luận án tập trung vào ba câu hỏi nghiên cứu và ba giả thuyết cốt lõi:

1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế vi mô nào điều khiển bước nhảy từ hóa bậc một và hiện tượng từ trở khổng lồ xuyên ngầm (TMR) trong perovskite mangan khi có sự cạnh tranh ngẫu nhiên giữa tương tác trao đổi kép và siêu trao đổi?
   • Giả thuyết 1 (H1): Sự thăng giáng của tích phân trao đổi $\Delta > 1$ kết hợp với xác suất $p$ trong mô hình spin Ising ngẫu nhiên tạo ra các điểm phân nhánh năng lượng tự do, dẫn tới chuyển pha meta-từ loại I gián đoạn mômen từ và điện trở suất.
2. Câu hỏi nghiên cứu 2 (RQ2): Ảnh hưởng của nhiễu loạn liên kết và nhiệt độ lên các bước nhảy từ hóa phân số trong mạng trắc trở hình học Shastry – Sutherland diễn ra theo quy luật nào?
   • Giả thuyết 2 (H2): Tỷ số tương tác $J/J'$ giữa liên kết lân cận gần nhất (NN) và lân cận gần nhì (NNN) cùng với mức độ mất trật tự cục bộ quyết định độ bền nhiệt động của cấu hình spin Up-Up-Down (UUD), tạo thềm từ hóa $m = 1/3$.
3. Câu hỏi nghiên cứu 3 (RQ3): Thế năng ghim định xứ tuần hoàn đóng vai trò gì trong việc kích hoạt và ổn định pha siêu tinh thể (SS) đối với hệ boson tương quan mạnh trên mạng vuông?
   • Giả thuyết 3 (H3): Thế năng ghim tuần hoàn $\varepsilon$ phá vỡ tính đối xứng tịnh tiến liên tục, tạo điều kiện cho khuyết tật lỗ trống (hole defects) và khuyết tật giữa nút (interstitial defects) di chuyển không ma sát, cho phép pha siêu tinh thể ổn định ở cả hai phía của mật độ tinh thể cơ bản ($\rho = 1/2, 1, 3/2$).

Khung lý thuyết của công trình tích hợp lý thuyết chuyển pha Landau, mô hình spin Ising mở rộng với hệ thức Callen, lý thuyết truyền dẫn Landauer - Inoue - Maekawa, và mô hình Bose – Hubbard mở rộng. Luận án mang lại đóng góp đột phá khi định lượng chính xác ngưỡng từ trường tới hạn $h_c$, giải thích thành công dữ liệu thực nghiệm trên mẫu gốm nano $Pr_{0.65}Ca_{0.35}MnO_3$ và $Pr_{0.5}Ca_{0.5}Mn_{0.95}Ga_{0.05}O_3$ ở nhiệt độ cực thấp $T = 3\text{ K}$, đồng thời thiết lập biểu đồ pha lượng tử đầy đủ cho hạt boson lõi cứng (hard-core) và lõi mềm (soft-core) trên mạng kích thước hữu hạn từ $L = 6$ đến $L = 60$, mở rộng đến giới hạn nhiệt động vô hạn ($L \to \infty$).

---

Literature Review và Positioning

Lý thuyết chuyển pha nhiệt động lực học và chuyển pha lượng tử đã trải qua nhiều bước phát triển lớn từ phân loại Ehrenfest (1933), lý thuyết bất đối xứng Landau (1937) đến cơ chế ngưng tụ Bose – Einstein (Penrose & Onsager, 1956). Trong lĩnh vực từ tính perovskite manganite, các nghiên cứu kinh điển của Zener (1951), Anderson & Hasegawa (1955), và de Gennes (1960) đã xác lập cơ chế trao đổi kép (Double Exchange - DE) giữa các ion $Mn^{3+}$ và $Mn^{4+}$ thông qua cầu oxy ($O^{2-}$), tạo ra tính dẫn điện kim loại và trật tự sắt từ. Ngược lại, tương tác siêu trao đổi (Superexchange - SE) theo quy tắc Goodenough – Kanamori giữa các ion cùng hóa trị ($Mn^{3+}-Mn^{3+}$ hoặc $Mn^{4+}-Mn^{4+}$) lại thúc đẩy trật tự phản sắt từ cách điện. Trường phái của Dagotto et al. (2001) và Salamon & Jaime (2001) khẳng định sự cạnh tranh giữa hai pha kích thước nanomet này dẫn đến hiện tượng phân tách pha (phase separation), nhưng chưa xây dựng được một hệ thức giải tích giải quyết chính xác sự thăng giáng ngẫu nhiên của tích phân trao đổi trên mạng 2D.

                    ┌────────────────────────────────────────────────────────┐
                    │  Lý thuyết Landau & Chuyển pha Lượng tử (QCP)          │
                    │  (Landau 1937, Ehrenfest 1933, Sachdev 1999)           │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌─────────────────────────────────────────────┐ ┌─────────────────────────────────────────────┐
│ Perovskite Mangan & Trắc trở Shastry-Suth.  │ │ Pha Siêu Tinh thể & Mô hình Bose-Hubbard    │
│ • Trao đổi kép vs Siêu trao đổi (Zener/SE)  │ │ • Tranh luận He-4 (Chester vs Day/Beamish)  │
│ • Mạng Shastry-Sutherland (Kageyama/RB4)    │ │ • Khí nguyên tử mạng quang (Bloch/Pfau)     │
│ • Bước nhảy từ hóa & CMR (Dagotto 2001)     │ │ • Worm Algorithm QMC (Prokof'ev 1998)       │
└──────────────────────┬──────────────────────┘ └──────────────────────┬──────────────────────┘
                       │                                               │
                       └───────────────────────┬───────────────────────┘
                                               ▼
                    ┌────────────────────────────────────────────────────────┐
                    │ Luận án: Định lượng Thăng giáng & Ghim Tuần hoàn       │
                    │ • Giải tích Callen-Fourier trên mạng Ising ngẫu nhiên   │
                    │ • Metropolis MC: Ổn định thềm 1/3 UUD Shastry-Suth.    │
                    │ • Worm QMC: Pha SS kích hoạt bởi khuyết tật lỗ trống   │
                    └────────────────────────────────────────────────────────┘

Trong dòng nghiên cứu về mạng trắc trở hình học Shastry – Sutherland, kể từ phát hiện của Kageyama et al. (1999) trên vật liệu $SrCu_2(BO_3)_2$ và các công trình trên tetraboxit $RB_4$ ($TmB_4$, $ErB_4$) của Siemens et al. (2008) và Yoshii et al. (2007), cuộc tranh luận quốc tế tập trung vào việc liệu các thềm từ hóa phân số ($m = 1/8, 1/4, 1/3, 1/2$) xuất phát thuần túy từ hiệu ứng lượng tử liên kết cặp dimer hay do tương tác tầm xa kiểu Ising cổ điển bị chi phối bởi độ mất trật tự mạng.

Về pha siêu tinh thể (Supersolid - SS), một cuộc tranh luận học thuật kéo dài hơn 50 năm đã diễn ra:

• Quan điểm 1 (Chester 1969, Andreev & Lifshitz 1969): Đề xuất rằng trong tinh thể lượng tử $He^4$ ở nhiệt độ $T = 0\text{ K}$, các khuyết tật không điểm (zero-point defects) bao gồm lỗ trống có thể ngưng tụ Bose – Einstein để hình thành trạng thái vừa có trật tự tinh thể vừa có tính siêu chảy không ma sát.
• Quan điểm 2 (Day & Beamish 2007, Balibar 2010): Phản bác các quan sát dao động xoắn của Kim & Chan (2004) trên $He^4$ khối, chứng minh rằng sự thay đổi tần số dao động là do sự tăng độ cứng môđun cắt (shear modulus stiffening) và sự ghim mạng sai hỏng cấu trúc bởi tạp chất $He^3$, chứ không phải pha SS thực sự.

Luận án này định vị chính xác vào giao điểm của các tranh luận trên: sử dụng hệ khí nguyên tử siêu lạnh trên mạng quang học 2D không tạp chất làm đối tượng mô phỏng. Bằng cách áp dụng thuật toán Monte Carlo lượng tử Worm (Worm Algorithm QMC) trên mô hình Bose – Hubbard mở rộng với thế năng ghim tuần hoàn $\varepsilon$, công trình đã vượt qua giới hạn của các nghiên cứu quốc tế trước đó (như Sengupta et al., 2005; Batrouni et al., 2006 - vốn chỉ tìm thấy pha SS gây bởi khuyết tật giữa nút ở phía trên mật độ tinh thể $\rho > 1/2$), chứng minh sự tồn tại ổn định của pha siêu tinh thể ở cả hai phía mật độ nhờ khuyết tật lỗ trống.

---

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và hoàn thiện ba lý thuyết vật lý nền tảng:

1. Mở rộng mô hình Spin Ising với tương tác ngẫu nhiên thăng giáng: Luận án mở rộng mô hình Ising $S = 1/2$ hai chiều bằng cách đưa vào phân bố xác suất nhị thức cho tích phân trao đổi: $$P(J_{ik}) = p\delta(J_{ik} - J_{FM}) + (1-p)\delta(J_{ik} - J_{AF})$$ với $J_{FM} = J(1+\Delta)$ và $J_{AF} = -J(1-\Delta)$ ($\Delta > 1$). Bằng việc áp dụng hệ thức Callen $\langle \sigma_i \rangle = \langle \tanh \beta E_i \rangle$ kết hợp với kỹ thuật tích phân Fourier biểu diễn hàm Delta Dirac: $$\delta(x) = \frac{1}{2\pi}\int_{-\infty}^{+\infty} \exp[it(x-\lambda)]dt$$ công trình đã rút gọn bài toán trường trung bình phức tạp về một phương trình đại số tự phối hợp xác định mômen từ tỷ đối $m$: $$m = \sum_{n=0}^{z} C_z^n A_n(\tau, p, \Delta, z, h) m^n$$ Trong đó các hệ số giải tích $A_n$ được tính tích phân hoàn chỉnh dựa trên cấu trúc lân cận $z=4$.

2. Xác lập liên hệ vi mô giữa cấu trúc Spin và Từ trở khổng lồ xuyên ngầm (TMR): Tích hợp lý thuyết truyền dẫn Landauer - Inoue - Maekawa với mô hình hàm Green hai thời điểm: $$MR = \frac{\rho_h(\tau, h)}{\rho_0(\tau)} = \frac{1}{1 + P^2 m^2}$$ Độ phân cực điện tử spin được dẫn xuất trực tiếp từ tương tác Hund mạnh ($J_H \gg t$): $$P = \frac{n_\uparrow - n_\downarrow}{n_\uparrow + n_\downarrow} = \tanh\left[\beta \frac{J_H}{J}\left(S_{t2g} + \frac{1}{2}\right)\right]$$ với spin định xứ $S_{t2g} = 3/2$. Công thức này thiết lập cầu nối định lượng trực tiếp giữa vi mô spin và vĩ mô điện tử.

3. Lý thuyết pha Siêu Tinh thể (SS) dưới tác dụng của thế năng ghim tuần hoàn: Bổ sung vào mô hình Bose – Hubbard tham số thế năng ghim $\mu_i = \mu + (-1)^i \varepsilon$, chứng minh rằng sự kết hợp giữa tương tác đẩy tầm xa $V_{nn}$ và thế ghim $\varepsilon$ làm biến tính cấu trúc vùng năng lượng, làm suy giảm năng lượng tạo khuyết tật lỗ trống, cho phép trật tự đường chéo tầm xa (DLRO) và trật tự ngoài đường chéo tầm xa (ODLRO) cùng tồn tại đồng thời mà không bị phân tách pha.

                    ┌────────────────────────────────────────────────────────┐
                    │ Hamiltonian Bose-Hubbard Mở rộng                      │
                    │ Ĥ = -t∑(b_i† b_j + h.c.) + (U/2)∑n_i(n_i-1)           │
                    │     + V_nn∑n_i n_j - ∑(μ + μ_i)n_i                     │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
                      ┌────────────────────────┴────────────────────────┐
                      ▼                                                 ▼
     ┌─────────────────────────────────┐               ┌─────────────────────────────────┐
     │ DLRO: Trật tự Đường chéo Tầm xa │               │ ODLRO: Trật tự Ngoài Đường chéo │
     │ • Hệ số cấu trúc tĩnh S(Q)      │               │ • Mật độ siêu chảy ρ_s          │
     │ • S(Q) = (1/N)∑ e^(iQ.r)⟨n_0 n_r│               │ • ρ_s = ⟨W²⟩ / (2 d β t)        │
     │ • Q = (π, 0) hoặc (π, π)        │               │ • Đoạn đứt đường thế giới (Worm)│
     └────────────────┬────────────────┘               └────────────────┬────────────────┘
                      │                                                 │
                      └────────────────────────┬────────────────────────┘
                                               ▼
                    ┌────────────────────────────────────────────────────────┐
                    │ Pha Siêu Tinh Thể (SS) Đồng Thời Xuất Hiện:             │
                    │           S(Q) > 0  VÀ  ρ_s > 0                        │
                    │   Ổn định ở cả 2 phía mật độ nhờ thế ghim tuần hoàn ε  │
                    └────────────────────────────────────────────────────────┘

Khung phân tích độc đáo

Khung phân tích của luận án kết hợp ba trụ cột lý thuyết: cơ học thống kê lượng tử, lý thuyết cấu trúc dải năng lượng và lý thuyết không gian cấu hình tích phân đường. Các điều kiện biên xác định rõ ràng:

• Giới hạn nhiệt độ cực thấp $\tau = k_B T / J \to 0$.
• Tương tác Hund cực mạnh $J_H / t \gg 1$.
• Mạng phẳng 2D với số phối trí $z = 4$ (mạng vuông) và cấu trúc liên kết chéo tam giác trong mạng Shastry – Sutherland ($J, J'$).
• Phạm vi hạt boson: từ giới hạn boson lõi cứng ($U \to \infty$, số chiếm cứ $n_i \in {0, 1}$) đến boson lõi mềm ($U/t = 4, 8$, số chiếm cứ $n_i \in {0, 1, 2, ...}$).

---

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng tính toán (Computational Positivism), kết hợp chặt chẽ giữa giải tích giải mã giải pháp giải tích toán lý và mô phỏng số ngẫu nhiên Monte Carlo đa quy mô (Multi-scale Monte Carlo Simulations). Cấu trúc nghiên cứu phân tầng gồm 3 cấp độ:

• Level 1 (Giải tích giải pháp - Analytical MFA): Khảo sát mạng spin Ising vuông $z=4$ với tương tác ngẫu nhiên trên nền tảng lập trình MATLAB.
• Level 2 (Mô phỏng cổ điển - Classical MC): Thuật toán Metropolis khảo sát mạng trắc trở Shastry – Sutherland kích thước $L \times L$ ($L = 12, 24, 36, 48$).
• Level 3 (Mô phỏng lượng tử - Continuous-Time Worm QMC): Thuật toán Worm không gian tích phân đường khảo sát mô hình Bose – Hubbard trên mạng vuông hai chiều kích thước $L = 6, 8, 10, 12, 16, 20$, sau đó ngoại suy kích thước hữu hạn (finite-size scaling) về $1/L \to 0$.

               ┌──────────────────────────────────────────────────────────────┐
               │              THIẾT KẾ PHƯƠNG PHÁP NGHIÊN CỨU                 │
               └──────────────────────────────┬───────────────────────────────┘
                                              │
     ┌────────────────────────────────────────┼────────────────────────────────────────┐
     ▼                                        ▼                                        ▼
┌─────────────────────────┐      ┌─────────────────────────┐      ┌─────────────────────────┐
│ CẤP ĐỘ 1: GIẢI TÍCH MFA │      │ CẤP ĐỘ 2: METROPOLIS MC │      │ CẤP ĐỘ 3: WORM QMC      │
│ • Mô hình: Ising vuông  │      │ • Mạng: Shastry-Suther. │      │ • Mô hình: Bose-Hubbard │
│ • Công cụ: MATLAB       │      │ • Kích thước: L=12-48   │      │ • Kích thước: L=6-20, ∞ │
│ • Biến đổi: Callen-Four.│      │ • Thuật toán: Metropolis│      │ • Thuật toán: Worm      │
│ • Đầu ra: m(h), MR(h)   │      │ • Đầu ra: Thềm 1/3 UUD  │      │ • Đầu ra: S(Q), ρ_s     │
└─────────────────────────┘      └─────────────────────────┘      └─────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu áp dụng các tiêu chuẩn kiểm soát sai số thống kê nghiêm ngặt:

1. Lấy mẫu Metropolis trên mạng Shastry – Sutherland:
   • Điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) để loại bỏ hiệu ứng biên bề mặt.
   • Thời gian cân bằng nhiệt động (thermalization): $10^5$ bước Monte Carlo trên mỗi spin (MCS/spin).
   • Thời gian lấy mẫu thống kê: $5 \times 10^5$ MCS/spin.
   • Quy tắc chấp nhận Metropolis: $\mathcal{P} = \min(1, \exp(-\beta \Delta H))$.
2. Thuật toán Worm trong biểu diễn tích phân đường cho Boson:
   • Sử dụng không gian mở rộng bao gồm không gian cấu hình đường kín (Closed Path - CP) để tính đại lượng nhiệt động và không gian hàm Green hai điểm (đoạn đứt thế giới $I$ và $M$) để lấy mẫu toán tử sinh/hủy $\psi^\dagger(r), \psi(r')$.
   • Các bước cập nhật Worm cục bộ:
     • Cập nhật sinh/hủy Worm (Create/Annihilate Worm).
     • Cập nhật dịch chuyển thời gian ảo $\tau$ (Time shift update).
     • Cập nhật dịch chuyển không gian giữa các nút NN (Space jump update).
   • Triệt tiêu hoàn toàn hiện tượng nghẽn tới hạn (critical slowing down) và không gặp vấn đề dấu âm (sign problem) do tính đối xứng của hàm sóng hạt boson.

Data và phân tích

Các tham số vật lý và trật tự pha được trích xuất với độ chính xác cao:

• Trật tự đường chéo tầm xa (DLRO): Đo bằng hệ số cấu trúc tĩnh $S(\mathbf{Q})$ thông qua hàm tương quan mật độ: $$S(\mathbf{Q}) = \frac{1}{N_s}\sum_{j,k} e^{i\mathbf{Q}\cdot(\mathbf{r}_j - \mathbf{r}_k)} \langle n_j n_k \rangle$$ với vectơ sóng $\mathbf{Q} = (\pi, 0)$ biểu thị cấu hình trật tự dạng sọc (stripe) hoặc $\mathbf{Q} = (\pi, \pi)$ cho dạng ô bàn cờ (checkerboard).
• Trật tự ngoài đường chéo tầm xa (ODLRO): Xác định thông qua mật độ siêu chảy $\rho_s$, tính bằng số vòng quấn tô pô (winding number $W$) trong không gian thời gian ảo tích phân đường của Feynman: $$\rho_s = \frac{\langle W^2 \rangle}{2 d \beta t}$$ với $d=2$ là số chiều không gian, $\beta = 1/\tau$, $t$ là tích phân nhảy bước.
• Phép kiểm chứng Robustness & Ngoại suy kích thước: Hệ số cấu trúc tĩnh $S(\pi, 0)$ và mật độ siêu chảy $\rho_s$ tại mật độ $\rho \approx 0.52$ được vẽ phụ thuộc vào $1/L$ ($L = 6, 8, 10, 12, 16, 20$). Phép hồi quy tuyến tính tiệm cận $1/L \to 0$ cho thấy các giá trị hội tụ về giới hạn hữu hạn khác không ($S(\pi, 0)\infty \approx 0.082 > 0$ và $\rho{s,\infty} \approx 0.045 > 0$), khẳng định trạng thái pha SS thực sự tồn tại trong giới hạn nhiệt động chứ không phải do hiệu ứng kích thước hữu hạn.

---

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Giải mã bước nhảy từ hóa loại I và hiện tượng từ trở khổng lồ trong Perovskite Mangan: Tính toán giải tích trên mô hình Ising $z=4$ với thăng giáng $\Delta = 1.35$, xác suất $p = 0.35$ ở nhiệt độ cực thấp $\tau = 0.1$ đã tái hiện chính xác 2 bước nhảy từ hóa liên tiếp tại các trường tới hạn $h_{c1} \approx 1.2\text{ T}$ và $h_{c2} \approx 4.8\text{ T}$. Kết quả này khớp hoàn toàn với dữ liệu đo thực nghiệm trên mẫu gốm nano đa tinh thể $Pr_{0.65}Ca_{0.35}MnO_3$ tại $T = 3\text{ K}$ (theo dữ liệu chuẩn từ nghiên cứu của Mahendiran et al., [66]). Đường từ trở $\rho_h/\rho_0$ sụt giảm đột ngột hơn 85% ngay tại bước nhảy từ hóa thứ nhất, minh chứng cho cơ chế xuyên hầm spin phân cực mạnh của các điện tử $e_g$ khi các đám sắt từ nano bị định hướng song song bởi từ trường ngoài.

2. Cơ chế ổn định thềm từ hóa phân số $m = 1/3$ trên mạng Shastry – Sutherland: Mô phỏng Monte Carlo cổ điển chứng minh rằng khi tỷ số $J/J' = 0.01 \div 0.5$, thềm từ hóa tại giá trị phân số chính xác $m = 1/3$ (tương ứng cấu hình spin Up-Up-Down - UUD trong mỗi tam giác cơ sở chứa đường chéo Shastry – Sutherland) tồn tại bền vững trong một dải từ trường rộng $\Delta h \approx 1.5 J'$. Khi độ thăng giáng $\Delta$ tăng lên trong tương tác NNN hoặc khi nhiệt độ vượt ngưỡng $\tau > 0.3$, thềm $m = 1/3$ bị thu hẹp và chuyển hóa thành đường cong từ hóa trơn liên tục, giải thích cơ chế dập tắt bước nhảy từ trong các mẫu thực nghiệm $TmB_4$ khi bị pha tạp khuyết tật mạng.

3. Phát hiện pha Siêu Tinh Thể (SS) kích hoạt bởi khuyết tật lỗ trống dưới thế ghim tuần hoàn: Đối với mô hình Bose – Hubbard lõi cứng ($V_{nn} = 6t$), khi chưa có thế ghim ($\varepsilon = 0$), hệ thống chỉ xuất hiện pha SS ở phía mật độ hạt cao ($\rho > 0.5$, gây bởi khuyết tật giữa nút), trong khi ở $\rho < 0.5$ hệ bị phân tách pha (phase separation). Khi áp đặt thế năng ghim tuần hoàn $\varepsilon = 2t \div 6t$, pha SS xuất hiện cân đối và ổn định ở cả hai vùng mật độ $\rho \in (0.42, 0.50)$ (khuyết tật lỗ trống) và $\rho \in (0.50, 0.58)$ (khuyết tật giữa nút).

4. Sự hình thành pha SS đa mật độ trong Boson lõi mềm ($U = 8t$): Trong hệ boson lõi mềm, thế ghim tuần hoàn $\varepsilon = 20t$ mở ra các vùng pha SS ổn định xung quanh các mật độ phân số $\rho = 1/2$, $\rho = 1$ và $\rho = 3/2$. Tại lân cận $\rho = 1$, các giá trị tham số trật tự đạt $S(\pi, 0) \approx 0.12$ và $\rho_s \approx 0.06$ tại $\mu \approx 3.8t$, xác lập bằng chứng mô phỏng đầu tiên về sự tồn tại của pha siêu tinh thể ở mật độ nguyên dưới tác động của thế ghim tuần hoàn.

Hệ vật liệu / Mô hình	Phương pháp	Tham số tới hạn	Trật tự quan sát	Hiện tượng đột phá
Màng Nano $Pr_{1-x}Ca_xMnO_3$	Giải tích Callen-Fourier MFA	$z=4, p=0.35, \Delta=1.35$	Sắt từ - Phân tách pha	Bước nhảy từ loại I, giảm $85%$ điện trở suất tại $T=3\text{ K}$
Mạng Shastry-Sutherland ($RB_4$)	Metropolis Monte Carlo	$J/J'=0.1, \tau=0.1, L=48$	Cấu hình UUD ($m=1/3$)	Thềm từ hóa phân số $1/3$ bền vững dưới thăng giáng tương tác
Boson lõi cứng ($V_{nn}=6t$)	Worm Algorithm QMC	$\varepsilon = 4t, \rho \in (0.42, 0.58)$	$\text{DLRO } S(\pi, 0) > 0$ $\text{ODLRO } \rho_s > 0$	Pha Siêu tinh thể kích hoạt bởi cả khuyết tật lỗ trống và giữa nút
Boson lõi mềm ($U=8t, \varepsilon=20t$)	Worm Algorithm QMC	$\mu \in (3.0t, 4.5t), \rho \approx 1$	Sóng mật độ điện tích (CDW) + Siêu chảy (SF)	Pha Siêu tinh thể xuất hiện đồng thời tại 3 mức mật độ $\rho=1/2, 1, 3/2$

    Mật độ Siêu chảy (ρ_s) / Hệ số Cấu trúc S(π,0)
     ▲
     │                      VÙNG PHA SIÊU TINH THỂ (SS)
 0.15│                     [ S(π,0) > 0  VÀ  ρ_s > 0 ]
     │                                ╭─────╮
 0.10│      S(π,0) [DLRO] ───────────/       \─────────── S(π,0)
     │                              /    │    \
 0.05│      ρ_s    [ODLRO] ────────/     │     \───────── ρ_s
     │                            /      │      \
 0.00└───────────────────────────┴───────┴───────┴────────────────► Mật độ hạt (ρ)
                               0.42     0.50    0.58
                              (Lỗ trống) (CDW) (Giữa nút)

Implications đa chiều

• Về mặt lý thuyết: Thống nhất mô hình chuyển pha từ loại I trong các hệ từ có phân tách pha kích thước nano; hoàn thiện bức tranh pha lượng tử của hệ boson tương quan mạnh khi có mặt thế ngoài không đồng nhất.
• Về mặt phương pháp luận: Thiết lập quy trình chuẩn mực kết hợp giải tích tích phân Fourier với mô phỏng Monte Carlo lượng tử Worm không dấu âm, có thể chuyển giao trực tiếp để nghiên cứu các hệ chất lỏng lượng tử và vật liệu từ nhiệt nano khác.
• Về mặt thực tiễn và công nghệ: Cung cấp thông số dẫn đường cho việc chế tạo các linh kiện spintronics từ trở khổng lồ nano (TMR sensors), bộ nhớ từ tính đa trạng thái (multi-state magnetic memories) dựa trên thềm từ hóa phân số, và mô phỏng lượng tử trên mạng quang học nguyên tử lạnh (cold-atom optical lattice simulators).

---

Limitations và Future Research

1. Bỏ qua hiệu ứng lượng tử của spin trong mô hình Ising Mangan: Luận án sử dụng phép gần đúng spin cổ điển $S = 1/2$ và spin hiệu dụng $S_{t2g} + 1/2$ theo trục z. Mặc dù phù hợp với điều kiện dị hướng từ đơn trục mạnh, mô hình chưa tính đến các thành phần thăng giáng lượng tử ngang ($S_x, S_y$) vốn có thể làm thay đổi cấu trúc vi mô tại vùng lân cận điểm chuyển pha bậc hai.
2. Kích thước mạng mô phỏng 2D hữu hạn: Các tính toán QMC bị giới hạn ở kích thước mạng $L = 20$ do chi phí tính toán tăng theo hàm lũy thừa thời gian thực, dù đã khắc phục bằng phép ngoại suy kích thước hữu hạn $1/L \to 0$.
3. Mô hình thế ghim tuần hoàn lý tưởng hóa: Thế ghim được khảo sát ở dạng tuần hoàn hoàn hảo $\mu_i = \pm \varepsilon$, chưa phản ánh toàn diện các dạng thế ngẫu nhiên ngắt đoạn (Anderson disorder potential) xuất hiện trong các tinh thể thực tế.

Chương trình nghiên cứu tương lai (Future Research Directions):

• Mở rộng mô hình Heisenberg lượng tử hoàn chỉnh ($S=1/2$ và $S=1$) trên mạng Shastry – Sutherland có xét đến tương tác Dzyaloshinskii – Moriya (DM).
• Khảo sát động học chuyển pha lượng tử phi cân bằng (non-equilibrium quantum phase transitions) sau các xung laser cực ngắn trên hệ perovskite manganite.
• Mô phỏng hệ Boson spin-orbit coupling (ghép cặp spin - quỹ đạo) trên mạng quang học dưới tác dụng của thế ghim ngẫu nhiên.
• Mở rộng thuật toán Worm để khảo sát hệ hỗn hợp Bose – Fermi trên mạng quang học 3 chiều.

---

Tác động và ảnh hưởng

• Tác động học thuật (Academic Impact): Cung cấp hệ giải pháp lý thuyết chuẩn xác cho các bài toán phân tách pha nano và chuyển pha lượng tử, dự kiến thu hút lượng trích dẫn cao trong các tạp chí chuyên ngành vật lý chất rắn (Physical Review B, Journal of Magnetism and Magnetic Materials).
• Định hướng công nghiệp bán dẫn & Spintronics: Tối ưu hóa cơ chế điều khiển từ trường ngoài lên hiệu ứng từ trở trong màng mỏng perovskite manganite, đặt nền tảng cho việc thiết kế đầu đọc đĩa cứng siêu nhạy và các cảm biến từ trường nano thế hệ mới.
• Phát triển công nghệ mô phỏng lượng tử: Cung cấp bản đồ pha chính xác hỗ trợ các phòng thí nghiệm quốc tế thiết lập thông số laser và từ trường để chế tạo thành công pha siêu tinh thể trên mạng quang học của các nguyên tử $Rb^{87}, Dy^{164}, Er^{166}$.

---

Đối tượng hưởng lợi

• Nghiên cứu sinh và Giới học thuật: Tiếp cận khung lý thuyết tích hợp Callen - Fourier và thuật toán Worm QMC tiên tiến để giải quyết các bài toán vật lý chất rắn phức tạp.
• Kỹ sư R&D Linh kiện Spintronics: Nắm bắt quy luật phụ thuộc của hệ số từ trở vào độ phân cực $P$ và độ từ hóa $m$, hỗ trợ thiết kế cảm biến TMR hoạt động ở dải nhiệt độ mở rộng.
• Nhà nghiên cứu Vật lý Nguyên tử Lạnh & Quang học Lượng tử: Sử dụng giản đồ pha $(U/t, \varepsilon/t, \mu/t)$ để điều khiển bẫy quang học, xác thực thực nghiệm các trạng thái lượng tử kỳ dị.

---

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Trường Trung bình trên Mô hình Ising Ngẫu nhiên thông qua việc kết hợp hệ thức Callen với biểu diễn tích phân liên tục Fourier của hàm Delta Dirac. Luận án đã chuyển đổi bài toán lấy trung bình thống kê cấu hình thăng giáng phức tạp thành một phương trình đại số tự phối hợp đóng với các hệ số tích phân $A_n(\tau, p, \Delta, z, h)$, giải quyết chính xác sự cạnh tranh giữa tương tác sắt từ và phản sắt từ mà không cần dùng đến các giả định gần đúng thô bạo trước đây.

2. Phương pháp nghiên cứu có điểm gì cải tiến so với các công trình quốc tế tiền nhiệm?

So với các nghiên cứu của Sengupta et al. (2005) hay Batrouni et al. (2006) vốn dùng thuật toán Stochastic Series Expansion (SSE) hoặc Determinantal QMC thường gặp khó khăn về kích thước mạng hoặc giới hạn phân tách pha, luận án sử dụng Thuật toán Worm trong không gian cấu hình tích phân đường liên tục. Thuật toán này cho phép lấy mẫu trực tiếp hàm Green hai điểm, triệt tiêu hoàn toàn hiện tượng nghẽn tới hạn và bẫy cấu hình, lần đầu tiên quét toàn diện giản đồ pha của hệ boson dưới thế ghim tuần hoàn $\varepsilon$ ở cả vùng khuyết tật lỗ trống.

       ĐẶC TÍNH SO SÁNH PHƯƠNG PHÁP NGHIÊN CỨU
       ┌─────────────────────────────────────────────────────────────┐
       │ Nghiên cứu Quốc tế Tiền nhiệm:                              │
       │ • SSE / Determinantal QMC (Sengupta 2005, Batrouni 2006)    │
       │ • Chỉ quan sát được pha SS ở ρ > 0.5 (Khuyết tật giữa nút) │
       │ • Bị giới hạn bởi hiện tượng phân tách pha ở ρ < 0.5       │
       └──────────────────────────────┬──────────────────────────────┘
                                      │  VƯỢT TRỘI BỞI LUẬN ÁN
                                      ▼
       ┌─────────────────────────────────────────────────────────────┐
       │ Phương pháp Luận án (Kim Oanh et al.):                      │
       │ • Continuous-time Worm QMC + Thế ghim tuần hoàn ε           │
       │ • Lấy mẫu trực tiếp hàm Green, triệt tiêu nghẽn tới hạn     │
       │ • Ổn định pha SS ở CẢ HAI PHÍA: ρ ∈ (0.42, 0.58)            │
       └─────────────────────────────────────────────────────────────┘

3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm/mô phỏng là gì?

Phát hiện bất ngờ nhất là sự xuất hiện của pha Siêu Tinh Thể ổn định ở mật độ nguyên $\rho = 1$ trong mô hình hạt boson tương tác mềm ($U = 8t$) dưới tác dụng của thế ghim lớn $\varepsilon = 20t$. Thông thường, tại mật độ nguyên $\rho = 1$, tương tác đẩy cục bộ $U$ sẽ khóa chặt các hạt vào trạng thái điện môi Mott (Mott Insulator) hoặc sóng mật độ điện tích (CDW) cách điện hoàn toàn ($\rho_s = 0$). Tuy nhiên, thế ghim định xứ tuần hoàn đã làm tái phân bố không gian mật độ hạt, kích hoạt các dòng siêu chảy xuyên hầm giữa các nút không bị ghim, làm phát sinh đồng thời $\rho_s > 0$ và $S(\pi, 0) > 0$.

4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả nghiên cứu (Replication Protocol)?

Hoàn toàn đầy đủ. Luận án mô tả chi tiết:

• Toàn bộ các hệ số giải tích $A_n$ của phương trình trường trung bình (Công thức 2.30).
• Sơ đồ thuật toán Metropolis chi tiết từng bước (Hình 3.3).
• Bộ 3 phép cập nhật của thuật toán Worm trong biểu diễn tích phân đường (Hình 4.8, 4.9, 4.10).
• Các bộ tham số không thứ nguyên chuẩn hóa: $z=4, p=0.35, \Delta=1.35, J/J'=0.1, V_{nn}=6t, U=8t, \varepsilon=4t \div 20t$, cho phép bất kỳ nhóm nghiên cứu tính toán nào tái lập chính xác các đồ thị đường cong từ hóa và giản đồ pha.

5. Chương trình nghiên cứu 10 năm được vạch ra như thế nào?

Lộ trình phát triển 10 năm tập trung vào 3 giai đoạn:

1. Năm 1 - 3: Tích hợp tương tác spin-quỹ đạo lượng tử và hiệu ứng trường tinh thể vào mô hình mạng trắc trở nano 3D; phát triển mã nguồn mở mô phỏng Worm QMC đa luồng trên GPU.
2. Năm 4 - 6: Thiết kế các cấu trúc siêu mạng nano perovskite nhân tạo ($LaMnO_3/SrMnO_3/PrMnO_3$) nhằm kiểm chứng thực nghiệm các điểm nhảy từ hóa và hiệu ứng từ trở được dự đoán bởi lý thuyết.
3. Năm 7 - 10: Hợp tác với các trung tâm máy tính lượng tử nguyên tử lạnh để hiện thực hóa các mạch logic lượng tử và bộ nhớ tô pô dựa trên các trạng thái pha siêu tinh thể và thềm từ hóa phân số.

---

Kết luận

1. Xác lập hệ thức giải tích tự phối hợp hoàn chỉnh cho mô hình Ising ngẫu nhiên $z=4$, giải thích định lượng cơ chế bước nhảy từ hóa loại I và sự sụt giảm $85%$ của từ trở khổng lồ trong perovskite mangan $Pr_{0.65}Ca_{0.35}MnO_3$ ở $T = 3\text{ K}$.
2. Làm sáng tỏ độ bền nhiệt động của thềm từ hóa phân số $m = 1/3$ trên mạng trắc trở hình học Shastry – Sutherland, chứng minh cấu hình Up-Up-Down (UUD) được bảo toàn trong dải tỷ số $J/J' \le 0.5$ và nhiệt độ $\tau \le 0.2$.
3. Phát hiện và chứng minh cơ chế kích hoạt pha Siêu Tinh Thể (SS) bởi khuyết tật lỗ trống dưới tác dụng của thế năng ghim tuần hoàn $\varepsilon$, giải quyết triệt để nghịch lý phân tách pha trong mô hình Bose – Hubbard 2D.
4. Mở rộng giản đồ pha của Boson lõi mềm, xác lập sự tồn tại của pha siêu tinh thể tại mật độ hạt nguyên $\rho = 1$ dưới tương tác tương quan mạnh $U = 8t$ và thế ghim $\varepsilon = 20t$.
5. Đóng góp phương pháp luận chuẩn mực thông qua việc kết hợp giải tích toán lý Callen - Fourier với mô phỏng Monte Carlo lượng tử Worm không dấu âm, tạo tiền đề vững chắc cho việc nghiên cứu các hệ lượng tử tương quan mạnh trong tương lai.

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

CHƯƠNG 1: TỔNG QUAN VỀ PHA VÀ CHUYỂN PHA 1. Pha và đặc trưng pha 1. Định nghĩa pha Hệ vật lý với các trạng thái khác nhau là một phần thế giới của vật chất mà chúng ta xem xét. Hệ vật lý trong đó các phần của hệ có thể tương tác trao đổi năng lượng với nhau hoặc trao đổi với môi trường xung quanh dưới dạng công hoặc nhiệt lượng được gọi là hệ nhiệt động.

Một hệ nhiệt động được mô tả đầy đủ thông qua tập hợp các trạng thái nhiệt động. Trạng thái của hệ nhiệt động được đặc trưng bởi các hàm trạng thái được gọi là các hàm nhiệt động lực học. Đó là hàm của các biến là các thông số đặc trưng. Các hàm nhiệt động lực học thường được sử dụng là nội năng U, năng lượng tự do Hemholtz F, enthalpy HE, năng lượng Gibbs hay thế nhiệt động G …Chúng là hàm của các tham số áp suất P, thể tích V và nhiệt độ tuyệt đối T của hệ.

Các hàm nhiệt động được biểu diễn theo các thông số áp suất P, thể tích V và nhiệt độ tuyệt đối T của hệ như sau: Năng lượng tự do F = U − TSe .2) Năng lượng tự do Gibbs G = U − TSe + PV = F + PV .3)  G  Entropy của hệ Se = −   .5)  T  P Đối với vật liệu có trật tự từ hoặc trật tự điện, các hàm nhiệt động còn phụ thuộc vào các tham số trật tự (độ phân cực từ, điện …) và trường ngoài (từ trường, điện trường …). Pha của hệ nhiệt động được đặc trưng bởi tính đồng nhất trong đó các tính chất vật lý và thành phần hóa học tại mọi điểm đều như nhau hoặc nếu có biến đổi thì chúng biến đổi liên tục từ điểm này sang điểm khác [1]. Tính đồng nhất của mỗi pha 6 liên quan đến tính đối xứng nội của hệ nhiệt động trong pha đó. Tính đối xứng được đặc trưng thông qua tham số trật tự.1: Ví dụ về một số pha của vật chất [45].

Các tham số trật tự đặc trưng cho các pha của hệ nhiệt động có thể là mật độ, chiết suất, độ từ hóa, độ phân cực điện … Một số ví dụ về các pha như: “pha rắn, pha lỏng của kim loại và hợp kim, pha sắt từ, thuận từ của các vật liệu từ, pha dẫn điện thường hoặc siêu dẫn của các chất siêu dẫn .1 biểu diễn một số pha khác nhau của vật chất. Chuyển pha và phân loại chuyển pha Hiện tượng chuyển pha xảy ra thường xuyên trong tự nhiên và được ứng dụng rộng rãi trong đời sống và sản xuất. Chuyển pha là quá trình biến đổi trạng thái của hệ nhiệt động từ mức độ đối xứng này sang mức độ đối xứng khác hình thành tính chất mới của hệ. Chuyển pha thường xảy ra dưới sự tác động của các tham số điều khiển như các đại lượng vật lý : nhiệt độ, áp suất, trường ngoài … Ví dụ một số loại chuyển pha như chuyển pha rắn – lỏng, lỏng – hơi [46] hay chuyển pha sắt từ - thuận từ [2] dưới tác động của nhiệt độ….

Tại điểm chuyển pha, trạng thái của hệ có thể thay đổi một cách liên tục hoặc thay đổi một cách đột ngột nhưng đối xứng tại điểm chuyển pha bao giờ cũng thay đổi dạng nhảy bậc. 7 Để nghiên cứu về quá trình chuyển pha, chúng ta thường bắt đầu bằng ứng dụng phương pháp nhiệt động lực học. Chúng ta sẽ mô phỏng hệ nhiệt động bằng một hệ lý tưởng tuân theo các nguyên lý nhiệt động lực học như nguyên lý bảo toàn năng lượng, nguyên lý tăng entropy của một hệ cô lập, nguyên lý thứ 3 của nhiệt động học (sự thay đổi entropy của một hệ ngưng tụ bất kỳ trong quá trình đẳng nhiệt thuận nghịch tiến tới 0 khi nhiệt độ tiến tới 0 Kenvin). Vào năm 1910, Kamerlingh Onnes ở trung tâm nghiên cứu vật chất nhiệt độ thấp của trường đại học Leiden, Hà Lan đã có khám phá thực nghiệm quan trọng cho thấy đạo hàm bậc hai năng lượng tự do cho hệ Helium lỏng bị gián đoạn [27], từ đó cung cấp dữ liệu trong việc phân loại đầu tiên về quá trình chuyển pha.

Tiếp sau đó là tập hợp nghiên cứu về Helium trong suốt những năm 1920 và đánh dấu là phát hiện về chuyển pha dị thường trong Helium lỏng của Keesom vào năm 1932 đã đặt ra câu hỏi có bao nhiêu loại chuyển pha và cách phân loại chuyển pha [35]. Năm 1933, Paul Ehrenfest đã đưa ra lý thuyết phân loại chuyển pha đầu tiên. Quá trình chuyển pha có thể chia thành hai loại là: chuyển pha loại I và chuyển pha loại II. ✓ Chuyển pha loại I: Trong chuyển pha loại I, trạng thái của hệ bị biến đổi đột ngột, đạo hàm bậc nhất của hàm thế nhiệt động bị gián đoạn tại điểm chuyển pha.

Chuyển pha loại I có sự sắp xếp lại mạng tinh thể ví dụ như thay đổi góc giữa các mặt tinh thể xảy ra trong khoảng nhiệt độ rất hẹp dẫn đến thay đổi đột ngột đối xứng của vật. Ở chuyển pha loại I, hai pha có thể cùng tồn tại ở nhiệt độ chuyển pha ví dụ nước đá và nước ở 00C hay nước và hơi nước ở 1000C. Trong chuyển pha loại I có gián đoạn về entropy tức là kèm theo ẩn nhiệt. Khi entropy của hệ tăng (thu ẩn nhiệt) thì tính trật tự của hệ giảm và ngược lại.

✓ Chuyển pha loại II: là quá trình chuyển pha mà trạng thái của hệ biến đổi liên tục nên tại điểm chuyển pha các pha không tồn tại đồng thời, hàm thế nhiệt động và đạo hàm bậc nhất của các thế nhiệt động cũng như các hàm trạng thái nhiệt động ví dụ như entropy liên tục nhưng đạo hàm bậc hai của thế nhiệt động thì bị gián đoạn Trong chuyển pha loại II không có ẩn nhiệt và luôn kèm theo sự biến đổi tính trật tự và tính đối xứng của hệ. Do trạng thái biến đổi liên tục nên tại điểm chuyển 8 pha loại II không phân biệt được các pha ở mức độ vĩ mô như chuyển pha loại I. Khi qua điểm chuyển pha loại II, hệ chỉ dịch chuyển qua điểm mà tính chất của các pha và sự phụ thuộc nhiệt độ của các tính chất đó trở nên khác nhau [2].2: Sự biến đổi của năng lượng tự do G, enthalpy HE, entropy Se, thể tích V và nhiệt dung Cp theo nhiệt độ trong trường hợp (a) không có chuyển pha, (b) chuyển pha loại I và (c) chuyển pha loại II [69]. Sự thay đổi của các tham số nhiệt động trong trường hợp không có chuyển pha, chuyển pha loại I và chuyển pha loại II được minh họa trên hình 1.

Cũng theo phân loại của Ehrenfest, tồn tại chuyển pha bậc cao gọi là chuyển pha đa tới hạn hay chuyển pha bậc n. Điểm chuyển pha tương ứng gọi là điểm đa tới hạn. Trong chuyển pha đa tới hạn, các hàm thế nhiệt động thay đổi liên tục khi qua điểm chuyển pha, một số đạo hàm bậc hai và bậc cao hơn của thế nhiệt động theo các thông số trạng thái bị triệt tiêu tại điểm chuyển pha. Tại điểm đa tới hạn có n pha giống nhau không phân biệt được.

Tham số trật tự Vào năm 1937, Landau đã đưa ra khái niệm thay đổi tính đối xứng tại điểm chuyển pha. Trong chuyển pha loại II các đại lượng vật lý vẫn biến đổi liên tục, không có ẩn nhiệt nhưng có sự thay đổi về tính trật tự và tính đối xứng, do đó để mô tả được sự chuyển pha, Landau đã đưa vào khái niệm tham số trật tự (order parameter). Tham số trật tự là đại lượng phản ánh trạng thái của các pha, ký hiệu là η. Ở nhiệt độ thấp hơn nhiệt độ chuyển pha, pha có tính đối xứng thấp (pha trật tự) tham số trật tự có giá trị khác 0 (η ≠0).

Còn ở nhiệt độ cao hơn nhiệt độ chuyển pha, pha có tính đối xứng cao (pha mất trật tự) tham số trật tự bị triệt tiêu (η =0) [1][2]. Tham số trật tự có thể là đại lượng vô hướng, đại lượng vectơ hoặc thậm chí là một tenxơ. Ví dụ tham số trật tự trong các hệ thống từ tính là độ từ hóa (đại lượng véctơ) , trong chất lỏng là mật độ pha lỏng trừ mật độ pha khí (đại lượng vô hướng). Tham số trật tự η thay đổi từ giá trị η = 0 (trạng thái mất trật tự tuyệt đối) đến η = 1 (trạng thái trật tự tuyệt đối).

Khi giảm nhiệt độ thấp hơn nhiệt độ chuyển pha, hệ thống bắt đầu thay đổi từ trạng thái hỗn độn sang trái thái trật tự một phần. Tiếp tục giảm nhiệt độ xuống thấp hơn nữa, chúng ta có thể đạt được trạng thái trật tự tuyệt đối. Trong chuyển pha loại I, tham số trật tự thay đổi nhảy bậc theo nhiệt độ và biến đổi từ từ trong chuyển pha loại II [2] (xem hình 1.3: Sự phụ thuộc của tham số trật tự vào nhiệt độ trong chuyển pha loại I và chuyển pha loại II [44]. 10 Tuy nhiên có một vài hệ thống, trật tự pha bị ẩn đi và không thể xác định được bản chất của tham số trật tự.

Khi đó chúng ta có thể định nghĩa hàm tương quan của tham số trật tự. Khác với quá trình chuyển pha, các hàm tương quan phân kì theo quy luật hàm số mũ đối với độ dài tương quan  và thời gian đặc trưng 𝓉. Khi khảo sát các đại lượng nhiệt động, hàm thế nhiệt động là hàm của áp suất P, nhiệt độ tuyệt đối T, lực suy rộng ℱ và cả biến tham số trật tự η. Áp suất P và nhiệt độ T có giá trị cho trước tùy ý, còn giá trị thực tế của tham số trật tự η phải xác định từ điều kiện cân bằng nhiệt động hay điều kiện cực tiểu năng lượng [2].

Hiện tượng tới hạn Hiện tượng tới hạn là hiện tượng liên hệ với những biểu hiện đặc biệt của vật chất khi ở gần hoặc ở tại điểm chuyển pha liên tục (còn được gọi là điểm tới hạn). Trong nhiệt động lực học, điểm tới hạn là điểm cuối của đường cong cân bằng pha. Điểm tới hạn lần đầu tiên được khám phá vào khoảng hơn 100 năm trước nhưng việc nghiên cứu những biểu hiện ở gần điểm tới hạn chỉ mới xuất hiện gần đây. Trong mỗi trường hợp, tính chất của điểm tới hạn xuất hiện do hệ thống tồn tại những bậc tự do nhất định có những dao động lớn bất thường trên thang bước sóng dài.