chuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang họcchuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang học chuyển pha kim loại – điện môi trong một số hệ tương quan mạnh trên mạng quang học
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI ––––––––––––– TRẦN THỊ THU TRANG CHUYỂN PHA KIM LOẠI - ĐIỆN MÔI TRONG MỘT SỐ HỆ TƯƠNG QUAN MẠNH TRÊN MẠNG QUANG HỌC Chuyên ngành: Vật lý lý thuyết vật lý tốn Mã số: 9.44.01.03 TĨM TẮT LUẬN ÁN TIẾN SĨ VẬT LÝ Hà Nội – 2023 Cơng trình hồn thành tại: Khoa Vật lý, Trường Đại học Sư Phạm Hà Nội; Viện Vật lý, Viện hàn lâm khoa học công nghệ Việt Nam Người hướng dẫn khoa học: PGS.TS Hoàng Anh Tuấn PGS.TS Lê Đức Ánh Phản biện 1: PGS.TS Phan Văn Nhâm Trường Đại học Duy Tân Phản biện 2: PGS.TS Đỗ Vân Nam Trường Đại học Phenikaa Phản biện 3: PGS.TS Phạm Văn Hải Trường Đại học Sư phạm Hà Nội Luận án bảo vệ trước Hội đồng chấm luận án cấp Trường họp Trường Đại học Sư phạm Hà Nội vào hồi … … ngày … tháng… năm 2023 Có thể tìm hiểu luận án thư viện: - Thư viện Quốc Gia, Hà Nội Thư viện Trường Đại học Sư phạm Hà Nội MỞ ĐẦU Lí chọn đề tài Chuyển pha kim loại - điện môi (Metal - Insulator Transition: MIT) vấn đề nghiên cứu vật lý hệ cô đặc Rất nhiều nghiên cứu lý thuyết lẫn thực nghiệm thực đạt nhiều kết tiến nhằm hiểu rõ vật liệu với vùng lấp đầy phần chất điện mơi chất điện mơi trở thành kim loại thông số điều khiển thay đổi Trên phương diện lý thuyết, Mott đưa lý giải trạng thái điện môi hình thành tương quan điện tử điện tử, từ gọi điện mơi Mott Hai số mơ hình thường sử dụng để mô tả hệ điện tử tương quan mơ hình Hubbard thơng thường (the standard Hubbard model) mơ hình Falicov- Kimball Mơ hình Hubbard mơ hình đơn giản hạn chế giải xác trường hợp hệ chiều hệ có số chiều vơ Mơ hình Falicov – Kimball (FKM) xem mơ hình Hubbard giản lược hạt với hướng spin có tham số nhảy nút 0, tức chúng không chuyển động Tương tự mơ hình Hubbard, MIT hệ lấp đầy nửa xảy FKM U thay đổi Điều khác biệt hai mơ hình pha kim loại mơ hình Hubbard mơ tả chất lỏng Fermi pha kim loại FKM chất lỏng không Fermi Một kết hợp tự nhiên hai mô hình nói mơ hình Hubbard bất đối xứng (Asymmetric Hubbard Model: AHM), loại hạt với spin khác có tham số nhảy nút khác Đa phần cơng trình lý thuyết AHM tập trung thiết lập giản đồ pha hệ chiều cho trường hợp tương tác hút lẫn tương tác đẩy Các trạng thái tìm thấy đa dạng, bao gồm điện mơi Mott, sóng trật tự điện tích siêu dẫn Tuy nhiên, cần phải thấy rằng: vấn đề phân tích mơ hình kiểu Hubbard thiếu phương pháp giải tích đáng tin cậy cho trường hợp U lớn Các phương pháp tính số mơ Monte - Carlo hay chéo hóa xác áp dụng tốt cho hệ có kích thước nhỏ lại hao tốn nhiều thời lượng tài ngun tính tốn Các gần giải tích khơng địi hỏi nhiều thời lượng tính tốn phương pháp lại có miền áp dụng hạn chế riêng nhiều trường hợp kết thu từ gần bộc lộ chất vật lý mơ hình hệ gần toán học sử dụng Chính vậy, cần sử dụng nhiều phương pháp khác để nghiên cứu, tiếp cận làm rõ giản đồ pha AHM Ngoài AHM, chúng tơi cịn nghiên cứu MIT mơ hình tổng qt đồng thời có tham số nhảy nút khơng đối xứng tương tác phụ thuộc vào nút mạng Những mơ hình nghiên cứu luận án thiết lập mạng quang học Mục tiêu nghiên cứu Làm rõ điều kiện xuất MIT đặc điểm pha điện môi số hệ tương quan mạnh, cụ thể: • Mơ hình Hubbard với tương tác phụ thuộc vào nút (HMSDI), • Mơ hình Hubbard bất đối xứng (AHM), • Mơ hình Hubbard bất đối xứng với tương tác phụ thuộc vào nút (AHMSDI) Đối tượng khách thể nghiên cứu Khách thể nghiên cứu: số hệ tương quan mạnh mô tả mơ hình HMSDI, AHM, AHMSDI Đối tượng nghiên cứu: Điều kiện xảy MMIT hệ tương quan mạnh đặc điểm chuyển pha Phạm vi nghiên cứu Một số hệ tương quan mạnh mạng quang học trạng thái khơng từ tính lấp đầy nửa Phương pháp nghiên cứu Phương pháp tính giải tích: • Kết hợp hình thức luận hàm Green CPA dẫn phương trình cho phép nghiên cứu MIT HMSDI • Sử dụng DMFT phương pháp phương trình chuyển động hàm Green (EMM) để nghiên cứu MIT AHM • Sử dụng lý thuyết trường trung bình động nút (2S-DMFT) DMFT kết hợp với EMM để nghiên cứu MIT AHMSDI Phương pháp tính số: Sử dụng phần mềm tính tốn (Fortran) để tính số mật độ trạng thái, giản đồ pha, số cư trú đôi làm sở đánh giá điều kiện chuyển pha đặc điểm pha Mott Ý nghĩa khoa học thực tiễn Các kết nghiên cứu ảnh hưởng cân khối lượng phụ thuộc vào nút tương tác xen kẽ lên điều kiện chuyển pha kim loại - điện môi chất pha điện môi Mott hỗn hợp khí fermion hai thành phần hệ quang học với tương tác không đồng Trong trường hợp riêng hệ có cân khối lượng tương tác xen kẽ không đồng nhất, kết thu phù hợp tốt với kết nghiên cứu thu trước Các mơ hình lý thuyết nghiên cứu luận án thực hố mơ lượng tử thơng qua mạng quang học, kiểm tra so sánh chúng với kết mô lượng tử để hiểu rõ vai trò pha điện mơi Mott định hướng việc tìm kiếm chúng vật liệu tương quan tương ứng Những điểm luận án Luận án thiết lập giản đồ pha, điều kiện xảy chuyển pha làm rõ chất pha điện môi Mott thông qua số cư trú đôi hệ tương quan mạnh trạng thái lấp đầy nửa, không từ tính Sử dụng 2S-DMFT, chúng tơi nghiên cứu AHMSDI, trạng thái lấp đầy nửa, nhiệt độ T=0K Chúng thu giản đồ pha hệ mặt giải tích, vùng kim loại bị giảm cân khối lượng tăng lên | 𝑈𝐴 𝑈𝐵 | = (𝑡↑ + 𝑡↓ + √𝑡↑2 + 𝑡↓2 + 14𝑡↑ 𝑡↓ ) Chúng thuộc tính trạng thái hệ: cân khối lượng lớn, fermion nhẹ tái chuẩn hóa nhiều fermion nặng, chúng mạng với tương tác nút nhỏ cường độ tương tác tới hạn giảm cân khối lượng ngày tăng Kết đáng tin cậy tìm lại kết riêng trường hợp tới hạn, cho thấy 2S-DMFT phương pháp tốt để nghiên cứu MIT hệ tương quan mạnh Bố cục luận án Về bố cục, luận án chúng tơi ngồi phần mở đầu, phần kết luận tài liệu tham khảo cấu trúc thành chương với nội dung cụ thể sau: Chương 1: Tổng quan chuyển pha kim loại – điện môi Mott Chương 2: Chuyển pha kim loại – điện mơi mơ hình Hubbard với thơng số bổ sung Chương 3: Chuyển pha kim loại – điện mơi mơ hình Hubbard bất đối xứng với tương tác phụ thuộc vào nút Chương TỔNG QUAN VỀ CHUYỂN PHA KIM LOẠI - ĐIỆN MÔI MOTT 1.1 Điện mơi Mott Theo lí thuyết vùng lượng, số lấp đầy định tính chất dẫn điện vật rắn trạng thái lấp đầy nửa, hay số electron ô sở lẻ ta có vùng khơng lấp đầy, chắn vật rắn kim loại Tuy nhiên, điều không Hạn chế lý thuyết vùng dựa tranh hạt, bỏ qua tương tác đẩy electron tinh thể Như vậy, tương tác Coulomb electron tinh thể yếu lí thuyết vùng mơ tả tốt, song với hệ tương quan mạnh - mà tương tác electron bỏ qua tranh khơng cịn đáp ứng Do nghiên cứu sâu Mott MIT thúc đẩy tương tác electron – electron, nhiều loại điện mơi hình thành loại tương tác gọi "điện mơi Mott" 1.2 Mơ hình Hubbard chuyển pha kim loại - điện môi Mott Một cột mốc cho nghiên cứu lí thuyết MIT dựa mơ hình đơn giản cho hệ fermion - Mơ hình Hubbard - mơ hình thành cơng nghiên cứu vấn đề J.Hubbard đưa từ năm 1960 Mơ hình quan tâm đến electron dải đơn với Hamiltonian đơn giản 𝑈 + 𝐻 = ∑ 𝑡𝑖𝑗 (𝑎𝑖,𝜎 𝑎𝑗,𝜎 + ℎ 𝑐 ) + ∑ 𝑛𝑖,𝜎 𝑛𝑖,−𝜎 , (1.1) ,𝜎 đó: 𝑖,𝜎 + 𝑎𝑖,𝜎 (𝑎𝑖,𝜎 ) + 𝑎𝑖,𝜎 𝑎𝑖,𝜎 tốn tử sinh (hủy) electron có spin σ nút i; 𝑛𝑖,𝜎 = tốn tử số hạt có spin σ nút i; 𝑡𝑖𝑗 tích phân nhảy nút, đặc trưng cho tính linh động electron; U tương tác Coulomb nút, định tính định xứ electron Hai thơng số cốt lõi mơ hình Hubbard cường độ tương quan electron U/t số lấp đầy n MIT mơ hình Hubbard mở rộng gần nghiên cứu với nhiều phương pháp khác chưa có lý thuyết hoàn thiện 1.3 Mạng quang học Mạng quang học tinh thể ánh sáng tạo thành giao thoa chùm laser tạo hiệu dụng tuần hồn bẫy ngun tử trung hồ qua dịch chuyển Stark đủ lạnh Trong mạng quang, nguyên tử bị bẫy vị trí có cực tiểu giống mạng tinh thể thực Tinh thể thực phức tạp với nhiều tương tác cạnh tranh, trật tự dao động mạng nên việc tính tốn tượng quan sát từ thực nghiệm khó khăn Trong đó, mạng quang cung cấp mơ hình tinh thể lí tưởng, khơng sai hỏng, khơng dao động tương tác điều chỉnh cách tinh tế tạo mơ hình kiểm tra cho nghiên cứu lí thuyết vật lí tinh thể rắn Có nhiều phương pháp để kiểm sốt các thơng số hệ ngun tử siêu lạnh mạng tinh thể quang học: Hình học kích thước mạng; Phonon; Chui hầm (tham số nhảy nút - t); Tương tác nút – U; Các tương tác lân cận tương tác tầm xa; Mạng quang phụ thuộc spin; Tương tác đa hạt (plaquette); Thế tương tác; Nhiệt độ; Phụ thuộc thời gian Với thành tựu kỹ thuật làm lạnh laser, mạng tinh thể quang học nguyên tử trung hoà siêu lạnh thiết lập chúng mơ mơ hình Hubbard Mơ mơ hình Hubbard cho phép dễ dàng kiểm soát điều chỉnh thơng số mơ hình 1.4 Một số phương pháp nghiên cứu chuyển pha kim loại - điện môi Mott 1.4.1 Phương pháp CPA Bản chất CPA là: • Thay hệ ngẫu nhiên mơ tả hàm Green G hệ tuần hoàn hiệu dụng với hàm Green 𝐺𝑝 cho ⟨𝐺⟩ = 𝐺𝑝 • Hệ tuần hồn hiệu dụng xây dựng đảm bảo yêu cầu tự hợp, đại lượng vật lí đo hệ ngẫu nhiên phải có thăng giáng xung quanh giá trị tương ứng hệ hiệu dụng (⟨𝑇⟩ = 0) CPA phép gần đơn giản, áp dụng tốt trường hợp độ rộng miền hẹp hay mật độ tạp thấp mà thăng giáng khơng có nhiều ảnh hưởng đến tính chất vật lí hệ Áp dụng CPA giải gần mơ hình Hubbard đơn giải giá trị tới hạn 𝑈 = 𝑈𝐶 = cho thấy CPA phương pháp dễ sử dụng với tính tốn khơng q phức tạp kết thu có sai khác khơng lớn so với số phương pháp phổ biến khác (Sử dụng gần phân tán ngẫu nhiên (RDA) Noack tìm kết 𝑈𝐶 ≈ 1) 1.4.2 Phương pháp DMFT Trong DMFT, chất vật lí hệ nhiều hạt giữ lại mơ hình tạp, nút mạng (tạp) lấy với liên kết với phần cịn lại, phần mạng cịn lại với hốc trống thay bể hạt đóng vai trị trường trung bình động Sự khác biệt tranh CPA bể có tương tác thơng qua lai hố với hốc Trong DMFT có hai gần áp dụng: thứ nhất, nghiệm giả thiết bất biến tịnh tiến đồng nhất; thứ hai, lượng riêng giả thiết định xứ Những gần xác trường hợp vơ hạn chiều, song cho nhiều kết gần tốt trường hợp hữu hạn chiều, chí ba chiều thăng giáng khơng gian nhỏ Bằng việc thay mơ hình mạng nhiều hạt thành mơ hình tạp, số bậc tự giảm đáng kể, vấn đề đơn giản với DMFT Bên cạnh đó, mơ hình hạt nghiên cứu từ lâu tất phương pháp giải mô hình tạp Anderson sử dụng để giải phương trình DMFT 1.4.3 Phương pháp DMFT hai nút Lý thuyết trường trung bình động hai nút Potthoff đề xuất vào năm 2001, cách đơn giản hố mơ hình DMFT xét giới hạn nhỏ thực cách ánh xạ mơ hình mạng tương quan vào mơ hình tạp liên kết với bể không tương tác chứa nút, trường hợp bể đơn giản Tuy phép ánh xạ gần song mơ hình tạp nút lại giải xác 2S-DMFT cung cấp cách tính tốn đơn giản, nhanh chóng mà khơng làm xáo trộn tiếp cận trường trung bình tới mơ hình mạng tương quan Một số kết tính số cho chuyển pha Mott mơ hình Hubbard đơn tạp cho so sánh kết thu phương pháp 10 cho ta biết khả tồn hạt tải trạng thái Hình 2.1, hàm mật độ trạng thái mạng A (B) với UB = UA = UA = Trong hai trường hợp, đường mật độ trạng thái mạng B xuất đỉnh nhọn lân cận mức Fermi (ω = 0) nên trạng thái hệ kim loại Hình 2.1 Hàm mật độ trạng thái mạng A (B) với 𝑈𝐵 = giá trị khác 𝑈𝐴 Hình 2.2 hàm mật độ trạng thái mạng A(B) với 𝑈𝐵 = 2UA Khi 𝑈𝐴 = −0.6 mật độ trạng thái cho hai mạng mức Fermi khác cho biết trạng thái hệ kim loại Nhưng 𝑈𝐴 = −1.5, xuất khe quanh mức Fermi đồ thị DOS tương ứng hai mạng cho thấy xuất pha điện môi So sánh trường hợp xảy chuyển pha kim loại – điện môi ta thấy chuyển pha xảy |UA| lớn |UB| nhỏ ngược lại Hình 2.2 Hàm mật độ trạng thái mạng A (B) với UB = 2UA Hình 2.3: Giản đồ pha HMSDI nhiệt độ T = 0K lấp đầy nửa 11 Khảo sát kĩ thay đổi trạng thái điện môi – kim loại 𝑈𝐴 , 𝑈𝐵 thay đổi ta thu điểm chuyển pha Hình 2.3 Khi tương tác yếu vật liệu trạng thái kim loại bình thường Tập hợp điểm chuyển pha thỏa mãn 𝑈𝐴 𝑈𝐵 = ±𝑊 (2.11) Biểu thức (2.11) kết chính, quan trọng mơ tả điều kiện xảy MIT mơ hình Hubbard với tương tác xen kẽ sử dụng CPA mà chúng tơi tìm Hình 2.4: Hàm mật độ trạng thái mạng A (B) mức Fermi theo 𝑈𝐴 với 𝑈𝐵 /𝑈𝐴 = 1.0(𝑎), 2.0(𝑏), −0.6(𝑐), −2.0(𝑑) Hình 2.5: Số cư trú đôi Dα phụ thuộc vào 𝑈𝐴 với giá trị 𝑈𝐵 /𝑈𝐴 = 1.0(𝑎), 2.0(𝑏), −0.6(𝑐), −2.0(𝑑) Hình 2.4 biểu diễn mật độ trạng thái mức Fermi theo 𝑈𝐵 ứng với giá trị khác 𝑈𝐴 Ta thấy ngoại trừ 𝑈𝐴 = 𝑈𝐵 = vùng lân cận, 𝜌𝛼 (0) có giá trị lớn mạng có tương tác nút nhỏ Đồ thị DOS đường cong trơn chuyển pha kim loại – điện môi hệ nhiệt độ T=0K chuyển pha liên tục, kết thu phương pháp DMFT NGR Các trạng thái Mott làm rõ phân tích số cư trú đơi 𝐷𝛼 = 〈𝑛𝛼↑ 𝑛𝛼↓ 〉 (Hình 2.5) trình bày bảng sau: Vùng chuyển pha Mạng A Mạng B Mott i Chuyển pha Mott Chuyển pha Mott Mott ii Chuyển pha kết cặp Chuyển pha Mott Mott iii Chuyển pha kết cặp Chuyển pha kết cặp Mott iv Chuyển pha Mott Chuyển pha kết cặp 12 So sánh với kết thu từ lí thuyết DMFT, kết phù hợp tốt công bố Chúng tin CPA phương pháp đơn giản, cho phép tìm kết giải tích, giải thích tính chất vật lí cốt yếu nhiệt độ thấp cho kết định tính MIT hệ với tương tác xen kẽ 2.2 Mơ hình Hubbard bất đối xứng 2.2.1 Mơ hình hình thức luận Trong AHM, loại spin có tích phân nhảy nút hóa học khác Hamiltonian mơ hình + H = ∑ 𝑡𝑖𝑗 (𝑎𝑖,𝜎 𝑎𝑗,𝜎 + ℎ 𝑐 ) − ∑ 𝜇𝜎 𝑛𝑖,𝜎 + 𝑈 ∑ 𝑛𝑖,↑ 𝑛𝑖,↓ , (2.13) ,𝜎 𝑖,𝜎 Tham số bất đối xứng xác định 𝑟 = 𝑖 𝑡↓ 𝑡↑ với hai giới hạn 𝑟 = tương ứng với mơ hình Falikov-Kimball 𝑟 = tương ứng với mơ hình Hubbard thơng thường Chúng ta cần lưu ý AHM sử dụng để mô tả hỗn hợp fermion hai thành phần tải mạng tinh thể quang học Sử dụng DMFT kết hợp với phương pháp phương trình chuyển động, tách phương trình chuyển động hàm Green bậc thứ hai, chúng tơi tìm biểu thức giải tích cho điều kiện tới hạn 𝑈𝐶 = [ ( 𝑡↑2 + 𝑡↓2 + √{ 𝑡↑4 + 𝑡↓4 + 14 𝑡↑2 𝑡↓2 })] (2.44) Biểu thức cho UC Stasyuk cộng thu cách sử dụng kỹ thuật chiếu sở tốn tử fermion Hubbard Ở đây, chúng tơi tái hình thành cách đơn giản 2.2.2 Chuyển pha kim loại - điện mơi mơ hình Hubbard bất đối xứng Hình 2.7 cho thấy mật độ trạng thái cho loại spin ba giá trị tương tác Coulomb nút, trình chuyển pha Mott 13 hệ xảy U = 1.22D Hình 2.8 mơ tả mật độ trạng thái mức Fermi Ta thấy rằng, hàm mật độ trạng thái ρσ(0) cho hai loại spin đồng thời bị triệt tiêu vùng kết cặp mạnh Trong trường hợp 𝑟 = 0.4, cách sử dụng phép ngoại suy spline đơn giản từ miền liệu 𝑈 < 1.1𝐷, thu 𝑈𝐶 ≈ 1.22𝐷 Hình 2.9 Lặp lại với nhiều giá trị khác r, thu hàm tương tác tới hạn 𝑈𝐶 phụ thuộc vào r thể Hình 2.10 gần giống với kết giải tích phương trình (2.44) tồn phạm vi r Hình 2.7: DOS cho spin up spin down với r = 0.4 giá trị tương tác nút Hình 2.8: DOS mức Fermi phụ thuộc vào tương tác nút với giá trị r Hình 2.9: Tổng DOS mức Fermi ρ(0) = ρ↑(0)+ρ↓(0) phụ thuộc vào đẩy nút với r = 0.4 Hình 2.10: Thế tới hạn UC phụ thuộc vào r (kết tính số (đường chấm xanh) kết giải tích (2.44) (đường liền đỏ)) theo EMM so sánh với kết tính theo DMFT [85] L-DMFT [15] 14 Kết tính số cư trú đơi vẽ Hình 2.11 với giá trị khác U r Trong trường hợp không tương tác (U = 0), số cư trú đơi 0.25 nhanh chóng giảm U tăng Trong chương này, sử dụng EMM cơng cụ giải tốn tạp cho DMFT để khảo sát MIT AHM lấp đầy nửa Kỹ thuật thực trực tiếp Hình 2.11: Số cư trú đơi < 𝑛↑ 𝑛↓ > phụ thuộc vào tương tác nút theo EMM trục tần số thực, cho so sánh với DMFT [85] hiệu mặt tính tốn Các kết chúng tơi cơng bố, so sánh cho thấy phù hợp tốt với kết thu cách sử dụng phương pháp chéo hóa xác kỹ thuật Monte Carlo lượng tử Kết cho thấy, tiếp cận phương trình chuyển động phương pháp giải tạp đơn giản đáng tin cậy để nghiên cứu MIT AHM Chương CHUYỂN PHA KIM LOẠI – ĐIỆN MÔI TRONG MƠ HÌNH HUBBARD BẤT ĐỐI XỨNG VỚI TƯƠNG TÁC PHỤ THUỘC VÀO NÚT 3.1 Mơ hình Chúng tơi nghiên cứu hệ fermion hai thành phần không cân khối lượng mạng tinh thể quang học với tương tác điều biến theo không gian, mơ tả mơ hình Hubbard bất đối xứng với tương tác phụ thuộc vào nút sau mạng tinh thể hai cực tạo hai mạng lồng A, B xếp cho quanh nút A nút B ngược lại Chúng sử dụng DMFT, Hamiltonian ánh xạ vào mơ hình đơn tạp tự hợp, cho bởi: 15 + + ∗ 𝐻𝛼𝑖𝑚𝑝 = ∑ 𝜀𝛼𝑘𝜎 𝑐𝑘𝜎 𝑐𝑘𝜎 + ∑( 𝑉𝛼𝑘𝜎 𝑐𝑘𝜎 𝑑𝜎 + 𝑉𝛼𝑘𝜎 𝑑𝜎+ 𝑐𝑘𝜎 ) − ∑ 𝜇𝜎 𝑛𝑑𝜎 𝑘𝜎 𝑘𝜎 𝜎 + 𝑈𝛼 [𝑛𝑑↑ 𝑛𝑑↓ + (𝑛𝑑↑ + 𝑛𝑑↓ )], (3.2) Tính tốn ta được: 𝑈𝛼 − Δ𝛼𝜎 (𝜔) (3.7) Ở đây, 𝐺0𝛼𝜎 (𝜔) hàm Green khơng tương tác mơ hình tạp hiệu dụng mạng 𝛼 𝐺0−1𝛼𝜎 (𝜔) = 𝜔 + 𝜇𝜎 + 3.2 Tiếp cận EMM Để tính tốn hàm Green mơ hình đơn tạp Anderson, chúng tơi sử dụng EMM, tách phương trình chuyển động mơ hình đơn tạp Anderson (3.2) thành bậc thứ hai, thu giá trị gần sau cho hàm Green tạp chất 𝑈𝛼 Δ𝛼𝜎 (𝜔) −1 (𝜔) 𝐺0𝛼𝜎 + −1 𝐺0𝛼𝜎 (𝜔) − 𝑈𝛼 − Δ𝛼𝜎 (𝜔) 1 + 𝑈 Δ (𝜔) −1 𝐺0𝛼𝜎 (𝜔) − 𝑈𝛼 + −1 𝛼 𝛼𝜎 𝐺0𝛼𝜎 (𝜔) − Δ𝛼𝜎 (𝜔) 𝐺𝛼𝜎 (𝜔) = (3.26) Để nghiên cứu MIT hệ có tương tác xen kẽ, chúng tơi tính tốn DOS cho mạng 𝜌𝛼𝜎 (𝜔) = − π ℑGασ (ω), DOS mức Fermi ρασ(0) số cư trú đôi 𝐷𝛼 = 〈 𝑛𝛼↑ 𝑛𝛼↓ 〉 Sau đó, chúng tơi xây dựng giản đồ pha cho pha đồng 𝑇 = 0𝐾 Chúng vẽ biểu đồ mật độ trạng thái phụ thuộc spin cho mạng với 𝑈 = 1.5𝐷 𝛾 = 0.8 số giá trị khác 𝑟 (Hình 3.1) tương ứng với trạng thái kim loại, điện môi MIT DOS mức Fermi thuộc tính dẫn hệ thống tính tốn Hình 3.2 Ta thấy hai 𝜌𝛼𝜎 (0) đồng thời triệt tiêu điểm chuyển pha Mott Để mô tả ảnh hưởng tương tác điều biến không gian đến ổn định trạng thái kim loại bình thường Khi U cố định, tham số điều biến không gian γ nhỏ hệ cân khối lượng dễ chuyển từ pha kim loại sang pha điện mơi Mott 16 Hình 3.1 Mật độ trạng thái phụ thuộc spin mạng U = 1.5D, tham số điều biến không gian γ = 0.8 giá trị khác tham số bất đối xứng r Hình 3.3: Mật độ trạng thái phụ thuộc spin mạng U = 1.5D, r = 0.8 giá trị khác tham số điều biến không gian γ Hình 3.2 DOS mức Fermi phụ thuộc vào tham số bất đối xứng U = 1.5D với nhiều giá trị khác tham số điều biến khơng gian γ Hình 3.4: DOS mức Fermi phụ thuộc vào tham số điều biến không gian U = 1.5D với nhiều giá trị khác tham số bất đối xứng r Trong Hình 3.3, vẽ biểu đồ DOS phụ thuộc spin trạng thái cho mạng với 𝑈 = 1.5𝐷 𝑟 = 0.8 số giá trị khác γ Để tìm giá trị tới hạn tham số điều biến không gian 𝛾𝐶 , DOS mức Fermi hàm tham số điều biến không gian 𝛾 cho 𝑈 = 1.5𝐷 với số giá trị khác 𝑟 (Hình 3.4) Ta thấy 17 hai 𝜌𝛼𝜎 (0) đồng thời triệt tiêu điểm chuyển pha Mott Khi U cố định, tham số cân khối lượng 𝑟 nhỏ, tính bất đối xứng tương tác phải lớn để thúc đẩy hệ chuyển từ pha kim loại sang pha Mott Hình 3.5: Số cư trú đôi < n↑n↓ > phụ thuộc vào U với γ = 0.5 giá trị cố định khác r Hình 3.7: Thế tương tác tới hạn UC mơ hình lấp đầy nửa phụ thuộc vào tham số bất đối xứng Hình 3.6: Tham số điều biến không gian tới hạn mô hình lấp đầy nửa phụ thuộc vào tham số bất đối xứng Hình 3.8: Giản đồ pha trạng thái mơ hình lấp đầy nửa phụ thuộc vào tương tác mạng Các kết tính số cư trú đơi vẽ Hình 3.5 cho 𝛾 = 0.5 với giá trị khác 𝑟 Như biết, trường hợp không tương tác (𝑈 = 0), số cư trú đôi 0.25 nhanh chóng giảm 𝑈 tăng Trạng thái kim loại đặc trưng giảm tuyến tính 18 số cư trú đơi gia tăng tương tác 𝑈, trong vùng điện môi, giá trị lớn tương tác, số cư trú đôi nhỏ phụ thuộc yếu vào 𝑈 Ở giá trị nhỏ 𝑟, số cư trú đôi giảm nhanh giá trị tương tác tới hạn bị giảm xuống Sự phụ thuộc tham số điều biến không gian tới hạn 𝛾𝐶 (𝑟, 𝑈) cho mơ hình lấp đầy nửa dạng hàm tham số cân khối lượng 𝑟 giá trị khác tương tác 𝑈 = 1.0𝐷, 1.5𝐷 2.0𝐷 (Hình 3.6) Thế 𝑈 tăng vùng kim loại giảm tương quan mạnh Trong Hình 3.7, tương tác tới hạn 𝑈𝐶 cho mơ hình lấp đầy nửa dạng hàm tham số cân khối lượng r giá trị khác tham số điều biến khơng gian γ = 1.0,0.8 0.6 Dễ dàng tìm lại kết trường hợp giới hạn với tương tác đồng nhất, 𝑈𝐶 = ( 𝑡↑2 + 𝑡↓2 + √𝑡↑4 + 𝑡↓4 + 14 𝑡↑2 𝑡↓2 ) Các kết tóm tắt thành giản đồ pha Hình 3.8 Đối với trường hợp cân khối lượng (𝑟 = 1), kết phù hợp tốt với kết thu từ DMFT với phương pháp NGR Biểu thức cho giá trị tới hạn thiết lập |𝑈𝐴 𝑈𝐵 | = ( 𝑡↑2 + 𝑡↓2 + √𝑡↑4 + 𝑡↓4 + 14 𝑡↑2 𝑡↓2 ) (3.17) Biểu thức (3.17) kết chúng tơi Trong trường hợp mơ hình Hubbard khơng đối xứng thơng thường UA = UB, tái tạo kết biết trước 3.3 Tiếp cận 2S-DMFT Sử dụng 2S-DMFT, vận dụng vào mơ hình AHM thu mơ hình Anderson hai nút sau: 𝐻𝛼2−𝑠𝑖𝑡𝑒 = ∑ 𝜀𝑐𝛼𝜎 𝑐𝜎+ 𝑐𝜎 + ∑ 𝑉𝛼𝜎 (𝑐𝜎+ 𝑑𝜎 + 𝑑𝜎+ 𝑐𝜎 ) + ∑ 𝜀𝑑𝛼𝜎 𝑑𝜎+ 𝑑𝜎 𝜎 𝜎 + 𝑈𝛼 𝑛𝑑↑ 𝑛𝑑↓ , 𝜎 (3.20) 19 Chúng tơi thực tính tốn đến số hạng bậc hai 𝑉𝜎 Các lượng riêng hàm sóng riêng tương ứng điện tử dẫn đến bốn phương trình tuyến tính cho 𝑉𝛼𝜎 với (𝛼 = 𝐴, 𝐵; 𝜎 =↑, ↓): 𝑉𝛼𝜎 = 2𝑡𝜎 ( 2𝑉𝛼𝜎 + 𝑉𝛼 𝜎 ) |𝑈𝛼 | (3.30) từ thu biểu thức sau cho tương tác tới hạn: |𝑈𝐴 𝑈𝐵 | = (𝑡↑ + 𝑡↓ + √{𝑡↑2 + 𝑡↓2 + 14𝑡↑ 𝑡↓ }) (3.32) Phương trình (3.32) kết giải tích thu chương Trong Hình 3.9, chúng tơi trình bày giản đồ pha dạng hàm tương tác UA UB giá trị khác tham số bất đối xứng r nhiệt độ không Đối với hệ thống cân khối lượng (r = 1), kết thu theo 2S-DMFT phương trình (3.32) phù hợp với kết thu từ DMFT với phương pháp nhóm tái chuẩn hóa số (NRG) Khi cân khối lượng tăng lên (r giảm), vùng kim loại bị giảm, D t↑ cố định, khác biệt khối lượng trần lớn hệ dễ định xứ Hình 3.9: Giản đồ pha trạng thái cho trường hợp hệ lấp đầy nửa hàm số tương tác UA/D UB/D với giá trị khác r Hình 3.10: Mật độ trạng thái T = mạng A(B) với r = 0.4,UB = −2UA giá trị khác UA/D 20 Để xác nhận kết thu từ phương trình (3.32), chúng tơi tính tốn mật độ trạng thái (DOS) cho mạng 𝜌𝛼 (𝜔) = − Σσ ℑGασ (ω) DOS mức Fermi 𝜌𝛼𝜎 (0) Hình 3.10 cho thấy DOS π mạng với r = 0.4,UB = −2UA giá trị khác UA Hình 3.11 mơ tả DOS mức Fermi cho thành phần mạng thông số mô hình giống Chúng tơi kết luận chuyển pha kim loại - điện môi thực xảy hệ với tương tác phụ thuộc vào nút Ngoài ra, vùng kim loại phía ranh giới pha xây dựng Hình 3.9 chất lỏng Fermi ngoại trừ giới hạn Falicov-Kimball r = Trong Hình 3.12, vẽ đồ thị trọng số chuẩn hạt cho mơ hình Qua trường hợp, chúng tơi nhận thấy fermion nhẹ tái chuẩn hóa hạt nặng, tái chuẩn hóa mạnh xuất mạng với tương tác nút lớn Cả cân khối lượng điều biến không gian tương tác ảnh hưởng đến ổn định trạng thái kim loại thông thường, ảnh hưởng chúng khác Hình 3.11: Mật độ trạng thái mức Fermi phụ thuộc vào 𝑈𝐴 Hình 3.12: Trọng số chuẩn hạt mức Fermi phụ thuộc vào 𝑈𝐴 trường hợp 21 Kết tính số cư trú đơi cho biết (Hình 3.13), tương tác đẩy tương tác hút xen kẽ hệ, trạng thái điện môi Mott ổn định vùng liên kết mạnh, nơi trạng thái chiếm đơn lẻ tồn mạng với tương tác đẩy, trạng thái trống Hình 3.13: Số cư trú đôi dA(dB) cho chiếm đôi xuất mạng A(B) hệ có UB = −2UA mạng phụ thuộc vào UA trường hợp khác r = 1.0,0.4,0.05 Trong chương này, nghiên cứu hệ fermion hai thành phần không cân khối lượng mạng tinh thể quang học với tương tác phụ thuộc vào nút cách sử dụng DMFT kết hợp với EMM 2S-DMFT Với cách thứ nhất, tính số phụ thuộc spin mật độ trạng thái mức Fermi số cư trú đôi giúp kiểm chứng thơng qua thực nghiệm đặc tính vật lý hệ Với cách tiếp cận thứ hai, thuộc tính trạng thái hệ fermion thảo luận từ đặc điểm phụ thuộc spin trọng số chuẩn hạt mức Fermi số cư trú đôi cho mạng hàm cường độ tương tác nút giá trị khác cân khối lượng Kết chúng tơi tìm được công bố, cho thấy phù hợp tốt với kết nghiên cứu thí nghiệm Fermion cực lạnh xác minh thông qua quan sát thực nghiệm bao trùm kết biết trường hợp đặc biệt giới hạn cân khối lượng (r = 1) trường hợp tương tác đồng (UA = UB) 22 KẾT LUẬN Luận án nghiên cứu MMIT mơ hình Hubbard với số thơng số thay đổi ba phương pháp thường dùng hệ tương quan mạnh: CPA – phương pháp dễ sử dụng với tính tốn khơng q phức tạp thu kết tương đối phù hợp, đặc biệt mặt định tính; DMFT 2S-DMFT mang lại kết khả quan cho chuyển pha Mott pha chất lỏng Fermi mơ hình Hubbard vùng với cố gắng tính tốn Sử dụng CPA, đặc điểm chuyển pha Mott mơ hình Hubbard lấp đầy nửa với tương tác xen kẽ làm rõ: thu giản đồ pha pha kim loại điện môi nhiệt độ T = 0K, chất trạng thái điện môi Mott thông qua số cư trú đôi, xác định MIT hệ chuyển pha liên tục Đối sánh với DMFT, nghiên cứu cho kết phù hợp, đặc biệt mặt định tính, cịn số hạn chế CPA cách tiếp cận đơn giản thích hợp sử dụng để nghiên cứu giai đoạn khởi đầu Khảo sát MMIT AHM lấp đầy nửa, sử dụng DMFT kết hợp với EMM để giải toán tạp cho DMFT, tìm hàm phụ thuộc tương tác tới hạn hệ tăng theo mức độ bất đối xứng Chúng tơi tính tốn mật độ trạng thái mức Fermi số cư trú đôi, làm rõ đặc điểm chuyển pha Các kết chúng tơi công bố, so sánh cho thấy phù hợp tốt với nghiên cứu sử dụng ED kỹ thuật Monte Carlo lượng tử EMM cách giải tạp đơn giản, dễ khai triển Phương pháp có ưu điểm lớn thực tính tốn đơn giản, đưa biểu thức giải tích cách rõ ràng cho điều kiện chuyển pha, kết trường hợp riêng tương đương gần Hubbard III không xác định điểm tới hạn qua triệt tiêu số cư trú đôi Cuối cùng, nghiên cứu hệ bất đối xứng với tương tác phụ thuộc vào nút theo hai cách tiếp cận sử dụng 2S-DMFT sử dụng DMFT kết hợp với EMM Trong trạng thái lấp đầy nửa, nhiệt độ T = 0K, thu giản đồ pha hệ mặt giải tích, vùng kim loại bị giảm cân khối lượng tăng lên Phân tích đặc điểm trọng số chuẩn hạt mức Fermi, mật độ trạng thái mức Fermi 23 số cư trú đôi cho mạng phụ thuộc vào tương tác nút tương ứng với số giá trị xác định tham số cân khối lượng cho thấy thuộc tính trạng thái hệ: Khi cân khối lượng lớn, fermion nhẹ tái chuẩn hóa nhiều fermion nặng, chúng mạng với tương tác nút nhỏ cường độ tương tác tới hạn giảm cân khối lượng ngày tăng; Trong vùng liên kết mạnh hệ fermion hai thành phần với tương tác đẩy tương tác hút xen kẽ, chuyển pha Mott chuyển pha kết cặp xảy đồng thời mạng tương ứng chuyển pha nhiệt độ T = 0K liên tục Biểu thức giải tích cho điều kiện chuyển pha tương tác tới hạn biểu diễn hàm tham số cịn lại mơ hình, theo cách tiếp cận sử dụng EMM: |𝑈𝐴 𝑈𝐵 | = ( 𝑡↑2 + 𝑡↓2 + √𝑡↑4 + 𝑡↓4 + 14 𝑡↑2 𝑡↓2 ) theo cách tiếp cận 2S-DMFT: | 𝑈𝐴 𝑈𝐵 | = (𝑡↑ + 𝑡↓ + √𝑡↑2 + 𝑡↓2 + 14𝑡↑ 𝑡↓ ) Các phương pháp gần khác áp dụng nghiên cứu toán MMIT hệ tương quan mạnh với số biến thể mơ hình Hubbard cho kết khác nhau, phù hợp tốn cụ thể Các kỹ thuật tính tốn mà chúng tơi sử dụng khơng q phức tạp hay yêu cầu tài nguyên máy tính đắt đỏ cho kết tốt xét phương diện, so sánh tương quan với nghiên cứu trước phương pháp tính tốn phức tạp cho kết tốt Quan trọng nhất, lựa chọn 2S-DMFT phương pháp tốt, cho kết đáng tin cậy nghiên cứu chuyển pha Mott AHMSDI Kết nghiên cứu mơ hình Hubbard vừa có tương tác phụ thuộc nút lại vừa có tham số nhảy nút khác mới, bao trùm trường hợp riêng có tương tác phụ thuộc nút có tham số nhảy nút khác Trong luận án này, nghiên cứu MMIT giới hạn hệ với trạng thái khơng từ tính nhiệt độ T = 0K Các trạng thái trật tự từ tính, siêu lỏng, trật tự điện tích hay spin chưa xem xét chúng tơi tiếp tục nghiên cứu tương lai DANH MỤC CƠNG TRÌNH KHOA HỌC THUỘC LUẬN ÁN ĐÃ CƠNG BỐ Các cơng bố tạp chí chuyên ngành thuộc danh mục ISI: [1] Le Duc Anh, Tran Thi Thu Trang, Hoang Anh Tuan (2013), “Coherent potential approximation study of the Mott transition in optical lattice system with site-dependent interactions”, The European Physical Journal B, 86(12), 503 [2] Hoang Anh Tuan, Tran Thi Thu Trang, Le Duc Anh (2016), “Mott transition in the asymmetric hubbard model at half-filling: Equation of motion approach”, Journal of the Korean Physical Society,68(2), pp.238-242 [3] Hoang Anh Tuan, Nguyen Thi Hai Yen, Tran Thi Thu Trang, Le Duc Anh (2016), "Two-component Fermions in Optical Lattice with Spatially Alternating Interactions”, Journal of the Physical Society of Japan,85(10), 104702 [4] Le Duc Anh, Tran Thi Thu Trang, Hoang Anh Tuan, Nguyen Toan Thang, Tran Minh Tien (2018), “Mass-imbalanced Hubbard model in optical lattice with site-dependent interactions”, Physica B: Condensed Matter, 532, pp.204-209