ĐẠI HỌC HUẾ TRƯỜNG ĐẠI HỌC KHOA HỌC NGUYỄN THÀNH ĐƯỢC TÍNH TỐN CÂN BẰNG LỎNG-HƠI CỦA Ar, N2, Cl2, CO BẰNG PHƯƠNG PHÁP HĨA LƯỢNG TỬ VÀ MƠ PHỎNG TỒN CỤC MONTE CARLO Chun ngành: Hóa lý thuyết Hóa lý Mã số: 944.01.19 LUẬN ÁN TIẾN SĨ HĨA HỌC Người hướng dẫn khoa học: PGS.TS PHẠM VĂN TẤT PGS.TS TRẦN DƯƠNG HUẾ, NĂM 2020 LỜI CAM ĐOAN Tôi xin cam đoan nghiên cứu riêng Tôi hướng dẫn PGS.TS Phạm Văn Tất PGS.TS Trần Dương Các số liệu kết nêu luận án trung thực chưa công bố luận văn, luận án khoa học khác Luận án thực Trường Đại học Khoa Học – Đại học Huế Chưa có kết nghiên cứu tương tự công bố hình thức trước thực luận án Một phần kết cơng trình cơng bố trên: Tạp Chí Khoa học Cơng nghệ - Viện Hàn Lâm Khoa học Công nghệ, Tạp chí Hóa học - Viện Hàn Lâm Khoa học Cơng nghệ, Tạp chí SmartScience, Tạp chí Đại học Huế, Tạp chí Khoa học Cơng nghệ - Trường Đại học Khoa học – Đại học Huế Nghiên cứu sinh Nguyễn Thành Được LỜI CẢM ƠN Để hoàn thành luận án trước hết Tơi xin bày tỏ lòng biết ơn sâu sắc đến PGS.TS Phạm Văn Tất- Trường Đại học Hoa Sen; PGS TS Trần DươngTrường Đại học Sư Phạm Huế giao đề tài, hướng dẫn trực tiếp truyền đạt kinh nghiệm kiến thức quý báu, tận tình dẫn, động viên tạo điều kiện thuận lợi giúp đỡ Tơi hồn thành luận án Tôi xin gửi lời cảm ơn chân thành đến GS.TS Trần Thái Hòa, PGS.TS Hồng Thái Long, PGS.TS Đinh Quang Khiếu, PGS.TS Nguyễn Thị Ái Nhung, TS Lê Thị Hòa - Khoa Hóa, Trường Đại học Khoa Học Huế Các Thầy Cô giúp đỡ, động viên suốt q trình Tơi học tập trường Tơi xin gửi lời cảm ơn Thầy Cơ Khoa Hóa, Thầy Cơ Phòng Sau đại học Ban Giám hiệu Trường Đại học Khoa Học Huế cho phép tạo thuận lợi cho tơi hồn thành luận án Tôi xin gửi lời cảm ơn chân thành đến gia đình tất bạn bè đồng nghiệp đặc biệt anh TS Nguyễn Tiến Dũng, chị Trần Thanh Nhung em Lê Văn Phi Long động viên, giúp đỡ cho tơi hồn thành luận án Nghiên cứu sinh Nguyễn Thành Được MỤC LỤC DANH MỤC CÁC CHỮ VIẾT TẮT GEMC COSMO COSMOSAC MP CCSD EOS PR-EOS Gibbs Ensemble Monte Carlo Conductor-like Screening Model Conductor-like Screening Tồn cục Gibbs Monte Carlo Mơ hình sàng lọc vật dẫn Model Segment Activity Hệ số hoạt độ phân đoạn Coefficient Möller-Plesset Coupled Cluster Single Double Equation of State Peng Robinson - Equation of Mưller-Plesset Nhóm cặp đơn đơi Phương trình trạng thái Phương trình trạng thái State D-EOS Deiters - Equation of State ANN Artificial Neural Network RMSE Root Mean Square Error MARE Mean Average Relative Error ARE Average Relative Error Cal Exp Reco Calculation Experiment Recover Peng – Robinson Phương trình trạng thái Deiters Mạng thần kinh nhân tạo Căn bậc hai sai số bình phương trung bình Giá trị trung bình ARE,% Giá trị tuyệt đối sai số tương đối Tính toán Thực nghiệm Phục hồi DANH MỤC CÁC BẢNG Bảng 3.1 Năng lượng tương tác ab initio Cl2-Cl2 cho L H tính hàm sở CCSD(T)/aug-cc-pVmZ 49 Bảng 3.2 Năng lượng tương tác ab initio Cl2-Cl2 cho T X tính hàm CCSD(T)/aug-cc-pVmZ 51 Bảng 3.3 Tối ưu hóa độ dài liên kết nitơ hàm khác 55 Bảng 3.4 Tối ưu hóa độ dài liên kết CO hàm khác 55 Bảng 3.5 Tối ưu hóa tham số hàm (2.3) cho đime Ar-Ar .56 Bảng 3.6 Tối ưu hóa tham số hàm (2.3) cho tương tác N2-N2 .58 Bảng 3.7 Tối ưu hóa tham số hàm (2.4) cho tương tác N2-N2 .59 Bảng 3.8 Các tham số tối ưu hàm (2.4) cho tương tác Cl2-Cl2 59 Bảng 3.9 Các tham số tối ưu hàm (2.5) cho tương tác Cl2-Cl2 60 Bảng 3.10 Tối ưu hóa tham số hàm (2.3) cho CO-CO 62 Bảng 3.11 Các tham số tối ưu hàm (2.4) cho tương tác CO-CO 64 Bảng 3.12 Hệ số virial bậc hai B2(T) Cl2 thu từ hàm (2.4), (2.5), D-EOS (1.34) liệu thực nghiệm 68 Bảng 3.13 Dữ liệu ban đầu tính chất tới hạn Pc, Vc, Tc, TL TU, hệ số a, b c .73 Bảng 3.14 Ma trận tương quan tham số tới hạn Pc, Vc, Tc, TL, TU hợp chất 74 Bảng 3.15 Dữ liệu thành phần tương ứng với tính chất tới hạn PC, VC, TC, TL TU, hệ số a, b c chuyển đổi thành logarit 77 Bảng 3.16 Kết hệ số tính tốn (Cal.) so với giá trị gốc từ thực nghiệm (Exp.) Các hệ số tính tốn từ mạng I(5)-HL(6)-O(3) phục hồi trở lại (Reco.) từ giá trị logarit 78 Bảng 3.17 Các giá trị nhiệt động nitơ GEMC-NVT sử dụng hàm (2.4) liệu thực nghiệm (TN) 83 Bảng 3.18 Sai số MARE,% giá trị nhiệt động nitơ sử dụng hàm (2.4) cho mô GEMC-NVT 84 Bảng 3.19 Các giá trị nhiệt động nitơ sử dụng phương trình (1.32) liệu thực nghiệm (TN) 84 Bảng 3.20 Sai số MARE,% giá trị nhiệt động nitơ sử dụng phương trình (1.32 85 Bảng 3.21 Các giá trị nhiệt động cacbon monoxit GEMC-NVT sử dụng hàm (2.4) liệu TN 86 Bảng 3.22 Sai số MARE,% giá trị nhiệt động cacbon monoxit sử dụng hàm (2.4) cho mô GEMC-NVT 87 Bảng 3.23 Các giá trị nhiệt động cacbon monoxit sử dụng phương trình (1.32) liệu TN 87 Bảng 3.24 Sai số MARE,% giá trị nhiệt động cacbon monoxit sử dụng phương trình (1.32) .88 Bảng 3.25 Các tính chất tới hạn nitơ cacbon monoxit từ kết mô GEMC-NVT; hàm (2.3) (2.4); phương trình (1.32) Thực nghiệm (TN) 88 Bảng 3.26 Sai số ARE,% giá trị nhiệt động N2 CO 89 Bảng 3.27 Hệ số virial bậc hai B2(T) CO-CO thu từ hàm (2.4) tính tốn mức lý thuyết CCSD(T)/aug-cc-pVQZ .91 Bảng 3.28 So sánh giá trị mơ hình RMS, MRDp MDy 94 DANH MỤC CÁC HÌNH Hình 1.1 a) Mơ tả q trình solvat hóa lý tưởng; b) So sánh q trình solvat hóa lý tưởng Klamt bình thường Hình 1.2 Sơ đồ hệ mô tả mô GEMC 22 Hình 1.3 Các tiến trình diễn q trình mơ GEMC 23 Hình 1.4 Bộ não người cấu trúc nơron sinh học .27 Hình 1.5 Mạng thần kinh nhân tạo 28 Hình 2.1 Sơ đồ tiến trình cơng việc thực .33 Hình 2.2 Cấu hình đime bốn kiểu tương tác đặc biệt .36 Hình 2.3 Sơ đồ xây dựng hàm tương tác phân tử 38 Hình 2.4 Sơ đồ tính hệ số virial bậc hai B2(T) .40 Hình 2.5 Cấu trúc mạng thần kinh ANN- PCA I(5)-HL(6)-O(3) 42 Hình 2.6 Sơ đồ mô GEMC-NVT 43 Hình 2.7 Sơ đồ mơ hình COSMO 45 Hình 3.1 Bề mặt ab initio đime Ar-Ar .47 Hình 3.2 Bề mặt ab initio cấu hình đặc biệt đime N2-N2 48 Hình 3.3 Bề mặt ab initio cấu hình L H đime Cl2-Cl2 50 Hình 3.4 Bề mặt ab initio cấu hình T X đime Cl2-Cl2 52 Hình 3.5 Bề mặt 15 cấu hình khác Cl2-Cl2 .53 Hình 3.6 Bề mặt ab initio cấu hình đặc biệt đime CO-CO .54 Hình 3.7 Mối quan hệ tương quan lượng ab initio với lượng khớp từ hàm (2.3) Ar-Ar 57 Hình 3.8 Tương quan lượng ab initio lượng khớp clo theo a) hàm (2.4); b) hàm (2.5) 61 Hình 3.9 Mối quan hệ tương quan lượng ab initio với lượng khớp từ hàm (2.3) CO-CO 63 Hình 3.10 Tương quan lượng ab initio khớp ab initio theo hàm (2.4) CO-CO 65 Hình 3.11 Hệ số virial bậc hai B2(T) argon 66 Hình 3.12 Hệ số virial bậc hai đime N2-N2 67 Hình 3.13 Hệ số virial bậc hai hệ Cl2 tính phương trình (1.34) so sánh với hàm (2.4)-(a), (2.5)-(b) liệu thực nghiệm .71 Hình 3.14 Hệ số virial bậc hai đime CO-CO 72 Hình 3.15 Đồ thị sườn dốc lựa chọn (a) đồ thị kép (b) có lợi để xác định số lượng thành phần thích hợp khơng gian 75 Hình 3.16 Hệ số virial bậc hai cho chất khí a) argon; b) khí N 2; c) khí CO; d) khí Cl2 79 Hình 3.17 Sự phụ thuộc hàm phân bố g(rN-N) g(rM-M) vào nhiệt độ q trình mơ GEMC-NVT N2 .81 Hình 3.18 Sự phụ thuộc hàm phân bố g(rC-C), g(rO-O), g(rC-O) g(rM) vào nhiệt độ q trình mơ GEMC-NVT CO 82 M Hình 3.19 Giản đồ cân lỏng – nitơ 83 Hình 3.20 Giản đồ cân lỏng-hơi cacbon monoxit 86 Hình 3.21 Tương quan bề mặt điện tích σm với lượng ab initio 90 Hình 3.22 a) Bề mặt điện tích b) Áp suất hệ CO-CO xác định từ kết tính COSMO .92 Hình 3.23 Các giá trị sigma phân tử đơn CO Cl2 xác định từ tính tốn COSMO 93 Hình 3.24 Giản đồ cân lỏng-hơi P-x-y hỗn hợp CO (1) Cl2 (2) .94 Hình 3.25 Giản đồ x-y xây dựng từ hệ số hoạt độ điều kiện đẳng nhiệt 95 Hình 3.26 Giản đồ cân lỏng T-x-y hệ bậc hai CO (1)-Cl2 (2) 96 Hình 3.27 Giản đồ x-y xây dựng từ hệ số hoạt độ điều kiện đẳng áp .97 ĐẶT VẤN ĐỀ Ngày nay, với phát triển mạnh mẽ cơng nghệ chế tạo máy tính với đời hệ máy tính có tốc độ tính tốn cực nhanh hỗ trợ cho hướng nghiên cứu khoa học tính tốn, có ngành mơ máy tính Từ mơ hình hệ thống thực chủ yếu dạng phương trình, hệ phương trình tốn học, máy tính với khả tính tốn nhanh độ xác cao hỗ trợ cho việc giải mơ hình Các kết tính tốn so sánh với kết thực nghiệm để rút kết luận tính đắn mơ hình [2] Phương pháp mơ tồn cục Gibbs Ensemble Monte Carlo (GEMC) khơng giống phương pháp mơ xác khác chỗ phương pháp mơ ngẫu nhiên, phương pháp không giải phương trình tường minh mà nhận kết mô hạt riêng lẻ ghi lại tham số trạng thái trung bình chúng Trạng thái trung bình hạt vật lý rút từ trạng thái trung bình hạt mơ [59] Trên giới, kỹ thuật mô GEMC phát triển đưa vào nghiên cứu từ sớm, tác giả nhóm tác giả đạt thành tựu tiêu biểu lĩnh vực Michels nghiên cứu hệ số virial bậc hai hiđro deuterium năm 1960 [12], Dymond Smith nghiên cứu hệ số virial khí tinh khiết hỗn hợp năm 1980 [47], Panagiotopoulos nghiên cứu tính chất cân pha chất lỏng năm 1987 [13], Spyriouni nghiên cứu cân pha lỏng n – hexadecan tinh khiết nhiệt độ cao năm 1998 [43], Martin Siepmann nghiên cứu chuyển tiếp cân pha họ ankan tính tốn lượng tự Gibbs năm 1998 1999 [59], Leonhard Deiters nghiên cứu việc dùng ab initio 10 r/Å 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.6 3.8 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 α 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Tọa độ góc β φ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Năng lượng ab initio E-pVDZ E-pVTZ 11.201 9.349 10.186 8.501 9.272 7.738 8.448 7.050 7.705 6.430 7.035 5.871 6.430 5.366 5.883 4.909 5.388 4.496 4.939 4.121 4.533 3.782 4.164 3.474 3.829 3.195 3.524 2.941 3.247 2.709 2.995 2.498 2.764 2.306 10214.445 8962.072 5265.437 4559.201 2551.138 2175.327 1112.371 927.123 383.845 303.167 40.169 13.405 -101.990 -103.467 -143.640 -135.250 -138.976 -128.501 -116.076 -107.041 -88.748 -82.418 -63.210 -59.728 -41.791 -40.803 -24.929 -25.915 -12.233 -14.668 -3.018 -6.444 3.435 -0.609 7.776 3.400 10.541 6.045 12.156 7.690 136 E-pV23Z 8.585 7.806 7.106 6.474 5.905 5.391 4.927 4.508 4.128 3.784 3.473 3.190 2.933 2.700 2.488 2.294 2.117 8434.642 4261.699 2017.030 849.118 269.209 2.143 -104.093 -131.732 -124.109 -103.254 -79.764 -58.266 -40.382 -26.318 -15.675 -7.863 -2.283 1.589 4.183 5.841 r/Å 7.6 7.8 8.2 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.4 3.6 3.8 α 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Tọa độ góc β φ 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 90 90 90 90 Năng lượng ab initio E-pVDZ E-pVTZ 12.951 8.616 13.172 9.034 13.003 9.100 12.576 8.929 11.989 8.604 11.309 8.183 10.586 7.710 9.851 7.213 9.127 6.712 8.429 6.222 7.766 5.750 7.142 5.302 6.560 4.881 6.020 4.489 5.522 4.125 5.063 3.788 4.643 3.479 4.257 3.194 3.905 2.933 3.583 2.694 3.289 2.476 3.021 2.276 2.776 2.093 2.553 1.926 2.349 1.774 2.163 1.634 1.993 1.507 1.838 1.390 1.696 1.284 1.567 1.186 1.448 1.097 1.339 1.015 1.240 0.940 18977.529 16872.085 10214.445 8962.072 5265.437 4559.201 2551.138 2175.327 137 E-pV23Z 6.822 7.322 7.487 7.422 7.205 6.892 6.522 6.123 5.715 5.310 4.918 4.543 4.189 3.857 3.549 3.263 2.999 2.756 2.532 2.328 2.140 1.969 1.811 1.668 1.536 1.416 1.306 1.206 1.114 1.029 0.952 0.881 0.816 15986.030 8434.642 4261.699 2017.030 r/Å 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 α 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Tọa độ góc β φ 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Năng lượng ab initio E-pVDZ E-pVTZ 1112.371 927.123 383.845 303.167 40.169 13.405 -101.990 -103.467 -143.640 -135.250 -138.976 -128.501 -116.076 -107.041 -88.748 -82.418 -63.210 -59.728 -41.791 -40.803 -24.929 -25.915 -12.233 -14.668 -3.018 -6.444 3.435 -0.609 7.776 3.400 10.541 6.045 12.156 7.690 12.951 8.616 13.172 9.034 13.003 9.100 12.576 8.929 11.989 8.604 11.309 8.183 10.586 7.710 9.851 7.213 9.127 6.712 8.429 6.222 7.766 5.750 7.142 5.302 6.560 4.881 6.020 4.489 5.522 4.125 5.063 3.788 4.643 3.479 4.257 3.194 3.905 2.933 3.583 2.694 138 E-pV23Z 849.118 269.209 2.143 -104.093 -131.732 -124.109 -103.254 -79.764 -58.266 -40.382 -26.318 -15.675 -7.863 -2.283 1.589 4.183 5.841 6.822 7.322 7.487 7.422 7.205 6.892 6.522 6.123 5.715 5.310 4.918 4.543 4.189 3.857 3.549 3.263 2.999 2.756 2.532 2.328 r/Å 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.6 3.8 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 α 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Tọa độ góc β φ 90 90 90 90 90 90 90 90 90 90 90 90 90 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 Năng lượng ab initio E-pVDZ E-pVTZ 3.289 2.476 3.021 2.276 2.776 2.093 2.553 1.926 2.349 1.774 2.163 1.634 1.993 1.507 1.838 1.390 1.696 1.284 1.567 1.186 1.448 1.097 1.339 1.015 1.240 0.940 12358.691 12913.491 5755.764 6074.569 2538.961 2724.933 1020.173 1124.466 336.920 388.282 54.110 70.619 -44.240 -50.370 -63.026 -83.281 -51.582 -80.019 -31.964 -64.495 -13.176 -47.062 1.780 -31.694 12.485 -19.495 19.494 -10.385 23.619 -3.870 25.635 0.619 26.183 3.592 25.754 5.462 24.708 6.549 23.300 7.090 21.706 7.260 20.046 7.181 18.397 6.944 16.806 6.609 139 E-pV23Z 2.140 1.969 1.811 1.668 1.536 1.416 1.306 1.206 1.114 1.029 0.952 0.881 0.816 14231.093 6829.270 3166.099 1372.922 511.236 110.456 -64.756 -131.580 -147.967 -142.250 -128.029 -111.623 -95.796 -81.611 -69.339 -58.912 -50.129 -42.756 -36.572 -31.379 -27.008 -23.320 -20.198 -17.548 r/Å 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.4 3.6 3.8 4.2 4.4 4.6 4.8 α 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 45 45 45 45 45 45 45 Tọa độ góc β φ 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 45 45 45 45 45 45 45 45 Năng lượng ab initio E-pVDZ E-pVTZ 15.303 6.220 13.903 5.807 12.612 5.388 11.429 4.978 10.352 4.584 9.374 4.211 8.490 3.862 7.691 3.537 6.971 3.237 6.322 2.962 5.737 2.709 5.211 2.478 4.737 2.267 4.310 2.075 3.925 1.899 3.578 1.740 3.264 1.595 2.982 1.463 2.726 1.342 2.495 1.233 2.286 1.133 2.096 1.042 1.924 0.960 1.768 0.884 1.626 0.815 1.497 0.752 1.379 0.695 1.272 0.642 1.175 0.594 13002.559 11539.581 6940.445 6063.563 3584.021 3074.832 1775.822 1487.014 831.666 670.249 358.173 268.022 134.228 82.659 38.162 6.763 140 E-pV23Z -15.290 -13.362 -11.709 -10.288 -9.062 -8.003 -7.085 -6.286 -5.590 -4.982 -4.449 -3.981 -3.570 -3.207 -2.886 -2.602 -2.349 -2.125 -1.926 -1.747 -1.588 -1.445 -1.317 -1.202 -1.099 -1.006 -0.921 -0.845 -0.776 10924.178 5694.372 2860.373 1365.373 602.273 230.064 60.944 -6.465 r/Å 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 α 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Tọa độ góc β φ 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Năng lượng ab initio E-pVDZ E-pVTZ 4.554 -16.691 -0.808 -17.123 4.848 -9.116 13.324 0.537 21.080 9.010 26.889 15.407 30.587 19.703 32.450 22.212 32.883 23.337 32.284 23.455 30.990 22.881 29.264 21.860 27.304 20.573 25.248 19.152 23.191 17.685 21.198 16.234 19.307 14.837 17.538 13.518 15.903 12.287 14.403 11.152 13.035 10.111 11.793 9.162 10.670 8.301 9.655 7.521 8.740 6.816 7.916 6.179 7.175 5.606 6.508 5.089 5.908 4.623 5.368 4.203 4.882 3.825 4.445 3.484 4.050 3.177 3.694 2.899 3.373 2.648 3.083 2.421 2.821 2.216 141 E-pV23Z -25.647 -24.005 -15.007 -4.855 3.924 10.574 15.126 17.911 19.330 19.752 19.482 18.758 17.756 16.601 15.383 14.159 12.970 11.838 10.778 9.795 8.891 8.065 7.313 6.630 6.013 5.455 4.951 4.497 4.087 3.718 3.384 3.084 2.813 2.567 2.346 2.146 1.965 r/Å 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.4 3.6 3.8 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 8.8 α 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Tọa độ góc β φ 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Năng lượng ab initio E-pVDZ E-pVTZ 2.584 2.031 2.369 1.862 2.174 1.710 1.997 1.571 1.837 1.445 1.691 1.331 1.558 1.226 1.437 1.131 1.326 1.044 10679.643 9331.602 5619.936 4812.331 2810.456 2342.675 1304.495 1040.515 530.706 384.423 156.034 75.403 -8.382 -53.680 -66.913 -93.942 -75.645 -93.661 -63.693 -77.497 -45.450 -57.381 -27.492 -38.585 -12.393 -23.027 -0.785 -11.036 7.560 -2.262 13.187 3.868 16.692 7.943 18.616 10.478 19.405 11.896 19.406 12.524 18.884 12.609 18.030 12.331 16.982 11.823 15.839 11.179 14.665 10.463 13.506 9.720 12.389 8.980 11.332 8.263 142 E-pV23Z 1.800 1.651 1.516 1.394 1.282 1.181 1.088 1.004 0.927 8764.113 4472.292 2145.717 929.377 322.843 41.460 -72.754 -105.329 -101.255 -83.315 -62.406 -43.252 -27.496 -15.338 -6.380 -0.036 4.280 7.075 8.758 9.650 9.989 9.953 9.672 9.236 8.711 8.142 7.560 6.985 r/Å 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.2 3.4 3.6 3.8 4.2 4.4 4.6 4.8 5.2 α 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Tọa độ góc β φ 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 Năng lượng ab initio E-pVDZ E-pVTZ 10.343 7.580 9.428 6.939 8.585 6.343 7.813 5.791 7.109 5.285 6.468 4.821 5.887 4.397 5.359 4.012 4.882 3.661 4.449 3.342 4.058 3.053 3.703 2.790 3.383 2.552 3.092 2.335 2.829 2.139 2.591 1.961 2.375 1.799 2.178 1.652 2.000 1.518 1.838 1.396 1.691 1.285 1.557 1.184 1.435 1.092 1.323 1.008 1.222 0.931 1.129 0.860 12579.669 10766.114 6765.434 5615.482 3347.091 2658.095 1449.436 1053.707 457.997 238.492 -18.974 -137.311 -217.313 -280.040 -273.085 -306.678 -261.822 -281.048 -223.709 -236.428 -178.958 -189.109 143 E-pV23Z 6.430 5.903 5.410 4.950 4.526 4.136 3.778 3.451 3.153 2.882 2.635 2.410 2.206 2.021 1.852 1.699 1.559 1.433 1.317 1.212 1.116 1.029 0.949 0.876 0.810 0.749 10002.865 5131.307 2367.974 887.100 146.105 -187.099 -306.424 -320.812 -289.146 -241.797 -193.405 r/Å 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 α 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 Tọa độ góc β φ 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 45 90 Năng lượng ab initio E-pVDZ E-pVTZ -136.872 -146.256 -101.036 -110.340 -72.201 -81.574 -49.840 -59.203 -32.958 -42.160 -20.481 -29.373 -11.433 -19.895 -4.993 -12.945 -0.504 -7.899 2.545 -4.276 4.544 -1.707 5.787 0.086 6.489 1.310 6.813 2.121 6.873 2.631 6.757 2.926 6.524 3.068 6.219 3.103 5.872 3.063 5.507 2.973 5.137 2.851 4.773 2.710 4.422 2.558 4.087 2.403 3.771 2.248 3.476 2.096 3.201 1.950 2.946 1.812 2.711 1.681 2.495 1.557 2.295 1.442 2.113 1.335 1.945 1.236 1.791 1.143 1.650 1.058 1.521 0.979 1.403 0.907 144 E-pV23Z -150.231 -114.281 -85.543 -63.166 -46.052 -33.130 -23.469 -16.301 -11.019 -7.152 -4.341 -2.316 -0.871 0.145 0.845 1.314 1.614 1.792 1.881 1.908 1.890 1.843 1.775 1.695 1.607 1.516 1.425 1.335 1.247 1.164 1.084 1.008 0.938 0.871 0.809 0.752 0.698 r/Å 12.8 13 13.2 13.4 13.6 13.8 14 2.6 2.8 3.2 3.4 3.6 3.8 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 α 45 45 45 45 45 45 45 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Tọa độ góc β φ 45 90 45 90 45 90 45 90 45 90 45 90 45 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Năng lượng ab initio E-pVDZ E-pVTZ 1.295 0.840 1.196 0.778 1.105 0.721 1.022 0.668 0.946 0.620 0.876 0.576 0.811 0.534 13384.562 11682.851 6439.967 5379.846 2915.200 2258.083 1195.869 787.470 395.809 139.590 48.632 -114.612 -83.692 -189.863 -119.303 -190.074 -115.044 -163.511 -97.160 -131.287 -77.017 -101.703 -58.983 -77.291 -44.295 -58.174 -32.895 -43.621 -24.283 -32.707 -17.875 -24.581 -13.146 -18.547 -9.670 -14.063 -7.117 -10.720 -5.241 -8.219 -3.861 -6.337 -2.842 -4.914 -2.088 -3.832 -1.529 -3.003 -1.113 -2.366 -0.804 -1.872 -0.573 -1.488 -0.400 -1.188 -0.272 -0.952 -0.176 -0.765 145 E-pV23Z 0.648 0.602 0.559 0.520 0.483 0.450 0.418 10966.970 4933.586 1981.412 615.520 31.720 -183.338 -234.567 -219.879 -183.929 -145.668 -112.107 -85.006 -64.022 -48.138 -36.252 -27.402 -20.817 -15.908 -12.233 -9.467 -7.375 -5.782 -4.561 -3.620 -2.889 -2.318 -1.870 -1.516 -1.235 -1.011 r/Å 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.2 3.4 3.6 3.8 4.2 4.4 4.6 α 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Tọa độ góc β φ 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Năng lượng ab initio E-pVDZ E-pVTZ -0.105 -0.617 -0.053 -0.499 -0.015 -0.405 0.013 -0.329 0.033 -0.268 0.047 -0.219 0.056 -0.179 0.062 -0.146 0.065 -0.120 0.067 -0.098 0.067 -0.080 0.066 -0.066 0.064 -0.054 0.062 -0.044 0.059 -0.036 0.057 -0.029 0.054 -0.024 0.051 -0.019 0.048 -0.015 0.046 -0.012 0.043 -0.009 0.040 -0.007 0.038 -0.006 0.036 -0.004 0.033 -0.003 0.031 -0.002 0.029 -0.001 0.027 0.000 2078.882 1537.063 680.598 340.229 54.323 -159.616 -193.471 -328.648 -264.628 -350.879 -259.523 -315.334 -226.238 -262.998 -186.585 -211.301 -149.565 -166.563 146 E-pV23Z -0.830 -0.685 -0.568 -0.472 -0.394 -0.329 -0.277 -0.233 -0.197 -0.167 -0.142 -0.121 -0.103 -0.088 -0.075 -0.065 -0.056 -0.048 -0.042 -0.036 -0.031 -0.027 -0.024 -0.020 -0.018 -0.016 -0.014 -0.012 1308.889 196.914 -249.680 -385.552 -387.192 -338.840 -278.488 -221.722 -173.731 r/Å 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 α 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Tọa độ góc β φ 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 Năng lượng ab initio E-pVDZ E-pVTZ -118.187 -130.151 -92.825 -101.442 -72.831 -79.174 -57.269 -62.032 -45.222 -48.863 -35.904 -38.731 -28.681 -30.907 -23.061 -24.834 -18.665 -20.092 -15.205 -16.365 -12.466 -13.416 -10.284 -11.067 -8.533 -9.184 -7.120 -7.665 -5.973 -6.431 -5.036 -5.423 -4.267 -4.595 -3.632 -3.911 -3.104 -3.344 -2.664 -2.870 -2.296 -2.473 -1.985 -2.139 -1.723 -1.857 -1.501 -1.617 -1.311 -1.413 -1.149 -1.238 -1.010 -1.089 -0.890 -0.960 -0.786 -0.848 -0.697 -0.751 -0.619 -0.667 -0.551 -0.594 -0.492 -0.530 -0.440 -0.474 -0.394 -0.425 -0.354 -0.382 -0.318 -0.343 147 E-pV23Z -135.198 -105.076 -81.847 -64.037 -50.393 -39.917 -31.839 -25.575 -20.687 -16.847 -13.810 -11.392 -9.453 -7.889 -6.619 -5.581 -4.729 -4.026 -3.441 -2.954 -2.546 -2.202 -1.911 -1.665 -1.454 -1.275 -1.120 -0.988 -0.873 -0.773 -0.687 -0.612 -0.546 -0.488 -0.438 -0.393 -0.354 r/Å 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 3.4 3.6 3.8 4.2 4.4 4.6 4.8 5.2 5.4 5.6 5.8 6.2 6.4 6.6 6.8 7.2 7.4 7.6 7.8 8.2 8.4 8.6 α 90 90 90 90 90 90 90 90 90 90 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 Tọa độ góc β φ 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 Năng lượng ab initio E-pVDZ E-pVTZ -0.287 -0.310 -0.259 -0.280 -0.234 -0.253 -0.212 -0.229 -0.193 -0.208 -0.175 -0.189 -0.160 -0.172 -0.145 -0.157 -0.133 -0.143 -0.121 -0.131 13364.493 11950.331 5552.373 4846.914 2149.141 1770.293 717.043 496.746 149.637 11.802 -49.568 -141.015 -99.687 -162.982 -94.973 -140.069 -74.292 -107.062 -52.484 -76.634 -34.330 -52.321 -20.657 -34.184 -10.927 -21.181 -4.272 -12.105 0.123 -5.905 2.915 -1.758 4.595 0.948 5.520 2.655 5.940 3.674 6.030 4.227 5.909 4.466 5.660 4.498 5.337 4.397 4.979 4.214 4.608 3.983 4.241 3.728 3.887 3.464 148 E-pV23Z -0.319 -0.288 -0.260 -0.236 -0.214 -0.195 -0.177 -0.162 -0.148 -0.135 11355.262 4549.788 1610.717 403.982 -46.225 -179.516 -189.637 -159.065 -120.867 -86.806 -59.896 -39.874 -25.489 -15.392 -8.431 -3.712 -0.574 1.461 2.733 3.480 3.869 4.019 4.011 3.901 3.728 3.519 3.292 r/Å 8.8 9.2 9.4 9.6 9.8 10 10.2 10.4 10.6 10.8 11 11.2 11.4 11.6 11.8 12 12.2 12.4 12.6 12.8 13 13.2 13.4 13.6 13.8 14 α 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 Tọa độ góc β φ 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 Năng lượng ab initio E-pVDZ E-pVTZ 3.553 3.202 3.240 2.949 2.952 2.708 2.686 2.482 2.443 2.272 2.222 2.077 2.020 1.898 1.838 1.734 1.673 1.584 1.523 1.448 1.388 1.323 1.265 1.210 1.154 1.107 1.054 1.013 0.963 0.928 0.881 0.850 0.807 0.780 0.739 0.716 0.678 0.658 0.623 0.605 0.572 0.557 0.526 0.513 0.484 0.472 0.446 0.436 0.412 0.402 0.380 0.372 0.351 0.344 149 E-pV23Z 3.061 2.832 2.611 2.401 2.203 2.020 1.850 1.693 1.550 1.418 1.298 1.188 1.088 0.997 0.914 0.839 0.770 0.707 0.650 0.598 0.551 0.508 0.468 0.432 0.399 0.369 0.341 Bảng Input sử dụng cho phép tính tốn lượng tử #T CCSD(T)/aug-cc-pVmZ MaxDisk=8192MW SCF=Tight Test A–A A–A #T CCSD(T)/aug-cc-pVmZ MaxDisk=8192MW nosymm SCF=Tight Test A–A A – ABq #T CCSD(T)/aug-cc-pVmZ MaxDisk=8192MW nosymm SCF=Tight Test A – ABq A–A 150 ... Tính tốn cân lỏng- hơi Ar, N2, Cl2, CO phương pháp hóa lượng tử mơ tồn cục Monte Carlo 11 Mục tiêu luận án Tính tốn hệ số virial bậc hai xác định giá trị nhiệt động cân lỏng- hơi cho chất Ar,. .. lỏng- hơi cho chất Ar, N2, Cl2, CO phương pháp hóa lượng tử kỹ thuật mơ tồn cục Monte Carlo Ý nghĩa khoa học luận án Luận án đưa hướng nghiên cứu – tính tốn giá trị nhiệt động cân lỏng- hơi áp suất tới... mơ hình COSMO COSMO phương pháp tính tốn để xác định tương tác tĩnh điện phân tử có dung môi Phương pháp thường sử dụng hóa học tính tốn để mơ hình hóa hiệu ứng solvat hóa [29] Mơ hình COSMO Klamt