Đề tài Các lý thuyết hiện đại về nhiệt phát quang

CÁC LÝ THUYẾT HIỆN ĐẠI VỀ NHIỆT PHÁT QUANG Thái Ngọc Ánh ∗ Mục lục 1 Mở đầu 2 2 Các bẫy và các trạng thái tái hợp. 2 3 Tương tác động học. 4 4 Phân bố các bẫy liên tục. 8 5 Kết luận 14 ∗Cao học Vật lý - Đại học Khoa học 1 Tóm tắt. Trong tiểu luận này tôi cố gắng trình bày khá chi tiết lý thuyết về các bẫy và các trạng thái tái hợp, tương tác động học, phân bố các bẫy liên tục. 1 Mở đầu Nhiệt phát quang (TL) là sự phát ra bức xạ từ các chất cách điện hoặc bán dẫn khi vật liệu được nung nóng sau khi được chiếu xạ ở nhiệt độ thấp (nhiệt độ phòng, nhiệt độ ni tơ lỏng, ...). TL có nhiều ứng dụng trong việc xác định các sai hỏng, khuyết tật của tinh thể; dùng liều kế TL; tính tuổi của khoáng vật và cổ vật. Nhiệt phát quang và ứng dụng là học phần không thể thiếu cho các học viên Cao học chuyên ngành Quang - Quang phổ để tiến hành thực nghiện. Lý thuyết của TL rất phong phú và đa dạng và đang từng ngày phát triển ở trên thế giới. Việc ứng dụng TL vào các quá trình như đã nói ở trên rất hiệu quả và chính xác. Để hiểu hơn môn học tôi nghiên cứu phần "Các lý thuyết hiện đại về nhiệt phát quang" để làm đề tài tiểu luận. Rất mong sự góp ý của quý thầy cô giáo và các bạn đọc để bài viết hoàn thiện hơn. 2 Các bẫy và các trạng thái tái hợp. Sự khảo sát trước đó đã diễn tả sự lệ thuộc kiểu một bẫy, một tâm đơn giãn, từ đó phát triển được biểu thức diễn tả đường TL cơ bản đối với vật liệu thực - tức là dạng của các đỉnh TL và sự phụ thuộc của chúng vào nồng độ điện tích bị bẫy, nhiệt độ và độ sâu các bẫy. Tuy nhiên, mặc dù sự có ích của các kiểu này và các khảo sát, sự kiện còn lại thì không có vật liệu thật như vậy, vật liệu mà chỉ có một bẫy và một tâm tái hợp. Kết quả quan trọng nhất của kiểu này là bàn về n = h, tức là số bẫy điện tử bằng số bẫy lỗ trống. Đây chính là điều kiện trung hoà về điện và thêm giả thuyết chuẩn cân bằng và nc0 = 0 (cùng dẫn đến nc ' 0). Trong đa số vật liệu thực, không tránh khỏi sự hiện diện của các bẫy lân cận của sự cân bằng nhiệt của điện tích bẫy là lớn hơn tín hiệu mà TL ghi được. Nói cách khác, tại mọi nhiệt độ; có thể tồn tại các bẫy sâu với mật độ nào đó trong suất quá trình làm sách các bẫy nông và ghi tín hiệu TL. Kết quả của điều này là n 6= m. 2 Thật vậy, vẫn còn điều kiện nc0 = 0, biểu thức điều kiện cân bằng nhiệt trở thành: n+ h = m (1) với h là nồng độ của các bẫy điện tử ở các bẫy sâu hơn. Các bẫy sâu được xem như tháo ra nhờ nhiệt và nhiều tác giả thích quan tâm đến điều này trong phân tích quá trình TL. Phương trình tốc độ bây giờ thêm vào số hạng dh dt = (H − h)υσh (2) và kết quả là dnc dt = dm dt − dn dt − dh dt (3) với H là tổng số các bẫy có sẳn, sâu, tách do nhiệt và σn là tiết diện bắt của các bẫy sâu. Biểu thức dn dt và dm dt đã có ở mục trước. Viết các phương trình này trong trường hợp chúng ta quan tâm đến các vùng sâu, các bẩy tách nhiệt có thể bắt điện tử giải thoát từ các bẫy nông tại năng lượng Et. Trong con đường này, các bẫy được gọi là đã tương tác với các bẫy sâu cạnh tranh với các vị trí tái hợp đối với sự bắt các điện tử được giải phóng từ các bẫy nông. Trong điều kiện đặc biệt h ' H thì dh dt ' 0 và ta có thể viết theo Kelly và Bra¨unlich, tương đương với phương trình tổng quát với trường hợp hai bẫy ITL = ns exp{−Et kT }(n+ h)σmn [(N − n)σn + (n+ h)σmn] (4) Với giới hạn này, áp dụng giả thuyết tái bắt chậm (tức là (N − n)σn ¿ (n + h)σmn) từ (4) đưa thẳng đến phương trình TL bậc 1, Randall - Wilkins. Trái lại, trường hợp tái bắt nhanh (tức là (N − n)σn À (n+ h)σmn) với n ¿ N , ta được ITL = ns exp{−Et kT }(n+ h)σmn Nσn (5) hay, cho σn = σmn, ta thấy rằng (4) trở thành ITL = ns exp{−Et kT }(n+ h) (N + h) (6) Phương trình (5) và (6) cả hai có thể biểu diễn dạng ITL = s ′n(n+ h) exp {−Et kT } (7) với s′ = sσmn Nσn , hay s′ = s N+h . Khai triển (7) ta được ITL = s ′nh exp {−Et kT } + s′n2 exp {−Et kT } (8) 3 và, như đã chỉ ra bởi tác giả Chen, điều này trong như là một sự pha trộn của động học bậc một và động học bậc hai. Rõ ràng, nếu h ¿ n, thì phương trình đưa về dạng bậc hai trong khi nếu h À n thì phương trình đưa về dạng bậc nhất. Với lý do này mà đã có số trộn bậc động học được đưa ra bởi Chen. Nghiệm của phương trình (8) là ITL = s′h2α exp{(hs′ β ) ∫ T T0 exp{−Et kθ }dθ} exp{−Et kT }[ exp{(hs′ β ) ∫ T T0 exp{−Et kθ }dθ} − α ]2 (9) với α = n0 n0+h . Dạng của phương trình (9) đồng nhất với dạng phương trình Randall - Wilkins khi α → 0, và dạng nó giống phương trình Garlick - Gibson khi α → 1. Từ những điều đã nghiên cứu trước rõ ràng rằng hình dạng, vị trí, kích cở (theo nhiệt độ) và dáng điệu (như là một hàm của sự tập trung đầy bẫy và tốc độ nhiệt)có thể được gói gon trong một phương trình cơ sở, phụ thuộc vào các giả thuyết ban đầu được dùng. Trong mỗi trường hợp, đĩnh TL được diễn tả bằng bốn thông số cơ bản n0, E, s và b (hay α) và phương trình tốc độ phức tạp có thể được rút ra hoặc là dạng động học bậc nhất, bậc hai hoặc là dạng trung gian (sự pha tron giữa các bậc) bằng cách áp dụng nhiều giả thuyết. Có lẽ hai giả thuyết quan trọng là dnc dt = 0 (chuẩn cân bằng (QE)) và h ' H (không có tương tác động học). Opanowicz so sánh phương trình (4) với biểu thức bậc tổng quát ta được nbs′ = nsγ (10) với γ = (n+ h)σmn [(N − n)σn + (n+ h)σmn (11) đối với trường hợp h = H. Thêm vào điều kiện n0 = N và s ′ = sn1−b0 . Opanowicz phát triển sự phụ thuộc vào nhiệt độ của thông số động học b(T ) b = ln(γn N ) ln( n N ) (12) Vì γ và n N phụ thuộc mạnh vào nhiệt độ nên b cũng vậy, trái với giả thuyết thường dùng thông số này là hằng số. Hệ quả kéo theo là b là thông số hình dạng, không có ý nghĩa vật lý. 3 Tương tác động học. Cách chung nhất của việc viết các phương trình tốc độ để diễn tả dòng các điện tử đi vào và đi ra khỏi vùng dẫn đối với hệ thống gồm nhiều bẫy và tâm tái hợp. Đối với 4 trường hợp tập hợp các bẫy điện tử rời rạc cho chỉ số i = 1 đến u và tập hợp các bẫy lỗ trống (tâm tái hợp) cho j = 1 đến υ ta có thể viết lại một cách hoàn chỉnh phương trình tốc độ như sau. Cho i = 1 đến u dni dt = −nisiexp {−Eti kT } + nc(Ni − ni)Ani (13) Cho j = 1 đến υ dmj dt = −ncmjAmnj (14) với Ani = υnσni và Amnj = υnσmnj. Tốc độ thay đổi theo thời gian của nồng độ điện tử tự do có thể được viết dnc dt = u∑ i=1 nisi exp {−Eti kT } − nc ( υ∑ j=1 mjAmnj + u∑ i=1 (Ni − ni)Ani ) (15) và vì chỉ có sự giải phóng điện tử bị bẩy do nhiệt được đề cập nên ta có dnυ dt = 0. Để phân tích tập hợp các phương trình này ta có thể tiến hành theo nhiều cách. Một trong các cách được phát triển bởi Levy người mà giữ gần chuẩn cân bằng và phát triển giống như các phương trình động học tổng quát đối với trường hợp này phức tạp hơn, tức là ITL = υ∑ j=1 E mjAmnj R + U (16) với E = u∑ i=1 nisi exp { −Eti kT } , (17) R = υ∑ j=1 mjAmnj (18) U = u∑ i=1 (Ni − ni)Ani (19) Trong cách viết (16) ta đã giả thuyết tất cả các quá trình tái hợp điện tử - lỗ trống đều phát xạ và đống góp vào tín hiệu TL. Nếu điều này không thể thì chỉ một phần quá trình tái hợp được dùng. Ngoài ra phương trình cũng giả thuyết rằng tất cả photon phát xạ đều được phát hiện với khả năng như nhau. Ngoài ra mỗi số tái hợp phải được nhân với một hiệu suất ζi < 1. Một ví dụ của tập hợp của các đường công nhiệt phát quang tích phân (glow curve) (GC), như được tính bởi Levy dùng phương trình (16) với u = 3 và υ = 1, được biểu diễn ở hình 1. Hình này minh hoạ các dạng đường thế có thể chấp nhận được tương tác 5 Hình 1: Tính toán các đường TL dùng tương tác đối với ba bẫy điện tử và một tâm tái hợp. Các đường liền đạt được với nồng độ điện tích bẫy ban đầu trong mỗi bẫy thay đổi từ 1× 1020 đến 5× 1021m−3 và σni σmnj . Cũng diễn tả đường TL bậc nhất (÷) và đường TL bậc hai (∗) đối với mật độ điện tích cư trú bẫy ban đầu 5× 1021m−3. hệ thống và những sự lệch sạch có thể đạt được liên quan với đường chuẩn bậc nhất hay bậc hai. Đối với số liệu được minh hoạ ta có σni σmnj = 0.1 và Nt = 10 22m−3. Đối với toàn bộ số liệu này có thể thấy rằng đường TL tương tác là giống dạng bậc hai cho trường hợp các giá trị ban đầu thấp của nồng độ lấp đầy trên các bẫy nhưng trở thành giống dạng bậc nhất hơn khi các bẫy được lấp đầy. Thay thế cách tiếp cận là không giả thuyết chuẩn cân bằng. Mặc dù ta sẽ không bàn về sự phân nhánh của tính gần đúng chuẩn cân bằng, nó có ích tại đây để nói rõ ra các tính chất của tương tác động học với sự hạn chế này được loại bỏ. Ngoài gần đúng chuẩn cân bằng các phương trình tốc độ (13) - (15)không thể làm đứt quảng và vì vậy không thể giải nghiệm bằng phân tích như đã làm ở trên. Thay vì phải xem xét nhiều dạng đường GC này ta có thể dùng phép phân tích số đối với các phương trình tốc độ. Một trong các áp dụng này được tiến hành bởi Bull, người mà dùng phép tính số để giải các phương trình tốc độ đối với nhiều bẫy điện tử và tâm tái hợp, dùng lại u = 3 và υ = 1. Số liệu 6 Hình 2: (A) Đường TL tương tác đối với trường hợp tái bắt chậm, và tổng của ba đường TL bậc nhất (o). (B) Đường TL tương tác được tính với sự bằng nhau của quá trình tái bắt và tái hợp. Tỉ trọng của các bẫy lấp đầy giãm từ (a) 1.0; (b) 0.3 đến (c) 0.1. Tổng của ba đường TL bậc hai được diễn tả bằng đường (o) của ví dụ này được biểu diễn ở hình 2. Ở đây ta thấy rằng khi quá trình tái hợp trội hơn (hình a) đường GC tương tác có thể được trình bày như là tổng của ba đường bậc nhất. Đây là điều được mong đợi từ trước, nếu quá trình tái bắt là chậm, tương tác nhỏ (bẫy đi vào các bẫy khác) sẽ được đối chiếu với quá trình tái hợp. Đối với trường hợp khi mà quá trình tái bắt và tái hợp như nhau, tuy nhiên đường GC tương tác là khác nhau đạt được từ tổng của các đường bậc hai. Hình 2 a, n0i = Nt = 10 15m−3, và σni σmnj = 1, đối với hình 2 b σni σmnj = 1, Ni = 10 15m−3 và tỉ số n0i Ni là không trong quá trình bắt. Kết luận đạt được từ các phân tích trên là trong một hệ thống mà quá trình tái hợp chiếm ưu thế hơn quá trình tái bắt thì đường TL GC có thể được diễn tả chính xác bằng sự chồng chập của các đỉnh bậc một kiểu Randall - Wilkin. Tuy nhiên, đối với các hệ mà quá trình tái bắt giống tái hợp thì chắc chắn dạng GC sẽ diễn tả sự sai lệch quan trọng từ sự chồng chập của các quá trình bậc hai. Trong sự phân tích này các dữ liệu dựa trên các phương trình đơn giản chẳng hạn bậc nhất, bậc hai, bậc trung gian hoặc hỗn hợp, có thể phát hiện không chính xác. 7 4 Phân bố các bẫy liên tục. Vấn đề đã bàn ở mục trước giả sử rằng các độ sâu bẫy kết hợp với các trạng thái định xứ là đơn trị theo năng lượng. Trong khi đó điều này có thể hi vọng đúng cho các vật liệu có tính đơn tinh thể cao, các vật liệu có các sai hỏng đặc biệt các vật liệu vô định hình hay thuỷ tinh độ sâu bẫy kết hợp với các sai hỏng đặc biệt sẽ trải ra trên phạm vi nhiều giá trị. Trong các vật liệu sau này phổ của các sai hỏng lân cận làm tăng tín hiệu TL có thể diễn tả sự khác nhau một cách ngẫu nhiên trong sự phụ thuộc gốc gần nhất và các phụ thuộc độ dài gần nhất. Kết quả là năng lượng phát xạ (với các bẫy sâu) phụ thuộc vào nồng độ. Đối với các vật liệu vô định hình nó có độ rộng vùng cấm không rộng hơn độ rộng vùng cấm của chất rắn kết tinh. Thay vì xem xét khe di động ở bên trong có thể tồn tại một số có hạn các trạng thái ở gần mức Fermi. Cho tạp chất của mạng tinh thể thay đổi có thể hi vọng vật liệu sai hỏng hay vật liệu vô định hình, nhiều loại phân bố năng lượng đã được diễn tả trong các trạng thái lân cận. Bao gồm một trong các dạng sau: (1)Phân bố đều N(E EB) và N(EA < E < EB) = Nt = sự đông đúc bẫy đồng đều (tính theo đơn vị m3eV −1) giữa EA và EB. (2) Phân bố tuyến tính dạng N(E) = Nm [ E − EA EB − EA ] , (20) với N(E) = 0 tại EA và Nm là nồng độ bẫy tích nạp cực đại tại năng lượng EB. (3)Phân bố theo hàm e mũ N(E) = Nc exp {−Et kTc } (21) với Ne là một hằng số (m 3eV −1) và Tc thông số nhiệt độ cho phân bố. (4) Phân bố Gaussian N(E) = N0 exp{−a(E − E0)2}, (22) với N0 là mật độ lớn nhất tại năng lượng E0 và a là hằng số. TL từ phân bố của các bẫy sâu đã xảy ra với nhiều vật liệu cả vô định hình và kết tinh chẳng hạn ZnIn2S4, polistilen, ZnSiO4 : Mn, ZnS, đioxit silíc và thạch anh.Tái hợp giữa các mật độ trạng thái đã được đề cập bởi Fowler và Randall và Wilkins, bàn về 8 tính lân quang và TL. Hornyak và các đồng nghiệp đã bàn về đặc trưng của đường tắt lân quang và TL GC như là giả thuyết của Gaussian hay phân bố bẫy như nhau. Theo các tác giả sau này ta sẽ xem xét các phương trình tốc độ và TL GC kết quả từ sự lọc lựa của các phân bố năng lượng. Biểu thức chung cho tốc độ thay đổi của điện tử tự do trong vùng dẫn đã được cho từ trước. Ta quan tâm ở đây chỉ là nồng độ của các bẫy năng lượng và vì vậy ta thay N(E) cho N(Et). Hơn nữa, ta chỉ quan tâm giải thoát do nhiệt của các điện tử bị bẫy tái hợp với các bẫy lỗ trống, và ta viết lại tích phân ∫ Ec EDp dE → ∫∞ 0 dEt (được gọi là điện tử bị bẫy, được định nghĩa, chỉ tồn tại giữa EDp và Ec. Chúng ta có thể viết được phương trình như sau dnc dt = ∞∫ 0 pn(Et)N(Et)f(Et)dEt − ncAn ∞∫ 0 N(Et)(1− f(Et))dEt − ncmAmn (23) Trong cách viết (22) ta có thể giả thuyết rằng tiết diện bắt không phụ thuộc vào năng lượng. Quan tâm đến mật độ các điện tử bị bẫy n(Et) = f(Et)N(Et) (24) Quan tâm sâu hơn ta phân ra các phần dEt của phân bố N(Et), đây là một phần rất nhỏ chứa n(Et)dEt điện tử tại thời điểm t. Toàn bộ số điện tử bị bẫy tại thời điểm t là n = ∞∫ 0 f(Et)N(Et)dEt (25) và toàn bộ các bẫy có sẳn là N = ∞∫ 0 N(Et)dEt (26) tại thời điểm t = 0 n0 = ∞∫ 0 f0(Et)N(Et)dEt (27) Quan sát rằng tại t = 0, n0 = m0. Vì vậy ta cũng có thể viết dn(Et) dt dEt = N(Et) df(Et) dt dEt = −pn(Et)f(Et)N(Et)dEt + ncAn[1− f(Et)]N(Et)dEt (28) 9 từ đây ta thấy rằng df(Et) dt = −pn(Et)f(Et) + ncAn[1− f(Et)] (29) Lấy tích phân phương trình (29) ngay lập tức ta được sự phụ thuộc vào thời gian của hàm lấp đầy f(Et) = f0(Et) exp{−pn(Et)}+ An exp{−pn(Et)} t∫ 0 nc(t ′) exp{pn(Et)}dt′ (30) Bằng cách đưa vào giả thuyết chuẩn cân bằng, nó tương ứng giả thuyết rằng n(t′) là hàm biến thiên rất chậm và có thể xem như hằng số. Trong trường hợp này phương trình (30) trở thành f(Et) ' f0(Et) exp{−pn(Et)}+ ncAn pn(Et) [1− exp{−pn(Et)}] (31) Bằng cách giả thuyết ncAn ¿ pn(Et) ta có thể ước lược f(Et) ' f0(Et) exp{−pn(Et)} (32) Tiến hành với giả thuyết chuẩn cân bằng, có thể nhận được phương trình tương tự phương trình bậc tổng quát đối với hệ này, tức là ITL = ∞∫ 0 pn(Et)N(Et)f(Et)dEt [ mσmn (N − n)σn +mσmn ] (33) với N và n được định nghĩa như ở phương trình (25) và (26). Việc đưa vào tính gần đúng tái bắt chậm (bậc nhất) đưa đến phương trình tương tự dạng Randall - Wilkins, do vậy ITL = ∞∫ 0 N(Et)f0(Et)s. exp { −Et kT } exp − ( s β ) T∫ T0 exp{−Et kθ }dθ  dEt (34) với β là hệ số gia nhiệt. Để tiến hành sâu hơn ta phải đưa vào các hàm phân bố rõ ràng N(Et). Hornyak và Chen đề cập đến phân bố bẫy như nhau, dạng (34) dễ dàng trở thành ITL = n0s ∆E EB∫ EA exp { −Et kT } exp − sβ T∫ T0 exp{−Et kθ }dθ  dEt (35) với n0 là nồng độ toàn phần của các điện tử bị bắt vào thời điểm t = 0, và ∆E = EB−EA là độ rộng phân bố. Dùng chuổi tiện cận để tính tích phân theo nhiệt độ T∫ exp { −Et kθ } dθ ' T exp { −Et kT }∑ α=1 ( kT Et )α (−1)α−1α! (36) 10 phương trình (35) có thể gần bằng ITL ' n0s ∆E EB∫ EA exp { −Et kT − skT 2 βEt exp{−Et kT } [ 1− 2kT Et + 6k2T 2 E2t − · · · ]} dEt (37) sau khi dùng phép gần đúng skT 2 βEt [ 1− 2kT Et + 6k2T 2 E2t − · · · ] ' skT 2 βE0 [ 1− 2kT E0 + 6k2T 2 E20 − · · · ] = γ (38) phương trình (37) có thể lấy tích phân cho ITL = n0s ∆E kT γ [ exp { −γ exp ( −EB kT )} − exp { −γ exp ( −EA kT )}] (39) với E0 = EA+EB 2 là năng lượng trung tâm của phân bố. Kết quả gần đúng này là tốt nếu ∆E là nhỏ (≤ 0.1eV ). Chẳng hạn đường GC được tính dùng (39) được biểu diễn ở hình 3. Có thể thấy rằng độ rộng nồng độ, độ rộng của đỉnh TL. Đối với ∆E nhỏ dạng đỉnh chủ yếu là đường bậc nhất đối với năng lượng E = 0.7eV . Tuy nhiên, độ rộng mật độ tăng thì đỉnh trở nên rộng và nhỏ hơn nhận dạng đỉnh bậc hai, nhưng lưu ý rằng phải cận thận trong việc phân biệt hai loại đỉnh này. Phương pháp để phân biệt một phân bố và đỉnh bậc hai là nồng độ điện tử bị bẫy ban đầu n0. Các tác giả khác xử lí kiểu giống của phân bố bằng cách chia nhỏ phân bố thành 50 đến 100 khoảng nhỏ, để tính GC cho mỗi khoảng ta dùng phương trình chuẩn Randall - Wilkins và rồi thêm vào kết quả để đạt được đường GC cần thiết. Cách khảo sát gần đúng này thời gian phá huỷ nhanh hơn dùng gần đúng (39), nhưng đúng trên khoảng rộng của giá trị ∆E. Các hàm phân bố khác N(Nt) là khó phân tích. Với các hàm này có một phương pháp tương tự như trên gọi là tổng nhân GC để đạt được đường GC cần thiết. Chẳng hạn dạng đường có thể dùng cho phương pháp này là phương pháp Gaussian và được biểu diễn ở hình 4 Ban đầu của việc phân tích trên được diễn đạt rằng ta quan tâm chỉ trường hợp phân bố theo các năng lượng bẫy và năng lượng các tâm tái hợp được giả thuyết gián đoạn và tiết diện bắt chỉ quan tâm các giá trị đơn, σn (và σmn) cả hai không phụ thuộc vào năng lượng. Một phân bố theo năng lượng tâm tái hợp không cần thiết để trình bày lại các phương trình vì tất cả nó đưa đến phân bố rộng của năng lượng phát xạ, chẳng hạn độ rộng phát xạ TL có thể đạt kết quả. Trong ý nghĩa này ta có thể đề cập đến thông số m trong các phương trình trên cho đến toàn bộ các tâm tái hợp có sẳn theo trường hợp 11 Hình 3: Các đĩnh của đường TL tích phân đạt được bằng cách giả sử phân bố đều của các bẫy cho trường hợp E0 = 0.7eV, s = 1× 1012s−1, β = 10K/s, và ∆E = 0.1, 0.5 và 0.01eV . Các đĩnh được chuẩn hoá cho cường độ chuẫn bằng 1 đối với đĩnh ∆E = 0.01eV các phương trình chính nó còn lại như nhau. Tuy nhiên, phân bố của thiết diện tái hợp m(σmn) (hay m(Amn)) có thể đưa ra theo dạng mong muốn vì các tâm này với xác suất tái hợp lớn nhất sẽ tái hợp phát xạ bậc nhất, và phát xạ TL với phổ phát xạ tương ứng sẽ xuất hiện tại nhiệt độ thấp hơn phát xạ tương ứng với thiết diện nhỏ hơn. Ở đây, nếu TL được ghi thông qua bộ phận phát xạ lộc lựa, thì đỉnh TL sẽ xuất hiện tại các nhiệt độ khác nhỏ không đáng kể, theo bộ lọc đã dùng. Phân bố theo thiết diện bẫy đối với bẫy n(σn) hay n(An), tuy nhiên tiến hành phân bố theo hệ số tần số n(s). Phân bố theo s đã được tiến hành với một số tác giả. Rudlof đã làm tương tự dạng đường TL dùng phương pháp như đã dùng ở trên nhưng đối với phân bố Gaussian theo s. Giả thuyết nếu phân bố Γ theo s và áp dụng điều này đối với trường hợp động học bậc một và giá trị đơn của Et có thể đạt được các phương trình động học trung gian. Xem xét hàm phân bố theo tần số g(ln s). Hàm này được chuẩn hoá như sau +∞∫ −∞ g(ln s)d(ln s) = 1 và là một hàm theo thời gian. Tại thời điểm t = 0 ta chỉ rõ giá trị g0(lnS), giả thuyết rằng hệ số gia nhiệt β là hằng số, thì tại thời điểm t ta có thể viết, 12 giả thiết động học bậc một n = n0 ∞∫ 0 g0(ln s) exp − sβ T∫ T0 exp{−Et kθ }dθ  dss (40) từ đó ta có ITL = n0 exp{−Et kT } ∞∫ 0 g(ln s) exp − sβ T∫ T0 exp{−Et kθ }dθ  dss (41) Phương trình bậc động học trung gian có thể được viết ITL = s ′′n0 exp{−Et kT } 1 + ((b− 1)s′′ β ) T∫ T0 exp{−Et kθ }dθ − b b−1 Hàm phân bố g0(ln s) có đỉnh động học bậc một trong giống đỉnh động học trung gian g0(ln s) = (a0s) 1 b−1 Γ( 1 b−1) exp{−a0s} (42) với a0 = 1 (b− 1)s′nb−10 (43) Phương trình (42) là hàm phân bố Γ theo ln s và trở thành hàm đenta khi b = 1. Hàm (42)thay đổi theo nhiệt độ trong suốt quá trình TL đọc ra, nhưng nó quay lại như dạng hàm Γ, nhưng thông số a thay cho thông số a0 a = a0 + 1 β T∫ T0 exp{−Et kθ }dθ (44) Chú ý phân bố là hàm của ln s rộng hơn s. Tuy nhiên như được chỉ ra bởiChristodoulides, hàm Γ theo ln s cũng là hàm Γ theo s. Phân bố theo s phải được phát triển từ phân bố tương tự theo σn. Tuy vậy, quá trình vật lý mà cho phân bố Γ tăng lên thì không biết. Mặc dù nó được bàn luận trong mục trước rằng, trong suốt quá trình ghi TL, hàm phân bố theo s chắc chắn là một hàm Γ, sự giống nhau này không thể lý giải phân bố theo Et. Nếu độ cư trú của các trạng thái khả dụng là như nhau tại thời điểm ban đầu của chu kì nhiệt, nó không như nhau trong suốt các nguồn nhiệt. Tái bắt điện tử từ bẫy Et thấp đến bẫy Et cao hơn có thể xảy ra, với kết quả này sự thay đổi động học của hàm phân bố xảy ra trong suốt quá trình nhiệt. Nói cách khác, tương tác động học trở nên quan trọng - nhưng thời gian tương tác ở giữa các trạng thái ngoài phân bố, tuy nhiên nếu tái bắt yếu thì cho đường động học bậc một - kết quả là nhỏ và thay đổi của phân bố trong quá trình nhiệt không cần xét đến. 13 5 Kết luận Cùng với các phương pháp khác để nghiên cứu tính chất vật rắn thì phổ TL cho chúng ta các thông tin cần thiết về vật liệu. Đây cũng coi như một chìa khoá quan trọng để xác định tính chất, nghiên cứu vật rắn và vật liệu hiện nay. Để dùng được TL thì các vật liệu xem xét phải có tính chất TL. Cùng với các phổ khác nữa thì phổ PL hay đường cong GC cho ta một cái nhìn khá chi tiết về vật liệu mà chúng ta nghiên cứu. Do thời gian nghiên cứu ngắn nên tiểu luận không tránh được một số sai sót xin chân thành cảm ơn quý thầy cô và bạn đọc để tiểu luận hoàn thiện hơn. Xin chân thành cảm ơn TS Nguyễn Mạnh Sơn, người đã trực tiếp hướng dẫn chúng tôi học phần này. Xin chân thành cảm ơn học viên Bùi Tiến Đạt đã có nhiều đống góp giúp cho tôi hoàn thành bài viết này. Huế, tháng 10 năm 2007 Học Viên: Thái Ngọc Ánh Chuyên ngành: Quang học Tài liệu [1] Đặng Mộng Lân, Ngô Quốc Quỳnh, " Từ điển Vật lý Anh - Việt", Nhà Xuất bản khoa học và kỹ thuật, Hà Nội, 1976 [2] Reuven Chen, Stephen W. S. McKeever, "Theory of Thermoluminescence and Related Phenomend" 14 Hình 4: Sự so sánh hình dạng các đường GC đạt được từ dạng bạc nhất, sự phân bố độ sâu các bẫy đều (a) và sự phân bố độ sâu bẫy Gaussian (b). Các sự phân bố được dùng để biểu diễn (1); sự thay đổi độ sâu bẫy Et(T ) theo nhiệt độ thông qua đĩnh GC. 15