Luận án: Cân bằng lỏng-hơi Ar, N2, Cl2, CO bằng Hóa lượng tử và Monte Carlo

Các phương pháp ab initio khởi đầu từ các nguyên lý cơ bản của cơ học lượng tử. Chúng không sử dụng dữ liệu thực nghiệm nào ngoài các hằng số vật lý. Phương pháp này cho phép tính toán năng lượng tương tác chính xác cao. Ví dụ bao gồm phương pháp Hartree-Fock (HF), phương pháp nhiễu loạn Møller-Plesset (MP2), và lý thuyết nhóm cặp (Coupled Cluster). Lý thuyết hàm mật độ (DFT) cũng được áp dụng rộng rãi. DFT tập trung vào mật độ electron thay vì hàm sóng phức tạp. Phương pháp DFT cung cấp sự cân bằng tốt giữa độ chính xác và chi phí tính toán. Nó hiệu quả khi xử lý các hệ lớn hơn. Các tính toán này tạo ra các bề mặt thế năng tương tác. Bề mặt thế năng là cơ sở để mô tả lực giữa các phân tử. Sự lựa chọn phương pháp cụ thể phụ thuộc vào độ chính xác mong muốn và tài nguyên tính toán sẵn có. Các kỹ thuật này là công cụ không thể thiếu trong hóa học tính toán hiện đại.

Tương tác liên phân tử quyết định các tính chất vật lý của vật chất. Chúng bao gồm lực Van der Waals, tương tác tĩnh điện và liên kết hydro. Đối với cân bằng lỏng-hơi, độ mạnh của các tương tác này rất quan trọng. Chúng ảnh hưởng trực tiếp đến nhiệt hóa hơi và áp suất hơi bão hòa. Các tính toán hóa lượng tử định lượng chính xác năng lượng tương tác này. Chúng giúp xác định các tham số cho các hàm thế tương tác. Ví dụ, hàm thế Lennard-Jones hoặc Morse thường được sử dụng. Độ chính xác của các tham số này ảnh hưởng trực tiếp đến kết quả mô phỏng. Nghiên cứu sâu về tương tác liên phân tử cung cấp cái nhìn sâu sắc về cấu trúc chất lỏng. Điều này cũng giải thích các đặc tính như độ nén, độ nhớt và sức căng bề mặt. Hiểu rõ các tương tác này là chìa khóa để dự đoán hành vi pha của vật liệu mới.

Kỹ thuật Gibbs Ensemble Monte Carlo (GEMC) được thiết kế đặc biệt để nghiên cứu cân bằng pha. Nó tạo ra hai hộp mô phỏng riêng biệt, đại diện cho pha lỏng và pha hơi. Các hạt và thể tích có thể trao đổi giữa hai hộp. Điều này đảm bảo thế hóa học và áp suất bằng nhau trong cả hai pha. Các thao tác chính trong GEMC bao gồm di chuyển hạt, thay đổi thể tích hộp và chuyển hạt giữa các hộp. Mỗi thao tác được chấp nhận hoặc từ chối theo tiêu chí Monte Carlo. Mục tiêu là đạt được trạng thái cân bằng pha động. Ưu điểm chính của GEMC là tránh được việc mô phỏng giao diện pha. Giao diện pha thường khó xử lý và cần nhiều tài nguyên tính toán. GEMC cung cấp cách tiếp cận hiệu quả để thu được dữ liệu cân bằng lỏng-hơi. Phương pháp này được áp dụng rộng rãi để tính toán các đường cong cân bằng pha, điểm tới hạn và các tính chất khác.

Mô phỏng GEMC cho phép tính toán trực tiếp nhiều tính chất nhiệt động. Các tính chất này bao gồm mật độ của cả pha lỏng và pha hơi. Áp suất hơi bão hòa cũng được xác định từ các phép đo. Nhiệt hóa hơi, một đại lượng năng lượng quan trọng, cũng có thể được tính toán. Nhiệt hóa hơi phản ánh năng lượng cần thiết để chuyển một chất từ trạng thái lỏng sang trạng thái hơi. Việc tính toán này dựa trên sự khác biệt năng lượng trung bình giữa hai pha. Kết quả mô phỏng cung cấp giản đồ cân bằng lỏng-hơi. Giản đồ này biểu diễn sự phụ thuộc của áp suất hơi vào nhiệt độ. Chúng cũng cho phép xác định các tham số quan trọng như điểm tới hạn. Độ chính xác của các tính toán này phụ thuộc vào chất lượng của hàm thế tương tác được sử dụng. Các tính chất nhiệt động thu được rất hữu ích cho thiết kế quy trình công nghiệp và nghiên cứu cơ bản.

Quá trình xây dựng hàm thế bắt đầu bằng việc tính toán năng lượng tương tác ab initio. Các cặp phân tử (đime) được đặt ở nhiều cấu hình và khoảng cách khác nhau. Mỗi cấu hình tương ứng với một điểm trên bề mặt thế năng. Phương pháp hóa lượng tử cấp cao, như MP2 hoặc CCSD(T), được sử dụng để có độ chính xác cao. Dữ liệu năng lượng tương tác thu được sau đó được sử dụng để khớp với một hàm thế phân tích. Các hàm thế phổ biến bao gồm Lennard-Jones (LJ), Morse, hoặc các dạng phức tạp hơn. Hàm thế Lennard-Jones mô tả tốt tương tác Van der Waals. Nó bao gồm cả lực hút và lực đẩy. Việc xác định các tham số của hàm thế (ví dụ: epsilon, sigma trong LJ) là rất quan trọng. Các thuật toán tối ưu hóa được sử dụng để tìm bộ tham số tốt nhất. Chúng giảm thiểu sai số giữa năng lượng tính toán lượng tử và năng lượng từ hàm thế phân tích. Quá trình này đảm bảo hàm thế đại diện chính xác cho tương tác thực tế.

Sau khi xây dựng, hàm thế tương tác cần được đánh giá độ chính xác. Điều này bao gồm việc so sánh năng lượng do hàm thế dự đoán với các điểm năng lượng ab initio gốc. Các chỉ số thống kê như sai số bình phương trung bình (RMSE) được sử dụng. Hàm thế cũng được kiểm tra khả năng tái tạo các tính chất vật lý khác. Ví dụ, hệ số virial bậc hai là một đại lượng nhiệt động nhạy cảm với hàm thế. Việc so sánh hệ số virial tính toán từ hàm thế với dữ liệu thực nghiệm hoặc phương trình trạng thái giúp đánh giá độ tin cậy. Nếu hàm thế không đủ chính xác, quá trình tối ưu hóa tham số cần được lặp lại. Mục tiêu là đạt được sự phù hợp tốt nhất có thể với dữ liệu hóa lượng tử. Đánh giá kỹ lưỡng đảm bảo hàm thế đủ mạnh để sử dụng trong mô phỏng phân tử lớn hơn. Nó là yếu tố quyết định chất lượng của các kết quả mô phỏng về sau.

Áp suất hơi bão hòa là áp suất của pha hơi khi nó cân bằng với pha lỏng ở một nhiệt độ nhất định. Đại lượng này rất quan trọng trong nhiều ứng dụng công nghiệp. Trong mô phỏng GEMC, áp suất hơi bão hòa được xác định trực tiếp từ áp suất trung bình trong pha hơi. Áp suất được tính toán bằng phương pháp virial hoặc từ các phương trình trạng thái. Các giá trị áp suất hơi được tính toán ở nhiều nhiệt độ khác nhau. Điều này cho phép xây dựng đường cong áp suất hơi bão hòa. Đường cong này có thể được so sánh trực tiếp với dữ liệu thực nghiệm. Độ chính xác của các giá trị này phụ thuộc vào chất lượng của hàm thế và thời gian mô phỏng. Kết quả xác định áp suất hơi bão hòa cung cấp thông tin quý giá. Nó giúp hiểu rõ hơn về tính bay hơi của các chất. Điều này hỗ trợ việc thiết kế và vận hành các hệ thống bay hơi, chưng cất.

Nhiệt hóa hơi là năng lượng cần thiết để chuyển một mol chất lỏng thành pha hơi ở nhiệt độ và áp suất không đổi. Đây là một đại lượng nhiệt động quan trọng. Nó phản ánh độ mạnh của tương tác liên phân tử trong pha lỏng. Giá trị nhiệt hóa hơi cao hơn cho thấy tương tác liên phân tử mạnh hơn. Trong mô phỏng GEMC, nhiệt hóa hơi được tính từ sự khác biệt năng lượng nội trung bình của hai pha. Nó được xác định từ năng lượng tiềm năng trung bình và các hiệu chỉnh khác. Kết quả nhiệt hóa hơi cung cấp cái nhìn định lượng về quá trình chuyển pha. Nó giúp đánh giá mức độ ổn định của pha lỏng. Các giá trị này cũng được sử dụng để kiểm tra tính nhất quán của mô hình. So sánh với dữ liệu thực nghiệm cho thấy mức độ chính xác của các dự đoán. Phân tích này là yếu tố then chốt để hiểu sâu sắc về cân bằng lỏng-hơi.

Mô hình COSMO (COnductor-like Screening MOdel) là một phương pháp phổ biến. Nó dùng để mô tả ảnh hưởng của dung môi lên các phân tử. COSMO xử lý dung môi như một môi trường điện môi liên tục. Nó tạo ra một vỏ bọc dẫn điện bao quanh phân tử chất tan. Điều này cho phép tính toán năng lượng solvat hóa một cách hiệu quả. Việc tích hợp COSMO vào tính toán cân bằng lỏng-hơi mở rộng khả năng. Nó giúp mô phỏng các hệ thống trong dung dịch hoặc pha lỏng đặc. Mô hình này đặc biệt hữu ích khi tương tác giữa chất tan và dung môi là đáng kể. Nó cải thiện độ chính xác của các tính toán nhiệt động trong môi trường phức tạp. COSMO cung cấp cái nhìn chân thực hơn về hành vi của các chất trong pha lỏng. Điều này quan trọng cho các ứng dụng trong hóa dược và hóa sinh.

Mạng thần kinh nhân tạo (ANN) là một công cụ mạnh mẽ trong lĩnh vực trí tuệ nhân tạo. Nó có khả năng học hỏi từ dữ liệu và nhận dạng các mẫu phức tạp. Trong hóa học tính toán, ANN có thể được sử dụng để tối ưu hóa hàm thế. Nó cũng có thể dự đoán các tính chất nhiệt động nhanh chóng. ANN có thể học mối quan hệ giữa cấu trúc phân tử và năng lượng tương tác. Sau khi được huấn luyện trên một tập dữ liệu hóa lượng tử lớn, ANN có thể dự đoán năng lượng cho các cấu hình mới. Điều này nhanh hơn nhiều so với tính toán lượng tử trực tiếp. Ứng dụng ANN giúp tăng tốc quá trình mô phỏng phân tử. Nó mở ra cánh cửa cho việc nghiên cứu các hệ thống lớn hơn và phức tạp hơn. ANN cũng hỗ trợ việc xác định các tham số hàm thế một cách hiệu quả. Đây là một bước tiến quan trọng trong việc tích hợp học máy vào nghiên cứu khoa học vật liệu.

Hệ số virial bậc hai (B2) là một đại lượng nhiệt động quan trọng. Nó mô tả tương tác giữa hai phân tử trong pha khí loãng. B2 cung cấp thông tin về độ lệch của khí thực so với khí lý tưởng. Hệ số virial bậc hai được tính toán từ các hàm thế tương tác đã xây dựng. Các giá trị này được so sánh với dữ liệu thực nghiệm. Chúng cũng được so sánh với các giá trị từ phương trình trạng thái đã biết. Sự phù hợp tốt giữa các giá trị tính toán và thực nghiệm khẳng định độ chính xác của hàm thế. Nó cũng chứng minh hiệu quả của phương pháp lượng tử được sử dụng. Phân tích này là một bước kiểm tra quan trọng cho toàn bộ quá trình mô hình hóa. Nó giúp đánh giá mức độ tin cậy của các kết quả cân bằng pha khác. Hệ số virial bậc hai là một chỉ số nhạy cảm. Nó phản ánh chất lượng của mô tả tương tác liên phân tử.

Giản đồ cân bằng lỏng-hơi là biểu đồ minh họa sự phụ thuộc của áp suất hơi vào nhiệt độ. Chúng cũng thể hiện mật độ của pha lỏng và pha hơi. Các giản đồ này đã được xây dựng cho Ar, N2, Cl2, CO từ mô phỏng GEMC. Các kết quả mô phỏng cho thấy sự phù hợp tốt với dữ liệu thực nghiệm có sẵn. Điều này chứng tỏ khả năng dự đoán của phương pháp. Giản đồ cung cấp thông tin về điểm tới hạn và đường cong bão hòa. Đường cong bão hòa cho thấy nhiệt độ và áp suất mà tại đó hai pha tồn tại cân bằng. Phân tích các giản đồ này giúp hiểu rõ hơn về hành vi pha của từng chất. Sự khác biệt trong giản đồ giữa các chất phản ánh sự khác biệt trong tương tác liên phân tử của chúng. Các giản đồ này là tài liệu tham khảo quý giá cho các ứng dụng kỹ thuật và nghiên cứu khoa học.