So sánh các hàm tiềm năng đơn giản trong mô phỏng nước lỏng
Tóm tắt
Các mô phỏng Monte Carlo cổ điển đã được thực hiện cho nước lỏng trong tập hợp NPT ở nhiệt độ 25 °C và áp suất 1 atm, sử dụng sáu hàm tiềm năng giữa các phân tử đơn giản cho dimmer nước: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P và TIP4P. Các so sánh được thực hiện với dữ liệu nhiệt động lực học và cấu trúc thực nghiệm, bao gồm cả kết quả nhiễu xạ neutron gần đây của Thiessen và Narten. Các mật độ và năng lượng tiềm năng được tính toán tương đối phù hợp với thực nghiệm, ngoại trừ mô hình BF ban đầu, mô hình này đưa ra ước lượng quá cao 18% về mật độ và kết quả cấu trúc kém. Các tiềm năng TIPS2 và TIP4P cung cấp các hàm cấu trúc phân tử oxy-oxy có sự phù hợp tốt với các kết quả nhiễu xạ neutron. Tuy nhiên, sự phù hợp với các hàm cấu trúc phân tử OH và HH thực nghiệm kém hơn; mặc dù vậy, các kết quả tính toán cho các hàm này là tương tự đối với tất cả các hàm tiềm năng. Do đó, sự sai khác có thể do các hạng mục sửa chữa cần thiết trong quá trình xử lý dữ liệu neutron hoặc do một hiệu ứng bị bỏ qua đồng nhất trong các tính toán. Các so sánh cũng được thực hiện cho các hệ số khuếch tán tự phát thu được từ các mô phỏng động lực học phân tử. Tổng thể, các mô hình SPC, ST2, TIPS2 và TIP4P cung cấp các mô tả cấu trúc và nhiệt động lực học hợp lý cho nước lỏng và sẽ hữu ích trong các mô phỏng dung dịch nước. Sự đơn giản của các hàm SPC, TIPS2 và TIP4P cũng rất hấp dẫn từ quan điểm tính toán.