Trong công trình này chúng tôi đưa ra các kết quả mô phỏng cân bằng lỏng hơi của flo lỏng bằng kỹ thuật mô phỏng Monte Carlo toàn cục sử dụng thế cặp tương tác phân tử 5 vị trí mới ab initio đã được chúng tôi đưa ra từ các tính toán cơ học lượng tử đối với dime F2-F2. Các kết quả mô phỏng bao gồm tỷ trọng, áp suất hơi và entanpi hóa hơi đã được so sánh với các số liệu thực nghiệm và với những số liệu lấy từ các tài liệu.
Tạp chí KHOA HỌC ĐHSP TP HCM Phạm Văn Tất MÔ PHỎNG MONTE CARLO CÂN BẰNG LỎNG HƠI CỦA FLO SỬ DỤNG CÁC THẾ TƯƠNG TÁC PHÂN TỬ AB INITIO MỚI Phạm Văn Tất* Giới thiệu Mô máy tính trở thành cơng cụ bắt buộc cần thiết cho nghiên cứu chất lỏng hỗn hợp chất lỏng [1] Một cố gắng Nasrabad Deiters dự đoán cân pha lỏng áp suất cao hỗn hợp khí [2,3] từ mơ toàn cục sử dụng tương tác phân tử Leonhard Deiters sử dụng Morse vị trí để mơ tả cặp nitơ [4] dự đoán áp suất tỷ trọng cân pha mơ Monte Carlo tồn cục[5] Trong cơng trình chúng tơi đưa kết mô cân lỏng flo lỏng kỹ thuật mơ Monte Carlo tồn cục sử dụng cặp tương tác phân tử vị trí ab initio đưa từ tính tốn học lượng tử dime F2-F2 Các kết mô bao gồm tỷ trọng, áp suất entanpi hóa so sánh với số liệu thực nghiệm với số liệu lấy từ tài liệu Phương pháp tính tốn 2.1 Phương pháp mơ Mơ GEMC-NPT sử dụng để kiểm tra xác cặp Mô khảo sát điều kiện đẳng áp 1,0 MPa 10,0 MPa khoảng nhiệt độ từ 90,0 K đến 270,0 K Mô GEMC-NVT tiến hành để nhận tỷ trọng lỏng, áp suất khoảng 60,0 K đến 140,0 K với thay đổi 10,0 K Các cặp vị trí flo dime dùng cho hai trường hợp mô phỏng: 5 u1 Deij ((1 e i 1 i 1 * ij ( rij qi q j C ij ) 1) f1a (rij ) nn f ( rij ) n6,8,10 r 4 rij ij ij ) (1) TS – Trường ĐH Đà Lạt 63 Tạp chí KHOA HỌC ĐHSP TP HCM 5 u Deij ((1 e i 1 i 1 ij Số 16 năm 2009 qi q j Cnij ) 1) f1b (rij ) f (rij ) n n6,8,10,12 r 4 rij ij ( rij ij ) f 1a (rij ) (1 e 2 ( ij rij ) 10 ) 15 f 1b ( rij ) e r ij ij k 0 (2) ( ij rij ) k k! Tổng số hạt N = 512 dùng hai trường hợp mô GEMC với điều kiện biên chuẩn Đối với việc chạy mô GEMC-NVT, cân hai pha thiết lập cần 10 đến 2.106 vòng lặp Tất chuyển động tiến hành cách ngẫu nhiên với xác suất xác định Số liệu mô đưa lần ứng với 1000 vòng lặp Mơ bắt đầu với tỷ trọng hai pha Hệ mô cân với khoảng 106 vòng lặp Bán kính giới hạn sử dụng 8,5 Å flo 2.2 Tính chất cân pha Các tính chất tới hạn: nhiệt độ Tc/K, tỷ trọng c / g.cm-3 thể tích Vc/ cm3 mol-1 flo lỏng tinh khiết rút từ phương pháp bình phương tối thiểu dựa vào tỷ trọng cân pha sử dụng phương trình quan hệ: 1 c A(Tc T ) 1 B (T Tc ) (3) l v tỷ trọng lỏng tỷ trọng hơi, hệ số mũ tới hạn ( = 0,325) A B số hiệu chỉnh Áp suất tới hạn Pc/ MPa tính phương trình Antoine Áp suất nhiệt hóa vapH tính phương trình Clausius-Clapeyron Kết thảo luận 3.1 Tính chất cấu trúc Hàm phân bố cặp vị trí tương tác rút từ hai kỹ thuật mô GEMCNVT –NPT flo lỏng tinh khiết mơ tả Hình Hình Các cặp ab initio Eq Eq flo sử dụng hai kỹ thuật mô tương ứng 64 Tạp chí KHOA HỌC ĐHSP TP HCM Phạm Văn Tất T=90K T=120K T=150K T=180K T=210K T=270K 4 2 1 T=60 K T=80 K T=90 K T=100 K T=110 K T=120 K T=130 K T=140 K g(F-F) g(F-F) 4.5 5.0 5.5 6.0 4.5 r/Å 5.0 a) 5.5 r/Å b) Hình 1: Sự phụ thuộc nhiệt độ g F-F P = 1,0 MPa; a) mô GEMC-NPT b) mô GEMC-NVT, hai trường hợp mô sử dụng Eq.1 Tất đỉnh Hình hàm tương quan vị trí nằm khoảng 3,98 Å 4,88 Å Các đỉnh thứ hai nằm 5,256 Å 5,717 Å 3.2 Tính chất cân pha Các kết mơ đưa Bảng Bảng Đường cong cân lỏng flo lỏng mô tả Hình Các số liệu thực nghiệm [8,9], giá trị tính tốn từ hàm trạng thái Deiters [7] tính tốn từ LennardJones [6] đưa đồ thị cân pha 150 T/K 140 130 120 110 100 90 80 70 60 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 -3 Hình Giản đồ cân pha; với: : số liệu thực nghiệm [8,9]; /g cmo: LennardJones; : phương trình trạng thái Deiters D1-EOS [7]; , *: cặp Eq Eq 65 6.0 Tạp chí KHOA HỌC ĐHSP TP HCM Số 16 năm 2009 Hvap/kJ mol p /MPa -1 0 60 70 80 90 100 110 120 130 140 150 60 70 80 90 100 110 120 130 T/K b) Hình a) Áp suất hơi, b) entanpi hóa hơi; xem giải thích Hình Các áp suất khác biệt so với thực nghiệm khoảng 1,4% 5,8% Những khác biệt nằm khoảng không chắn nguồn liệu thực nghiệm [8,9] Bảng Các tính chất tới hạn nhận từ mơ GEMC-NVT sử dụng cặp Eq Eq 2; D1-EOS: phương trình Deiters; LJ: LennardJones; Exp.: giá trị thực nghiệm Tc / K c/ g.cm-3 Pc /MPa Vc/ cm3mol-1 Tài liệu Eq 146,41 0,592 4,911 64,207 Cơng trình Eq 147,65 0,565 5,380 67,298 Cơng trình LJ 143,63 0,567 5,039 66,996 [6] D1-EOS 144,16 0,568 5,050 66,954 [7] Exp 144,30 0,574 5,215 66,200 [8] Exp 144,12 0,571 5,172 66,545 [9] Các tính chất tới hạn flo lỏng tinh khiết nhận từ tỷ trọng cân pha lỏng phương pháp bình phương tối thiểu dựa vào mối quan hệ (3), đưa Bảng Áp suất tới hạn flo phù hợp tốt với số liệu thực nghiệm 66 150 T/K a) Phương pháp 140 Tạp chí KHOA HỌC ĐHSP TP HCM Phạm Văn Tất Bảng Entanpi hóa vapH, entropi hóa vapS nhiệt độ sơi Tb P = 101,3 kPa tính tốn từ áp suất mô Phương pháp vapH/ kJ mol-1 vapS/ kJ.mol1 K-1 Tb/ K Tài liệu Eq 6,805 0,0791 85,978 Cơng trình Eq 7,160 0,0821 87,181 Cơng trình LJ 7,646 0,0871 87,768 [6] D1-EOS 7,131 0,0821 86,835 [7] Exp 6,921 0,0809 85,597 [8] Exp 6,941 0,0811 85,570 [9] Những khác biệt kết dự đốn số liệu thực nghiệm khơng có ý nghĩa Kết luận Cân pha lỏng tính chất nhiệt động học flo lỏng tính tốn thành cơng chương trình mơ GEMC-NVT GEMC-NPT phát triển sử dụng cặp tương tác phân tử ab initio Các kết mô cho thấy phù hợp tốt với số liệu thực nghiệm TÀI LIỆU THAM KHẢO [1] D R Lide, 2002, Handbook of Chemistry and Physics, CRC Press, 82nd Edition., Boca Raton A E Nasrabad and U K Deiters, J Chem Phys, 2003, 119, 947-952 [2] E Nasrabad, R Laghaei, and U K Deiters, J Chem Phys, 2004,121, 6423 [3] K Leonhard and U K Deiters, Mol Phys, 2002, 100, 2571-2585 [4] K M de Reuck, 1990 Fluorine international thermodynamic Tables of the Fluid State, vol-11 IUPAC Chemical Data series No 36, Oxford [5] M P Allen and D J Tildesley, 1991, Computer Simulation of Liquids., Clarendon Press, Oxford [6] U K Deiters homepage: http://thermoc.uni-koeln.de/index.html A Z Panagiotopoulos, Mol Phys, 1987, 61, 813-826 [7] Z Panagiotopoulos homepage: http://kea.princeton.edu/ppe/index.html 67 Tạp chí KHOA HỌC ĐHSP TP HCM Số 16 năm 2009 Tóm tắt Mơ Monte Carlo cân lỏng flo sử dụng tương tác phân tử ab initio Cân lỏng flo lỏng dự đốn kỹ thuật mơ Monte Carlo (GEMC) sử dụng hai cặp tương tác phân tử vị trí ab initio Các cặp ab initio xây dựng từ tính tốn sử dụng mức lý thuyết CCSD(T) với tập sở Dunning aug-cc-pVmZ (m = 2, 3) [8] Giản đồ pha, tính chất tới hạn, tính chất nhiệt động, áp suất tỷ trọng pha nhận phù hợp tốt với số liệu thực nghiệm Abstract Monte Carlo simulations of vapor – liquid equilibria of fluorine using new ab initio intermolecular interaction potentials The vapor-liquid equilibria of pure fluid fluorine were predicted by Gibbs ensemble Monte Carlo (GEMC) simulation techniques using our two new 5-site intermolecular pair potentials ab initio The ab initio pair potentials were established from coupled-cluster calculations, using the CCSD(T) level of theory and Dunning's correlation consistent basis sets aug-cc-pVmZ (m =2, 3) [8] The coexistence phase diagram, critical properties, thermodynamic properties, vapor pressures and orthobaric densities based on them are found to be in very good agreement with experimental data 68 ... tắt Mô Monte Carlo cân lỏng flo sử dụng tương tác phân tử ab initio Cân lỏng flo lỏng dự đốn kỹ thuật mơ Monte Carlo (GEMC) sử dụng hai cặp tương tác phân tử vị trí ab initio chúng tơi Các cặp ab. .. Tính chất cấu trúc Hàm phân bố cặp vị trí tương tác rút từ hai kỹ thuật mô GEMCNVT –NPT flo lỏng tinh khiết mô tả Hình Hình Các cặp ab initio Eq Eq flo sử dụng hai kỹ thuật mô tương ứng 64 Tạp chí... nghĩa Kết luận Cân pha lỏng tính chất nhiệt động học flo lỏng tính tốn thành cơng chương trình mơ GEMC-NVT GEMC-NPT phát triển sử dụng cặp tương tác phân tử ab initio Các kết mô cho thấy phù