Giáo trình hóa tin học

TRÂN V Ĩ N H Q GIÁO TRÌNH Hóa NHÀ XUẤT BẢN ĐẠI HỌC s PHẠM PGS TS TRẦN VĨNH Q GIÁO TRÌNH HỐ TIN HỌC ■ (Các tốn nhiệt động, thống kê lí thuyết phản ứng hố học) (Tái lần thứ nhất, có sủa chữa) DẠI HỤCTHÁỊ_ NGUYỀN TRUNG TẤiM H Ọ CLIỆt' NHÀ XUẤT BẢN ĐẠI HỌC s PHẠM Mục lục Lời nói đẩu Chương Mở đầu NGÔN NGỮ LẬP TRÌNH PASCAL VÀ CÁC THUẬT TỐN TÌM NGHIỆM Cơ BÀN DÙNG TRONG HOÁ HỌC 0.1 Ngôn ngữ Turbo Pascal 7.0 Khởi động Tu rbo Pascal Làm việc m ôi trường Turbo P ascal Kí hiệu đại lượng Turbo P a sca l 10 Biểu thức hàm mẫu 11 Cấu trúc chương trinh Turbo Pascal 13 Các toán tử 15 Biến đếm , m ảng, chu trìn h 18 Giải thích 20 Case of 21 10 Các chương trình : PROCEDURE (Thủ tục) FUNCTION (hàm) 21 11 Kiểu liệu có cấu trúc (R E C O R D ) 27 12 Kiểu liệu có cấu trúc: tệp (F IL E ) 28 0.2 Các phương pháp tìm nghiệm phương trình đại sơ' 30 Phương pháp dây c u n g 30 Phương pháp tiếp tuyến - Phương pháp N e w to n 32 Phương pháp kết h ợ p 35 Phương pháp chia đôi khoảng chứa n g h iệ m .37 Phương pháp lặp .40 Phương pháp bước lặp giảm dần 41 Tài liệu tham khảo chương MỞ đầu 4 Chương CÁC PHƯƠNG PHÁP NHIỆT ĐỘNG XÁC ĐỊNH HẰNG SÔ CÂN BẰNG HOÁ HỌC 45 / Phương p háp phương trình N e rnst Các b i tập ứng dụng phương pháp phương trinh N e rn st 60 Các b i toán tính n g ợ c 63 1.2 Phương p háp Tem kin - Schw arsm an C ác tập ứng dụng phương pháp Temkin - S ch w a rsm a n 96 Các b i toán ngư ợ c 9 I.3 Phương pháp th ế p h i s a o 101 Các tập ứng dụng phương pháp thê phi s a o ??6 Các toán tinh ngư ợ c Tài liệu tham khảo chương 121 Chương PHƯƠNG PHÁP HỌC THỐNG KẺ XÁC ĐỊNH HẰNG SỐ CÂN BẰNG HOÁ HỌ C 123 II.1 Sơ lược li thuyết Co học thống k ê 123 Hàm phân bố hè đoan n h iê t 125 Hàm phân bô’ hệ đẳng n h iệ t 125 11.2 Áp dụng H thuyết Co học thống kê tinh hàng số cản bàng hoá h ọ c 128 Hoá theo Cơ học thống k ê 129 Tinh sô’ cân b ằ n g 130 11.3 Tinh tổng trạng thái (tổng thống k ê ) 132 Đôi với đơn nguyên tử .133 Tổng thống kê cùa phân tử hai nguyên tử 136 Đối với phàn tử nhiéu nguyên tử (từ ba nguyên từ trờ lê n ) 149 II.4 ứng dụng: Tinh hàng só cân bàng phản ứng chất k h í 151 Tinh sơ' cân phản ứng phân li nguyên tử 151 Tính sơ cân cùa phàn ứng không phản li ng uyên t .199 Các tập chương 207 C ác toán tinh n g ợ c .210 Các loàn loại k h c 272 Tài liệ u th a m k h ả o c h n g 213 Chương NGHIÊN cứu TÍNH CHẤT HỐ LÍ HỆ PHÀN ỨNG CÂN BẰNG NHIỆT ĐỘ CAO .215 III Li thuyết hệ cân b àng 215 Ba định luật chi phối hệ phàn ứng cân b ằ n g .215 Nhặn x é t 219 111.2 Phương pháp N ew ton 220 Khai triển T a y lo r 220 Giải toán nghiên cứu cân nhièt đô c a o 223 111.3 Phương pháp Newton cải tiến 231 Thiết ỉảp hệ phương trinh xác đinh nồng đô cấu tử 231 NÔI dung cùa phương phap Newton cài tiên 232 Bài tốn cụ thể chương trình giải theo phương pháp Newton cải tiến 235 B ài to n cụ th ể g iả i b ằ n g hai phương ph áp N ew to n v N e w to n c ải tiế n đ ể s o s n h 239 III Cân bàng ion dung d ịc h 253 Các tập ứng dụng chương 261 Tài liệu tham khảo chương 264 Chương LÍ THUYẾT HỒI QUY 265 IV N ội dung phương p h p 265 N guyên lí tổng binh phương độ lệch cực tiể u 265 Phương pháp xác định hệ sơ' phương trình 265 IV Phép khử G a u ss 268 Phương pháp tổng q u t 268 Sơ đồ G auss bó g ọ n 271 IV Các toán hoá học sử dụng li thuyết hồi q u y 273 Tinh hệ s ố phương trình IgKp = F(T) xuất phát từ liệu thực n g h iệ m .273 Chương trình 274 C ác ứng dụng khác lí thuyết hồi q u y 276 Tài liệu tham khảo chương 286 Chương PHƯƠNG PHÁP GIẢI HỆ PHƯƠNG TRÌNH VI PHÂN ĐỘNG HOÁ HỌC XÁC ĐỊNH Cơ CHẾ PHÀN ỨNG HOÁ HỌC 287 V.1 Phương pháp tích phân s ố hệ phương trình vi phân mơ tả ca c h ế phản ứng hoá học phức tạp 287 Nội dung phương p h p 289 Đ ánh giá độ xác phương pháp g iả i 294 V.2 Tinh tốn động học q trình clo hố etilen phân tích độ chinh xác thuật to n 296 Hệ phản ứng hoá học c h ế 296 Hệ phương trình vi phân m õ tả chê phản ứng clo hoá e tile n 297 Phương pháp giải hệ phương trình vi phân mô tả ch ế phản ứng clo hoá e tile n 301 Đ ánh giá độ xá c kết tích phân s ố hệ phương trình (5 ) 303 Giải gấn phương trình vi phân thường 307 Chương trình, kết cá c đố thị phép tích phân s ố hệ (5 ) 309 B ài tập ứng dụng chương 316 Tài liệu tham khảo chương 322 Chương ĐỘNG HỌC PHẢN ỨNG ĐƠN PHÂN TỬ LÍ THUYẾT RRKM 323 Mỏ đầu 323 VI Co sỏ li thuyết R R K M 323 Mơ hình lí thuyết R R K M 323 Biểu thức S k 1(E.đE tỗE )/ k 325 Biểu thức ka( E * ) 326 Biểu thức ku„, the o lí thuyết R R K M 330 Giới hạn áp suất c a o 330 Giới hạn áp suất th ấ p 332 Những thừa sô thống kê 332 Số trạng t h i 333 Mật độ trạng thái lượng tử N (E ) 333 VI.2 Các phương pháp tinh s ố trạng thái lượng tử .334 Sự tinh toán trực tiếp trạng thái dao đ ộ n g 334 Phép gẩn VVhitten Rabinovitch 337 VI.3 Ap dụng li thuyết RRKM nghiên cứu động học phản ứng đơn phân tử 339 Các đại lượng cẩn tính nghiên cứu phàn ứng đơn phân tử 339 Áp dụng: Tính s ố tốc độ phản ứng đồng phân hoá m etyl is o x ia n u a 343 C ác chương trinh m áy tính áp dụng cho tốn nghiên cứu động học phản ứng phân hoá nhiệt metyl is o x ia n u a 348 Bài tập 355 Tài liệu tham khảo .358 Lời nói đẩu Hố tin học - mơn khoa học liên ngành ứng dụng phương pháp tính tốn, xử lí thơng tin vào việc nghiên cứu vấn đề hoá học, từ đời thúc đẩy phát triển mạnh mẽ lĩnh vực môn khoa học, lĩnh vực nghiên cứu giảng dạy lí thuyết Giáo trình Hố tin học (Các tốn nhiệt động, thống kê lí thuyết phản ứng hố học) trình bày bước phát triển phần lí thuyết hệ vĩ mơ hố học cách ứng dụng tin học vào trợ giúp lí thuyết vượt qua trở ngại tiến Phần cốt lõi, tảng giáo trình áp dụng Cơ học thống kê vào vấn đề lí thuyết hố học (lí thuyết cân hố học lí thuyết phản ứng hố học) Giáo trình tồng kết nội dung nghiên cứu tập lớn, luận văn sinh viên, học viên cao học, Nghiên cứu sinh (lĩnh vực lí thuyết hệ vĩ mơ) năm vừa qua ỏ khoa Hoá học, trường Đại học Sư phạm Hà Nội Nội dung giáo trình xếp theo trình tự phát triển vấn đề bắt đầu từ: + Ngơn ngữ lập trình Pascal thuật tốn tìm nghiệm dùng hố học - chương Mở đẩu, tiếp đến chương: + Các phương pháp Nhiệt động xác định số cân hoá học chương + Phương pháp Cơ học thống kê xác định số cân hoá học chương + Nghiên cứu hệ cân nhiệt độ cao - chương + Lí thuyết hổi quy - chương + Phương pháp giải hệ phương trình vi phân động hoá học nghiên cứu chế phản ứng - chương Và kết thúc chương 6: Động học phản ứng đơn phân tử - lí thuyết RRKM Giáo trình trình bày lí thuyết hố học cách đơn giản hệ thống chặt chê Người đọc sinh viên học viên cao học Hoá học nắm lí thuyết hố học nói chung, đặc biệt, Cơ học thống kê hoá học cách chắn vận dụng nghiên cứu ỏ bậc học tập cao hơn, theo dõi cách chương giáo trinh Ngôn ngữ lập trình sách ngơn ngữ Pascal, ngơn ngữ mà sinh viên Hoá học học Các chương giáo trinh có tập để người đọc có Ihể vận dụng vấn đề lí thuyết vào giải toán cụ thể tất mơn Hố học Các tập với nội dung sô' lượng phong phú giải cẩn thận nghiêm túc Tuy vậy, giáo trình đưa đáp số ngắn gọn để người đọc kiểm tra tính tốn Phần chương trình giải tập đưa vào đĩa CD kèm theo sách (nếu có đủ điều kiện làm điều đó) Một Bài tặp Hố tin học kèm theo giáo trình bổ ích thời gian tới Tác giả chân thành cảm ơn Thạc sĩ Ngô Tuấn Cường sửa chữa đánh máy thảo cách nghiêm túc, Thạc sĩ Nguyễn Văn Thạc sĩ Nguyễn Đình Độ ưng thuận để tác giả sử dụng chương chủ yếu luận văn Cao học làm nội dung để viết Chương Mở đầu, Chương Chương Những đóng góp Thạc sĩ giúp tác giả hoàn thành sách điều kiện thời gian vật chất eo hẹp Cuối cùng, tác giả chúc bạn đọc thành công việc đọc sách ứng dụng vào việc nghiên cứu Tác giả thiết tha mong bạn đọc góp ý thẳng thắn cho thiếu sót khơng thể tránh khỏi nội dung hình thức sách để phục vụ bạn cách đắc lực tương lai Xin chân thành cảm ơn T Á C G IẢ Chương Mở đẩu NGƠN NGỮ LẬP TRÌNH PASCAL VÀ CÁC THUẬT TỐN TÌM NGHIỆM Cơ BẢN DÙNG TRONG HỐ HỌC 0.1 NGƠN NGỮ TURBO PASCAL 7.0 Khởi động T u rb o Pascal 7.0 Đế dù ng Turbo P ascal (TP) phải có n h ấ t tệp: TURBO.EXE TURBO.TPL H tệp lưu trữ đĩa mềm đĩa cứng Đế khởi động TP, bạn vào thư mục có chứa tệp TURBO.EXE sa u dù ng lệnh: TURBO.EXE J S au m àn h ìn h soạn thảo TP x u ất hiện, bạn có th ể thực việc soạn th ả o chạy chương trình Tổ hợp phím Alt+X: Ra khỏi chương tr ìn h TP Làm việc mơi trường Turbo Pascal Soạn thảo TP: Nói chung việc soạn th ả o TP giông n h hệ soạn thảo khác (Word, Notepad ) Phím ( E n te r J): d ù n g để tạo dịng Bốn phím m ũi tên (< -í- > ị) : dùng đê dịch chuyển trỏ Các phím PgDn, PgUp dùng để lậ t trang Phím Home: đế đưa trỏ vê' đầu dịng Phím End: đưa trỏ cuối dịng Phím In sert: đê th a y đôi chê độ chèn đè Phím Delete: để xố kí tự vị trí trỏ Phím («- Backspace): đế xố kí tự bên tr i trỏ Tổ hợp phím Ctrl+K B: đế đ ánh dấu đầu khối Tơ hợp phím C trl+K K: để đ n h dấu cuối khối Tơ hợp phím C trl+K C: đề chép khôi đ n h dấu Tổ hợp phím Ctrl+K Y: đẻ xố khơi đ n h dấu Tố hợp phím Ctrl+K V: để chuyển khối chép tới vị trí trỏ - Một sơ p h ím chức năng: F l: Hướng dẫn F2: Ghi tệp lên đĩa F3: Mỏ tệp F5, F : Chuyển cửa sổ mở F9: Soát lỗi ngữ pháp Tổ hợp phím Ctrl+F9: Chạy thử chương trình Tơ hợp phím Ctrl+Pause: Thốt khỏi chương trình chạy thử Tố’ hợp phím AU+F3: T hốt khỏi sổ thời Kí hiệu dại lượng Turbo Pascal í ề K i hiệu C hữ cái: TP dù ng chữ tiếng Anh viết hoa: A B c .z viết thường: a b c z chữ sô 0, .9 Dâu: + - */ : Các phép tốn sơ học: cộng, trừ, n h ân , chia = > > = < < = : Các phép toán q uan hệ: bằng, lớn hơn, không nhỏ nhỏ không lỏn khác : Phép gán : Chấm th ập phân, dấu hết chương trình : : Kêt thúc tốn tử mơ tả : Mơ tả : D ấu nháy dùng cho h n g v ăn Từkhoá : AND DO GOTO FUXCTION\ 3.2 Các đại lượng H ằng (const) Có hai loại hãng: h ả n g sô h ằn g v ăn bàn Các h ằ n g sô biểu diễn d ạng nguyên thực Ví dụ: 1076 123 h ằ n g nguyên 10 Qua bảng này, ta th ấ y chênh lệch E , từ miền áp s u ấ t cao đến miền áp su ấ t th ấ p vào khoảng 2kcal/mol, giá tr ị chấp nhận E„ 38,4kcal/mol Trên sỏ thông số trên, áp d ụng lí t h u y ế t RRKM, ta tính số tốc độ p h ả n ứng đồng p hân hoá metyl isoxianua theo phương pháp tính trực tiếp phương pháp tính gần W itten - Rabinovitch hai nhiệt độ: 230, 4°c 199, 4°c Năng lượng tới hạn E0 tính nhiệt độ 230°c dùng chung cho n h iệ t độ khác Kết cuối thông số n h iệ t động, h ằ n g số tốc độ phản ứng hai n h iệ t độ tr ê n kê bảng sau: N h iệ t d ộ 9 ,4 °c N h iệ t đ ộ 230,4° C ,2 9 ' 58,1 0' 3,816 2.408 4,330 2,563 2,419 2,663 Tổng thống kê dao động phức tính the o phương pháp trực tiế p (Q23) 2,179 2,429 1,797 2,062 1,365 1,587 S.(A) 6,465 7,077 S,(A*) 4,635 5,091 AS, 0,353 0,226 Hằng số tốc độ ỏ áp suất cao the o thực nghiệm (k „,) ,5 '5 9,25.1(T5 Hằng số tốc độ ỏ áp suất cao tính the o phương pháp 7,82.10'5 9.68.10 7.86.10"5 9.71.10 7.08.1Ũ'5 8.86.10'5 C ác th ô n g s ô ' z Tổng thống kê dao động phân tử ban đầu (Q2) Tổng thống kê dao động phức tính the o phương pháp giải tích (Q 21) ^'ổng thống kê dao động phức tính the o phương pháp V Vhitten-R abinovitch (Q22) giải tích (k l2 ) Hằng số tốc độ áp suất cao tinh the o phương pháp trực tiếp (k 13) Hằng số tốc độ áp suất cao tính the o phương pháp W h itte n -R a b in o v itc h (k t4) Năng lượng tới hạn (E 0) Các kết n ày cho thấy, giá trị k 37,845 kcal/m o l tín h theo phương pháp khác n h a u k h ằ p h ù hợp với n h a u với giá trị thực nghiệm 34Ỉ Sự phụ thuộc h ằng số tốc độ p h ả n ứng (kunj) vào p s u ấ t (p) (hay phụ thuộc lg(kuni / k r ) vào lgp n h iệ t độ 199 4°c n h sau: Igp lg(k„™/k.) thực nghiêm Ig (K J k x ) lg (K ,„ /M ‘ heo pp W hitten tình trực - Rabinovitch tiếp Igp thực nghiệm tinh trực tiếp Ig ( iO M theo pp W hitten - Rabinovitch -0,7932 -2 ,872 -2,262 -2 ,190 0,8603 -0 ,723 -0 ,7 -0 ,809 -0,7645 -1.9281 -2,262 -2,161 0.9576 -0,655 -0 7 -0.6392 -0,7212 -1.8633 -2,191 -2 ,119 1.0000 -0,619 -0 ,6 7 -0,691 -0,6144 -1.8508 -2 ,085 -2,0141 1,1038 -0 ,578 -0.5791 -0.5469 -0,5272 -1 ,742 -2 000 -1,9287 1,1959 -0 540 -0,5221 -0,493 -0,4989 -1.705 1,9728 -1,9011 1,3945 -0 ,3 9 -0 ,4 2 -0.390 -0,4908 -1 ,714 -1 ,964 -1,8931 1,5527 -0 ,3 -0,337 -0,319 -0,4473 -1,653 -1,9221 -1 ,850 1,8591 -0 ,2 -0 ,2 -0,2101 -0 ,169 -0,3487 -1,556 -1,8257 -1,7551 2,0065 -0 ,159 -0 ,1 -0,3045 -1,560 -1,872 -1,712 2,0795 -0,128 - Í9 -0 ,1 -0,2815 -1,565 -1 ,7 -1,690 2,1909 -0 130 -0 ,1 -0 ,1 7 -0,2692 -1 ,504 -1 ,748 -1 ,678 2,3181 -0 095 -0 ,1 0 -0 ,104 -0,2328 -1 ,511 -1,7131 -1 ,643 2,5276 -0,054 -0 -0 ,0 -0,1884 -1,423 -1 ,6 -1,6011 2,7559 -0,0311 -0 -0 ,0 -0,1549 -1,398 -1 ,6 -1,569 2,9768 -0 ,043 -0 ,0 -0 035 -0,0942 -1,3841 -1 ,579 -1.5117 3,0580 -0,011 -0 ,0 -0 ,0 -0,0315 -1,3401 -1 ,519 -1,452 3.1761 -0 ,014 -0 ,0 -0 ,025 0.0128 -1 ,264 1,4778 -1 ,411 3,2504 -0 ,008 -0 ,0 4 -0,0231 0,0569 -1 ,307 -1 ,436 -1,3707 3,3518 -0 ,0 -0.0211 -0 ,0 0 0,0934 -1 ,251 -1 ,402 -1,3371 3,4606 -0,0031 -0 ,0 -0 ,0 0,1492 -1 ,2 -1,350 -1.2862 3,6998 -0 ,005 -0 ,0 -0 ,0 0,3075 -1 ,1 ,3 -1 ,205 -1,1451 3.7767 -0 ,000 -0 ,0 -0 0,4200 -1,0070 -1 ,103 -1,0128 3,8516 -0,000 -0,0121 -0 1 0,6990 -0 ,8 -0,8711 -0 ,886 3.9445 -0 ,012 -0 1 -0 010 0,7168 -0 754 -0 ,856 -0.8234 5.0000 0.0000 -0 ,007 - 0075 346 H inh 6.1 Đ ó th ị b iế u diễn s ự p h ụ th u ộ c lg (k un/ k j vào Igp n h iệ t đ ộ 199,4°c the o p h ng p h p tin h trụ c tiếp H ình 6.2 Đ õ th ị b iế u diễn s ự p h ụ th u ộ c lg (k un/ k j vào Igp n h iệ t đ ộ 199,4°c the o p h o n g p h p W h itte n - R a b in o v itc h : cá c đ iể m thực nghiệm : đường tính từ li th u y ế t R RKM H ình Đ ó th ị b iế u d iễ n s ự p h ụ th u ộ c lg (k u„ / k r ) vào Ig p n h iệ t đ ộ 199,4°c the o cà h a i p h o n g pháp 1: C c đ iể m thực n gh iệ m so sánh 2: Đ ường kết tính th e o phương p há p tính trực tiếp Đ ường kế t q u tin h th e o phư ong p há p kinh n g h iê m W h itte n -R a b in o v itc h 347 N h ậ n xét: Từ kết đồ thị th u cho th ấ y kêt q tính theo hai phương pháp phù hợp tương đơi tốt với thực nghiệm , đặc biệt ỏ miền áp s u ấ t cao Đôi với miên áp s u ấ t th ấ p k ết thu cịn sai sơ" Một ngun n h â n ta giả thiêt hiệu su ấ t va chạm \ = 0,6 miền áp suất Điều kh ông áp s u ấ t nhỏ (cỡ 0.1 torr), m ậ t độ p h â n tử th ấ p va chạm hoi Bên cạnh đó, lựa chọn mơ hình p hàn ứng từ rút tỉ sơ tổng thống kê quay đoạn nhiệt tr ìn h thực nghiệm đo đạc đêu mắc phải sai sô n h ấ t định Nói chu ng kết tìm phù hợp k h tốt với th ụ c nghiệm Điều cho th ấ y việc áp dụng lí thu y ế t RRKM đê nghiên cứu dộng học p hản ứng đồng p h ân hố nói riêng p h ả n ứng tro n g pha khí nói chung th u n hữ n g k ết k h ả quan Giá trị đ n g tin cậy kết n ày cho phép ta k h ẳ n g định r ằ n g phương pháp tôt đê nghiên cứu động học p h ả n ứng đơn p h â n tử Các chương trìn h m áy tính áp dụng cho to án ngh iên cứu động học phản ứng dồng phân hoá n h iệt m etyl /so xian ua 3.1 Chương trình tinh E0 PROGRAM TinhEo; uses crt.p rin ter: Const h= 9.537E -14: C=3E+10;R=0.001987{Kcal/mo]}; N A =6.023E+23;Rl=6.2362E+4;{cm3/torr}k=1.3805E-16; Var TBE.TBEc,EO.Evc:real; iiinteger: q:array[1 10] of integer: v:array[1 1000] of real: b:text: BEGIN Clrscr; Assign(b,'C :\w it e n doc'); Rewrite(b); 348 TK:=503; {TINH CHO PHAN TU HOA} TBE:=0; v[l]:=2998;q[l]:=3; v[2]:=2161;q[2]:=l; v[3]:=1443;q[3]:=3; v[4]:=1041;q[4]:=2; v[5]:=945;q[5]:=l; v[6]:=270; q[6]:=2; For i:=l to TBE:=TBE+q[i]*hc*v[i]/(exp(hc*v[i]/RT)-l); {TINH CHO PHUC HO AT DONG} TBEc:=0; v[l]:=2998;q[l]:=3; v[2]:=1990;q[2]: = l; v[3]:=1443;q[3]:=3; v[4]:=1041;q[4]:=2; v[5]:=600;q[5]:=l; v[6]:=270;q[6]:=l; For i:=l to TBEc:=TBEc+q[i]*hc*v[i]/(exp(hc*v[i]/RT)-l); EO:=Evc+TBE-TBEc-RT; Writeln('Nang luong toi h an EO - ,E0:2:3,' Kcal/mol'); Writeln(b,'Nang luong toi h a n EO =',E0:2:3,' Kcal/mol'); close(b); Readln; END 3.2 Chương trinh tính hàng sô'tốc độ theo phương pháp Whitten-Rabinovitch PROGRAM TinhTheoWR; uses c rt,printer; Const h=9.537E-14; C=3E+10;R=0.001987(iCcal/mol}; N A =6.023E+23;Rl=6.2362E+4;{cm3/torr(k=1.3805E-16; Var {Danh cho c h a t h o a t hoa} Sv,TBE,Q2,Ts,EZ,E,Ep,w ,beta,TB,TB2,M.T.T2,dw:real; {Danh cho phuc hoat dong} 349 Tsp,EZp,Epp,Tp.w p,betap,TBp,TB2p.Mp,T2p:real: Q2cl ,Q2c2.TBEc.Svp,dwp:real: ¡Dung chung} DelSv.TL,la,Z,hc,EO.Evc.kvcl ,kvc2:real; kuni,T K ,p,L nPE a.L nN E.RT ,d,K L Ir,L p:real; i,n,s,j,ST,l,STp:integer; q:arrav[1 10] of integer: v,N E vp,ka.P E P E l:array[1 1000] of real; b.text; BEGIN Clrscr; A ssign(b.'C:\w iten.doc'); Rewrite(b); la:=0.6: KL:=20.05:d:=4.5E-8:Evc:=38.4; Lp:=4.012: Ir:=3: hc:=h*c; Repeat Begin write('T :=');Readln(TK); RT:=R*TK;Z:=(NA*sqr(d)/Rl)*sqrt(16*pi*NA*k/KL)/sqrt(TK); WritelnCTINH CHO PHAN TU HOAT HOA’); M :=l: T2:=0: T s:= l; ST:=0:T:=0; Q2:=1;TBE:=0; v[l]:=2998;q[l]:=3: v[2]:=2161:q[2]:=l: v[3]:=1443;q[3]:=3; v[4]:=1041;q[4]:=2; v[5]:=945; q[5]: = l; v[6]:=270: q[6]:=2; For i:=l to Begin T:=T+v[i]*q[i]; T2:=T2+q[i]*v[i]*v[i]; M:=M*exp(q[i]*ln(hc*v[i])):ST:=ST+q[i]: Q2:=Q2*exp((-q[i])*ln(l-exp( hc*v[i]/RT))); TBE:=TBE+q[i]*hc*v[i]/(exp(hc*v[i]/RT)-l); End; Sv:= R * (ln (Q 2)+T B E /R T )* 10 00 : 350 Writeln('TINH CHO PHUC HOAT DONG'); Mp:=l; Tp:=0; T2p:=0; T sp:= l; STp:=0; Q 2cl:=l;TBEc:=0; v[l]:=2998;q[l]:=3; v[3]:=1443;q[3]:=3; v[5]:=600;q[5]:=l; v[2]:=1990;q[2]:=l; v[4]:=1041;q[4]:=2; v[6]:=270;q[6]:=l; For i:=l to Begin Q2cl:=Q2cl*exp((-q[i])*ln(l-exp(-hc*v[i]/(RT)))); TBEc:=TBEc+q[i]*hc*v[i]/(exp(hc*v[i]/RT)-l); Tp:=Tp+v[i]*q[i]; T2p:=T2p+q[i]*v[i]*v[i]; Mp:=Mp*exp(q[i]*ln(hc*v[i])); STp:=STp+q[i]; End; EO:=Evc+TBE-TBEc-RT; Svp:=R*(ln(Q2cl)+TBEc/RT)*1000; DelSv:=Svp-Sv+1000*R*ln(Ir); I******************************************************** I WritelnCTINH s o TRANG THAI CHO PTU HOAT HOA'); Writeln(b,'TINH s o TRANG THAI CHO PTU HOAT HOA’); For 1:=1 to ST Ts:=Ts*l; EZ:=T*hc/2;TB:=T/ST;TB2:=T2/ST; beta:=((ST-l)*TB2)/(ST*TB*TB); For j:=0 to 799 Begin E:=(j*0.05+0.025); Ep:=(E+E0)/EZ; Begin If (Ep