Phổ nhiệt phát quang (TL) là những đường cong phức tạp không tuân theo những phân bố thông thường mà theo mô hình bậc một, hai, tổng quát, một bẫy một tâm tái hợp hoặc mô hình trộn. Bài báo này đưa ra phương pháp sử dụng phần mềm Python để mô phỏng và làm khớp đường cong thực nghiệm của phổ TL theo các mô hình khác nhau. Phương pháp mô phỏng phổ TL dựa vào thông số bẫy hoặc thông số đỉnh phổ theo các phương trình động học. Việc xử lí và phân tích phổ TL tìm ra được các thông số đặc trưng của phổ TL của vật liệu như năng lượng bẫy, tần số thoát và thời gian sống ở bẫy. Kết quả cho thấy, phổ TL được mô phỏng và làm khớp phù hợp nhất với mô hình bậc tổng quát. Hệ số khớp phổ nhỏ cho thấy phổ TL mô phỏng và thực nghiệm là tương đồng nhau