Luận văn Phương pháp toán tử cho bài toán exciton hai chiều

Phụlục 9:Chuẩn hóa bộhàm sóng cơsởcho bài toán exciton hai chiều Trước hết, ta chọn bộhàm sóng của dao ñộng tử ñiều hòa (vì hàm này chắc chắn là nghiệm riêng của các toán tửtrung hòa nên sẽlà nghiệm riêng của 0ˆH )

pdf81 trang | Chia sẻ: maiphuongtl | Lượt xem: 2015 | Lượt tải: 0download
Bạn đang xem trước 20 trang tài liệu Luận văn Phương pháp toán tử cho bài toán exciton hai chiều, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
+ ( ) ( ) ( ) 1 11 1 ˆ ˆ ˆ0 1 0 1 1 ! 1 ! n n a n a n a n n n n + ++ + +    = = + = + +  +  . Tương tự, ta cũng thấy rằng mỗi toán tử sinh có tác dụng “sinh” (tăng) lên một bậc của vector trạng thái. Như vậy cứ có bao nhiêu toán tử sinh tác dụng lên vector trạng thái thì sẽ sinh thêm bấy nhiêu bậc của nó. 5. Chứng minh sự liên hợp của ˆ ˆ+a,a , 1 , 1 ˆ 1 , ˆ 1 , n j n j n a j j n j j j a n j j n j − + − = − = = − = δ δ ˆ ˆn a j j a n+⇒ = . Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 41 Nhận xét: Từ các tính chất (3, 4, 5) ở trên ta thấy rằng: nếu như tác dụng một toán tử chứa cùng số toán tử sinh và toán tử hủy lên một vector trạng thái, thì sẽ không làm vector này thay đổi bậc, và ta gọi các toán tử như thế là toán tử “trung hòa”; ngược lại nếu toán tử chứa số toán tử sinh – hủy khác nhau thì sẽ làm thay đổi bậc của vector trạng thái. Đây là một tính chất rất quan trọng trong các tính toán đại số khi sử dụng biểu diễn toán tử và cũng chính là yếu tố để ta tiến hành việc tách toán tử Hamilton của hệ thành hai thành phần: trung hòa và nhiễu loạn. Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 42 Phụ lục 2. Dạng chuẩn (normal) của một số toán tử trong luận văn Dạng chuẩn (normal) của một toán tử được định nghĩa là dạng đã được biến đổi sao cho toán tử hủy luôn về phía bên phải của biểu thức, toán tử sinh luôn về phía bên trái của biểu thức. aˆ+ trái aˆ phải. Mục đích của việc đưa các biểu thức toán tử về dạng chuẩn là giúp cho việc tính toán trong các bài toán chứa nhiều loại toán tử được dễ dàng hơn rất nhiều. Thực vậy, khi biểu biễn tất cả trạng thái qua trạng thái cơ bản 0( )ω thì lợi dụng tính chất ˆ 0( ) 0a =ω và ˆ 0( ) 0b ω = , chúng ta sẽ biểu diễn tất cả trạng thái còn lại qua biểu thức chỉ còn một loại toán tử sinh tác dụng. A. Trường hợp các toán tử sinh, hủy với số mũ lũy thừa Trường hợp này ta chỉ cần áp dụng các tính chất của giao hoán tử trên là có thể đưa về dạng chuẩn. Ví dụ: Đưa toán tử ( )22ˆ ˆa a+ về dạng chuẩn ta thực hiện như sau: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 22 2 2 2 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ1 ˆ ˆ ˆ ˆ ˆ ˆ1 1 1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ1 1 2 ˆ ˆ ˆ ˆ ˆ ˆ2 3 1 ˆ ˆ ˆ ˆ ˆ ˆ2 3 ˆ ˆ ˆ ˆ2 4 . a a a aa a a a a a aa aa aa a a a a a a a a a a a aa a a a a a a a a a a a a a a a a a + + + + + + + + + + + + + + + + + + + + + + + = = + = + = + + + + = + + + + = + + + = + + + = + + Các phép biến đổi trên thường được áp dụng khi các biểu thức toán tử có dạng như các đa thức. Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 43 B. Trường hợp hàm e mũ của các toán tử sinh, hủy Đối với dạng hàm e mũ thì khi vận dụng phép biến đổi như trên sẽ gặp khó khăn. Vì các toán tử sinh hủy trên mũ khi khai triển để đưa về dạng chuẩn sẽ có bậc lũy thừa rất cao. Nên ta phải áp dụng phương pháp biến đổi khác như dưới đây. Ví dụ: ( )ˆ ˆt a ae + + Vì ta có hệ thức giao hoán ˆ ˆ, 1a a+  =  nên từ đây các toán tử ˆ ˆ,a a + và số 1 tạo thành một đại số kín. Như vậy ta có thể viết: ( ) ( )ˆ ˆ ˆ ˆ( ) ( ) ( )t a a f t a g t a h te e e e F t+ ++ = = . (A2.1) và tiến hành tìm các hàm số ( ), ( ), ( )f t g t h t theo các bước sau: Bước một: Lấy đạo hàm hai vế của (2.1) theo biến số t ta có: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )ˆ ˆˆ ˆ ˆ ˆ' ' 't a aa a e f t a F t g t aF t h t F t+ ++ ++ = + + . (A2.2) Định nghĩa hàm nghịch đảo của ( )F t là ( )1F t− sao cho ( ) ( )1. 1F t F t− = ta có: ( ) ˆ ˆ1 ( ) ( ) ( )h t g t a f t aF t e e e +− − − −= . (A2.3) Nhân hai vế (2.2) cho ( )1F t− và thu gọn các số hạng ta được: ( ) ( ) ( )ˆ ˆ( ) ( )ˆ ˆ ˆ ˆ' ' 'f t a f t aa a f t a g t e ae h t+ ++ + −+ = + + (A2.4) Bước hai: Sử dụng công thức quen thuộc (phụ lục 1): ˆ ˆ 1 1 2! 3! ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ A A e B e = B+ A,B + A, A,B + A, A, A,B +...−                    cùng với hệ thức giao hoán của ˆ ˆ,a a+ ta có: ( ) ( )ˆ ˆ( ) ( )ˆ ˆ ˆ ˆ ˆ, ...f t a f t ae ae a f t a a a f t+ +− + = + + = −  . Thay vào (2.4), ta có: ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ˆ ˆ ˆ ˆ' ' ' ˆ ˆ' ' ' ' . a a f t a g t a f t h t f t a g t a h t g t f t + + + + = + − + = + + − (A2.5) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 44 Bước ba: Đồng nhất hai vế của (2.5) và chọn điều kiện biên Đồng nhất hai vế, ta có hệ phương trình: ( ) ( ) ( ) ( ) ( ) ' 1, ' 1, ' ' 0. f t g t h t g t f t =  =  − = Giải hệ này ta có: ( ) ( ) ( ) 1 2 2 1 3 , , . 2 f t t c g t t c th t c t c        = + = + = + + Dựa vào biểu thức (2.1), ta có điều kiện khi t = 0 thì: f(t) = g(t) = h(t)= 0. Suy ra: c1= c2 = c3 = 0. Như vậy dạng chuẩn của ( )ˆ ˆt a ae + + là: ( ) 2ˆ ˆ ˆ ˆ / 2t a a ta ta te e e e+ ++ = . (A2.6) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 45 Phụ lục 3: Yếu tố ma trận cho toán tử Hamilton của dao động tử phi điều hòa A. Tính các yếu tố ma trận của toán tử Hamilton (phương pháp giải tích) 3.1 Tính các yếu tố ma trận: Theo [1] ta có: 1 1 1( ) . ( ) ( ) 2n n n H n H Hξ ξ ξ ξ − += + . Khi đó yếu tố ma trận H được tính: 2 (0)* (0) (0)* 2 (0) 0 2 1 1 ˆ ( ) 2 2nn n n n n dH H dx x dx dx +∞ +∞ −∞ −∞ = Ψ Ψ = Ψ − + Ψ∫ ∫ 2 2 2 2 2 (0)* 2 2 1 1( 1) exp exp 2 2 2 2 n n n x x n n n n d x d x dA e x e dx dx dx dx +∞ − − −∞       = − Ψ − +              ∫ 2 2 2 2 2 2 2 2 2 1 2 1 2 (0)* 2 2 1 2 2 1 2 exp exp exp 2 2 21 1( 1) exp 2 2 2 exp exp 2 2 n n n x x x n n n n n x n n nn n x x n n x d x d x d e x e x e dx dx dx x dA x e dxx d x d x e e dx dx + − − − ++∞ − + + −∞ − − + +         + +                 = − Ψ − +           + +           ∫ dx  ( ) 2 (0)* 1 2 1 1 ˆ2 exp 2 2 2 1 11 2 2 n n n n n A x n n n n A xH x dx A n n +∞ + + −∞   = − − + + Ψ −    = − + + = + ∫ 3.2 Tính yếu tố ma trận V Từ 1 1 1( ) . ( ) ( ) 2n n n H n H Hξ ξ ξ ξ − += + , ta tính được: 2 1 1 2 2 1 1 1 . ( 1) . , 2 2 4n n n n n n H n H H n H n n H Hξ ξ ξ − + − +   = + = + + − +    3 2 2 2 1 3 1 3 1 1 3 3 1( ) ( 1) . ( 1)( 2) ( 1) 2 4 2 4 8n n n n n n n n H n H n n H H n H n n n H n H Hξ ξ ξ ξ ξ − + − − + +   = + + − + = + − − + + +    4 2 2 4 4 2 3 1 1(2 2 1) (2 3) ( 1)( 2)( 3) ( 1)(2 1) 4 4 16n n n n n n H n n H n H H n n n H n n n Hξ + + − −= + + + + + + − − − + − − Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 46 Tính: ( ) ( ) 2 2 *(0) (0) *(0) 4 *(0) 22 2 4 4 2 ˆ( ) ( ) exp 2 3 1 1 . 2 2 1 (2 3) ( 1)( 2)( 3) ( 1)(2 1) 4 4 16 m n m n n x m n n n n n mn x x V dx x A x H x dx e n n H n H H n n n n H n n n H dx V λ λ λ +∞ +∞ −∞ −∞ +∞ − + + − − −∞   = Ψ Ψ = Ψ − =      = Ψ + + + + + + − − − + − −    ∫ ∫ ∫ ( )2 2 4 4 2 2 4 4 2 3 2 2 1 (2 3) ( 1)( 2)( 3) ( 1)(2 1) 4 4 16 n n n n mn n n n n n n n n A A A A n n n n n n n n n n A A A A λ δ δ δ δ δ+ + − − + + − −   = + + + + + + − − − + − −    Khi đó: ( ) 2 *(0) (0) *(0) 4 , 4 4 4 ˆ( ) ( ) exp ( 4)( 3)( 2)( 1) 2 4n n n n n n n xV x V dx x A x H x dx n n n nλλ λ +∞ +∞ + + + −∞ −∞   = Ψ Ψ = Ψ − = + + + +    ∫ ∫ ( ) 2 *(0) (0) *(0) , 2 2 2 ˆ( ) ( ) exp (2 3) ( 2)( 1) 2 2n n n n n n n xV x V dx x A H x dx n n nλλ λ +∞ +∞ + + + −∞ −∞   = Ψ Ψ = Ψ − = + + +    ∫ ∫ B. Tính các yếu tố ma trận của toán tử Hamilton( OM) Ta có: ( ) ( ) ( ) ( ) ( ) ( ) 2 2 2 22 4 4 3 24 3 2 2 1 1 3 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 1 2 2 1 4 4 4 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ4 4 6 6 . 4 H a a a a a a a a a a a a a a a a ω ω λ ω ω ω λ ω + + + + + + + +                 + − = + + + + + + + + + + + + Ta tách ˆH thành hai phần: 0ˆ ˆ ˆH H V= + , ( ) ( )2 20 21 3ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 1 2 2 14 4H a a a a a aω λω ω+ + +   += + + + + , ( )( ) ( ) ( ) ( )2 2 4 3 22 3 2 421ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ4 6 4 6 .4 4V a a a a a a a a a aω λω ω+ + + + +   −= + + + + + + + Ta tính các phần tử ma trận khác không của ˆH : ( ) ( ) ( ) ( ) ( ) 2 2 0 2 2 2 2 1 3 ˆ ˆ ˆ ˆ ˆ ˆ2 1 2 2 1 4 4 1 32 1 2 2 1 , 4 4 nn nn H H n a a a a a a n n n n ω λ ω ω ω λ ω ω + + +     + = = + + + + + = + + + + ( ) ( ) ( )( ) ( ) ( ) 2 2 3 2 , 2 2 2 2 2 2 1 ˆ ˆ ˆ ˆ4 6 2 4 4 2 !1 14 6 2 1 2 3 , 4 4 !4 2 n nV n a a a a n n n n n n n ω λ ω ω ω λ ω λ ω ωω ω + + − = + + + +   − − = + + + + = + +        Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 47 ( )4 , 4 2 2 4 ! ˆ 4 4 4 !n n n V n a n n λ λ ω ω+ + = + = ; các phần tử ma trận khác được tính dựa vào tính đối xứng: nm mnV V= . Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 48 Phụ lục 4: Phương trình Schrödinger cho bài toán exciton hai chiều Theo cơ học cổ điển, năng lượng của hệ gồm hai hạt tương tác ( ) 2 2 1 2 1 22 2 p pE U r m m = + + . trong đó r là khoảng cách giữa hai hạt, một cách tương ứng Hamiltonain của hệ bằng: ( ) 2 2 2 2 1 2 1 2 ˆ 2 2 H U r m m = − ∇ − ∇ +  . Gọi r1 và r2 là các bán kính vector của hạt 1 và hạt 2, r là bán kính vector từ hạt 1 sang hạt 2, R là bán kính vector của tâm bán kính G. Chúng ta có các hệ thức: 1 1 2 2 2 1 1 2 , m r m r r r r R m m + = − = + . Chiếu hai biểu thức này xuống trục x ta có: 1 1 2 2 2 1 1 2 1 2 ; m x m x x x x X m m m m = − = + + + . Theo hệ thức trên ta có: 1 1 1 1 1 2 x X m x x x x x x m m X ∂ ∂ ∂ ∂ ∂ ∂ ∂ = + = − + ∂ ∂ ∂ ∂ ∂ ∂ + ∂ ; 2 2 22 2 2 2 1 1 1 2 2 2 1 1 1 2 1 2 1 2 2m m m x x x m m X x m m x X m m X      ∂ ∂ ∂ ∂ ∂ ∂ ∂ = = − + = − +     ∂ ∂ ∂ + ∂ ∂ + ∂ ∂ + ∂      . Tương tự: 2 2 22 2 2 2 1 2 2 2 2 2 2 2 1 2 1 2 1 2 2m m m x x x m m X x m m x X m m X      ∂ ∂ ∂ ∂ ∂ ∂ ∂ = = − + = − +     ∂ ∂ ∂ + ∂ ∂ + ∂ ∂ + ∂      . Từ hai biểu thức trên ta có: 2 2 2 2 2 2 2 2 1 1 2 2 1 2 2 1 1 1 1 1 1 m x m x m m x m m X  ∂ ∂ ∂ ∂ + = + + ∂ ∂ ∂ + ∂  . Tiến hành trên hai trục còn lại ta thu được Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 49 2 2 2 2 1 2 1 2 2 1 1 1 1 1 r G m m m mµ ∇ + ∇ = ∇ + ∇ + , trong đó 1 2 1 2 m m m m µ = + gọi là khối lượng rút gọn. Khi đó toán tử Hamitonain có dạng: ( ) 2 2 2 2 2 1 ˆ 2 2( )r GH U rm mµ   = − ∇ − ∇ Ψ +  +    (A4.1) * Dạng không thứ nguyên của phương trình thứ (2.9) ( ) ( ) 2 2 2 2 r r r Ze r E r r ψ ψ µ   − ∇ − =     . Hamilton có dạng: 2 2 2 ˆ 2 r ZeH rµ = − ∇ − . đặt . , .x ya x a yρ ρ= = và E bε= ; 2 2 2 2 2 2 2 2 2 1 ; x y x x a a r x y x x a ρ ρ ρ ρ ∂ ∂ ∂ ∂ = ⇒ = = + = + ∂ ∂ ∂ ∂ ; 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 2 2x y x y a aZe e Z bb a a µ µψ εψ ψ εψ µ ρ ρ ρ ρ ρ ρ       ∂ ∂ ∂ ∂    − + − = ⇒ − + − =         ∂ ∂ ∂ ∂             . Ta đặt: 2 2 2 2 0 . 11e ea a r µ µ = ⇒ = =   ; 2 2 4 2 2 2 .1b a eb a µ µ µ = ⇒ = =    . Khi đó ta thu được phương trình Schrödinger không thứ nguyên sau: 2 2 2 2 1 2 x y Z ψ εψ ρ ρ ρ   ∂ ∂  − + − =   ∂ ∂    . (A4.2) Để thuận tiện ta có thể viết lại phương trình Schrödinger không thứ nguyên có dạng: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 50 2 2 2 2 2 2 1 2 Z x y x y ψ εψ   ∂ ∂  − + − =  ∂ ∂ +    . (A4.3) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 51 Phụ lục 5: Hamilton cho bài toán exciton hai chiều Chuyển từ tọa độ vuông góc sang hệ tọa độ cực: os sin x rc y r ϕ ϕ =  = Chuyển từ tọa độ cực sang hệ tọa độ vuông góc: 2 2 arctg r x y y x ϕ  = +   =  . 1. Toán tử Hamilton ˆH Để chuyển toán tử 2 2 2 2x y ∂ ∂∆ = + ∂ ∂ sang hệ tọa độ cực, ta tiến hành ta tiến hành chuyển các đạo hàm theo tọa độ vuông góc sang tọa độ cực. Lấy ví dụ là đạo hàm theo biến x. Theo quy tắc đạo hàm của hàm số hợp: r r x x x ϕ ϕ ∂ ∂ ∂ ∂ ∂ = + ∂ ∂ ∂ ∂ ∂ (A5.1) Muốn thực hiện phép chuyển này ta phải tính các đạo hàm riêng: ,r x x ϕ∂ ∂ ∂ ∂ , os 1 sin r x c x r x r ϕ ϕ ϕ ∂ = =∂ ∂ = −  ∂ , (A5.2) Thay (5.2) vào (5.1) ta được: 2 2 2 2 2 2 2 2 2 1 os sin r 1 1 os sin os sin r r r 1 2 2 1 os sin sin os sin os sin r rr c x r r c c x x r x r c c c r r r r ϕ ϕ ϕ ϕϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ∂ ∂ ∂ = − ∂ ∂ ∂    ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ = − + −   ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂    ∂ ∂ ∂ ∂ ∂ ∂ = + − + + ∂ ∂ ∂ ∂∂ ∂ Tương tự tính cho y ∂ ∂ , với các đạo hàm riêng: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 52 sin ; 1 os r y c y r ϕ ϕ ϕ ∂ =∂ ∂ =  ∂ Ta có: 2 2 2 2 2 2 2 2 2 2 1 1 sin os sin os r r r 1 1 1 sin sin os os sin os r r r 2 2 1 1 sin sin os sin os os os r rr r c c y r y ry c c c r r r c c c c r rr r ϕϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ ϕ    ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂ = + + +   ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂∂        ∂ ∂ ∂ ∂ ∂ ∂ = + + +   ∂ ∂ ∂ ∂ ∂ ∂    ∂ ∂ ∂ ∂ ∂ ∂ = − + + + ∂ ∂ ∂ ∂∂ ∂ 2 Khi chuyển 2 2 2 2x y ∂ ∂∆ = + ∂ ∂ sang hệ tọa độ cực ta được: 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 = . rr r r r r rx y r rϕ ϕ ∂ ∂ ∂ ∂ ∂ ∂ ∂ ∂  + + + = + ∂ ∂ ∂∂ ∂ ∂ ∂ ∂  . Từ công thức của toán tử Laplace trong tọa độ cực, ta có công thức của toán tử Hamilton của electron. 2 2 2 1 1 ˆ ( ) 2 2 H r U r r r r r ϕ ∂ ∂ ∂  = − − + ∂ ∂ ∂  (A5.3) 2. Toán tử ˆxL ˆ ˆ ˆ x z yL yp zp i y z z y  ∂ ∂ = − = − − ∂ ∂   (A5.4) Thay các biểu thức ở phần trên vào (3) ta có: 1 sin sin os sin 1 1 os os sin sin os sin r sin y r c z r r c z rc c y r r θ ϕ θ θ θ ϕθ θ ϕ θ ϕ θ θ ϕ ∂ ∂ ∂  = − ∂ ∂ ∂   ∂ ∂ ∂ ∂ = + + ∂ ∂ ∂ ∂  Từ đó ta thu được: ˆ ˆ ˆ sin cot osx z yL yp zp i g cϕ θ ϕθ ϕ  ∂ ∂ = − = − + ∂ ∂   Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 53 3. Với toán tử ˆ ˆ,y zL L : ˆ ˆ y Z L i z y y z L i x y y x  ∂ ∂ = − − ∂ ∂   ∂ ∂ = − − ∂ ∂    Tương tự như toán tử ˆxL , ta cũng thay các đạo hàm riêng có được ở trên vào: ˆ os cot sin ˆ y Z L i c g L i ϕ θ ϕ θ ϕ ϕ  ∂ ∂ = − − ∂ ∂  ∂ = − ∂   4. Tìm riêng và trị riêng của toán tử ˆzL Phương trình hàm riêng- trị riêng của ˆzL : ˆ z z z uL u L u i L u ϕ ∂ = → − = ∂ . Vì hàm U chỉ phụ thuộc vào biến số ϕ nên ta thay đạo hàm riêng toàn phần thành đạo hàm toàn phần: ( ) . z z iL dui L u d u C e ϕ ϕ ϕ − = = , hệ số C được xác định từ điều kiện chuẩn hóa: 2 2 0 ( ) 1u d pi ϕ ϕ =∫ , ta được 1 . 2 C pi = Khi ϕ thay đổi một lượng 2pi thì hạt trở lại vị trí ban đầu. Do đó , để ( )u ϕ xác định đơn giá thì ( 2 ) 2 1 2 2 m= 0, 1, 2... z z z iL iL i L z z e e e L m L m ϕ ϕ pi pi ϕ pi pi + = = = = ± ± Vậy hàm riêng của toán tử ˆzL là : 1( ) . 2 imu e ϕϕ pi = , Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 54 và trị riêng là m= 0, 1, 2...zL m= ± ± Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 55 Phụ lục 6: Các toán tử sinh – hủy hai chiều ˆ ˆ( ) , ( ) , 2 2 ˆ ˆ( ) , ( ) ; 2 2 x x x x x x y y y y y y 1 1 a x a x x x 1 1b y b y y y ω ω ω ω ω ω ω ω ω ω ω ω + +    ∂ ∂ = + = −   ∂ ∂       ∂ ∂ = + = −      ∂ ∂    Suy ra ( ) ( )22ˆ ˆ ˆ ˆ 22 a a a a x x ωω + ++ + = → = ( ) ( )22ˆ ˆ ˆ ˆ 22 b b b b y y ωω + ++ + = → = ( ) ( ) ˆ ˆ 2 ˆ ˆ 2 a a x b b y ω ω + + ∂ = − ∂ ∂ = − ∂ Suy ra : ( ) ( ) ( ) ( ) ( )2 2 2 21 1 1ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ2 2 2x y a a a a b b b b a a b bx y + + + + + +∂ ∂+ = + − + + − = − + −∂ ∂ Ta có: ( ) ( ) ( ) 2 2 22 2 2 22 2 ˆ ˆ ˆ ˆ ˆ ˆ2 1 ; 2 2 ˆ ˆ ˆ ˆ2 1 . 2 x x y a a a a a a x b b b b y ω ω ω + + + + + ∂   = − = − − +  ∂ ∂   = − − +   ∂ ( ) ( ) ( ) ( ) 2 2 22 2 2 2 2 2 2 ˆ ˆ ˆ ˆ 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 2 2 ( ) ( ) 2 1 ˆ ˆ ˆ 2 a a b b x y a b a a b b a b M N M ω ω ω + + + + + + + ∂ ∂  + = − + −   ∂ ∂ = + − − − + + = − + Mặt khác: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 56 ( ) ( ) 22 2 22 2 1 ˆ ˆ ˆ ˆ2 1 ; 2 1 ˆ ˆ ˆ ˆ2 1 ; 2 x y x a a a a y b b b b ω ω + + + +   = + + +     = + + +    Suy ra: ( ) ( ) ( ) ( ) ( ) 22 2 2 222 2 ˆ ˆ ˆ ˆ 2 2 1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ 2 1 ˆ ˆ ˆ 2 b ba a x y a b aa bb a a b b a b M N M ω ω ω ω ++ + + + + + + + ++ + = +   = + + + + + + +    = + + ( )( ) ( )( ) ( )( ) ( )( ) ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ 2 yx z y x x y x y x y iL ix y b b a a a a b b y x i ab a b ab a b ωω ω ω ω ω ω ω ω ω + + + + + + + +   ∂ ∂  =− − = + − − + −     ∂ ∂     = − − + + −   Khi x yω ω ω= = , ta có: ( ) ( )2 2 222 22 2 ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ2 2 2 ;2 a b a b a a b bx y ω + + + +∂ ∂  + = + + + − − −  ∂ ∂ ( ) ( )222 2 2 21 ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ2 2 2 .2x y a b a b a a b bω + + + + + = + + + + + +   ( )ˆ ˆˆ ˆ ˆzL i ab a b+ += − . A. Để thuận tiện trong tính toán, ta sử dụng các toán tử: ( ) ( ) 22 22 ˆ ˆˆ ˆ ˆ ˆ ˆ2 1, , , ˆ ˆ ˆ ˆˆ ˆ ˆ2 1, , , x y N a a A a A a N b b B b B b + + + + + + = + = = = + = = ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ; ;x yN N N M A B M A B + + + = + = + = + trong đó từng bộ ba toán tử ˆ ˆˆ , ,xN A A+ , ˆ ˆ ˆ, ,yN B B+ , ˆ ˆ ˆ, ,M M N + tạo thành các đại số kín, thỏa mãn các hệ thức giao hoán: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 57 ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ , 2 , , 4 , , 4 , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ , 2 , , 4 , , 4 , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ , 2 , , 4 , , 4 . x x x y y y A A N A N A N A A B B N B N B N B B M M N M N M N M M + + + + + + + + +      = = =           = = =           = = =      các giao hoán tử khác bằng 0. Tính các giao hoán tử: ( ) ( ) ( ) ( ) ( ) 2 2 22 1 1 1 1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , , , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , , , ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 2 2 1 2 x A A a a a a a a a a aa a a a a a a a a a a a a a a a a aa a a N + + + + + + + + + + + + = = = = + + +        = = +                    = + + +        = + = + =     ( ) [ ]  [ ]  2 2 2 0 0 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, ,2 1 2 , 2 , , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 , , , , ˆ ˆ4 4 xA N a a a a a a a a a a a a a a aa a a a a a a a a a a a a a a A + + + + + + + + = =          = + = = +               = + + +       = = ( ) ( ) ( ) ( ) ( )( ) ( ) 2 2 2 2 0 2 2 ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, 2 1, 2 , 2 , , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 , , 4 4 xN A a a a a a a a a a a a a a a a a a a a a A + + + + + + + + + = + + + + + + +             = + = = +                        = + = =     0 0 ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ , , , , , , ˆ ˆ ˆ2 2 2x y M M A B A B A A A B B A B B N N N + + + + + + + = =           = + + = + + +           = + =   00 ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ , , , , , , ˆ ˆ ˆ4 4 4 x y x y x yM N A B N N A N A N B N B N A B M ==           = + + = + + +           = + =  0 0 ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ , , , , , , ˆ ˆ ˆ4 4 4 x y x x y yN M N N A B N A N B N A N B A B M + + + + + + + = = + + +           = + + = + + +           = + =   B. Chứng minh toán tử ˆH giao hoán với toán tử ˆzL Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 58 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, ( )( ) ( ) ( ) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ( ) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ( 1) ( 1) ( 1) ZL M i ab a b a a b b a a b b i ab a b i ab a a a a ab b b a b a bb b i a a a b a aa b b b b a + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +   = − + − + −  = − + − = + − − + − − ˆ ˆ ˆ ˆ( 1) ˆ ˆ ˆ ˆ 2 2 0 b b b a i a b a b + + + + + + +  +   = − =  ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, ( )( ) ( ) ( ) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ( ) ˆ ˆ ˆ ˆ 2 2 0 ZL M i ab a b aa bb aa bb i ab a b i ab bb aaab aaa b a aab i ab ab + + + + + + + +   = − + − + −  = − + −  = − + =  + + + + + + + + + + + + + + + + + + ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, ( )2( ) 2( ) ( ) ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ =2i(aa ab -a aab +b bab -bb ab +a aa b-a a ab+bb a b-bba b ˆ ˆ ˆ ˆ =2i(ab -ab ZL N i ab a b a a bb a a bb i ab a b + + + + + + + +  = − + − + −  + +ˆ ˆ ˆ ˆ+ab -ab )=0 Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 59 Phụ lục 7: Dạng chuẩn của toán tử { }ˆ ˆ ˆ ˆexp ( )S M M Nτ += − + + Do các toán tử ˆ ˆ ˆ, ,M M N+ tạo thành một đại số kín, ta có thể sử dụng công thức cho hai toán tử không giao hoán bất kỳ ˆ ˆ,X Y : ( ) ( ) ( ) ( )ˆ ˆ ˆ ˆ ˆ ˆ1exp exp [ , ] exp exp2X Y X Y X Y+ = − , (A7.1) ta có thể đưa toán tử ˆS về dưới dạng chuẩn như sau: ( ) ( ) ( ) ( )ˆ ˆ ˆ ˆ ˆ ˆexp ( ) exp ( ) exp ( ) ( ) exp ( )M M N f M g N h Mτ τ τ τ+ +− + + = . (A7.2) Ở đây chúng ta cần xác định các hàm số ( ), ( ), ( )f g hτ τ τ với điều kiện biên: (0) 0, (0) 0, (0) 0f g h= = = . (A7.3) Bước một: Lấy đạo hàm hai vế của (A7.2) theo τ : ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ˆ ˆ ˆ ˆ ˆ ˆ( )exp ( ) ˆ ˆ ˆ ˆ ˆ '( ) '( ) exp ( ) exp ( ) exp ( ) ˆ ˆ ˆ ˆ '( )exp ( ) exp ( ) exp ( ) . M M N M N M f M F g f M N g N h M h f M g N M h M τ τ τ τ τ τ τ τ τ τ τ + + + + + − + + − + + = = + + (A7.4) Nhân (A7.4) với toán tử ngược 1ˆS − : với toán tử 1ˆ ˆ. 1S S − = , ( ) ( ) ( ) ( ) ( ) ( ) ( ) ˆ ˆ ˆ ˆ ˆ ˆ ˆ( ) ( ) exp ( ) exp ( ) ˆ ˆ ˆ ˆ ˆ( ) exp ( ) exp ( ) exp ( ) exp ( ) M M N f M g t f M N f M h f M g N M g N f M τ τ τ τ τ τ τ τ + + + + + + ′ ′ − + + = + − ′+ − − (A7.5) Bước hai: Ta sử dụng công thức ( ) ( ) 1 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆexp exp , , , , , ,2! 3!X Y X Y X Y X X Y X X X Y         − = + + + +          … Tính lần lượt các thành phần của (A7.5): ( ) ( ) 21ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ( )exp ( ) exp ( ) ( ) ( ) , ( ) , , ...2!I g f M N f M g N f M N f M M Nτ τ τ τ τ τ+ + + + +     ′ ′= − = + + +       ( ) 21 ˆ ˆ ˆ ˆ( ) ( )( 4 ) ( ) , 4 ... 2! ˆ ˆ( ) 4 ( ) g N f M f M M g N f M τ τ τ τ τ + + + +   ′= + − + − +    ′= − (A7.6) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 60 Trước hết ta tính thành phần: ( ) ( ) 2 2 4 ( ) 1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆe x p ( ) ex p ( ) ( ) , ( ) , . .. 2 ! ˆ1 6 ˆ ˆ4 ( ) . .. 2 ! ˆ . g J g N M g N M g N M g N N M MM M g M e τ τ τ τ τ τ −     = − = + + +     = − + + = (A7 .7) Tiếp theo ta tính thành phần: ( ) ( ) ( ) ( ) ( ) ( ) ( ) 4 ( ) 4 ( ) 2 4 ( ) 2 ˆ ˆ ˆ ˆ ˆ( )exp ( ) exp ( ) exp ( ) exp ( ) ˆ ˆ ˆ( ) exp ( ) exp ( ) 1 ˆ ˆ ˆ ˆ ˆ ˆ( ) ( ) , ( ) , , ... 2! ˆ ˆ ˆ( ) 2 ( ) 4 ( ) g g g K h f M g N M g N f M e h f M M f M e h M f M M f M M M e h M Nf M f τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ + + − + + − + + + − + ′= − − ′= −      ′= + + +       ′= − + (A7.8) Thay vào trong biểu thức: ( ) ( ) ( )4 ( ) 2 4 ( ) 4 ( ) 4 ( ) 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ( ) ( ) 4 ( ) ( ) 2 ( ) 4 ( ) ˆ ˆ ˆ ˆ ˆ ˆ( ) ( ) 4 ( ) ( ) ( ) 2 ( ) ( ) 4 ( ) ( ) g g g g M M N f M g N f M e h M Nf M f f M g N g f M e h M e h Nf e h f M τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ τ + + + − + + + − − − + ′ ′ ′ − + + = + − + − + ′ ′ ′ ′ ′ ′= + − + − + Bước ba: Đồng nhất các hệ số trước các toán tử ˆ ˆ ˆ, ,M M N+ ta thu được hệ phương trình vi phân trên để xác định các hàm số ( ), ( ), ( )f g hτ τ τ : 4 ( ) 2 4 ( ) 4 ( ) 1 ( ) 4 ( ) ( ) 4 ( ) ( ) 1 2 ( ) ( ) ( ) 1 ( ) g g g f g f e h f e h f g e h τ τ τ τ τ τ τ τ τ τ τ τ − − − ′ ′ ′− = − +  ′ ′ − = − +  ′ − = (A7.9) . ( )2 1 2 1( ) 1 2 ( ) ( ) . 2 2 Cf f f C τ τ τ τ τ − + − ′→ = − + → = + áp dụng điều kiện biên 1(0) 0 (0) 0 1 ( ) 2 2 1 Cf f C f C τ τ τ − − = → = = → = → = + (A7.10) Thay vào trong hệ phương trình: 1( ) ln 2 1 2 g Cτ τ= − + + Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 61 ( ) ( ) 1 1 1 1 1( ) 1 2 2 1 2 2 2 1 2 1 h ττ τ τ τ   − = − = − =  + + +  Áp dụng điều kiện biên: ( ) ( ) 1 1 1 1 1( ) 1 2 2 1 2 2 2 1 2 1 h ττ τ τ τ   − = − = − =  + + +  (A7.11) ( ) ( ) 1 1 1 1 1( ) 1 2 2 1 2 2 2 1 2 1 h ττ τ τ τ   − = − = − =  + + +  (A7.12) Như vậy ta đã tìm được dạng chuẩn của toán tử ( ) 1ˆ ˆ ˆ ˆ ˆ ˆexp ( ) exp exp ln 2 1 ( ) exp .2 1 2 2 1M M N M N Mτ ττ ττ τ+ +− −     − + + = − +     + +      (A7.13) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 62 Phụ lục 8: Chuyển toán tử Hamilton qua biểu diễn toán tử sinh –hủy 2 22 2 ( ) 2 2 0 1 ˆ 2 t x yZ eH dt x y tpi +∞ − + ∂ ∂ = − + − ∂ ∂  ∫ ( ) ( ) 0 ˆ ˆ ˆ ˆ ˆ ˆ ˆexp 4 2 Z dt tH M M N M N M t ω ωpi +∞ + +−  = − + − − + +    ∫ ( ) ( ) 0 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆexp 4 dH M M N Z M N Mω ω τ τ pi τ +∞ + + = − + − − − + +  ∫ với toán tử có dạng hàm mũ ( ) ( )ˆ ˆ ˆ ˆexpS N M Mτ τ + = − + +  có thể đưa về dạng chuẩn như sau ( xem phụ lục 5): ( ) 1ˆ ˆ ˆ ˆexp exp ln 2 1 ( ) exp 2 1 2 2 1 S M N Mτ ττ τ τ τ +      = − − + −     + +      Khai triển S theo chuỗi Taylor ta được: ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) 1 ˆ ln 1 2 2 0 0 2 ˆ ˆ2 /2 /2 0 0 0 1 2 1 ˆ ˆ ˆ ! ! 1 2 1 1 1 1 ˆ ˆ ˆ ˆ 1 2 ! ! 1 2! 1 2 1 2 ˆ ˆ i j i N j i j i i j i ii j N N i i j i j S M e M i j M M M M i ji S S ττ τ τ τ τ ττ τ + ∞ ∞ − + + = = + ∞ ∞ ∞ + + = = = ≠ −  =  +  − −    = +   + +   + + = + ∑∑ ∑ ∑∑ Khi đó ta có thể tách Hamiltonian thành hai thành phần: ( ) ( ) ( ) 2 0 ˆ2 / 2 00 2 1 1 ˆ ˆ ˆ ˆ 4 1 2! 1 2 i i i N i dH N Z M M i ω ω τ τ pi ττ τ +∞ ∞ + =   = −  +  + ∑∫ , ( ) ( ) ( ) ( )( ) ˆ / 20 00 12 1 ˆ ˆ ˆ ˆ ˆ 4 ! ! 1 2 1 2 i j i j i j N i j i j dV M M Z M M i j ω ω τ τ pi ττ τ + ++∞ ∞ ∞ + + = = ≠ −   = − + −  +  + ∑∑∫ Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 63 Phụ lục 9: Chuẩn hóa bộ hàm sóng cơ sở cho bài toán exciton hai chiều Trước hết, ta chọn bộ hàm sóng của dao động tử điều hòa (vì hàm này chắc chắn là nghiệm riêng của các toán tử trung hòa nên sẽ là nghiệm riêng của 0ˆH ) ( ) ( ) ( )ˆˆ, 0 ,yx x y nn x y n n x yn n C a b ω ω + + = , trong đó ,x yn n là các số nguyên dương và 0 là trạng thái chân không được định nghĩa: ( ) ( )ˆˆ 0 , 0, 0 , 0x y x ya bω ω ω ω= = ; và điều kiện chuẩn hóa là 0 0 1= . Như vậy nghiệm riêng của phương trình Schrödinger ta sẽ viết dưới dạng tổ hợp tuyến tính của các vector sóng trên: ( ) ( ) ( ) , , ˆ ˆ, 0 ,yx x y x y x y x y nn n n x y n n x y n n n n C n n C a bψ ω ω+ += =∑ ∑ . Nhận xét: tổng số mũ của hai toán tử ˆˆ ,a b+ + là x yn n+ . (*) Mặt khác do toán tử ˆzL là đại lượng bảo toàn nên hàm riêng của phương trình Schrödinger phải đồng thời là hàm riêng của toán tử này: ˆ zL mψ ψ= , với ( )ˆ ˆˆ ˆ ˆzL i b a a b+ += − . Ta có: ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) , 1 11 1 , ˆ ˆ ˆˆ ˆ ˆ ˆ 0 , ˆ ˆ ˆ ˆ 0 , . yx x y y yx x x y nn z x y n n n nn n x y x y n n L Ci b a a b a b Ci n a b n a b ψ ω ω ω ω + + + + + − − ++ + + + = −   = −    ∑ ∑ Nhận xét: tổng số mũ của hai toán tử ˆˆ ,a b+ + vẫn là x yn n+ .(**) Như vậy, từ hai nhận xét (*) và (**), kết hợp với công thức khai triển nhị thức Newton, ta có thể chọn dạng của hàm sóng cơ sở như sau: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 64 ( )22 )ˆ ˆˆ ˆ, [( ) ( ] 0mkkmk m C a b a ib+ + + += + ± . Xét: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 1 1 2 2 3 1 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ = , , , ... , ˆˆ ˆ ˆ ˆ ˆ, , , k k k k k k k k k k k a M M a a M M a a M M M a M M a a M M M a M M M a M M M a M a M a a a b − − + + + + + + + + − − − − + + + + + + + + + + + + + +     = + = + +                + + + + +            = +     ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 1 1 1 ˆ2 ˆ ˆ , 0 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 ... 2 ˆ ˆ ˆ ˆ 2 . k k k k k k k a a M a M M a a M M a M a k M a + + + − − + + + + + + − + + +   =     =    → = + + + = + Tính các toán tử sau: ( ) ( ) ( ){ } ( ) ( ) ( ) ( ) ( ) 12 1 1 2 22 ˆ ˆ ˆ ˆ ˆ ˆ ˆ2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ4 2 4 ( 1) ; k k k k k k k a M a M a k M a M a k M a a k M k k M a − + + + + − − − + + + + + + = + = + + + − ( ) ( ) ( ) ( ) ( ) ( )21 1 22 2ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ4 2 4 ( 1)k k k k kb M M b k M b b k M k k M b− − −+ + + + + + += + + + − ; ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 21 2 22 2 1 1 ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ4 1 4 ( 1) ˆ ˆ ˆ ˆ ˆ2 4 ( 1) ; k k k k k k k M M M a b k M a a b b k k M a b M M k M N k k M − − + + + + + + + + − − + + +  → = + + + + + − +   = + + − ( ) ( ) ( ) ( ) ( ) ( )1ˆ ˆ ˆˆ ˆ ˆˆ ˆ ˆ2 ;k k ka ib M M a ib k M a ib−+ + + + +± = ± + ± ( ) ( ) ( ) ( )1 2ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ2k k ka a M M a a k M a−+ + + + + += + ; ( ) ( ) ( ) ( )21ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ2k k kb b M M b b k M b−+ + + + + += + ; ( ) ( )( ) ( ) ( ) ( ) ( ) ( )ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ2 1 2 1 4 4k k k k k kN M a a b b M M a a b b k M M N k M+ + + + + + + + + += + + = + + + = + Tính Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 65 ( ) ( ) ( ) ( ) ( ) ( )1ˆ ˆ ˆˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ, 0 0 2 0m m mk k kkm km kma k m C a M a ib C M a a ib kC M a a ib−+ + + + + + + + + += ± = ± + ± ( ) ( ) ( ) ( )ˆ ˆˆ ˆˆ ˆ ˆ ˆ ˆ, 0 0m mk kkm kma k m C a M a ib C M a a ib+ + + + + + + + += ± = ± ( ) ( ) ( ) ( )1ˆ ˆˆ ˆ ˆ ˆˆ ˆ, 0 0m mk kkm kmM k m C M M a ib C M a ib++ + + + + + + += ± = ± ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 1 0 1, ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, 0 2 0 4 ( 1) 0 , 4 1, . m m mk k k km km km km k m M k m C M M a ib kC M N a ib k k C M a ib Ck k m k m C − − + + + + + + + + + = − = ± + ± + − ± = + −  Tính : ( ) 1 1, , , ˆ ˆ , 4 1, 4 ( 1)( 2)...( 1)( )( 1)...( 1) , !( )! 4 , .( )!( )! j jkm k m j km k j m j km k j m CM k m k k m M k m C Ck k k k j m k m k m k j k j m C k m k C k j m k j m k j C − − − − = + − = − − − + + + − + − + − + = − − + − ( ) ( ) ( )( ) ( ) ( )( ){ } ( ) 1 , 1 ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ, 2 0 2 , 1 m mk k km km k m a ib k m C M a ib a ib k M a ib a ib Ck m k m C − + + + + + + + + − = ± + ± = + − ∓ ∓ ∓ Tính ( ) ( ) ( ) ( ) ( ) 1 , 1 , ˆ ˆ ˆ ˆ, 2 , 1 ... ! 2 , . ! l l km k m l km k m l C a ib k m k m a ib k m C k m C k m l Ck m l − − − = + − = + = − + − ∓ ∓ Từ điều kiện chuẩn hóa ta thu được như sau: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 66 ( ) ( ) ( )  2 0, 2 0, 0, 1 2 2 ˆ ˆ ˆ, , 0 , !( )! ˆ ˆ4 0 0, ! !( )! 4 2 ! 0 0 ! 2 !( )! 1 m k km mk km m mk km m m k m km m k k m C a ib M k m k m k C a ib m m C k m k C m C m C k m k C = + = + = + = = + = ∓ ∓ , 1 2 2 !( )!k m mk C k m k = + . Vậy ta thu được hàm sóng sau khi chuẩn hóa có dạng: ( ) ( )1 ˆˆ ˆ, 0 2 2 !( )! mk mk k m M a ib k m k + + + = ± + . Tính các tác dụng của các toán tử lên hàm sóng này: ( ) ( )( ) ( ) ˆ , 2 1, ; ˆ , 2 1 1 1, ; ˆ , 2 2 1 , . M k m k k m k m M k m k k m k m N k m k m k m + = + − = + + + + = + + Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 67 Phụ lục 10: Các thành phần ma trận cho bài toán exciton hai chiều Từ điều kiện chuẩn hoá ta thu được bộ hàm cơ sở cho nguyên tử hydro như sau: ( ) ( )1 ˆˆ ˆ, 0 2 2 !( )! mk mk k m M a ib k m k + + + = ± + . ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( )( ) ( ) ˆ , 2 2 1 , ; ˆ , 2 1, ; ( )! ! ˆ , 2 , ;( )! ! ˆ , 2 1 1 1, ; ! ! ˆ , 2 , . ! ! j j i i N k m k m k m M k m k k m k m k k m M k m k j m k j k m j M k m k k m k m k i k m i M k m k i m k k m + + = + + = + − + = − − + − = + + + + + + + = + + * Tính thành phần ma trận của toán tử ˆS Do toán tử ˆS đựơc tách làm hai thành phần, ta tiến hành tính lần lượt các thành phần ma trận của toán tử ˆS là 1 2ˆ ˆS ,S như sau: ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) 2 1 ˆ2 /2 0 2 2 2 2 1 0 2 2 2 1 0 1 1 ˆ ˆ ˆ , , 1 2! 1 2 ( )! ! ! !1 12 2 , 1 2 ( )! ! ! !! 1 2 ( )! !1 12 ( )! !! 1 2 i i i N i ik i i k i m i k i k m i S k m M M k m i k k m k k m k m k i k m i k i k m ii k k m k i k m ii τ τ τ τ τ τ τ τ ∞ + = − + + = + + =   −  =   +  +   + + −  =   + − + − − + −  +   +  = −   − + − +  ∑ ∑ ∑ ,k m Thành phần 2ˆS : Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 68 ( ) ( )( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 ˆ /2 0 0 2 2 1 0 0 1 1 ˆ ˆ ˆ , , ! ! 1 2 1 2 ( )! ! ! !1 12 2 , ! ! 1 2 ( )! ! ! !1 2 ( )! ! ! !1 2 ! ! ( i j i j N i j i j i jk j i k j m i j i j i j S k m M M k m i j k k m k i j k m i j k i j m i j k j k m j k j k m j k k m k i j k m i j i j k j τ τ τ τ τ τ τ + ∞ ∞ + = = ≠ + ∞ − + + = = ≠ + −  =  +  + + + − + + − −  + − + − + − − + −  + + + − + + − = − − = ∑∑ ∑∑ ( ) ( )( )( ) 2 10 0 1 ,)! ! 1 2 k k j i m i j i j k i j m k m j τ ∞ − + + + = = ≠ + − + − + ∑∑ Khi đó ta tính ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( )( ) 2 1 ˆ2 /2 0 2 2 2 1 0 1 2 2 2 1 0 1 1 ˆ ˆ ˆ , , , , 1 2! 1 2 ( )! !2 1 , ,( )! !! 1 2 ( )! !2 1 ( )! !! 1 2 i i i N i ik k m i ik k m i m k S k m m k s M M k m i k k m m k k m k i k m ii k k m k i k m ii τ τ τ τ τ τ τ ∞ + = + + = = + + =   −  = +   +  +  + − = − + − + + − = − + − + ∑ ∑ ∑  Tính ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( )( ) 2 ˆ /2 0 0 1 4 4 20 0 2 1 1 ˆ ˆ ˆ , , , , ! ! 1 2 1 2 ( )! ! ! !1 2 1 , , ! ! 1 2 ( )! ! 1 2 1 2 ! ! 1 2 i j i j N i j i j i jk s k k j mi j i j m k s S k m m k s M M k m i j k k m k i j k m i j m k s k i j m i j k j k m j i j τ τ τ τ τ τ τ τ + ∞ ∞ + = = ≠ ++ − + + = = ≠   −  + = +   +  +    + + − + + − −  = + + − + − + −  + − = + ∑∑ ∑∑ ( ) ( ) ( ) ( ) ( )( )( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ,2 2 1 0 0 2 2 1 ( )! ! ! ! 1 ( )! ! 1 2 ( )! ! ! !2 1 !( )! ( )! ! 1 2 i jk s k s i jk j m i j i j i sk s k s m i s k k m k i j k m i j k j k m j k k m k s k m s i i s k s i k m s i δ τ τ τ ++ − − + + = = ≠ − + + + + = + + − + + −   − + −  + + + + + − = − + − + + − + ∑∑ ∑ * Tính các phần tử ma trận: Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 69 ( ) ( )( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) , 0 1 0 2 2 2 1 00 2 2 2 1 0 0 2 ˆˆ ˆ , , , , 4 ( )! !22 12 2 1 4 ( )! !! 1 2 ( )! ! 22 12 1 2 2 2 ( )! !! 1 2 km km ik k i i k m i k m dH m k H k m m k N Z S k m k k mdk m Z k i k m ii k k m k m Z d k i k m ii ω ω τ pi τ τω ω τ pi τ τ τω ω τ pi τ +∞ +∞ + + = +∞ + + = = = − + − = + + − − + − + + − = + + − − + − + ∫ ∑∫ ∑ ∫ Đặt 2 2 d t dt ττ τ = → = ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) 2 , 2 2 1 0 0 2 2 12 0 ( )! ! 212 1 2 2 ( )! !! 1 2 ( )! !1 = 2 1 2 ( )! !! i k k k m i k k i k m i k k m H k m Z dt k i k m ii k k m k m Z I k i k m ii τω ω pi τ ω ω +∞ + + = + + = + − = + + − − + − + + + + − − + − ∑ ∫ ∑ với ( ) ( ) ( ) ( )2 2 1 00 , 1 (2 2 3)!! 2 1 !!2 (1 ) 2 ( 1)! q q p p p q p q Z t p q qq dtp t pI pi pi +∞ − > ≥ ∈ − − − − − = = + −∫ ( ) ( ) ( ) ( ) ( )( )( ) ( ) ( ) ( ) ( ) ( )( )( ) 2, , , , 0 , 1 , 1 ,2 1 0 00 2 ˆˆ ˆ ˆ , , , , 4 2 2 1 1 4 ( )! ! ! !2 1 2 1 ! ! 1 2 ( )! ! 1 2 kmkm k s m k s m k s k k s k i jk k s kk s m i j i j dH H m k s V k m m k s M M Z S k m k k m k k m k k m k s k m sdZ i j k s i k m s i ω ω τ pi τ ω δ δ ω τ τ δ pi ττ τ +∞ + + + + − + + ++∞ ∞ + ++ + + = = ≠ = = + = + − + − = − + + + + + + + + + −  −  + + − + + −  + ∫ ∑ ∑∫ ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) , 1 ,1 22 2 120 , 1 ,1 2 2 1 1 1 2 ( )! ! ! !1 2 !( )! ( )! ! 1 1 1 2 ( )! ! ! !1 ! ! ( )! ! i j s s i s k s k s m i s s s k m s k k m k k m k k m k s k m s t Z i i s k s i k m s i t k k m k k m k k m k s k m s Z I i i s k s i k m s i ω δ δ ω pi ω δ δ ω − − − +∞+ + + + = − + + + = − + + + + + + + + + − − − + − + + − + = − + + + + + + + + + − − + − + + − ∑ ∫ k s i s i s + − = ∑ ( ) ( ) ( ) ( ) ( ) ( ) ( )1 1 1, , , , , , , , ,k m k s m k s m k m k m k s mH H H k k s− − += = = − Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 70 ( )( ) ( ) ( ) ( ) ( ) ( ) 1 2 1 , 1 2 2 1 ( )! ! 1 ! 1 !11 1 2 ! 1 ! ( 1 )! 1 ! k i k k k m i k k m k k m H k k m Z I i i k i k m i ω ω + − + + + = + + + + = − + + + − − + − + + − ∑ ( ) ( ) ( ) ( ) ( ) ( )1 2 , 2 1 ( )! ! ! !1 ! ! ( )! !s k s i s k k s k m s i s k k m k s k m s H Z I i i s k s i k m s i ω > + − + + + + = + + + + = − − + − + + − ∑ * Xác định ( )f ω : Ta có: ( ) ( ) ( ) ( ) ( ) 0 2 2 12 0 ( )! !1 = 2 1 2 ( )! !! i k kk k m i k k k m H k m Z I k i k m ii ω ε ω + + = + = = + + − − + − ∑ ( ) ( ) ( ) ( ) ( ) 0 2 2 12 0 ( )! !1 10 2 1 2 ( )! !2 ! in k m i k k k mZk m I k i k m ii ε ω ω + + = +∂ = = + + − ∂ − + −∑ ( ) ( ) ( ) ( ) 2 2 2 12 0 !!1 .(2 1) ( )! !! i k k m i m kkZ I k m k i m k ii ω + + =  + =   + + − + −   ∑ Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 71 Phụ lục 11: Hàm Gamma Hàm Gamma được định nghĩa bằng đẳng thức tích phân sau đây: ( ) ( )1 0 0xe x dtαα α ∞ − −Γ = >∫ Tính tích phân từng phần, ta có: ( ) ( ) 1 100 0 0 01 x x x x xe x dx x d e x e e x dx e x dxα α α α αα α α∞ ∞ ∞ ∞− − − ∞ − − − −Γ + = − = − + =∫ ∫ ∫ ∫ hay ( ) ( )1α α αΓ + = Γ Khi 1α = và 1 2 α = , hàm ( )αΓ có giá trị: ( ) 11 1, 2 pi  Γ = Γ =    Từ đó ta xác định được giá trị của hàm ( )αΓ đối với các α nguyên và bán nguyên như sau: ( ) ( )1 !n nΓ = − ( )2 1 !!1 2 2n n n pi − Γ + =    trong đó n=1,2,3... Đối với các giá trị khác của α , hàm ( )αΓ có thể tìm trong các bảng riêng. Tích phân có dạng ( ) ( ) ( ) ( ) 2 2 00 , 1 11 2 2 2 (1 ) q q p p q p q Z p q qtq dtp t p I pi pi +∞ > ≥ ∈     − Γ − − Γ +    −     = = + Γ∫ Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 72 với 1 1 1 (2 2 3)!!1 2 2 2 p q p qp q p q pi − − − −   Γ − − = Γ − − + =        , ( )2 1 !!1 2 2q q q pi − Γ + =    , ( ) ( 1)!p pΓ = − khi đó tích phân có dạng: ( ) ( ) ( ) ( )2 2 1 00 , 1 (2 2 3)!! 2 1 !!2 (1 ) 2 ( 1)! q q p p p q p q Z t p q qq dtp t pI pi pi +∞ − > ≥ ∈ − − − − − = = + −∫ , với , 1, 2,3...p q = Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 73 Phụ lục 12: Lập trình fortran cho mức năng lượng của exciton hai chiều c Calculate Exciton2D excited energy PROGRAM MAIN integer i,m,k double precision w,Hmatrix * From the mininum of energy -> omega=Pi w=0.0021 m=6 k=0 CALL MAINSUB(w,k,m) END * MAIN subroutine calculate approximated energy SUBROUTINE MAINSUB(w,k,m) INTEGER Z,i,j,s,L,m,k DOUBLE PRECISION w,Hmatrix,E,C,H,tuso,mauso,temp,msum PARAMETER (smax=100,kmax=101) * Chu y thay kmax=smax+k DIMENSION E(0:smax),C(0:kmax,0:smax),H(0:kmax,0:kmax) Z=1 * Initialize matrix Hmatrix DO i=0,smax+k write(*,*) i DO j=0,i H(i,j)= Hmatrix(w,Z,m,i,j) H(j,i)=H(i,j) ENDDO ENDDO WRITE(*,*) 'Hmatrix done!' * Initialize C coefficient matrix DO i=0,smax DO j=0,smax+k C(j,i)=0.0 ENDDO C(k,i)=1.0 ENDDO WRITE(*,*) 'C matrix done!' Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 74 * Initialize E(s) DO i=0,smax E(i)=0.0 ENDDO E(0)=H(k,k) * Calculate E(s): Energy in s th approximation OPEN(10,FILE='Energy.dat',STATUS='Unknown') WRITE(10,*) 's=0','E=',E(0) DO s=1,smax * Calculate C(L,s) DO L=0,s+k IF (L.NE.k) THEN tuso=H(L,k) mauso=E(s-1)-H(L,L) DO i=0,k+s IF (i.NE.k.AND.i.NE.L) THEN tuso=tuso+C(i,s-1)*H(L,i) ENDIF ENDDO IF (mauso.EQ.(0.0)) STOP 'Error, division by zero' C(L,s)=tuso/mauso ENDIF ENDDO * Calculate E(s) from C(L,s) msum=0.0 DO L=0,k+s IF (L.NE.k) msum=msum+C(L,s)*H(L,k) ENDDO E(s)=H(0,0)+msum WRITE(10,*) 's=',s,'E=',E(s) ENDDO RETURN END * function calculate Hmatrix C Calculate elements of H matrix H(omega,Z,m,row,col) FUNCTION Hmatrix(w,Z,m,row,col) INTEGER Z,row,col,r,c,s,k,p,q,I,m1 DOUBLE PRECISION Hmatrix,w,H,msum,coeff1,coeff2 PARAMETER (pi=3.14159265358979323846) Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 75 IF (row<0.OR.col<0) STOP 'Error! Wrong matrix indexes!' * Note: H(row,col)=H(col,row) IF (row.LT.col) THEN r=col c=row ELSE r=row c=col ENDIF s=r-c k=c m1=ABS(m) IF (s.EQ.0) THEN * Calculate Hkk msum= 0.0 p=2*k+m1+1 DO i=0,k msum=msum+coeff1(k,i,k+m1,i,p,2*i) ENDDO H=w/2.0*p-Z*SQRT(w*pi)*msum ELSEIF (s.EQ.1) THEN * Calculate Hk,k+1 msum= 0.0 p=2*k+m1+2 DO i=1,k+1 msum=msum & +coeff2(k,k+1,k+1-i,m1)*coeff1(k,i-1,k+1,i,p,2*i-1) ENDDO H=(-1.0)*w/2.0*SQRT((k+1.0)*(k+m1+1.0))-Z*SQRT(w*pi)*msum ELSE * Calculate Hk,k+s,s>1 msum= 0.0 p=2*k+s+m1+1 DO i=s,k+s msum=msum & +coeff2(k,k+s,k+s-i,m1)*coeff1(k,i-s,k+s,i,p,2*i-s) ENDDO H=(-1.0)*Z*SQRT(w*pi)*msum ENDIF Hmatrix=H RETURN Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 76 END function calculate coefficient 1, 2 c function coeff1 FUNCTION coeff1(n1,k1,n2,k2,p,q) INTEGER n1,k1,n2,k2,p,q INTEGER i1,i2,i3,i4,i5,ii1,ii2,ii3,ii4,ii5 DOUBLE PRECISION coeff1,temp,kq,maxnum,minnum maxnum=1.79769D308 minnum=3.00000D-308 kq=1.0 ii1=1 ii2=1 ii3=1 ii4=1 ii5=1 11 IF ((ii1.LE.k1.OR.ii2.LE.k2.OR.ii3.LE.(p-q-1).OR.ii4.LE.q) & .OR.ii5.LE.(p-1)) THEN DO i1=ii1,k1 temp=kq*(n1-i1+1.0)/DBLE(i1) IF (temp.GT.maxnum.OR.temp.LT.minnum) EXIT kq=temp ENDDO ii1=i1 DO i2=ii2,k2 temp=kq*(n2-i2+1.0)/DBLE(i2) IF (temp.GT.maxnum.OR.temp.LT.minnum) EXIT kq=temp ENDDO ii2=i2 DO i3=ii3,(p-q-1) temp=kq*(2.0*i3-1.0) IF (temp.GT.maxnum.OR.temp.LT.minnum) EXIT kq=temp ENDDO ii3=i3 DO i4=ii4,q temp=kq*(2.0*i4-1.0) IF (temp.GT.maxnum.OR.temp.LT.minnum) EXIT kq=temp ENDDO ii4=i4 Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 77 DO i5=ii5,p-1 temp=kq/(2.0*i5) IF (temp.GT.maxnum.OR.temp.LT.minnum) EXIT kq=temp ENDDO ii5=i5 GOTO 11 ENDIF IF (q.GT.0.AND.MOD(q,2).NE.0) kq=kq*(-1.0) coeff1=kq RETURN END c function coeff2 FUNCTION coeff2(k1,k2,k3,m) INTEGER i,k1,k2,k3,m DOUBLE PRECISION coeff2,temp temp=1.0 DO i=1,m temp=temp*SQRT(DBLE((k1+i)*(k2+i)))/DBLE(k3+i) ENDDO coeff2=temp END Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 78 TÀI LIỆU THAM KHẢO Tiếng Việt [1]. Hoàng Dũng (1999), Nhập môn cơ học lượng tử- Tập I, Nhà xuất bản Giáo dục [2]. Vũ Văn Hùng (2004), Cơ học lượng tử, Nhà xuất bản Giáo dục. [3]. Lê Văn Hoàng (2004), Phương pháp đại số giải phương trình Schrödinger cho nguyên tử Hydro trong từ trường với cường độ bất kỳ, Đề tài Khoa học công nghệ cấp cơ sở CS.2004.23.59 [4]. Lê Văn Hoàng (2005), Phổ năng lượng trạng thái exciton của khí điện tử hai chiều tạo ra do hệ nhiều lớp GaAs/GaAsAl trong từ trường đều, Đề tài KHCN cấp bộ B2005.23.72. [5]. Đặng Quang Khang (2006), Cơ học lượng tử, Nhà xuất bản Khoa Học Kĩ Thuật. [6]. Nguyễn Hữu Mình (2007), Bài tập cơ học lượng tử tâp II, Nhà xuất bản Giáo dục. [7]. Hoàng Đỗ Ngọc Trầm (2008), Phương pháp toán tử giải phương trình Schodinger cho exciton hai chiều trong từ trường đều với cường độ bất kỳ, Luận văn Thạc sĩ. Khoa Vật lý trường Đại học KHTN Tp. Hồ Chí Minh Tiếng Anh [8]. Le Van Hoang, Hoang Do Ngoc Tram, Lu Thanh Trung (2005), Analytical Solution of 2D Exciton in a Magnetic Field, Communications in Physics, Supplement 2005, p.101-106 [9]. S. H. Patil (2008) The helium atom and isoelectronic ions in two dimensions, Eur. J. Phys.29(2008)517–525. Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 79 MỤC LỤC MỞ ĐẦU ............................................................................................................................... 1 Chương 1: GIỚI THIỆU PHƯƠNG PHÁP TOÁN TỬ (OM) QUA BÀI TOÁN DAO ĐỘNG TỬ PHI ĐIỀU HÒA .................................................................. 5 1.1 Sơ đồ Rayleigh- Schrödinger cho phương pháp nhiễu loạn dừng ............... 5 1.2 Phương pháp nhiễu loạn và dao động tử phi điều hòa.................................. 8 1.3 Phương pháp toán tử cho bài toán dao động tử phi điều hòa .................... 10 Chương 2: EXCITON - BÀI TOÁN EXCITON HAI CHIỀU ...................................... 17 2.1 Exciton .............................................................................................................. 17 2.1.1 Khái niệm exciton ............................................................................................ 17 2.1.2 Phân loại exciton .............................................................................................. 17 2.1.3 Tính chất của exciton....................................................................................... 18 2.2 Bài toán exciton hai chiều ............................................................................... 19 2.2.1 Phương trình Schrödinger cho exciton hai chiều ......................................... 19 2.2.2 Phương pháp giải tích cho bài toán exciton hai chiều.................................. 20 Chương 3: PHƯƠNG PHÁP TOÁN TỬ CHO BÀI TOÁN EXCITON HAI CHIỀU..................................................................................................... 25 3.1 Phương trình Schrödinger cho exciton hai chiều biểu diễn qua toán tử sinh hủy.......................................................................................... 25 3.2 Phương pháp toán tử giải bài toán exciton hai chiều ..................................... 28 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI ......................................................... 36 PHỤ LỤC............................................................................................................................ 37 Phụ lục 1: Các toán tử sinh – hủy một chiều........................................................... 37 Phụ lục 2: Dạng chuẩn (normal) của một số toán tử trong luận văn.................... 41 Phụ lục 3: Yếu tố ma trận của toán tử Hamilton của dao động tử Luận văn tốt nghiệp GVHD: Th.S Hoàng Đỗ Ngọc Trầm  2010 SVTH: Trương Mạnh Tuấn Trang 80 phi điều hòa............................................................................................ 44 Phụ lục 4: Phương trình Schrödinger cho bài toán exciton hai chiều. ................ 47 Phụ lục 5: Hamilton cho bài toán exciton hai chiều ............................................... 50 Phụ lục 6: Các toán tử sinh – hủy hai chiều ............................................................ 54 Phụ lục 7 : Dạng chuẩn của toán tử { }ˆ ˆ ˆ ˆexp ( )S M M Nτ += − + + ....................... 58 Phụ lục 8: Chuyển toán tử Hamilton qua biểu diễn toán tử sinh –hủy ............... 61 Phụ lục 9: Chuẩn hóa bộ hàm sóng cơ sở cho bài toán exciton hai chiều ............ 62 Phụ lục 10: Các thành phần ma trận cho bài toán exciton hai chiều ................... 66 Phụ lục 11: Hàm Gamma.......................................................................................... 70 Phụ lục 12: Lập trình Fortran cho năng lượng exciton hai chiều......................... 72 TÀI LIỆU THAM KHẢO ................................................................................................. 77

Các file đính kèm theo tài liệu này:

  • pdfluanvan-31102311.pdf