Giáo trình chuyên đề Vật lý Nano Phương pháp trường tự hợp Hartree Fock áp dụng cho hệ nhiều điện tử Biên tập bởi TS Nguyễn Hồng Quang www princexml com Prince Non commercial License This document was[.]
Giáo trình chuyên đề Vật lý Nano Phương pháp trường tự hợp Hartree Fock áp dụng cho hệ nhiều điện tử Biên tập bởi: TS Nguyễn Hồng Quang Ket-noi.com kho tai lieu mien phi Giáo trình chuyên đề Vật lý Nano Phương pháp trường tự hợp Hartree Fock áp dụng cho hệ nhiều điện tử Biên tập bởi: TS Nguyễn Hồng Quang Các tác giả: TS Nguyễn Hồng Quang Phiên trực tuyến: http://voer.edu.vn/c/ac2ee593 MỤC LỤC Phương pháp trường tự hợp Hartree - Fock áp dụng cho hệ nhiều điện tử Áp dụng cho hệ nhiều điện tử lỗ trống chấm lượng tử dạng đĩa với giam cầm Parabolic Phụ lục I Phụ lục II Tài liệu tham khảo Tham gia đóng góp 1/38 Ket-noi.com kho tai lieu mien phi Phương pháp trường tự hợp Hartree - Fock áp dụng cho hệ nhiều điện tử Mơ hình ( ) → Xét hệ ba chiều gồm N điện tử có khối lượng m đặt trường V r Hamiltonian hệ phụ thuộc vào tọa độ N hạt mà hạt có ba thành phần theo ba phương khác nên Hamiltonian hệ phụ thuộc vào 3N toạ độ Hamiltonian hệ viết dạng: ( ) ^ ^ → → H = H r 1, , r N Phương trình Schrodinger hệ có dạng ^ Hψ = Eψ Thực tế phương trình khơng phải phương trình mà hệ 3N phương trình vi phân, phương trình khơng thể giải giải tích xác, nên hệ phương trình khơng giải xác mà phải giải gần Một phương pháp gần thông dụng phương pháp Hartree - Fock Nội dung phương pháp chuyển việc nghiên cứu giải phương trình Schrodinger hệ nhiều điện tử (hệ phương trình nhiều biến) việc nghiên cứu phương trình Schrodinger đơn điện tử (phương trình biến) Phương trình Schrodinger hệ N điện tử trạng thái dừng có dạng ( ) ( ) ^ → → → → Hψ r 1, , r N = Eψ r 1, , r N , với Hamiltonian ^ H= N ∑ i=1 ^ Hi + N ∑ i≠j=1 e2 , ϵrij 2/38 ( ) ^ → 2∇ Hi = − ℏ 2mi + V r i toán tử Hamiltonian điện tử thứ i trường ( ) → V r Số hạng thứ hai Hamiltonian mô tả tương tác Coulomb tất điện tử, ϵ số điện môi, rij = | | → → r i − r j khoảng cách hạt i j Để đưa phương trình Schrodinger hệ N điện tử phương trình điện tử ta đưa vào khái niệm trường trung bình Hãy đơn cử lấy điện tử thứ i Điện tử tương tác với tất N − điện tử cịn lại, mơ tả điện tử cách xét chuyển động trường tạo tất điện tử lại Giả sử → thời điểm ta tạo vị trí điện tử thứ i r i trường giống ( ) trường tạo thành điện tử cịn lại Kí hiệu trường điện tử thứ i → → trường điện tử lại Ueff r i Ueff r i phải mô tả gần ( ) ( ) tác dụng trung bình tất điện tử lên điện tử Gần Hartree ( ) → Giả sử cách ta biết trường trung bình Ueff r i Khi tốn tử Hamiltonian hệ N điện tử viết lại dạng N ^ ^' H = ? Hi , i=1 ^' với Hi = ℏ 2m ∇i ( ) ( ) → → + V r i + Ueff r i toán tử Hamiltonian điện tử thứ i Toán tử Hamiltonian hệ viết thành tổng N Hamiltonian với Hamiltonian → thứ i phụ thuộc vào tọa độ r i hệ, hàm sóng hệ tìm dạng tích trực tiếp N hàm sóng ( ) ( )( ) ( ) → → → → → → ψ r 1, r 2, , r N = ψ r ψ r ψ r N , ( ) → ^' với ψ r i hàm riêng toán tử Hamiltonian Hi với trị riêng ϵni, ta có ( ) ( ) ^' → → Hiψni r i = ϵniψni r i , 3/38 Ket-noi.com kho tai lieu mien phi với lượng hệ N E= ∑ϵ ni i=1 → → Mật độ xác suất tìm thấy điện tử thứ vị trí r 1, điện tử thứ hai vị trí r → điện tử thứ N vị trí r N ( → → → |ψ r 1, r 2, , r N )| ( )| = |ψn1 → r 2 | ( ) | | ( )| → r ψn2 → ψnN r N | ( )| ψ → Mât độ xác suất tìm thấy điện tử thứ k vị trí r k →nk r k Mật độ điện tích | ( )| với k = 1, , N → → điện tử thứ k vị trí r k ek ψnk r k → Thế tương tác Coulomb điện tử thứ i vị trí r i với điện tử thứ k vị trí → r k Vik = với ∫ | ( )| d →r , → eiek ψnk r k k rik | | → → → d r k = dxkdykdzkrik = r i − r k ei = ek = − e Thế tương tác điện tử thứ i với tất in t k cũn li k lại (k i) Ueff ( ) ∑ → r i = k≠i Vik = ∑ k≠i ∫ | ( )| d →r e ψnk → r k rik k Đó biểu thức hiệu dụng phương pháp trường trung bình gần Hartree đưa năm 1928 4/38 Gần Hartree - Fock Trong phép gần Hartree chưa tính đến nguyên lí hệ hạt đồng Các điện tử có spin bán nguyên s = / nên chúng tuân theo thống kê Fermi Dirac chúng thỏa mãn nguyên lí loại trừ Pauli Trạng thái điện tử i đặc trưng tọa độ xi, yi, zi thành phần hình chiếu spin si lên phương OZ Đối với điện tử sz có trị riêng msℏ với ms = ± / Hàm sóng điện tử i hàm biến số tọa độ xi, yi, zi si kí hiệu biến số ξi(i = 1, , N) Để mô tả trạng thái điện tử có tính đến spin ta đưa vào hàm s sau () ( ) ℏ −ℏ σ =0 σ (2) ( ) = ℏ =1 σ1 −ℏ = 0, σ1 1 −2 −2 Khi ta có ∑σ * α (σ )σ β (σ ) = δαβ σ Nếu bỏ qua tương tác mômen từ điện tử với từ trường điện tử chuyển động theo quỹ đạo gây nên ta biểu diễn hàm sóng điện tử i dạng ( ) → ψk(ξi) = ψnk r i σα(σi), số k hàm ψk(ξi) kí hiệu trạng thái lượng tử (nk, a) Điều kiện trực giao chuẩn hóa hàm ψk(ξi) ∫ψ * k (ξi)ψk(ξi)dξi = ∫ ( ) ( )∑ σ (σ )σ (σ ) → → ψn*k r i ψnl r i α i β i σi = δkl ≡ δnknlδαβ Phương trình Schrodinger có dạng ^ Hψ(ξ1, , ξN) = Eψ(ξ1, , ξN) 5/38 Ket-noi.com kho tai lieu mien phi Để phù hợp với nguyên lý loại trừ Pauli hàm Ψ(ξ1, , ξN) phải hàm phản đối xứng có dạng định thức Slater ψ(ξ1, , ξN) = = N √ ! ∑ ( − 1) P [ψ N √ ! | ν ν (ξ1) ψkN(ξN)] k1 ν ψk1(ξ1) ψk1(ξN) ⋱ ψkN(ξ1) [ ψkN(ξN) | , ] ký hiệu Pν ψk1(ξ1) ψkN(ξN) hàm nhận từ hàm ψk1(ξ1) ψkN(ξN) cách hoán vị ν cặp biến số ξi ? ξk cho Khi hốn vị cặp số ξi ? ξk hay cặp trạng thái ki ? kj cho định thức đổi dấu Khi ξi = ξk hay ki = kk định thức (khơng tồn hàm sóng) điều thỏa mãn ngun lí loại trừ Pauli (không tồn hạt trạng thái lượng tử) ( ) → ^' ^' Thực tế Ueff r i Hi cịn chưa biết nên hàm Ψni(ξi) hàm riêng Hi ( ) → cịn chưa xác định Ta dùng ngun lí biến phân để xác định Ueff r i ( ) → Gọi ψ0 r E0 hàm sóng lượng hệ trạng thái hệ lượng ( ) ^ → tử với toán tử Hamiltonian H ψ0 r , E0 thỏa mãn phương trình Schrodinger ( ) ( ) ^ → → Hψ0 r = E0ψ0 r , ( ) → ψ r hàm sóng khơng phải trạng thái (ở trạng thái kích thích) Năng lượng trung bình hệ lượng tử trạng thái ψ ∫ ( ) ( ) ¯ → ^ → → E = ψ * r Hψ r d r , ¯ E0 lượng trạng thái (là nhỏ ) nên E = E0 nghĩa 6/38 ∫ ( ) ( ) → ^ → → r Hψ r d r ≥ E0 ψ* ( ) ( ) → → ¯ Ta thấy hàm ψ r gần với hàm riêng ψ0 r E gần ( ) → E0 nhiêu Ta chọn trước lớp hàm ψ r có dạng thích hợp ( ) → ¯ lớp hàm chọn hàm ψ r cho giá trị E nhỏ (gần E0 nhất) nghĩa ¯ lời giải gần toán E ứng với hàm ψ cho nhỏ nên ¯ δE = E − E0 Vậy nghiệm gần ψ0 phải thỏa mãn điều kiện → ∫ ( →r )Hψ^ ( →r )d →r = δE = δ ψ * nội dung ngun lí biến phân Năng lượng trung bình hệ N điện tử ∫ ¯ ^ E = ψ * (ξ1, , ξN)Hψ(ξ1, , ξN)dΓ với dΓ = dξ1dξ2 ξN Thay hàm sóng (10) vào ta có lượng trung bình hệ N điện tử ¯ E N = ∑ ∫ψ * k ^ (ξi)H0(ξi)ψk(ξi)dξi k=1 + − 2 N ∑ ∫ψ ' * k (ξi)ψl* (ξj)U(ξi, ξj)ψk(ξi)ψl(ξj)dξidξj * k (ξi)ψl* (ξj)U(ξi, ξj)ψk(ξj)ψl(ξi)dξidξj, k, l = N ∑ ∫ψ ' k, l = ( ) → thay ψk(ξi) = ψnk r i σα(σi) ý ∑σ σα(σi)σβ(σ) = δαβ ta có i 7/38 Ket-noi.com kho tai lieu mien phi ¯ E= ∑ ∫ψ ( →r )H^ ( →r )ψ ( →r )d →r N * nk i i nk i i k=1 + − N ∑ ∫ψ ' * nk k, l = ( →r )ψ ( →r )U( →r , →r )ψ ( →r )ψ ( →r )d →r d →r * nl i N ∑ ∫ψ ' k, l = ↑ ↑ * nk j i j nk i nl j i j ( →r )ψ ( →r )U( →r , →r )ψ ( →r )ψ ( →r )d →r d →r * nl i j i j nk j nl i i j Trong số hạng cuối ta lấy tổng ứng với cặp điện tử có spin định hướng song song chiều ( ↑ ↑ , ↓ ↓ ) ¯ ¯ Ta tính δE sau cho δE = ta có ∫ ( ) ( ) ( ) → → → → → → → → + ∑ ∫δψ ( r )ψ ( r )U( r , r )ψ ( r )ψ ( r )d r d r ¯ → ^ → → → δE = δψn*k r i H0 r i ψnk r i δ r i N ' * nk i * nl j i j nk i nl j i j l=1 N − ∑ ∫δψ ' l = 1↑ ↑ * nk ( →r )ψ ( →r )U( →r , →r )ψ ( →r )ψ ( →r )d →r d →r , i * nl j i j nk j nl i i j ¯ thừa số 1/2 hai tổng cuối E lấy biên phân theo δψn*k ta gặp hai lần lần theo tổng k lần theo tổng l Từ điều kiện chuẩn hóa hàm sóng ∫ ( ) ( ) → → → ψn*k r i ψnl r i d r i = δnk, nl, ta suy ∫ ( ) ( ) → → → δψn*k r i →nl r i d r i = 0.vớimọink, nl 8/38