Ảnh hưởng của trường bức xạ laser lên hấp thụ sóng điện từ yếu bởi điện tử giam cầm trong siêu mạng pha tạp tán xạ điện tử phonon âm

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN BÙI HỒNG PHƯỢNG ẢNH HƯỞNG CỦA TRƯỜNG BỨC XẠ LASER LÊN HẤP THỤ SÓNG ĐIỆN TỬ YẾU BỞI ĐIỆN TỬ GIAM CẦM TRONG SIÊU MẠNG PHA TẠP (TÁN XẠ ĐIỆN TỬ - PHONON ÂM) LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội- 2011 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Bùi Hồng Phượng ẢNH HƯỞNG CỦA TRƯỜNG BỨC XẠ LASER LÊN HẤP THỤ SÓNG ĐIỆN TỬ YẾU BỞI ĐIỆN TỬ GIAM CẦM TRONG SIÊU MẠNG PHA TẠP (TÁN XẠ ĐIỆN TỬ - PHONON ÂM) Chuyên ngành: Vật lý lý thuyết vật lý toán Mã số: 60 44 01 LUẬN VĂN THẠC SĨ KHOA HỌC Cán hướng dẫn : GS.TS Nguyễn Quang Báu Hà Nội - 2011 MỤC LỤC MỞ ĐẦU 1 Lí chọn đề tài Cấu trúc luận văn CHƢƠNG SIÊU MẠNG PHA TẠP VÀ BÀI TOÁN ẢNH HƢỞNG CỦA TRƢỜNG BỨC XẠ LASER LÊN HẤP THỤ SÓNG ĐIỆN TỪ YẾU BỞI ĐIỆN TỬ TRONG BÁN DẪN KHỐI 1.1 Tổng quan siêu mạng pha tạp .4 1.1.1 Siêu mạng pha tạp 1.1.2 Hàm sóng phổ lƣợng electron siêu mạng pha tạp 1.2 Hệ số hấp thụ sóng điện từ yếu điện tử bán dẫn khối có mặt trƣờng laser 1.2.1 Xây dựng phƣơng trình động lƣợng tử cho điện tử bán dẫn khối .5 1.2.2 Tính hệ số hấp thụ 14 CHƢƠNG PHƢƠNG TRÌNH ĐỘNG LƢỢNG TỬ VÀ HỆ SỐ HẤP THỤ SĨNG ĐIỆN TỪ YẾU BỞI ĐIỆN TỬ GIAM CẦM TRONG SIÊU MẠNG PHA TẠP KHI CÓ MẶT TRƢỜNG LASER (TRƢỜNG HỢP TÁN XẠ ĐIỆN TỬ-PHONON ÂM) 20 2.1 Hamiltonian hệ điện tử-phonon siêu mạng pha tạp .20 2.2 Xây dựng phƣơng trình động lƣợng tử cho điện tử siêu mạng pha tạp: 21 2.3 Tính hệ số hấp thụ sóng điện từ yếu điện tử giam cầm siêu mạng pha tạp kể đến ảnh hƣởng trƣờng xạ Laser 33 CHƢƠNG TÍNH TỐN SỐ VÀ KẾT LUẬN 50 3.1 Tính tốn số vẽ đồ thị cho hệ số hấp thụ cho trƣờng hợp siêu mạng pha tạp n-GaAs/p-GaAs 50 3.2 Thảo luận kết thu đƣợc: 53 KẾT LUẬN 54 TÀI LIỆU THAM KHẢO 55 PHỤ LỤC 57 66 DANH MỤC HÌNH VẼ Hình 3.1: Sự phụ thuộc hệ số hấp thụ vào nhiệt độ T Hình 3.2: Sự phụ thuộc hệ số hấp thụ vào cường độ sóng điện từ mạnh E01 Hình 3.3: Sự phụ thuộc hệ số hấp thụ vào lượng sóng điện từ mạnh (Laser) Hình 3.4: Sự phụ thuộc hệ số hấp thụ vào lượng sóng điện từ yếu Hình 3.5: Sự phụ thuộc hệ số hấp thụ vào độ rộng hố lượng tử L MỞ ĐẦU Lí chọn đề tài Khởi đầu từ thành công rực rỡ vật liệu bán dẫn với phát triển mạnh mẽ công nghệ nuôi tinh thể, người ta chế tạo nhiều cấu trúc nanô Song song với phát triển công nghệ chế tạo phát triển kỹ thuật đo hiệu ứng vật lý cấp độ vi mơ Có thể nói rằng, hai thập niên vừa qua cấu trúc tinh thể nanô ( màng mỏng, siêu mạng, hố lượng tử, dây lượng tử, chấm lượng tử, …) dần thay vật liệu bán dẫn khối kinh điển 1  4 Trong cấu trúc nanô vậy, chuyển động hạt dẫn bị giới hạn nghiêm ngặt dọc theo hướng tọa độ với vùng có kích thước đặc trưng vào cỡ bậc bước sóng De Broglie Trong cấu trúc nanơ này, tính chất vật lí điện tử thay đổi đáng kể, xuất số tính chất khác, gọi hiệu ứng kích thước Ở quy luật học lượng tử bắt đầu có hiệu lực, đặc trưng hệ điện tử hàm sóng phổ lượng bị biến đổi Phổ lượng bị gián đoạn dọc theo hướng tọa độ giới hạn Do tính chất quang, điện hệ biến đổi mở khả ứng dụng cho linh kiện điện tử làm việc theo ngun lí hồn tồn cơng nghệ đại, siêu nhỏ – đa – thông Hệ vật liệu thấp chiều cấu trúc hoàn toàn mới, khác hẳn với vật liệu trước chia hệ thấp chiều làm loại: hệ không chiều (0D), hệ chiều (1D), hệ hai chiều (2D) Đối với hệ hai chiều (2D), cụ thể siêu mạng pha tạp, phổ lượng điện tử trường hợp trở nên gián đoạn theo chiều siêu mạng pha tạp điện tử chuyển động tự theo hai chiều, chiều bị hạn chế Chính gián đoạn phổ lượng hạn chế chuyển động điện tử theo chiều lần lại ảnh hưởng lên tính chất phi tuyến hệ 5  20 Trong lĩnh vực nghiên cứu lý thuyết, cơng trình ảnh hưởng sóng điện từ mạnh lên sóng điện từ yếu bán dẫn khối nghiên cứu nhiều Thời gian gần có cơng trình nghiên cứu ảnh hưởng sóng điện từ laze lên hấp thụ phi tuyến sóng điện từ yếu điện tử giam cầm bán dẫn thấp chiều Tuy nhiên, siêu mạng pha tạp, ảnh hưởng trường xạ laze lên hấp thụ sóng điện từ yếu điện tử giam cầm đề tài mở Do đó, luận văn này, chúng tơi chọn vấn đề nghiên cứu “Ảnh hưởng trường xạ Laser lên hấp thụ sóng điện từ yếu điện tử giam cầm siêu mạng pha tạp (trường hợp tán xạ điện tử - phonon âm)” Về phương pháp nghiên cứu: Chúng ta sử dụng nhiều phương pháp lý thuyết khác để giải tốn vật lí kể phương pháp có ưu điểm riêng nên việc áp dụng chúng phụ thuộc vào tốn cụ thể Trong luận văn này, chúng tơi sử dụng phương pháp phương trình động lượng tử Từ Hamilton hệ biểu diễn lượng tử hóa lần hai ta xây dựng phương trình động lượng tử cho điện tử giam cầm, áp dụng phương trình động lượng tử để tính mật độ dịng hạt tải, từ suy biểu thức giải tích hệ số hấp thụ Đây phương pháp sử dụng rộng rãi nghiên cứu hệ bán dẫn thấp chiều, đạt hiệu cao cho kết có ý nghĩa khoa học định Đối tượng nghiên cứu luận văn cấu trúc bán dẫn thấp chiều thuộc hệ hai chiều, cụ thể siêu mạng pha tạp 2.Cấu trúc luận văn Ngoài phần mở đầu, kết luận, tài liệu tham khảo phụ lục, luận văn chia làm chương, mục, hình vẽ, tổng cộng 66 trang: Chƣơng I: Siêu mạng pha tạp toán ảnh hưởng trường xạ laser lên hấp thụ sóng điện từ yếu điện tử bán dẫn khối Chƣơng II: Phương trình động lượng tử hệ số hấp thụ sóng điện yếu từ điện tử giam cầm siêu mạng pha tạp có mặt trường xạ Laser (trường hợp tán xạ điện tử - phonon âm) Chƣơng III: Tính tốn số kết luận Trong đó, chương II chương III hai chương chứa đựng kết luận văn Kết luận quan trọng rút từ kết nghiên cứu luận văn là: số điều kiện thỏa mãn định liên quan đến nhiệt độ lượng sóng điện từ (tần số sóng điện từ), hệ số hấp thụ sóng điện từ yếu trở nên âm, tức hệ số hấp thụ trở thành hệ số gia tăng sóng điện từ yếu Điều mở khả gia tăng sóng điện từ yếu siêu mạng pha tạp có mặt sóng điện từ mạnh khác Đây điều mà bán dẫn khối xảy CHƢƠNG SIÊU MẠNG PHA TẠP VÀ BÀI TOÁN ẢNH HƢỞNG CỦA TRƢỜNG BỨC XẠ LASER LÊN HẤP THỤ SÓNG ĐIỆN TỪ YẾU BỞI ĐIỆN TỬ TRONG BÁN DẪN KHỐI Trong chương trình bày khái quát siêu mạng pha tạp (cấu trúc phổ lượng, hàm sóng điện từ) từ phương pháp phương trình động lượng tử đưa biểu thức giải tích cho hệ số hấp thụ sóng điện từ yếu điện tử bán dẫn khối chịu ảnh hưởng trường laser 1.1 Tổng quan siêu mạng pha tạp 1.1.1 Siêu mạng pha tạp Bán dẫn siêu mạng loại cấu trúc tuần hoàn nhân tạo gồm lớp bán dẫn thuộc hai loại khác có độ dày cỡ nanomet đặt Do cấu trúc tuần hoàn, bán dẫn siêu mạng, tuần hoàn mạng tinh thể, electron phải chịu tuần hoàn phụ siêu mạng tạo với chu kỳ lớn số mạng nhiều Thế phụ tạo nên khác biệt đáy vùng dẫn hai bán dẫn cấu trúc thành siêu mạng Trong bán dẫn siêu mạng, độ rộng lớp đủ hẹp để electron xuyên qua lớp mỏng nhau, coi siêu mạng tuần hoàn bổ sung vào mạng tinh thể Bán dẫn siêu mạng chia thành hai loại: bán dẫn siêu mạng pha tạp bán dẫn siêu mạng hợp phần Bán dẫn siêu mạng pha tạp có cấu tạo hố siêu mạng tạo thành từ hai lớp bán dẫn loại pha tạp khác 1.1.2 Hàm sóng phổ lƣợng electron siêu mạng pha tạp Hàm sóng bán dẫn siêu mạng pha tạp có dạng:     e  n,p  ei p r U n r   S0 j 1  n  z  jd  ik z jz Phổ lượng:   2 p   n,p   p  n    *  2m   Trong đó: n = 1,2,3 số lượng tử phổ lượng theo phương z    p  p   p z : vectơ xung lượng điện tử (chính xác vectơ sóng điện tử) S0: số chu kì siêu mạng d: chu kì siêu mạng e m* : điện tích khối lượng hiệu dụng điện tử  p  2e  nD : tần số plasma gây tạp chất donor  m* nD : nồng độ pha tạp  : số điện Từ ta thấy phổ lượng điện tử bị giam cầm siêu mạng pha tạp nhận giá trị lượng gián đoạn, không giống bán dẫn khối, phổ lượng liên tục tồn khơng gian Sự biến đổi phổ lượng gây khác biệt đáng kể tất tính chất điện tử siêu mạng pha tạp so với mẫu khối 1.2 Hệ số hấp thụ sóng điện từ yếu điện tử bán dẫn khối có mặt trƣờng laser 1.2.1 Xây dựng phƣơng trình động lƣợng tử cho điện tử bán dẫn khối Ta có Hamilton hệ điện tử - phonon bán dẫn khối là: H  H e  H ph  H e ph Với :   e     He     p  A(t )  a p a p  c   p H ph   q bqbq  q  H e ph   Cq ap  q ap bq  bq   q, p  (1.1) Phương trình động lượng tử cho điện tử có dạng: i np (t ) t   ap ap , Hˆ    (1.2) t Vế phải (1.2) có tương ứng ba số hạng với tốn tử Hamilton Ta tính số hạng Số hạng thứ nhất:     e        e          a ; a  p '  A ( t ) a a    p p    p  c A(t )  a p a p , a p 'a p '     p' p'    c   p ' p   t   e               p '  A(t )  a p a p ' p , p '  ap ap'ap ap '  ap'ap p ,p'  ap'ap ap 'ap   c   p'   e      p  A(t )  ap ap  ap ap  c           b b Số hạng thứ hai:  ap ap ;    q q q  q    toán tử a, b hai loại độc lập t chúng giao hốn với Số hạng thứ ba:       a ;  a  a b  b  a C  p p  q q    q p ' q p ' q, p '     Cq  q aa      p p q q t b      C ap ap ' p ,p'q  ap ' q ap p ,p' bq  bq    q q, p ' t  ap ap q bq t  ap q ap bq t  ap q ap bq t   Cq  Fp ,p q ,q (t )  Fp*q ,p , q (t )  Fp ,p q ,q (t )  Fp*,p q , q (t )     q Vậy phương trình (1.2) trở thành: i np (t ) t  Cq  Fp ,p q ,q (t )  Fp*q ,p , q (t )  Fp ,p q ,q (t )  Fp*,p q , q (t )     q   (t )  a a b Với Fp , p ,q p p q 2 (1.3) t  (t ) thơng qua phương trình: Để giải (1.3) ta cần tính Fp , p ,q i  (t ) Fp , p ,q t   ap a b ; H   p2 q  (1.4) t Do thi anpha - E01 x 10 16 T=150K T=350K 14 he so hap thu anpha 12 10 0 Bien song E01 10 12 V/m x 10 Hình 3.2 Sự phụ thuộc hệ số hấp thụ vào cƣờng độ sóng điện từ mạnh E01 Do thi anpha - nang luong song 25 T=80K T=400K he so hap thu anpha 20 15 10 -5 0.012 0.014 0.016 0.018 0.02 0.022 0.024 0.026 Nang luong song h*omega 0.028 0.03 eV Hình 3.3: Sự phụ thuộc hệ số hấp thụ vào lƣợng sóng điện từ mạnh (laser) 51 Do thi anpha - nang luong song 30 T=100K T=350K he so hap thu anpha 20 10 -10 -20 -30 0.01 0.02 0.03 0.04 Nang luong song h*omega 0.05 0.06 eV Hình 3.4: Sự phụ thuộc hệ số hấp thụ lƣợng sóng điện từ yếu Do thi anpha - L 0.018 T=200K T=350K 0.016 he so hap thu anpha 0.014 0.012 0.01 0.008 0.006 0.004 0.002 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Chieu dai sieu mang L 0.8 0.9 -7 x 10 m Hình 3.5: Sự phụ thuộc hệ số hấp thụ vào chiều dài siêu mạng 52 3.2 Thảo luận kết thu đƣợc: Nhìn vào kết tính số vẽ đồ thị hệ số hấp thụ phi tuyến sóng điện từ điện tử giam cầm siêu mạng pha tạp có từ trường ngồi, ta có số nhận xét sau: Hình 3.1 phụ thuộc hệ số hấp thụ vào nhiệt độ T hệ nhiệt độ tăng từ 100K tới 800K, trường hợp biên độ sóng E01  3,5.106V / m Đồ thị cho thấy hệ số hấp thụ tăng dần đạt giá trị cực đại, sau lại giảm dần nhiệt độ tiếp tục tăng Hình 3.2 phụ thuộc hệ số hấp thụ α vào cường độ sóng điện từ mạnh laser E01 hai nhiệt độ T  150K T  350K Từ đồ thị, ta nhận thấy rằng: α phụ thuộc phi tuyến vào cường độ E01 sóng điện từ mạnh, cường độ E01 tăng hệ số hấp thụ α tăng lên nhanh Hình 3.3 hình 3.4 phụ thuộc hệ số hấp thụ vào lượng sóng điện từ mạnh 1 lượng sóng điện từ yếu 2 Nhìn vào đồ thị ta thấy vùng nhiệt độ thấp, hệ số hấp thụ nhận giá trị âm Khi nhiệt độ cao hệ số hấp thụ nhận giá trị dương giảm nhanh lượng 1 tăng, tốc độ giảm hệ số hấp thụ chậm lượng 2 giảm Hình H3.5 biểu diễn phụ thuộc phi tuyến hệ số hấp thụ vào chiều dài siêu mạng L Từ đồ thị ta thấy hệ số hấp thụ đạt giá trị cực đại chiều dài siêu mạng cỡ 0,3.107 m Như vậy, đồ thị cho thấy phụ thuộc hệ số hấp thụ sóng điện từ yếu điện tử giam cầm có mặt xạ sóng điện từ Laser vào đại lượng kể nói chung khơng tuyến tính nhận giá trị âm 53 KẾT LUẬN Trên sở phương trình động lượng tử, tính tốn hệ số hấp thụ phi tuyến sóng điện từ yếu điện tử giam cầm siêu mạng pha tạp có mặt trường xạ Laser (trường hợp tán xạ điện tử - phonon âm) với kết là: Xuất phát từ Hamiltonian hệ điện tử - phonon, phương pháp phương trình động lượng tử xây dựng biểu thức hệ số hấp thụ sóng điện từ yếu điện tử giam cầm siêu mạng pha tạp có mặt trường xạ Laser (trường hợp tán xạ điện tử - phonon âm) Đã khảo sát tính số vẽ đồ thị biểu thức hệ số hấp thụ sóng điện từ yếu điện tử giam cầm siêu mạng pha tạp n-GaAs/p-GaAs phụ thuộc phi tuyến mạnh hệ số hấp thụ sóng điện từ vào tần số sóng điện từ yếu, biên độ tần số trường Laser nhiệt độ hệ kích thước lượng tử siêu mạng pha tạp Kết lý thú số điều kiện thỏa mãn định liên quan đến nhiệt độ lượng sóng điện từ (tần số sóng điện từ), hệ số hấp thụ sóng điện từ yếu trở nên âm, tức hệ số hấp thụ trở thành hệ số gia tăng sóng điện từ yếu Điều mở khả gia tăng sóng điện từ yếu siêu mạng pha tạp có mặt sóng điện từ mạnh khác Đây điều mà bán dẫn khối xảy 54 TÀI LIỆU THAM KHẢO Tài liệu tiếng việt: [1] Lê Đình (2008), Một số hiệu ứng cao tần tương tác electron – phonon dây lượng tử bán dẫn, Luận án tiến sĩ vật lý, ĐHKHTN, ĐHQG Hà Nội [2] Nguyễn Quang Báu (chủ biên), Đỗ Quốc Hùng, Lê Tuấn, “ Lý thuyết bán dẫn đại”, NXB ĐHQGHN (2011) [3] Nguyễn Quang Báu (chủ biên), Nguyễn Vũ Nhân, Phạm Văn Bền, “Vật lý bán dẫn thấp chiều”, NXB ĐHQGHN (2010) [4] Nguyễn Văn Hùng, “Lý thuyết chất rắn”, NXB ĐHQGHN (1999) Tài liệu tiếng anh: [5] A.P.Solin.Sov Phys Usp., 28, 972-993, (1985) [6] G.M.Shmelev et, Phys Stat Sol B (82), pp 391-1937, (1977) [7] G M Shmelev, L.A Chaikovkii and N Q Bau, Sov Phys Semicond 12,1932 (1978) [8] Le Dinh, Tran Cong Phong, Communications in physics, Vietnam Acedemy of Science and Technology, Suppl., 117.(2006) [9] L Esaki, R.Tsu IBM J.Res Develop., 14, (1970), 246-248 [10] Nguyễn Quang Báu and Trần Công Phong J.Phys Soc Jpn., 67, 3875, (1998) [11] Nguyen Quang Bau, Do Manh Hung, Nguyen Ngoc Bich, Jounal of the Korean Physical Society, Vol.54, No.2 -765.(2009) [12] Nguyen Quang Bau, VNU Journal of Science, Mathematics – physics, (XXII), 4.(2006) [13] N Nishguchi, Phys.Rev.B, 52, 5379 (1995) [14] P.Zhao, Phays Rev B,49,13589 (1994) [15] J.S Harris Jr, Int J Mod Phys B4, 1149 (1990) [16] R.A.Smith Wave Mechanics Crystal Solids, London, (1961) [17] R Dingle, Confined Carries Quantum states in ultra-thin Semiconductor Heterostructures, New York (1975) [18] R Tsu and L.Esaki, Appl Phys.Lett.22, 562 (1193) 55 [19] Tran Cong Phong, Jounal of the Korean Physical Society, Vol.42, No.5 -647.(2003) [20] V V Pavlovich and E M Epshtein, Sov Phys Semicond 11,809 (1977) 56 PHỤ LỤC Chƣơng trình tính số vẽ đồ thị cho hệ số hấp thụ sóng điện yếu điện tử giam cầm siêu mạng pha tạp n-GaAs/p-GaAs có mặt hai sóng Để làm đƣợc việc này, sử dụng ngôn ngữ lập trình Matlap 7.0 for Windows 1.Hệ số hấp thụ  phụ thuộc vào nhiệt độ T: clc;clear all;close all; T=linspace(100,800,150); %E02=linspace(1.5e7,1e9,1000); E01=3.5e6 E012=3.1e6; %T=245.1; L=90e-10; omeg1=2.1e12; omeg2=13.8e12;gama=0;k=1.3807e-23; L=85e-10; V=13.5*1.6021e-19; c=1;Xinf=10.9;%V chinh dien the ro=5320;vs=5370; e0=1; m0=9.109389e-31;h=1;%.05459e-34; m1=.067*m0;e=2.07*e0;n0=1e23; y1=y(omeg1,omeg2,T,L,E01,gama) y2=y(omeg1,omeg2,T,L,E012,gama) plot(T,y1,'.r');grid on;hold on plot(T,y2,'.g') title('Do thi anpha - T'); xlabel('Nhiet (K)'); ylabel('he so hap thu anpha'); 57 2.Hệ số hấp thụ  phụ thuộc vào cƣờng độ sóng điện từ E01: %tinh toan alpha phu thuoc vao bien dien truong E01 clc;clear all;close all; E01=linspace(7.1e6,10.1e7,100); %E02=linspace(1.5e7,1e9,1000); %omeg2=linspace(130e11,87.1e12,99); T1=150; T2=350 L=85e-10;omeg1=2.1e12; omeg2=13.8e12;gama=0;k=1.3807e-23; V=13.5*1.6021e-19;c=1;Xinf=10.9;%V chinh dien the ro=5320;vs=5370;Xif=10.9 e0=1;m0=9.109389e-31;h=1;%.05459e-34; m1=.067*m0;e=2.07*e0;n0=1e23; y1=y(omeg1,omeg2,T1,L,E01,gama) plot(E01,y1,'.b');grid on; hold on; y2=y(omeg1,omeg2,T2,L,E01,gama) plot(E01,y2,'-r') title('Do thi anpha - E01'); xlabel(' Bien song E01 V/m'); ylabel('he so hap thu anpha'); 3.Chƣơng trình tính số phụ thuộc hệ số hấp thụ  theo lƣợng sóng điện từ  1( trƣờng hợp phi tuyến) % Do thi alpha phu thuoc nang luong song h*omega1 clc;clear all;close all; %E01=1.5e6 58 %E02=2e7; %omeg1=linspace(23.8e12,47.1e12,99); T1=80; T2=400; %T3=350; E01=8.1e6; E02=2.1e7; omeg1=linspace(20e12,45e12,100); N1=omeg1; omeg2=2e11; L=85e-10;gama=10; k=1.3807e-23;V=13.5*1.6021e-19;c=1;Xinf=10.9 ro=5320;vs=5370; e0=1.6e-19;m0=9.109389e-31;h=1; m1=.067*m0;e=2.07*e0;n0=1e20; N1=omeg1*1.0544e-34/1.6*1e19; y1=y(omeg1,omeg2,T1,L,E01,gama) y2=y(omeg1,omeg2,T2,L,E01,gama) % y3=y(omeg1,omeg2,T3,L,E01,gama) plot(N1,y1,'.b');grid on;hold on; plot(N1,y2,'-r');grid on;hold on; %plot(N1,y3,'+g'); title('Do thi anpha - nang luong song 1'); xlabel(' Nang luong song h*omega eV'); ylabel('he so hap thu anpha'); 4.Chƣơng trình tính số phụ thuộc hệ số hấp thụ  theo lƣợng sóng điện từ  % Do thi alpha phu thuoc nang luong song h*omega2 clc;clear all;close all; 59 E01=5.5e7 %E02=2e7; %omeg2=linspace(13e12,88e12,99); omeg2=linspace(13e12,88e12,100); T1=100; T2=350; %T3=350.1; L=850e-10;omeg1=20e12; gama=0;k=1.3807e-23; V=13.5*1.6021e-19;c=1;Xinf=10.9;%V chinh dien the ro=5320;vs=5370;Xif=10.9 e0=1;m0=9.109389e-31;h=1;%.05459e-34; m1=.067*m0;e=2.07*e0;n0=1e23; %a11=2*pi^2*e^2*k.*T*V^2./(c*sqrt(Xinf)*ro*vs^2*m1^2.*omeg2.^3); N1=omeg2*1.0544e-34/1.6*1e19; y1=y(omeg1,omeg2,T1,L,E01,gama); plot(N1,y1,'.');grid on;hold on; y2=y(omeg1,omeg2,T2,L,E01,gama); plot(N1,y2,'-r');grid on;hold on; % y3=y(omeg1,omeg2,T3,L,E01,gama); % plot(N1,y3,'*g');grid on;hold on; title('Do thi anpha - nang luong song 2'); xlabel(' Nang luong song h*omega eV'); ylabel('he so hap thu anpha'); Chƣơng trình tính số phụ thuộc hệ số hấp thụ  theo chiều dài siêu mạng L: %tinh toan alpha phu thuoc vao chieu dai sieu mang L clc;clear all;close all; T1=200; 60 %T2=295; T3=350; omeg2=2e12; %L=90e-10; omeg1=13.8e12; gama=0;k=1.3807e-23;V=13.5*1.6021e-19;c=1;Xinf=10.9;%V chinh dien the ro=5320;vs=5370;Xif=10.9 e0=1;m0=9.109389e-31;h=1;%.05459e-34; m1=.067*m0;e=2.07*e0;n0=1e23; al=0; E01=1e6; L=linspace(50e-10,999.1e-10,200); y1=y(omeg1,omeg2,T1,L,E01,gama); plot(L,y1,'.r');grid on;hold on; %y2=y(omeg1,omeg2,T2,L,E01,gama); %plot(L,y2,'.b');grid on;hold on; y3=y(omeg1,omeg2,T3,L,E01,gama); plot(L,y3,'-g');grid on;hold on; title('Do thi anpha - L'); xlabel(' Chieu dai sieu mang L m'); ylabel('he so hap thu anpha'); Các chƣơng trình 6.1 function B= B(T,L,delta) n0=1e23; e0=1.60219e-19; e=2.07*e0;%kiem tra lai gia tri so cu the h=1.05e-34 ; k=1.3807e-23; 61 V0=1;%thay doi lai gia tri m1=(9.1095*1e-31).*0.066; n01=n0*(e*pi)^(3/2)*h^3./(V0*(m1*k.*T).^(3/2)); B=1/(2*pi)^3.*n01*m1/h^3.*sqrt(2*m1*k.*T*pi)*2.*pi./L.*(delta+1);%ham I thay doi end 6.2 function D= D(s,m,omeg1,omeg2,T,L,n,n1) h=1.05e-34 ; k=1.3807e-23;%thay doi lai gia tri m1=(9.1095*1e-31).*0.066; epxilonn=pi^2*h^2*n^2./(2*m1.*L.^2); epxilonn1=pi^2*h^2*n1^2./(2*m1.*L.^2); a= ksi(s,m,omeg1,omeg2,L,n,n1)./(2*k.*T); c=ksi(s,m,omeg1,omeg2,L,n,n1); b=besselk(1,abs(c)./(2*k.*T));%k1 D=pi.*exp(-a).*(4*m1^2.*c.^2/h^4).^(1/2).*b.*(exp((-epxilonn./(2*k.*T))-exp((1./(2*k.*T)).*(epxilonn1-c)))); End 6.3 function G = G(s,m,omeg1,omeg2,T,L,E01,gama,n,n1) h=1.05e-34 ; k=1.3807e-23;%thay doi lai gia tri e0=1.60219e-19; e=2.07*e0;%kiem tra lai gia tri so cu the m0=9.109389e-31; m1=.067*m0; a1=e.*E01./(m1.*omeg1.^2); a= ksi(s,m,omeg1,omeg2,L,n,n1)./(2*k.*T); c= ksi(s,m,omeg1,omeg2,L,n,n1); 62 b=besselk(3,abs(c)./(2*k.*T));%k3 epxilonn=pi^2*h^2*n^2./(2*m1.*L.^2); epxilonn1=pi^2*h^2*n1^2./(2*m1.*L.^2); %H=a1.^2*(pi/2+(pi/4)*cos(2*gama))*exp(a/(2*k*T)*((4*m1.^2*(a^2)/(h.^4)*b*(exp((-epxilon0(n)/(2*k*T)-exp((1/(2*k*T))*(epxilon(n1)*(-a); G=(a1.^4)*(3*pi/8+pi/4.*cos(2.*gama)).*exp(a).*((4*m1^2.*c.^2)/h^4).^(3/2).*b.*(exp(-epxilonn./(2*k.*T))-exp((1./(2*k.*T).*(epxilonn1-c)))); end 6.4 function H = H(s,m,omeg1,omeg2,T,L,E01,gama,n,n1) h=1.05e-34 ; k=1.3807e-23;%thay doi lai gia tri e0=1.60219e-19; e=2.07*e0;%kiem tra lai gia tri so cu the m0=9.109389e-31; m1=.067*m0; a1=e.*E01./(m1.*omeg1.^2); a= ksi(s,m,omeg1,omeg2,L,n,n1)./(2*k.*T); c= ksi(s,m,omeg1,omeg2,L,n,n1); b=besselk(2,abs(c)./(2*k.*T));%k2 epxilonn=pi^2*h^2*n^2./(2*m1.*L.^2); epxilonn1=pi^2*h^2*n1^2./(2*m1.*L.^2); %H=a1.^2*(pi/2+(pi/4)*cos(2*gama))*exp(a/(2*k*T)*((4*m1.^2*(a^2)/(h.^4)*b*(exp((-epxilon0(n)/(2*k*T)-exp((1/(2*k*T))*(epxilon(n1)*(-a); H=(a1.^2).*(pi/2+pi/4.*cos(2.*gama)).*exp(-a).*((4*m1^2.*c.^2)/h^4).*b.*(exp(epxilonn./(2*k.*T))-exp((-1./(2*k.*T).*(epxilonn1-c)))); end 63 6.5 function y=y(omeg1,omeg2,T,L,E01,gama) V=13.5*1.6021e-19;c=1;Xinf=10.9;%V chinh dien the ro=5320;vs=5370;Xif=10.9 e0=1;m0=9.109389e-31;h=1;%.05459e-34; m1=.067*m0;e=2.07*e0;n0=1e23; k=1.3807e-23; a11=2*pi^2*e^2*k.*T*V^2./(c*sqrt(Xinf)*ro*vs^2*m1^2.*omeg2.^3); y=0; for n=1:3 for n1=1:3 if n==n1 delta=1 else delta=0 end al=a11.*B(T,L,delta)/(2*pi)^2.*(D(0,1,omeg1,omeg2,T,L,n,n1)-D(0,2,omeg1,omeg2,T,L,n,n1)- 1/2.*(H(0,1,omeg1,omeg2,T,L,E01,gama,n,n1)-H(0,1,omeg1,omeg2,T,L,E01,gama,n,n1))+ 3/32.*(G(0,1,omeg1,omeg2,T,L,E01,gama,n,n1)-G(0,1,omeg1,omeg2,T,L,E01,gama,n,n1))+ 1/4.*(H(-1,1,omeg1,omeg2,T,L,E01,gama,n,n1)-H(-1,1,omeg1,omeg2,T,L,E01,gama,n,n1)+ H(1,1,omeg1,omeg2,T,L,E01,gama,n,n1)-H(1,1,omeg1,omeg2,T,L,E01,gama,n,n1))- 1/16.*(G(-1,1,omeg1,omeg2,T,L,E01,gama,n,n1)-G(-1,1,omeg1,omeg2,T,L,E01,gama,n,n1)+ G(1,1,omeg1,omeg2,T,L,E01,gama,n,n1)-G(1,1,omeg1,omeg2,T,L,E01,gama,n,n1))+ 64 1/64.*(G(-2,1,omeg1,omeg2,T,L,E01,gama,n,n1)-G(-2,1,omeg1,omeg2,T,L,E01,gama,n,n1)+ G(2,1,omeg1,omeg2,T,L,E01,gama,n,n1)-G(2,1,omeg1,omeg2,T,L,E01,gama,n,n1))); y=y+al; end end end 6.6 %xay dung ham ksi phu thuoc cac yeu to thay doi n,n' function Ksi = ksi(s,m,omeg1,omeg2,L,n,n1) % syms n1 n%n1 la n phay h=1.05e-34 ; w0=3.625e-4*1.60219e-19/h; m1=(9.1095*1e-31).*0.066; Ksi=(h^2*pi^2)./(2*m1.*L.^2)*(n1-n)+h*w0-s*h.*omeg1-m*h.*omeg2; end 65 ... cứu ảnh hưởng sóng điện từ laze lên hấp thụ phi tuyến sóng điện từ yếu điện tử giam cầm bán dẫn thấp chiều Tuy nhiên, siêu mạng pha tạp, ảnh hưởng trường xạ laze lên hấp thụ sóng điện từ yếu điện. .. NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Bùi Hồng Phượng ẢNH HƯỞNG CỦA TRƯỜNG BỨC XẠ LASER LÊN HẤP THỤ SÓNG ĐIỆN TỬ YẾU BỞI ĐIỆN TỬ GIAM CẦM TRONG SIÊU MẠNG PHA TẠP (TÁN XẠ ĐIỆN TỬ - PHONON ÂM) Chuyên... Siêu mạng pha tạp toán ảnh hưởng trường xạ laser lên hấp thụ sóng điện từ yếu điện tử bán dẫn khối Chƣơng II: Phương trình động lượng tử hệ số hấp thụ sóng điện yếu từ điện tử giam cầm siêu mạng