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ử

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 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 bản trực tuyến: http://voer.edu.vn/c/ac2ee593 MỤC LỤC 1. Phương pháp trường tự hợp Hartree - Fock áp dụng cho hệ nhiều điện tử 2. Áp dụng cho hệ nhiều điện tử và lỗ trống trong chấm lượng tử dạng đĩa với thế giam cầm Parabolic 3. Phụ lục I 4. Phụ lục II 5. Tài liệu tham khảo Tham gia đóng góp 1/38 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 một hệ ba chiều gồm N điện tử có khối lượng m đặt trong một trường V ( → r ) nào đó. Hamiltonian của hệ phụ thuộc vào tọa độ của N hạt mà mỗi hạt có ba thành phần theo ba phương khác nhau nên Hamiltonian của hệ phụ thuộc vào 3N toạ độ. Hamiltonian của hệ có thể viết dưới dạng: ^ H = ^ H ( → r 1 , , → r N ) . Phương trình Schrodinger của hệ có dạng ^ Hψ = Eψ. Thực tế thì phương trình trên không phải là một phương trình mà là một hệ 3N phương trình vi phân, mỗi phương trình không thể giải được giải tích chính xác, nên hệ phương trình trên cũng không giải chính xác được mà phải giải gần đúng. Một trong các phương pháp gần đúng thông dụng là phương pháp Hartree - Fock. Nội dung của phương pháp này là chuyển việc nghiên cứu giải phương trình Schrodinger của hệ nhiều điện tử (hệ phương trình nhiều biến) về việc nghiên cứu phương trình Schrodinger đơn điện tử (phương trình một biến). Phương trình Schrodinger của 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 ^ H i + 1 2 N ∑ i ≠ j = 1 e 2 ϵr ij , 2/38 trong đó ^ H i = − ℏ 2 ∇ i 2 2m + V ( → r i ) là toán tử Hamiltonian của điện tử thứ i trong trường V ( → r ) . Số hạng thứ hai của Hamiltonian mô tả tương tác Coulomb giữa tất cả các điện tử, ϵ là hằng số điện môi, r ij = | → r i − → r j | là khoảng cách giữa 2 hạt i và j. Để đưa phương trình Schrodinger của hệ N điện tử về phương trình của một điện tử ta đưa vào khái niệm trường trung bình. Hãy đơn cử lấy một điện tử thứ i nào đó. Điện tử này tương tác với tất cả N − 1 điện tử còn lại, và do đó có thể mô tả điện tử đó bằng cách xét chuyển động của nó ở trong trường được tạo ra bởi tất cả các điện tử còn lại. Giả sử tại mỗi thời điểm ta có thể tạo ra được ở vị trí của điện tử thứ i ( → r i ) một trường giống như trường được tạo thành bởi các điện tử còn lại. Kí hiệu trường thế của điện tử thứ i trong trường của các điện tử còn lại là U eff ( → r i ) . U eff ( → r i ) này sẽ phải mô tả gần đúng nhất tác dụng trung bình của tất cả các điện tử lên một điện tử nào đó. Gần đúng Hartree Giả sử bằng cách nào đó ta đã biết được trường thế trung bình này U eff ( → r i ) . Khi đó toán tử Hamiltonian của hệ N điện tử được viết lại dưới dạng ^ H = N ? i = 1 ^ H i ' , với ^ H i ' = ℏ 2 2m ∇ i 2 + V ( → r i ) + U eff ( → r i ) là toán tử Hamiltonian của một điện tử thứ i. Toán tử Hamiltonian của hệ được viết thành tổng của N Hamiltonian với Hamiltonian thứ i chỉ phụ thuộc vào tọa độ → r i của hệ, do đó hàm sóng của hệ có thể tìm được dưới dạng tích trực tiếp của N hàm sóng ψ ( → r 1 , → r 2 , , → r N ) = ψ ( → r 1 ) ψ ( → r 2 ) ψ ( → r N ) , với ψ ( → r i ) là hàm riêng của toán tử Hamiltonian ^ H i ' với trị riêng ϵ n i , ta có ^ H i ' ψ n i ( → r i ) = ϵ n i ψ n i ( → r i ) , 3/38 với năng lượng của hệ E = N ∑ i = 1 ϵ n i . Mật độ xác suất tìm thấy điện tử thứ nhất ở vị trí → r 1 , điện tử thứ hai ở vị trí → r 2 . . . điện tử thứ N ở vị trí → r N bằng | ψ ( → r 1 , → r 2 , , → r N ) | 2 = |ψ n 1 ( → r 1 ) | 2 | ψ n 2 ( → r 2 ) | 2 | ψ n N ( → r N ) | . 2 Mât độ xác suất tìm thấy điện tử thứ k ở vị trí ψ r k bằng | → n k ( → r k ) | 2 . Mật độ điện tích của điện tử thứ k ở vị trí → r k bằng e k | ψ n k ( → r k ) | 2 với k = 1, , N. Thế năng tương tác Coulomb giữa điện tử thứ i ở vị trí → r i với điện tử thứ k ở vị trí → r k là V ik = ∫ e i e k | ψ n k ( → r k ) | 2 r ik d → r k , với d → r k = dx k dy k dz k r ik = | → r i − → r k | e i = e k = − e. Thế năng tương tác giữa điện tử thứ i với tất cả các điện tử k còn lại k cßn l¹i ( k ≠ i ) bằng U eff ( → r i ) = ∑ k ≠ i V ik = ∑ k ≠ i ∫ e 2 | ψ n k ( → r k ) | 2 r ik d → r k . Đó là biểu thức thế năng hiệu dụng của phương pháp trường trung bình trong gần đúng Hartree đưa ra năm 1928. 4/38 Gần đúng Hartree - Fock Trong phép gần đúng Hartree ở trên chúng ta chưa tính đến nguyên lí hệ các hạt đồng nhất. Các điện tử có spin bán nguyên s = 1 / 2 nên chúng tuân theo thống kê Fermi - Dirac và chúng thỏa mãn nguyên lí loại trừ Pauli. Trạng thái của điện tử i được đặc trưng bởi 3 tọa độ x i , y i , z i và một thành phần nữa là hình chiếu của spin s i lên phương OZ. Đối với điện tử s z có trị riêng là m s ℏ với m s = ± 1 / 2. Hàm sóng của điện tử i là hàm của các biến số tọa độ x i , y i , z i và s i kí hiệu các biến số này là ξ i ( i = 1, , N ) . Để mô tả trạng thái của điện tử có tính đến spin ta đưa vào hàm s như sau σ 1 2 ( ℏ 2 ) = 1 σ − 1 2 ( ℏ 2 ) = 0 σ 1 2 ( −ℏ 2 ) = 0, σ − 1 2 ( −ℏ 2 ) = 1. Khi đó ta có ∑ σ σ α * ( σ ) σ β ( σ ) = δ αβ . Nếu bỏ qua tương tác giữa mômen từ của điện tử với từ trường do điện tử chuyển động theo quỹ đạo gây nên thì ta có thể biểu diễn hàm sóng của điện tử i dưới dạng ψ k ( ξ i ) = ψ n k ( → r i ) σ α ( σ i ) , chỉ số k ở hàm ψ k ( ξ i ) kí hiệu trạng thái lượng tử ( n k , a ) . Điều kiện trực giao và chuẩn hóa của hàm ψ k ( ξ i ) ∫ ψ k * ( ξ i ) ψ k ( ξ i ) dξ i = = ∫ ψ n k * ( → r i ) ψ n l ( → r i ) ∑ σ i σ α ( σ i ) σ β ( σ i ) δ kl ≡ δ n k n l δ αβ . Phương trình Schrodinger có dạng ^ Hψ ( ξ 1 , , ξ N ) = E ψ ( ξ 1 , , ξ N ) . 5/38 Để phù hợp với nguyên lý loại trừ Pauli hàm Ψ ( ξ 1 , , ξ N ) phải là hàm phản đối xứng và nó có dạng là định thức Slater. ψ ( ξ 1 , , ξ N ) = = 1 √ N! ∑ ν ( − 1 ) ν P ν [ ψ k 1 ( ξ 1 ) ψ k N ( ξ N ) ] 1 √ N! | ψ k 1 ( ξ 1 ) ψ k N ( ξ 1 ) ⋱ ψ k 1 ( ξ N ) ψ k N ( ξ N ) | , trong đó ký hiệu P ν [ ψ k 1 ( ξ 1 ) ψ k N ( ξ N ) ] là hàm nhận được từ hàm ψ k 1 ( ξ 1 ) ψ k N ( ξ N ) bằng cách hoán vị ν cặp biến số ξ i ? ξ k bất kì cho nhau. Khi hoán vị bất kì một cặp chỉ số ξ i ? ξ k hay một cặp trạng thái k i ? k j cho nhau thì định thức đổi dấu. Khi ξ i = ξ k hay k i = k k thì định thức bằng 0 (không tồn tại hàm sóng) điều này thỏa mãn nguyên lí loại trừ Pauli (không tồn tại hơn một hạt trên 1 trạng thái lượng tử). Thực tế thì thế U eff ( → r i ) trong ^ H i ' còn chưa biết nên hàm Ψ n i ( ξ i ) là hàm riêng của ^ H i ' vẫn còn chưa xác định. Ta dùng nguyên lí biến phân để xác định U eff ( → r i ) . Gọi ψ 0 ( → r ) và E 0 là hàm sóng và năng lượng của hệ ở trạng thái cơ bản của hệ lượng tử với toán tử Hamiltonian ^ H và ψ 0 ( → r ) , E 0 thỏa mãn phương trình Schrodinger ^ Hψ 0 ( → r ) = E 0 ψ 0 ( → r ) , ψ ( → r ) là hàm sóng không phải ở trạng thái cơ bản (ở trạng thái kích thích) Năng lượng trung bình của hệ lượng tử trong trạng thái ψ là ¯ E = ∫ ψ * ( → r ) ^ Hψ ( → r ) d → r , vì E 0 là năng lượng ở trạng thái cơ bản (là nhỏ nhất ) nên ¯ E = E 0 nghĩa là 6/38 ∫ ψ * ( → r ) ^ Hψ ( → r ) d → r ≥ E 0 . Ta thấy các hàm ψ ( → r ) càng gần với hàm riêng ψ 0 ( → r ) bao nhiêu thì ¯ E càng gần E 0 bấy nhiêu. Ta chọn trước một lớp hàm ψ ( → r ) nào đó có dạng thích hợp rồi trong lớp hàm này chọn một hàm ψ ( → r ) sao cho giá trị ¯ E là nhỏ nhất (gần E 0 nhất) nghĩa là lời giải gần đúng nhất của bài toán vì ¯ E ứng với hàm ψ đã cho là nhỏ nhất nên δE = ¯ E − E 0 → 0. Vậy nghiệm gần đúng ψ 0 nhất phải thỏa mãn điều kiện δE = δ ∫ ψ * ( → r ) ^ Hψ ( → r ) d → r = 0 đó là nội dung của nguyên lí biến phân. Năng lượng trung bình của hệ N điện tử ¯ E = ∫ ψ * ( ξ 1 , , ξ N ) ^ Hψ ( ξ 1 , , ξ N ) dΓ với dΓ = dξ 1 dξ 2 ξ N . Thay hàm sóng (10) vào ta có năng lượng trung bình của hệ N điện tử ¯ E = + − N ∑ k = 1 ∫ ψ k * ( ξ i ) ^ H 0 ( ξ i ) ψ k ( ξ i ) dξ i 1 2 N ∑ ' k, l = 1 ∫ ψ k * ( ξ i ) ψ l * ( ξ j ) U ( ξ i , ξ j ) ψ k ( ξ i ) ψ l ( ξ j ) dξ i dξ j 1 2 N ∑ ' k, l = 1 ∫ ψ k * ( ξ i ) ψ l * ( ξ j ) U ( ξ i , ξ j ) ψ k ( ξ j ) ψ l ( ξ i ) dξ i dξ j , thay ψ k ( ξ i ) = ψ n k ( → r i ) σ α ( σ i ) và chú ý ∑ σ i σ α ( σ i ) σ β ( σ ) = δ αβ ta có 7/38 ¯ E = N ∑ k = 1 ∫ ψ n k * ( → r i ) ^ H 0 ( → r i ) ψ n k ( → r i ) d → r i + 1 2 N ∑ ' k, l = 1 ∫ ψ n k * ( → r i ) ψ n l * ( → r j ) U ( → r i , → r j ) ψ n k ( → r i ) ψ n l ( → r j ) d → r i d → r j − 1 2 N ∑ ' k, l = 1↑ ↑ ∫ ψ n k * ( → r i ) ψ n l * ( → r j ) U ( → r i , → r j ) ψ n k ( → r j ) ψ n l ( → r i ) d → r i d → r j . Trong số hạng cuối ta chỉ lấy tổng ứng với các cặp điện tử có spin định hướng song song cùng chiều (↑ ↑, ↓ ↓). Ta tính δ ¯ E rồi sau đó cho δ ¯ E = 0 ta có δ ¯ E = ∫ δψ n k * ( → r i ) ^ H 0 ( → r i ) ψ n k ( → r i ) δ → r i + N ∑ ' l = 1 ∫ δψ n k * ( → r i ) ψ n l * ( → r j ) U ( → r i , → r j ) ψ n k ( → r i ) ψ n l ( → r j ) d → r i d → r j − N ∑ ' l = 1↑ ↑ ∫ δψ n k * ( → r i ) ψ n l * ( → r j ) U ( → r i , → r j ) ψ n k ( → r j ) ψ n l ( → r i ) d → r i d → r j , thừa số 1/2 trong hai tổng cuối của ¯ E sẽ mất đi vì khi lấy biên phân theo δψ n k * ta gặp hai lần một lần theo tổng k một lần theo tổng l. Từ điều kiện chuẩn hóa hàm sóng ∫ ψ n k * ( → r i ) ψ n l ( → r i ) d → r i = δ n k , n l , ta suy ra ∫ δψ n k * ( → r i ) → n l ( → r i ) d → r i = 0.vớimọin k , n l . 8/38 [...]... trỡnh [link] ri gii C th tip tc cho n khi ta tỡm c nghim gn ỳng tt nht (tc l giỏ tr ca hai nghim liờn tip lin nhau khỏc nhau khụng ỏng k Trng Ueff c tớnh nh trờn c gi l trng t hp 10/38 p dng cho h nhiu in t v l trng trong chm lng t dng a vi th giam cm Parabolic Phng trỡnh Hartree- Fock cho h nhiu in t v l trng Trong phn ny chỳng ta s m rng bi toỏn trong h N in t t do cho bi toỏn tng quỏt: h nhiu in... trỡnh Hartree- Fock di dng ma trn cho h cỏc in t v l trng H phng trỡnh ny cho phộp xỏc nh cỏc h s khai trin hm súng t hp ca in t v l trng Ci và Ck Gii h phng trỡnh ma trn ny bng phng phỏp chộo hoỏ ma trn ta tớnh c cỏc h s Ci và Ck v khi bit cỏc h s ny thỡ cú ngha l ta ó tỡm c hm súng ca h , Thay th biu thc t (23) - (26) ca f, (1) và f (1) vo cụng thc [link] ta tớnh v , thu c biu thc cui cựng cho cỏc... i ij i j i j nờn N i (1) | {h(1) + j=1 M | l=1 j (2) l(2) | ( ) ^ e2 1 P 12 r12 e2 r12 | | l(2) i j (2) N + j=1 e2 r12 }| i (1) = 0 Bng cỏch tng t chỳng ta nhn c phng trỡnh Hartree - Fock cho hm súng t hp ca in t v l trng trong h nhiu Exciton (N in t v M l trng) õy l phng trỡnh tng quỏt cn gii 20/38 | j (2) N N {h(1) + ( J j K j )+ j=1 J j K j )+ j=1 { M h (1) + { h l ( Jl... hng tng tỏc trao i hai ht trng thỏi cú spin song song (cho in t) Tng t cho l trng ' ' e2 * d r 2 Jl k (1) = (2)l (2) (1)-sốhạngt-ơngtácđẩyCoulomb r12 k l 21/38 r12 = r | 1 | r 2 = e2 * d r2 (2)k(2) (1) r12 l l = Kl k(1) e2 * d r2 ^ (2) P (2)k(1), r12 12 l l s hng tng tỏc trao i hai chun ht trng thỏi cú spin song song (cho l trng) Phng trỡnh (30) c vit li: = i i (1)i = 1,... ( ) U , d , (r ) (r )(r r ) r ( r ) N ' l = 1 ( ) nl r 2 U r 1, r 2 d r 2, ) vỡ U r 1, r 2 = nl r 1 nk 1 e2 4 r 1 r 2 1 | | * nl 2 nk 2 1 2 2 Phng trỡnh [link] l phng trỡnh Hartree - Fock cho phộp ta xỏc nh hm súng t hp trng thỏi nk trong ú Ueff r 1 l trng hiu dng c xỏc nh bi [link] ( ) 9/38 gii [link] ta chn nghim nk gn ỳng no ú ó bit (chng hn hm súng ca mt in t t do hay hm... trỡnh (30) c vit li: = i i (1)i = 1, , N f(1)i (1) = i i (1)i = 1, , N f (1)k (1) = k k (1)k = 1, , M f (1)k (1) = k k (1)k = 1, , M f(1)i (1) 2 Hỡnh thc lun Roothaan Trong hỡnh thc lun Hartree - Fock - Roothaan chỳng ta tỡm i và k di dng khai trin theo h hm c s no ú Trong trng hp chm lng t hai chiu vi th giam cm parabolic t trong t trng ngoi ta cú th chn h hm c s , ny l h hm riờng ca Hamiltonian... vi ch s ri sau ú tớnh tng t cho cỏc ch s cũn li) ta cú: f(1) C (1) i i = C (1), i C f (1) (1) i i = C (1) i * Nhõn hai v phng trỡnh trờn v phớa phi vi (1) ri ly tớch phõn theo ton b * khụng gian d r 1 (1) v s dng iu kin trc giao chun hoỏ (1) | (1) = ta thu c: C (1) | f (1) | (1) = C i i i C (1) | f (1) | (1) = C i i i Tng t cho l trng: k C k C ... j(2) ( | ) l (2) + r12 e2 r12 M l (2) l (2) l=1 | e2 r12 | e2 r12 M l=1 | j(2) , ^ trong ú kớ hiu P 12 l toỏn t trao i bin (c a vo tin cho vic tớnh toỏn) ^ P 12(1) (2) = (2) (1) * Ly bin phõn E theo i (1) (tớnh vi ch s , cỏc ch s khỏc tớnh tng t), sau ú cho E = 0, ta c |{ | ( h(1) + (2) j r12 j=1 N + j=1 j (2) | e2 r12 | j (2) M l=1 l(2) ) ^ e2 1 P 12 N E = i (1) | l (2) | l... phi cú dng phn i xng chỳng tha món nguyờn lớ loi tr Pauli ca h cỏc ht ng nht) Hm súng ca h cú dmg ( e1, , eN, h1, , hM ì ) | ( ), , ( )| = 1 e1 N eN | ( ), , ( )| 1 h1 M hM trong ú l bin s c trng cho c to v spin Vi nng lng ^ E = * e1, , eN, h1, , hM H e1, , eN, h1, , hM ( ) ( ) Cỏc hm () tha món in kin trc giao chun húa ( ) ( )d d * i ei j ej ei ej = eiej ij * k hk l hl dhkdhl = hkhl kl... 2 d r 1 d r 2 r12 | ( )| ( ) ( ) * e2 * k r 1 l r 1 l r 2 k r 2 d r 1 d r 2 r12 ( )( ) ( )| |i r 1 2 ( ) ( )( )( ) e2 2 k r 2 d r 1 d r 2 r12 | ( )| ( ) ( ) Trong ú kớ hiu ' l tng ng cho cỏc giỏ tr ca i j, k l Vit di dng khai trin theo N, N, M, M vi N, M, N, M l s in t v l trng cú spin lờn ( ) v spin xung ( ) ( N + N = N, M + M = M) 13/38 N E= i (1) | h(1) | i (1) i=1 + 1 2 . 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 Giáo trình chuyên đề Vật lý Nano - Phương pháp trường. tuyến: http://voer.edu.vn/c/ac2ee593 MỤC LỤC 1. Phương pháp trường tự hợp Hartree - Fock áp dụng cho hệ nhiều điện tử 2. Áp dụng cho hệ nhiều điện tử và lỗ trống trong chấm lượng tử dạng đĩa với thế giam cầm Parabolic 3 các phương pháp gần đúng thông dụng là phương pháp Hartree - Fock. Nội dung của phương pháp này là chuyển việc nghiên cứu giải phương trình Schrodinger của hệ nhiều điện tử (hệ phương trình nhiều