›

Khóa luận Đường tán sắc của exciton-Polariton hai chiều trong tương tác với phonon âm học

pdf 32 trang thiennha21 16/04/2022 4560
Download
Bạn đang xem 20 trang mẫu của tài liệu "Khóa luận Đường tán sắc của exciton-Polariton hai chiều trong tương tác với phonon âm học", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfkhoa_luan_duong_tan_sac_cua_exciton_polariton_hai_chieu_tron.pdf

Nội dung text: Khóa luận Đường tán sắc của exciton-Polariton hai chiều trong tương tác với phonon âm học

  1. BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TPHCM KHOA VẬT LÝ TRẦN DƯƠNG ANH TÀI ĐƯỜNG TÁN SẮC CỦA EXCITON-POLARITON HAI CHIỀU TRONG TƯƠNG TÁC VỚI PHONON ÂM HỌC KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC CHUYÊN NGÀNH: SƯ PHẠM VẬT LÝ TP. HỒ CHÍ MINH – 2018
  2. BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TPHCM KHOA VẬT LÝ TRẦN DƯƠNG ANH TÀI ĐƯỜNG TÁN SẮC CỦA EXCITON-POLARITON HAI CHIỀU TRONG TƯƠNG TÁC VỚI PHONON ÂM HỌC KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC CHUYÊN NGÀNH: SƯ PHẠM VẬT LÝ MÃ NGÀNH: 102 NGƯỜI HƯỚNG DẪN KHOA HỌC TS. PHẠM NGUYỄN THÀNH VINH TP. HỒ CHÍ MINH – 2018
  3. Lời cảm ơn Đầu tiên, tôi xin dành lời cảm ơn chân thành và sâu sắc nhất đến thầy hướng dẫn khoa học của tôi, TS. Phạm Nguyễn Thành Vinh. Trong quá trình học tập tại khoa Vật Lí, Trường Đại học Sư Phạm TPHCM, tôi may mắn được làm việc với thầy khi còn là một sinh viên năm nhất. Thầy đã kiên trì hướng dẫn và tận tình giúp đỡ khi tôi vừa bắt đầu thực hiện đề tài nghiên cứu đầu tiên, một điều hoàn toàn mới mẻ với một sinh viên năm nhất khi đó. Thầy không chỉ dạy tôi những kiến thức Vật Lí và kĩ năng cần thiết cho công việc nghiên cứu, trong quá trình làm việc dưới sự hướng dẫn của thầy, thầy còn dạy tôi nhiều bài học quý giá trong cuộc sống và luôn tạo điều kiện để tôi có thể phát triển bản thân một cách tốt nhất. Những bài học bổ ích ấy đã giúp tôi gặt hái được nhiều thành tích và có những trải nghiệm đáng nhớ trong suốt bốn năm đại học. Ngoài ra, tôi cũng học tập ở thầy về thái độ làm việc nghiêm túc, cách làm hiệu quả, và một số kĩ năng mềm. Suốt quãng thời gian thực hiện đề tài khoá luận tốt nghiệp, thầy luôn động viên, khích lệ tinh thần, giúp tôi vượt qua những khó khăn để hoàn thành khoá luận tốt nghiệp. Khoá luận tốt nghiệp này có thể sẽ không hoàn chỉnh nếu thiếu những nhận xét, góp ý của TS. Nguyễn Duy Vỹ, Viện Vật Liệu Tiên Tiến, Trường Đại học Tôn Đức Thắng và TS. Tomotake Yamakoshi, Viện Khoa học LASER, Trường Đại học Điện Tử–Viễn Thông (Institute for Laser Science, University of Electro– Communications). Những nhận xét phản biện này không chỉ góp phần đảm bảo tính chính xác về mặt khoa học cho khoá luận tốt nghiệp của tôi mà còn giúp tôi hiểu rõ hơn về bức tranh Vật Lí của đề tài mà tôi đang thực hiện. Ngoài ra, trong quá trình thảo luận với TS. Tomotake Yamakoshi, tôi học hỏi thêm về kĩ thuật lập trình với ngôn ngữ FORTRAN 77 và một số phương pháp toán lý mới. Tôi xin gửi lời cảm ơn đến các thầy cô trong khoa Vật Lí, Trường Đại học Sư Phạm TPHCM, những người đã tận tình giảng dạy, truyền đạt những kiến thức, i
  4. và kinh nghiệm quý giá trong bốn năm qua để tôi có thể hoàn thành khoá luận tốt nghiệp này và có được hành trang tốt nhất cho công việc trong tương lai của tôi. Ngoài ra, tôi cũng xin gửi lời cảm ơn riêng đến TS. Phan Thị Ngọc Loan, người đã dạy tôi học phần “Phương pháp nghiên cứu khoa học”, những bài giảng của cô đã tạo cho tôi cảm hứng với việc nghiên cứu Vật Lí và thầy cố vấn học tập, TS. Hoàng Văn Hưng, nhờ những buổi nói chuyện với thầy, tôi học hỏi thêm được nhiều điều bổ ích bên cạnh giải toả áp lực trong học tập và nghiên cứu. Tôi cũng cảm ơn những thành viên trong nhóm nghiên cứu của TS. Phạm Nguyễn Thành Vinh. Trong quá trình làm việc, tôi luôn nhận được sự hỗ trợ tận tình và động viên kịp thời từ các thành viên trong nhóm. Cùng với các thành viên trong nhóm, tôi có những hành trình đáng nhớ, đặc biệt là chuyến đi tham quan Vũng Tàu năm 2018 cùng với TS. Tomotake Yamakoshi. Trong bốn năm học tập tại trường Đại học Sư Phạm TPHCM, tôi may mắn được quen biết nhiều bạn bè cùng khoá và các anh chị khoá trên, những người luôn bên cạnh và giúp đỡ tôi những lúc tôi gặp những vấn đề khó giải quyết. Tôi trân trọng khoảng thời gian ôn tập cho những kì thi kết thúc học phần căng thẳng cùng với các bạn Hồ Hoàng Huy, Nguyễn Tấn Phú, Nguyễn Thành Nhân, Trương Ngô Bích Trâm. Tôi muốn gửi lời cảm ơn đến anh Trần Công Hiếu vì đã giúp đỡ tôi trong kì thi tuyển sinh đại học năm 2014, hỗ trợ tôi hoàn tất thủ tục nhập học và cung cấp tài liệu những học phần đại cương dành cho sinh viên năm nhất. Tôi cũng xin cảm ơn chị Hoàng Khánh Linh, chị Nguyễn Mai Khanh đã lắng nghe và cho tôi những lời khuyên để tôi vượt qua nhiều khó khăn trong lúc hoàn thành khoá luận tốt nghiệp này. Tôi sẽ không quên những lời khuyên về cách học Vật Lí và kinh nghiệm nghiên cứu được chia sẻ từ CN. Lê Đại Nam. Cuối cùng nhưng không kém phần quan trọng, tôi xin cảm ơn ba mẹ của tôi. Ba và mẹ tôi đã tạo mọi điều kiện để tôi có thể tập trung vào việc học tập suốt bốn năm qua và luôn ủng hộ những quyết định của tôi. Tôi không thể thành công như ngày hôm nay nếu không có sự hi sinh của ba, mẹ tôi. Tp.HCM, ngày 02 tháng 05, năm 2018 Sinh viên Trần Dương Anh Tài ii
  5. Mục lục Trang Danh sách hình vẽ ii Danh mục chữ viết tắt iii Mở đầu1 1 Cơ sở lý thuyết4 1.1 Toán tử sinh và huỷ trong cơ học lượng tử 4 1.2 Giả hạt Polariton 7 2 Định lý Floquet 12 3 Kết quả và thảo luận 15 3.1 Hướng tiếp cận sử dụng phương pháp gần đúng sóng quay . . . . . . 15 3.2 Hướng tiếp cận sử dụng định lý Floquet 19 Kết luận và hướng phát triển 22 Tài liệu tham khảo 23 i
  6. Danh sách hình vẽ Trang Hình 1.1: Hình vẽ thể hiện sự phụ thuộc năng lượng của LP và UP theo độ lệch năng lượng giữa exciton và photon. . . . . . . . 11 Hình 1.2: Hình vẽ thể hiện đường tán sắc năng lượng của LP và UP và sự phụ thuộc vào vector sóng song song của các hệ số Hopfield tương ứng với các trường hợp a) ∆ = 2g0, b) ∆ = 0, c) ∆ = −2g0 11 ii
  7. Danh mục chữ viết tắt Chữ viết tắt Tiếng Việt Tiếng Anh BEC Ngưng tụ Bose – Einstein Bose–Einstein Condensation HHG Sóng điều hoà bậc cao High Harmonic Generator LA Sóng âm học dọc Longtitudinal Acoustic LP Polariton nhánh dưới Lower Polariton Light Amplification by Stimulated LASER La–de Emisson Radiation RWA Phép gần đúng sóng quay Rotating Wave Approximation SAW Sóng âm học bề mặt Surface Acoustic Wave TA Sóng âm học ngang Tranverse Acoustic Phương trình Schrodinger¨ TDSE Time–Dependent Schrodinger¨ Equation phụ thuộc thời gian UP Polariton nhánh trên Upper Polariton iii
  8. Mở đầu Từ khi ra đời vào năm 1960, LASER (viết tắt của cụm từ Light Amplification by Stimulated Emisson Radiation) là một công cụ đắc lực giúp các nhà vật lý nghiên cứu cấu trúc của nguyên tử, phân tử thông qua các hiệu ứng phi tuyến như phát xạ sóng điều hòa bậc cao (HHG – High Harmonic Generation ) [1], quá trình ion hóa nguyên tử, phân tử [2,3] hay được dùng trong bẫy từ-quang (MOTs) [4, 5] để bẫy các nguyên tử cho các nghiên cứu sự biến đổi trạng thái của vật chất ở pha ngưng tụ Bose–Einstein (BEC) [6–8] từ đó giúp chúng ta hiểu thêm về thế giới tự nhiên. Với những tính chất đặc biệt như tính đơn sắc, kết hợp và có cường độ cao, LASER được ứng dụng rộng rãi trong nhiều lĩnh vực khác trong cuộc sống. Tuy nhiên, các LASER được dùng trong các phòng thí nghiệm hiện nay đều được tạo ra bằng cách tạo môi trường đảo mật độ sao cho electron trong các nguyên tử chủ yếu ở trạng thái kích thích, khi các electron này trở về trạng thái cơ bản, các photon phát ra phản xạ nhiều lần qua hệ cộng hưởng quang học và tạo thành LASER. Phương pháp này đòi hỏi phải tạo ra môi trường đảo mật độ, và giữ các electron ở trạng thái kích thích đủ lâu để có thể phát ra LASER, việc này tương đối khó khăn vì thời gian sống của electron ở trạng thái kích thích ngắn, vào cỡ 10−8s. Ngoài ra, ngưỡng năng lượng để xảy ra sự phát xạ này tương đối lớn, chúng ta phải cung cấp nhiều năng lượng để quá trình có thể xảy ra, do đó việc này gây tốn kém. Năm 1996, A. Imamoglu và cộng sự đã đưa khái niệm về một loại LASER hoàn toàn mới mà không cần đến môi trường đảo mật độ [9]. Các tác giả đã sử dụng các giả hạt polariton là sự kết hợp giữa exciton (cặp electron và lỗ trống) và photon trong cấu trúc tinh thể của chất bán dẫn được cấu hình sẵn. Các hạt polariton có spin nguyên do đó chúng có thể có cùng trạng thái lượng tử đơn (single quantum state), như các hạt boson trong pha BEC, . . . và phát ra những photon kết hợp và đơn sắc, đây chính là cơ sở để tạo nên polariton LASER. Quá trình 1
  9. phát xạ polariton LASER xảy ra ở nhiệt độ thấp khoảng 4K, được quan sát lần đầu tiên bởi L. S. Dang và các cộng sự vào năm 1998 [10]. Đến năm 2007, LASER polariton với bơm quang học lần đầu tiên được tạo ra ở nhiệt độ phòng [11]. Tuy Polartion LASER có nhiều ưu điểm như ngưỡng phát xạ thấp, không cần đến môi trường đảo mật độ, tần số của LASER có thể kiểm soát một cách dễ dàng bằng việc thay đổi tính chất của giếng lượng tử và vật liệu bán dẫn nhưng do công suất phát xạ ở nhiệt độ phòng vẫn còn rất nhỏ [11], nên chưa thể ứng dụng vào thực tế. Chúng tôi nhận thấy rằng việc tạo ra polariton LASER ở nhiệt độ phòng có cường độ cao có ý nghĩa vô cùng to lớn. Nhằm thực hiện điều này, chúng tôi thêm vào hệ các phonon âm học thông qua sóng âm học bề mặt (SAWs - Surface Ascoustic Waves). Khi các hạt polariton tương tác với các phonon âm học, đường tán sắc năng lượng bị thay đổi [12], do đó quá trình ngưng tụ BEC của các hạt polariton có thể bị thay đổi từ đó làm tăng nhiệt độ chuyển pha và cường độ của polariton LASER. Trong công trình [12], Ivanov và các cộng sự đã đưa ra phương trình đường tán sắc năng lượng của các hạt polariton khi có mặt sóng âm học (phương trình (4) [12]). Tuy nhiên, các tác giả chỉ trình bày kết quả cuối cùng mà thiếu đi quy trình toán học chặt chẽ để đưa ra kết quả này. Do đó, việc tìm ra một quy trình toán học phù hợp, chi tiết để dẫn dắt đến kết quả của Ivanov và cộng sự [12] là vô cùng cần thiết cho việc thực hiện các nghiên cứu tiếp theo của chúng tôi. Với những nhận xét nêu trên, chúng tôi thực hiện đề tài “Đường tán sắc của các exciton-polariton hai chiều trong tương tác với phonon âm học” cho khóa luận tốt nghiệp này nhằm đưa ra quy trình toán học chặt chẽ để đưa ra lại phương trình đường tán sắc năng lượng được nêu trong [12]. Khoá luận tốt nghiệp được trình bày thành ba chương, nội dung từng chương như sau: • Chương 1: Những tìm hiểu về hình thức luận lượng tử hoá lần hai trong cơ học lượng tử, quá trình phát xạ polariton LASER và sóng âm học bề mặt được trình bày trong chương này. • Chương 2: Chúng tôi trình bày các phương pháp tính toán được sử dụng để chéo hoá Hamiltonian mô tả sự tương tác giữa polariton và phonon âm học, cụ thể là định lý Floquet. 2
  10. • Chương 3: Các kết quả của khoá luận tốt nghiệp được trình bày trong chương này. Kết quả tính toán cho thấy khi ta sử dụng phương pháp gần đúng sóng quay để chéo hoá Hamiltonian thì kết quả thu được hoàn toàn khác với kết quả được đưa ra bởi Ivanov và các cộng sự [12]. Do đó chúng tôi sử dụng hướng tiếp cận khác, đó là sử dụng định lý Floquet kết hợp với phương trình Schr¨odingerphụ thuộc thời gian, để kiểm tra lại các kết quả trước đó của chúng tôi. Chúng tôi lại thu được kết quả hoàn toán khác hai kết quả trước đó. Đồng thời, chúng tôi cũng phát hiện ra rằng kết quả sau khi chéo hoá Hamiltonian trong công trình của Ivanov năm 2003 [12] hoàn toàn khác với kết quả của chính tác giả nó được công bố vào năm 2001 [13]. 3
  11. Chương 1 Cơ sở lý thuyết 1.1 Toán tử sinh và huỷ trong cơ học lượng tử Trong cơ học lượng tử, ngoài cách tiếp cận theo hướng giải tích, chúng ta có thể tiếp cận theo hướng đại số, sử dụng toán tử sinh và huỷ. Hướng tiếp cận này giúp ta tiết kiệm thời gian tính toán, thuận tiện trong việc tính toán các hệ nhiều hạt do đó nó cũng thường được sử dụng nhiều bởi các nhà vật lý. Trong phần này, chúng tôi trình bày tóm tắt về các tính toán với toán tử sinh và huỷ trong cơ học lượng tử. Đặt a và a† là hai toán tử tác động lên những trạng thái trong không gian Hilbert, và thoả mãn giao hoán tử [a, a†] = 1 (1.1) trong đó “1” kí hiệu cho toán tử đơn vị trong không gian Hilbert. Toán tử a và a† là các toán tử không tự liên hợp, hay nói cách khác các toán tử này không có tính chất Hermitic. Toán tử a† được gọi là toán tử sinh và a được gọi là toán tử huỷ. Ta gọi |αi là trạng thái được chọn sao cho vector riêng của toán tử Hermitic, a†a, có trị riêng là số thực α a†a|αi = α|αi. (1.2) Do đó α = hα|a†a|αi = ||a|αi||2 ≥ 0, (1.3) trong đó chúng tôi sử dụng tiên đề cơ bản của cơ học lượng tử, rằng “norm” của tất cả các trạng thái trong không gian Hilbert đều dương. Kết quả là trị riêng α của trạng thái riêng của a†a là một số thực không âm. 4
  12. Thêm vào đó, với các toán tử A, B, C, ta luôn có [AB, C] = A[B, C] + [A, C]B, (1.4) từ đây ra suy các giao hoán tử quan trọng của các toán tử sinh và huỷ † [ai ai, aj] = −ajδi,j, (1.5) † † † [ai ai, aj] = ajδi,j, (1.6) trong đó δi,j là Kronecker delta. Một trạng thái bất kỳ |ni được biễu diễn thông qua trạng thái cơ bản |0i, trạng thái bị huỷ bởi toán tử huỷ a|0i = 0, (1.7) qua biểu thức sau 1 |ni = √ (a†)n|0i, (1.8) n! và nội tích của nó thoả hm|ni = n!δm,n. (1.9) Tóm lại, các toán tử sinh và huỷ phải tuân theo các phương trình sau √ a†|ni = n + 1|n + 1i (1.10) √ a|ni = n|n − 1i (1.11) a†a|ni = n|ni (1.12) và do đó, các yếu tố ma trận lần lượt là √ hm|a|ni = nδm,n−1 (1.13) √ † hm|a |ni = n + 1δm,n+1 (1.14) Để minh hoạ việc vận dụng toán tử sinh và huỷ vào các bài toán vật lý lượng tử, chúng tôi trình bày lời giải bài toán dao động tử điều hoà một chiều bằng cách sử dụng toán tử sinh huỷ. Toán tử Hamilton mô tả một dao động tử điều hoà có dạng như sau pˆ2 1 Hˆ = + mω2xˆ2 (1.15) 2m 2 5
  13. d trong đó pˆ = −i là toán tử động lượng, xˆ là toán tử toạ độ, m và ω lần lượt là ~dx khối lượng và tần số của dao động tử điều hoà. Ta định nghĩa các toán tử sinh và toán tử huỷ cho dao động tử điều hoà như sau 1 rmω  d  a† = √ (mωx − ip) = x − ~ , (1.16) 2mω~ 2~ mω dx 1 rmω  d  a = √ (mωx + ip) = x + ~ . (1.17) 2mω~ 2~ mω dx Với định nghĩa toán tử sinh và huỷ như trên, các toán tử toạ độ và động lượng được biểu diễn thông qua các toán tử sinh và huỷ như sau r xˆ = ~ a + a† , (1.18) 2mω r mω pˆ = ~ a† − a i. (1.19) 2 Khi này, phương trình (1.15) được viết lại thành  1 Hˆ = ω a†a + . (1.20) ~ 2 Với Hamiltonian biểu diễn theo các toán tử sinh và huỷ, trị riêng năng lượng của dao động tử điều hoà được tìm phương trình Schr¨odingerdừng Hˆ |ni = En|ni  1 ⇔ ω a†a + |ni = E |ni ~ 2 n  1 ⇔ ω n + |ni = E |ni (1.21) ~ 2 n từ phương trình (1.21), ta suy ra trị riêng năng lượng của dao động tử điều hoà  1 E = ω n + , (1.22) n ~ 2 kết quả này tương tự như kết quả thu được khi tính theo phương pháp giải tích. Để đưa ra hàm sóng của dao động tử điều hoà, ta cần định nghĩa trạng thái chân không trước, trạng thái chân không của dao động tử điều hoà được định nghĩa bởi a|0i = 0, (1.23) hay  d  x + ~ ψ (x) = 0. (1.24) mω dx 0 6
  14. Phương trình (3.1) có nghiệm  mωx2  ψ0(x) = A exp − , (1.25) 2~ với A là hệ số chuẩn hoá được xác định từ điều kiện chuẩn hoá ∞ Z 2 |ψ0(x)| dx = 1. (1.26) −∞ Hàm sóng chuẩn hoá trạng thái chân không của dao động tử điều hoà có dạng r  2  4 mω mωx ψ0(x) = exp − . (1.27) π~ 2~ Một trạng thái bất kỳ |ni của dao động tử điều hoà được xây dựng bằng cách tác động toán tử sinh lên hàm sóng trạng thái chân không r  n r  2  1 mω ~ d 4 mω mωx ψn = √ x − exp − . (1.28) n! 2~ mω dx π~ 2~ 1.2 Giả hạt Polariton Polariton là một giả hạt trong chất rắn, được tạo ra bởi sự tương tác của ánh sáng (photon) và vật chất. Hàm sóng của polariton là sự chồng chập lượng tử từ hàm sóng mô tả photon và exciton, một sự kết cặp của một electron và một lỗ trống trong chất bán dẫn |ψi = X|ψxi + C|ψci, (1.29) trong đó ψc kí hiệu cho photon và ψx kí hiệu exciton. Các hệ số X, C được gọi là hệ số Hopfield đặc trưng cho tính chất của polariton. Chúng tôi giải thích cặn kẽ hệ số Hopfield ở bên dưới. Tiếp theo, chúng tôi trình bày chi tiết cách xây dựng đường tán sắc năng lượng của các hạt polariton. Hamiltonian đặc trưng cho tương tác giữa photon và exciton trong vi hốc quang học (optical microcavity) là tổng của Hamiltonian mô tả photon trong vi hốc Hˆcav, Hamiltonian biểu diễn cho exciton Hˆexc và thành phần Hamilotonian thể ˆ hiện sự kết cặp (coupling) giữa photon và exciton Hcoupling. Theo hình thức luận lượng tử hoá lần hai, Hamiltonian này có dạng như sau ˆ ˆ ˆ ˆ Hpol = Hcav + Hexc + Hcoupling X  c † xˆ†ˆ †ˆ ˆ†  = ~ωpaˆ aˆ + ~ωp b b + g0(ˆa b +a ˆb ) (1.30) p 7
  15. với aˆ† là toán tử sinh photon trong vi hốc quang học, ˆb† là toán tử sinh exciton, x c ωp và ωp lần lượt là tần số của photon và exciton, g0 là độ lớn lưỡng cực tương tác giữa exciton và photon và có độ lớn khác “không” với những mode có cùng vector sóng song song p. Lưu ý, ở đây chúng tôi dùng kí hiệu p để mô tả vector sóng song song, một số tài liệu khác có thể dùng kí hiệu pk hoặc kk. Hamiltonian trên có thể được chéo hoá thông qua phép biến đổi tuyến tính sau của các toán tử huỷ " # " #" # Pˆ XC ˆb = (1.31) Qˆ −CX aˆ trong đó các toán tử Pˆ và Qˆ là các toán tử huỷ của một giả hạt (quasiparticle) mới, giả hạt này được đặt tên là polariton và X, C là các hệ số đặc trưng cho phép biến đổi tuyến tính. Từ định nghĩa các toán tử huỷ Pˆ và Qˆ, ta suy ra định nghĩa các toán tử sinh Pˆ† và Qˆ† như sau " # " #" # Pˆ† X∗ C∗ ˆb† = . (1.32) Qˆ† −C∗ X∗ aˆ† Xét giao hoán tử [P,ˆ Pˆ†], ta có [P,ˆ Pˆ†] = [Xˆb + Ca,ˆ X∗ˆb† + C∗aˆ†], (1.33) sử dụng các tính chất của giao hoán tử của các toán tử sinh và huỷ cho bởi phương trình (1.1), phương trình (1.33) được viết lại thành [P,ˆ Pˆ†] = |X|2[ˆb, ˆb†] + |C|2[ˆa, aˆ†] = |X|2 + |C|2 (1.34) mặt khác, ta có [P,ˆ Pˆ†] = 1 nên ta suy ra được |X|2 + |C|2 = 1. (1.35) Ta cũng thu được kết quả tương tự nếu xét giao hoán tử [Q,ˆ Qˆ†]. Để biểu diễn Hamiltonian trong phương trình (1.30) theo các toán tử mới, ta biễu diễn các toán tử a,ˆ ˆb theo P,ˆ Qˆ " # " #−1 " # " #" # ˆb XC Pˆ CX Pˆ = = , (1.36) aˆ −CX Qˆ X −C Qˆ và aˆ†, ˆb† theo các toán tử Pˆ†, Qˆ† " # " #−1 " # " #" # ˆb† X∗ C∗ Pˆ C∗ X∗ Pˆ† = = . (1.37) aˆ† −C∗ X∗ Qˆ X∗ −C∗ Qˆ† 8
  16. Thay các phương trình (1.36) và (1.37) vào phương trình (1.30), và thu gọn ta được ˆ X  2 c 2 x ˆ† ˆ 2 c 2 x ˆ† ˆ Hpol = (|C| ~ωp + |X| ~ωp + 2g0XC)P P + (|X| ~ωp + |C| ~ωp − 2XCg0)Q Q p  c x 2 2  ˆ† ˆ ˆ† ˆ + XC(~ωp − ~ωp ) + g0(|X| − |C| ) (P Q + Q P ) (1.38) Theo công trình [14], Hamiltonian sau khi chéo hoá sẽ có dạng như sau ˆ X ˆ† ˆ ˆ† ˆ Hpol = ~ωLP P P + ~ωUP Q Q , (1.39) p trong đó tần số ωLP đặc trưng cho nhánh năng lượng dưới (lower eigenenergy) tương ứng với polariton dưới (lower polariton) trong khi tần số ωUP đặc trưng cho nhánh năng lượng trên (upper eigenenergy) tương ứng với polariton trên (upper polariton), do đó, hệ số của biểu thức thứ ba trong phương trình (1.38) phải thoả mãn c x 2 2 XC(~ωp − ~ωp ) + g0(|X| − |C| ) = 0. (1.40) Toán tử Pˆ†Qˆ hay Qˆ†Pˆ thể hiện việc thay đổi trạng thái của giả hạt thu được từ việc chéo hoá Hamiltonian cho bởi phương trình (1.30), trong các thí nghiệm, hiện tượng này chỉ xảy ra khi ta kích thích hệ rất “mạnh” và rất khó thực hiện, do đó trong những tính toán lý thuyết, số hạng này thường được bỏ qua, ta xem như không có sự nhảy trạng thái này. Điều đó giải thích ý nghĩa của phương trình (1.40). Phương trình (1.39) cho ta thấy rằng ta có thể biểu diễn Hamiltonian mô tả tương tác giữa các photon và exciton thông qua các toán tử là tổ hợp tuyến tính của các toán tử mô tả photon và exciton. Do đó, các hạt polariton là sự chồng chất của các hạt photon trong vi hốc và các hạt exciton trong cùng mặt phẳng vector sóng. Ta có thể thấy rằng các hệ số X, C trong phép biến đổi tuyến tính đặc trưng cho tỉ lệ pha trộn giữa photon và exciton. Bằng cách thay đổi tỷ lệ “trộn” photon và exciton, người ta có thể điều khiển tương tác giữa ánh sáng và vật chất. Trong thực nghiệm, điều này có thể là thực hiện được bằng cách dịch chuyển năng lượng của vi hốc (bơm thêm một số loại khí) hoặc exciton (bằng cách thay đổi nhiệt độ hoặc áp một điện trường ngoài). Từ phương trình (1.40), ta tìm được giá trị bình 9
  17. phương module của các hệ số X, C ! 1 ∆ |X|2 = 1 + , (1.41) 2 p 2 2 ∆ + 4g0 ! 1 ∆ |C|2 = 1 − , (1.42) 2 p 2 2 ∆ + 4g0 x c 2 2 với ∆ = ~ωp − ~ωp là độ lệch năng lượng giữa exciton và photon. Khi |X| < |C| , điều này có nghĩa photon trong vi hốc đóng góp nhiều hơn các exciton trong việc hình thành giả hạt polariton do đó, ta gọi các polariton này là photon-polariton và ngược lại, ta sẽ gọi chúng là exciton-polariton. Trong giới hạn của khoá luận tốt nghiệp này, chúng tôi chỉ khảo sát các exciton-polartion. Khi ∆ = 0, tương ứng 1 với các giá trị |X|2 = |C|2 = , lúc này các photon và exciton đóng góp bằng nhau 2 để tạo thành các polariton. Đồng nhất hệ số phương trình (1.38) và (1.39) ta thu được hệ thức biểu diễn đường tán sắc của các hạt polariton 1  q  E = ω = ωc + ωx + 4g2 + ( ωx − ωc)2 , (1.43) LP ~ LP 2 ~ p ~ p 0 ~ p ~ p 1  q  E = ω = ωc + ωx − 4g2 + ( ωx − ωc)2 . (1.44) UP ~ UP 2 ~ p ~ p 0 ~ p ~ p Các phương trình (1.43) và (1.44) là phương trình mô tả đường cong tán sắc của các hạt polariton tương ứng với mức năng lượng nhánh dưới (LP - Lower Polariton) và nhánh trên (UP – Upper Polariton). Khi các photon và exciton liên kết với x c nhau tại tần số cộng hưởng, ωp = ωp, năng lượng của LP và UP sai khác nhau một x c lượng có giá trị nhỏ nhất, ~ωp − ~ωp = 2g0, sự phân tách này tương tự như sự phân tách Rabi trong hệ hai mức năng lượng. Do đó, độ lớn lưỡng cực tương tác giữa exciton và photon thường được đặc trưng bởi tần số Rabi, Ωc, theo công thức [12] i Ω g = ~ c . (1.45) 0 2 Do sự kết cặp giữa exciton và photon, năng lượng của các polariton có xu hướng chống lại sự dịch chuyển năng lượng từ photon sang năng lượng exciton. Đây là một trong những dấu hiệu thể hiện sự kết cặp mạnh giữa photon và exciton và khi x c |~ωp − ~ωp|  g0, ta không thể phân biệt exciton và photon một cách rõ ràng. Điều này được thể hiện rõ trong hình 1.1. 10
  18. Hình 1.1: Hình vẽ thể hiện sự phụ thuộc năng lượng của LP và UP theo độ lệch năng lượng giữa exciton và photon [14]. (Lưu ý: kí hiệu Ω trong hình vẽ tương ứng là tần số Rabi Ωc mà chúng tôi trình bày ở trên.) Hình 1.2: Hình vẽ thể hiện đường tán sắc năng lượng của LP và UP và sự phụ thuộc vào vector sóng song song của các hệ số Hopfield tương ứng với các trường hợp a) ∆ = 2g0, b) ∆ = 0, c) ∆ = −2g0[14].(Lưu ý: kí hiệu kk trong hình vẽ tương ứng là vector sóng song song p mà chúng tôi trình bày ở trên.) 11
  19. Chương 2 Định lý Floquet Năm 1883, nhà toán học Floquet đề xuất cách giải phương trình vi phân có dạng dx = A(t)x, (2.1) dt với A(t) là một hàm số liên tục theo biến số t và tuần hoàn với chu kì T , A(t) = A(t + T ). (2.2) Ngày nay, phương pháp này được gọi là định lý Floquet. Từ phương trình (2.1), ta có thể thấy rằng nó tương đương về mặt toán học với phương trình Schr¨odingerphụ thuộc thời gian (TDSE) có Hamiltonian tuần hoàn theo thời gian với chu kì T H(t) = H(t + T ). (2.3) Do đó, chúng tôi sử dụng định lý Floquet để giải quyết bài toán được đặt ra trong khoá luận tốt nghiệp này theo đề xuất của TS. Tomotake Yamokoshi. Ngoài ra, định lý Floquet là công cụ toán học mạnh để giải quyết các bài toán trong các nghiên cứu về hệ lượng tử tuần hoàn do đảm bảo tính tuần hoàn của sự nhiễu loạn của tất cả các mức gần đúng và tránh được những biểu thức phụ thuộc tuyến tính hoặc không tuần hoàn theo thời gian. Ở đây, chúng tôi không trình bày chi tiết về định lý Floquet mà chỉ đưa ra cách áp dụng vào việc giải TDSE có Hamiltonian tuần hoàn theo thời gian, những tìm hiểu sâu hơn về định lý Floquet có thể tham khảo các tài liệu [15, 16]. Xét phương trình Schr¨odingerphụ thuộc thời gian (để đơn giản chúng tôi chỉ xét hệ một chiều không gian)  ∂  H(x, t) − i Ψ(x, t) = 0, (2.4) ~∂t 12
  20. trong đó H(x, t) = H0(x) + V (x, t). (2.5) Thành phần Hamiltonian không phụ thuộc thời gian H0(x) tương ứng với hàm riêng ψn(x) và trị riêng En và V (x, t) là thế năng tuần hoàn theo thời gian với chu kì T V (x, t) = V (x, t + T ). (2.6) Theo định lý Floquet, nghiệm của phương trình (2.4) có dạng như sau Ψ(x, t) = exp(−iεαt/~)Φα(x, t), (2.7) trong đó Φα(x, t), Floquet mode, tuần hoàn theo thời gian cùng chu kỳ với V (x, t) Φα(x, t) = Φα(x, t + T ). (2.8) Ở đây, εα là một thông số thực đặc trưng cho thành phần mũ. Do εα bằng bội số của ~ω với ω = 2π/T nên được gọi là giả năng lượng (quasienergy) tương tự như giả động lượng (quasimomentum) k, đặc trưng cho hàm sóng Bloch trong chất rắn. Thay phương trình (2.7), vào phương trình (2.4), ta suy ra phương trình để tìm trị riêng εα HF (x, t)Φα(x, t) = εαΦα(x, t). (2.9)  ∂  với H (x, t) = H(x, t) − i có tính chất Hermitic và được gọi là Hamiltonian F ~∂t Floquet. Từ phương trình (2.9), ta thấy rằng các Floquet mode |Φα0 (x, t)i = exp(inωt)|Φα(x, t)i ≡ Φαn(x, t), (2.10) có nghiệm tương tự như nghiệm của phương trình (2.4) với dịch chuyển trong giả năng lượng bởi phương trình εα = εα0 + n~ω = εαn, (2.11) với n ∈ Z. Điều này cho ta thấy rằng kí hiệu α đang được sử dụng liên quan đến các lớp nghiệm tương ứng với α0 = (α, n) với n ∈ Z. Do đó nên ta có thể rút gọn ω ω một giả năng lượng bất kỳ về vùng Brillouin thứ nhất, −~ < ε < ~ . Các vector 2 2 13
  21. riêng của toán tử HF (x, t) phải tuân theo điều kiện trực chuẩn trong không gian Hilbert, do đó T +∞ Z Z 1 ∗ hΦ 0 (t)|Φ (t)i = dt Φ 0 (t)Φ (t)dx = δ 0 0 = δ δ , (2.12) α β T α β α β α,β n,m 0 −∞ và tạo thành bộ đủ X X ∗ 0 Φαn(x, t)Φαn(y, t) = δ(x − y)δ(t − t ). (2.13) α n 14
  22. Chương 3 Kết quả và thảo luận 3.1 Hướng tiếp cận sử dụng phương pháp gần đúng sóng quay Như đã đề cập ở trên, chúng tôi kích thích hệ photon và exciton bằng sóng âm ac học bề mặt (SAWs – Surface Acoustic Waves) có vector sóng k, tần số Ωk = vsk, với vs là tốc độ truyền âm và cường độ Iac. Trong trường hợp này, Hamiltonian tổng cộng của hệ được cho bởi [12] ˆ ˆ ˆ ˆ ˆ H = Hcav + Hexc + Hcoupling + Hpump X i Ωc ac † ac =  ωca†a + ωxb†b + ~ (a†b − b†a ) + imx(b†b e−iΩk t − b b eiΩk t), ~ p p p ~ p p p 2 p p p p k p p−k p−k p p (3.1) q x dp/pe ph trong đó mk = mx−ac (k) N0 (k) là tham số liên kết giữa exciton và sóng âm học dp 1/2 được bơm vào hệ, với mx−ac = Dx[~k/(2ρvs)] là yếu tố ma trận thế năng biến pe 3/2 dạng (deformation potential matrix element), và mx−ac ∝ e14k là yếu tố ma trận áp điện (piezoelectron matrix element) lần lượt đặc trưng cho sóng âm học dọc (LA – Longtitudinal Acoustic) và sóng âm học ngang (TA – Tranverse Acoustic). Với, ρ là khối lượng riêng, Dx là thế năng biến dạng của tương tác giữa exciton và LA–phonon, e14 là tensor đặc trưng cho tương tác giữa TA–phonon và exciton. Số ph lượng phonon kết hợp được xác định bởi N0 = Iac/(~vsΩkac). Để biến đổi Hamiltonian phụ thuộc thời gian trong phương trình (3.1) thành Hamiltonian không phụ thuộc thời gian, chúng tôi sử dụng phương pháp gần đúng sóng quay – RWA theo đề xuất của Ivanov trong công trình [12]. Trong hệ qui 15
  23.  X † †  chiếu quay, thông qua phép biến đổi chính tắc S = exp it (vs · p)(bpbp + apap) , p ac 2 với vs = vsk/k do đó vs · p = vsk · p/k = Ωk (k · p)/k , trong đó vs và p lần lượt là tốc độ truyền âm và vector sóng song song mặt phẳng, Hamiltonian không phụ thuộc thời gian được suy ra từ Hamiltonian phụ thuộc thời gian bởi phương trình ∂S† Hˆ = SHSˆ † − iS . (3.2) TI ∂t Các toán tử trong phương trình (3.1) được biến đổi bằng cách áp dụng công thức Baker–Hausdorff (Baker–Hausdorff identity) 1 1 exp(A)B exp(−A) = B + [A, B] + [A, [A, B]] + + [A, [A, A[A, . . .[A, B] ]. (3.3) 2! n! | {z } n và các tính chất của giao toán tử của các toán tử sinh huỷ cho bởi phương trình (1.1), từ đó, ta suy ra các biểu thức sau −it(vs·q) bq → bqe , −it(vs·q) aq → aqe , † † † † aqbq − bqaq → aqbq − bqaq, (3.4) † † ac ac bqbq−k → bq, bq−kexp[it(Ωq − Ωq-k)], † † ac ac bq−kbq →, bq−kbqexp[it(Ωq-k − Ωq )], và ta có X i Ωc SHSˆ † =  ωxb†b + ωca†a + ~ (a†b − b†a )+ ~ p p p ~ p p p 2 p p p p p ac ac ac ac ac ac x † −iΩk t it(Ωp −Ωp−k) † iΩk t it(Ωp-k−Ωp )  + imk(bpbp−ke e − bp−kbpe e ) X h i Ωc = ωxb†b + ωca†a + ~ (a†b − b†a )+ ~ p p p ~ p p p 2 p p p p p ac ac ac ac ac ac i x † −i(Ωk −Ωp +Ωp−k)t † i(Ωk −Ωp +Ωp-k)t + imk(bpbp−ke − bp−kbpe ) . (3.5) † P † † P † † ∂S it (vs·p)(bpbp+apap) −it (vs·p)(bpbp+apap) X −iS = − ie p e p (−i) (v · p)(b†b + a†a ) ∂t s p p p p p X † † = − (vs · p)(bpbp + apap). (3.6) p 16
  24. Thay phương trình (3.5) và (3.6) vào phương trình (3.2), ta có  X i Ωc Hˆ = ( ωc − v · p)b†b + ( ωc − v · p)a†a + ~ (a†b − b†a )+ TI ~ p s p p ~ p s p p 2 p p p p p ac ac ac ac ac ac i x † −i(Ωk −Ωp +Ωp−k)t † i(Ωk −Ωp +Ωp-k)t +imk(bpbp−ke − bp−kbpe ) . (3.7) Theo định luật bảo toàn năng lượng, ta có ac ac ac Ωp−k + Ωk = Ωp , (3.8) từ đó, suy ra được Hamilonian không phụ thuộc thời gian có dạng được cho bởi biểu thức ˆ X  x † c † HTI = (~ωp − vs · p)bpbp + (~ωp − vs · p)apap p i Ω + ~ c (a†b − b†a ) + imx(b†b − b† b ) (3.9) 2 p p p p k p p−k p−k p Từ phương trình (3.9), ta thấy rằng năng lượng của photon và exciton bị dịch chuyển một lượng vs · p, để đơn giản ta đặt các biến số năng lượng mới cho exciton và photon x x ~ω˜p = ~ωp − vs · p, (3.10) c c ~ω˜p = ~ωp − vs · p, phương trình (3.9) được viết các biến số mới như sau X i Ωc Hˆ =  ω˜xb†b + ω˜ca†a + ~ (a†b − b†a ) + imx(b†b − b† b ). (3.11) TI ~ p p p ~ p p p 2 p p p p k p p−k p−k p p Một cách tương tự, chúng tôi chéo hoá HˆTI bằng phép biến đổi tuyến tính để suy ra phương trình đường tán sắc năng lượng của polariton khi có thêm phonon ˆ ˆ P † âm học. Hamiltonian HTI được chéo hoá thành Hc = ~ωcpcp bằng cách đặt p † ∗ † ∗ † cp = ubp + vap ⇒ cp = u bp + v ap. (3.12) Sử dụng các tính chất của các giao hoán tử [AB, C] = A[B, C] + [A, C]B và [A, BC] = B[A, C] + [A, B]C, chúng tôi tính toán các giao hoán tử sau X [cq, Hˆc/ ] = ωcpδpq = ωcq = ω(ubq + vaq), (3.13) ~ p ˆ ˆ [cq, HTI /~] =[ubq + vaq, HTI /~] = I0, (3.14) 17
  25. trong đó X iΩc X iΩc I =u [b , ω˜xb†b − b†a + imx(b†b − b† b )] + v [a , ω˜ca†a + a†b ] 0 q p p p 2 p p k p p−k p−k p q p p p 2 p p p p iΩc X iΩc =u[˜ωxb − a + imx(δ b − δ b )] + v(a ω˜c + b ) q q 2 q k pq p−k p−k,q p q q 2 q p iΩ iΩ =u[˜ωxb − c a + imx(b − b )] + v(˜ωca + c b ) q q 2 q k q−k q+k q q 2 q iΩ iΩ =(uω˜x + v c )b + (vω˜c − u c )a + uimx(b − b ). q 2 q q 2 q k q−k q+k Đồng nhất hệ số hai phương trình (3.13) và (3.14), ta có iΩ ωub =(uω˜c + v c )b + uimx(b − b ), (3.15) q q 2 q k q−k q+k iΩ ωv = − u c + vω˜c. (3.16) 2 q x x Với trường hợp k = 0 → mk = 0, mk ∝ k, ta có phương trình đường tán sắc thường của polariton  iΩc  " # " # ω − ω˜x − u 0 q 2  iΩ  = (3.17) c ω − ω˜c v 0 2 q 2 x Ωc ⇒ ω − ωq − c = 0. (3.18) 4(ω − ωq) Phương trình (3.18) cho nghiệm phù hợp với các tính toán tại mục 1.2, hai nghiệm phân biệt này tương ứng với phương trình tán sắc mô tả Polariton nhánh dưới (phương trình (1.43)) và nhánh trên (phương trình (1.44)). Tuy nhiên, phương trình chúng tôi đưa ra lại khác bậc so với phương trình mà Ivanov đã đưa ra vào năm 2003 [12]. Kết quả trong công trình năm 2003 như sau (ωΩ )2 4ω2|mx|2 ω2 − (˜ωx)2 − c − t k p ω2 − (˜ωc)2 2 p 2 x 2 (ωΩc) ω − (˜ωp+k) − 2 c 2 − Mp+2k ω − (˜ωp+k) 2 x 2 4ωt |mk| − 2 = 0, (3.19) 2 x 2 (ωΩc) ω − (˜ωp−k) − 2 c 2 − Mp−2k ω − (˜ωp−k) 2 x 2 4ωt |mk| với Mp±nk = 2 . Chúng tôi không thể giải 2 x 2 (ωΩc) /4 ω − (˜ωp±nk) − 2 c 2 − Mp±(n+1)k ω − (ωp±nk) thích được sự khác biệt này. Sau đó, chúng tôi phát hiện kết quả khác về phương 18
  26. trình đường tán sắc trong một công bố khác vào năm 2001 của chính tác giả Ivanov [13] 2 x Ωc/4 ω − ω˜p − c − Mp+k − Mp−k = 0, (3.20) ω − ω˜p x 2 |mk| trong đó Mp±nk = 2 , và kết quả này lại có bậc x Ωc/4 ω − ω˜p±nk − c − Mp±(n+1)k ω − ω˜p±nk phù hợp với kết quả tính toán của chúng tôi. Do đó, chúng tôi tiếp tục thực hiện các tính toán với hi vọng rằng có thể dẫn ra phương trình đường tán sắc trong tài liệu [13]. Tuy nhiên, chúng tôi không thể dẫn ra lại kết quả trong công trình [13] do gặp khó khăn trong việc xử lý hiệu số bp−k − bp+k. Do đó, chúng tôi đã tham khảo ý kiến từ TS. Tomotake Yamakoshi về bài toán này. TS. Tomotake Yamakoshi đề xuất dùng phương pháp khác dựa trên định lý Floquet để tiếp cận bài toán này do Hamiltonian tuần hoàn với chu kỳ T . Hướng tiếp cận dựa trên định lý Floquet chính là nội dung của phần sau. 3.2 Hướng tiếp cận sử dụng định lý Floquet Trong phần này, chúng tôi trình bày tính toán nhằm đưa ra phương trình đường tán sắc của polariton khi sóng âm học được đưa vào hệ. Chúng tôi dùng phương trình Schr¨odingerphụ thuộc thời gian để đưa ra phương trình đường tán sắc ∂Ψ i = Hˆ Ψ, (3.21) ~ ∂t để tiện theo dõi, chúng tôi viết lại Hamiltonian ˆ ˆ ˆ ˆ ˆ H = Hcav + Hexc + Hcoupling + Hpump X i Ωc ac † ac =  ωca†a + ωxb†b + ~ (a†b − b†a ) + imx(b†b e−iΩk t − b b eiΩk t). ~ p p p ~ p p p 2 p p p p k p p−k p−k p p (3.22) Theo định lý Floquet, hàm sóng của Ψ có dạng như sau +∞ ac −iωnt X −imΩk t n Ψ = e e φm, (3.23) m=−∞ 19
  27. n trong đó φm là hàm sóng mô tả polariton trên cơ sở |I1,I2, pi m X I φm = Cp |I1,I2, pi, (3.24) I,p với qui ước như sau exciton sẽ tương ứng với I = x, I1 = 1,I2 = 0, photon tương ứng với I = c, I1 = 0, I2 = 1. Các toán tử sinh và huỷ của exciton tác động lên cơ sở này theo qui tắc 0 0 bp I1,I2, p = δI11δI20δp p |0, 0, pi , (3.25) † 0 0 bp I1,I2, p = δI10δI21δp p |1, 0, pi , tương tự cho các toán tử sinh và huỷ của photon. Từ đó, ta dễ dàng suy ra các phương trình sau D 00 00 00 † 0 0 0E I ,I , p b bp I ,I , p = δp00pδ 00 δ 00 δp0pδ 0 δ 0 , 1 2 p 1 2 I1 1 I2 0 I11 I20 D 00 00 00 † 0 0 0E I ,I , p a ap I ,I , p = δp00pδ 00 δ 00 δp0pδ 0 δ 0 , 1 2 p 1 2 I1 0 I2 1 I10 I21 D 00 00 00 † 0 0 0E I1 ,I2 , p apbp I1,I2, p = δp00pδI00 0δI00 1δp0pδI0 1δI0 0, 1 2 1 2 (3.26) D 00 00 00 † 0 0 0E I ,I , p b ap I ,I , p = δp00pδ 00 δ 00 δp0pδ 0 δ 0 , 1 2 p 1 2 I1 1 I2 0 I10 I21 D 00 00 00 † 0 0 0E I1 ,I2 , p b bp I1,I2, p = δp00p−kδ 00 δ 00 δp0pδ 0 δ 0 , p−k I1 1 I2 0 I11 I20 D 00 00 00 † 0 0 0E I1 ,I2 , p bpbp−k I1,I2, p = δp00pδ 00 δ 00 δp0p−kδ 0 δ 0 . I1 1 I2 0 I11 I20 n Do đó, hàm sóng φm có dạng tường minh như sau n x c φm = Cp (n, m) |1, 0, pi + Cp(n, m) |0, 1, pi. (3.27) Tính đạo hàm phương trình (3.23), +∞ ∂Ψ X ac = e−iωnt (−imΩac − iω ).e−imΩk tφn (3.28) ∂t k n m m=−∞ và đưa vào phương trình Schr¨odingerphụ thuộc thời gian (3.21), ta có +∞ +∞ ac X ac −i(ωn+mΩk )t n X X  x † c † (~ωn + m~Ωk )e φm = ~ωp bpbp + ~ωpapap m=−∞ m=−∞ p  Ωc ac † ac ac +i~ (a†b − b†a ) + imx(b†b .e−iΩk t − b b eiΩk t) e−i(ωn+mΩk t)φn . (3.29) 2 p p p p k p p−k p−k p m 20
  28. 0 ac Nhân hai vế của phương trình (3.29) với h1, 0, p|e−i(ωn0 +m ~Ωk )t và lấy tích phân trên toàn miền thời gian, ta được phương trình i Ω ( ω + m Ωac)Cx(n; m) = ωxCx(n; m) − ~ c Cc(n; m) ~ n ~ k p ~ p p 2 p x  x x  + imk Cp−k(n; m − 1) − Cp+k(n; m + 1) , (3.30) 0 ac tương tự, khi nhân hai vế phương trình (3.29) với h0, 1, p|e−i(ωn0 +m ~Ωk )t và lấy tích phân trên toàn miền thời gian, ta được i Ω ( ω + m Ωac)Cc(n; m) = ωcCc(n; m) + ~ c Cx(n; m). (3.31) ~ n ~ k p ~ p p 2 p x c Từ phương trình (3.31), ta suy ra mối liên hệ giữa Cp (n, m) và Cp(n, m) c i~Ωc 1 x Cp(n, m) = ac c Cp (n, m). (3.32) 2 ~ωn + m~Ωk − ~ωp Kết hợp phương trình (3.30) và (3.32), ta suy ra hệ thức truy hồi của hệ số khai triển hàm sóng exciton  2  ac x (~Ωc) 1 x ~ωn + m~Ωk − ~ωp + ac c Cp (n, m) 4 ~ωn + m~Ωk − ~ωp x  x x  = imk Cp−k(n, m − 1) − Cp+k(n, m + 1) (3.33) Ta có thể thấy rằng, kết quả này khác biệt so với kết quả của Ivanov vào các năm 2001, 2003. Chúng tôi dự định thực hiện tính toán số để đánh giá các phương trình tán sắc của chúng tôi đưa ra và của Ivanov nhưng do thời gian thực hiện khoá luận có giới hạn nên chúng tôi chưa thực hiện được việc này. 21
  29. Kết luận và hướng phát triển Trong khoá luận tốt nghiệp này, chúng tôi thực hiện việc chéo hoá Hamiltonian mô tả tương tác giữa các polariton và các phonon âm học theo phương pháp RWA được đề xuất bởi Ivanov và các cộng sự vào năm 2003 để thu được phương trình đường tán sắc của các polariton khi có mặt các phonon âm học. Tuy nhiên, chúng tôi lại không thể thu được kết quả tương tự như phương trình (2) trong bài báo của tác giả Ivanov [12]. Để chắc chắn những tính toán của chúng tôi là chính xác, chúng tôi đề nghị TS. Nguyễn Duy Vỹ đưa ra những đánh giá về các kết quả chúng tôi. TS. Nguyễn Duy Vỹ kết luận rằng kết quả của chúng tôi là chính xác. Điều này đã gây khó khăn cho chúng tôi, và để giải quyết vấn đề này, chúng tôi nhờ đến sự giúp đỡ của TS. Tomotake Yamakoshi. Theo đề xuất của TS. Tomotake Yamakoshi, chúng tôi sử dụng định lý Floquet kết hợp với phương trình Schr¨odingerphụ thuộc thời gian và tính toán lại để dẫn ra phương trình đường tán sắc do phương pháp gần đúng sóng quay không phù hợp với quá trình tương tác đa photon (quan điểm của TS. Tomotake Yamakoshi). Với hướng tiếp cận này, chúng tôi thu được một phương trình mô tả đường tán sắc hoàn toàn khác so với kết quả trong công trình [12]. Sau đó, chúng tôi phát hiện ra một phương trình đường tán sắc khác cũng do chính tác giả Ivanov đưa ra vào năm 2001 [13]. Điều này đã làm chúng tôi rất bất ngờ, khi chéo hoá cùng một Hamiltonian lại cho nhiều kết quả khác nhau. Vì thời gian thực hiện khoá luận tốt nghiệp có giới hạn nên chúng tôi chưa thể đưa ra câu trả lời cuối cùng để khẳng định kết quả nào chính xác. Tuy nhiên, đây cũng là một kết quả thú vị đối với một khoá luận tốt nghiệp khi phát hiện ra sự “kì dị” này, đặc biệt là những kết quả này được công bố trên tạp chí hàng đầu trên thế giới trong lĩnh vực vật lý. Hướng phát triển tiếp theo của đề tài này là tham khảo ý kiến của một số chuyên gia để giải thích sự bất thường này. Ngoài ra, chúng tôi cũng sẽ tiến hành giải số để đánh giá các kết quả hiện tại. 22
  30. Tài liệu tham khảo [1] M. Lewenstein, P. Balcou, M. Y. Ivanov, et al., “Theory of high-harmonic generation by low-frequency laser fields”, Phys. Rev. A, vol. 49, no. 3, p. 2117, 1994. [2] V. Pham, O. Tolstikhin, and T. Morishita, “Molecular Siegert states in an electric field. II. Transverse momentum distribution of the ionized electrons”, Phys. Rev. A, vol. 89, no. 3, p. 033 426, 2014. [3] O. I. Tolstikhin, L. B. Madsen, and T. Morishita, “Weak-field asymptotic the- ory of tunneling ionization in many-electron atomic and molecular systems”, Phys. Rev. A, vol. 89, no. 1, p. 013 421, 2014. [4] E. Raab, M. Prentiss, A. Cable, et al., “Trapping of neutral sodium atoms with radiation pressure”, Phys. Rev. Lett., vol. 59, no. 23, p. 2631, 1987. [5] S. Chu, L. Hollberg, J. E. Bjorkholm, et al., “Three-dimensional viscous con- finement and cooling of atoms by resonance radiation pressure”, Phys. Rev. Lett., vol. 55, no. 1, p. 48, 1985. [6] M. Anderson, J. Ensher, M. Matthews, et al., “Observation of Bose-Einstein condensation in a dilute atomic vapor”, Science, vol. 269, no. 5221, pp. 198– 201, 1995. [7] T. Yamakoshi and S. Watanabe, “Wave-packet dynamics of noninteracting ultracold bosons in an amplitude-modulated parabolic optical lattice”, Phys. Rev. A, vol. 91, no. 6, p. 063 614, 2015. [8] S. Watanabe, S. Aizawa, and T. Yamakoshi, “Contrast oscillations of the Bose- Einstein-condensation-based atomic interferometer”, Phys. Rev. A, vol. 85, no. 4, p. 043 621, 2012. [9] A. Imamoglu, R. J. Ram, S. Pau, et al., “Nonequilibrium condensates and lasers without inversion: Exciton-polariton lasers”, Phys. Rev. A, vol. 53, no. 6, p. 4250, 1996. 23
  31. [10] L. S. Dang, D. Heger, R. André, et al., “Stimulation of polariton photolu- minescence in semiconductor microcavity”, Phys. Rev. Lett., vol. 81, no. 18, p. 3920, 1998. [11] S. Christopoulos, H. Von H. Baldassarri, A. J. D. Grundy, et al., “Room- temperature polariton lasing in semiconductor microcavities”, Phys. Rev. Lett., vol. 98, no. 12, p. 126 405, 2007. [12] A. L. Ivanov and P. B. Littlewood, “Resonant acousto-optics of microcavity polaritons”, Semicond. Sci. Technol., vol. 18, pp. 428–434, 2003. [13] A. L. Ivanov, “Acoustically Induced Stark Effect for Excitons in Intrinsic Semiconductors”, Phys. Rev. Lett., vol. 87, no. 13, p. 136 403, 2001. [14] D. Hui, H. Haug, and Y. Yamamoto, “Exciton-polariton Bose-Einstein con- densation”, Rev. Mod. Phys., vol. 82, pp. 1489–1537, 2010. [15] W. Magnus and S. Winkler, Hill’s equation. Courier Corporation, 2013. [16] E. L. Ince, “Ordinary Differential Equations Dover, New York, 1956”, First published by Longmans, Green and Co. in, 1926. 24
  32. Tp.HCM, ngày tháng 05, năm 2018 Xác nhận của Phản Biện Tp.HCM, ngày tháng 05, năm 2018 Xác nhận của Người Hướng Dẫn Khoa Học