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

Tổng quan về luận án

Trong khoa học nhiệt động học và hóa lý tính toán, việc xác định chính xác các đặc trưng cân bằng lỏng-hơi (Vapor-Liquid Equilibria - VLE) cùng các thông số trạng thái tới hạn của các đơn chất và hỗn hợp khí đóng vai trò nền tảng đối với thiết kế công nghệ hóa chất, công nghiệp luyện kim, xử lý môi trường và hóa dầu. Luận án tiến sĩ hóa học của nghiên cứu sinh Nguyễn Thành Được với đề tài "Tính toá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 toàn cục Monte Carlo" (chuyên ngành Hóa lý thuyết và Hóa lý, mã số 944.19, thực hiện tại Trường Đại học Khoa học – Đại học Huế dưới sự hướng dẫn của PGS. Phạm Văn Tất và PGS. Trần Dương) đã thiết lập một quy trình đa tầng tiên phong, kết nối chặt chẽ giữa tính toán cấu trúc điện tử lượng tử từ nguyên lý ban đầu (ab initio), xây dựng mặt thế tương tác liên phân tử phi thực nghiệm, tích phân nhiệt động thống kê và mô phỏng cơ học thống kê hiện đại.

Khoảng trống nghiên cứu (research gap) cốt lõi mà công trình hướng đến giải quyết là sự thiếu hụt các hàm thế tương tác cặp khả dĩ có độ chính xác cao xuất phát hoàn toàn từ lý thuyết lượng tử cấp cao (first-principles), cho phép mô tả chính xác cả vùng tương tác tầm ngắn (lực đẩy xen phủ đám mây electron), vùng cực tiểu thế năng và vùng phân tán tầm xa đối với các phân tử hai nguyên tử có tính bất đối xứng tĩnh điện như cacbon monoxit ($\text{CO}$), các phân tử có độ phân cực cao như clo ($\text{Cl}_2$), cũng như các khí công nghiệp thiết yếu như nitơ ($\text{N}_2$) và argon ($\text{Ar}$). Hơn nữa, việc đo đạc thực nghiệm các thông số nhiệt động của các chất khí độc hại hoặc ăn mòn mạnh như $\text{Cl}_2$ và $\text{CO}$ ở điều kiện nhiệt độ và áp suất cực hạn thường gặp rủi ro an toàn nghiêm trọng và chi phí đo đạc rất đắt đỏ.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập rõ ràng:

1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một mô hình thế năng tương tác cặp liên phân tử ab initio đa tâm (5-site model) có khả năng tái hiện chính xác bề mặt thế năng lượng tử của các đime $\text{Ar}-\text{Ar}$, $\text{N}_2-\text{N}_2$, $\text{Cl}_2-\text{Cl}_2$ và $\text{CO}-\text{CO}$ trên toàn bộ không gian cấu hình định hướng?
2. Câu hỏi nghiên cứu 2 (RQ2): Liệu các hiệu chỉnh lượng tử bậc nhất kết hợp với tích phân cầu phương bốn chiều Gauss-Legendre có thể nâng cao độ chuẩn xác của việc tính toán hệ số virial bậc hai $B_2(T)$ để kiểm chứng hàm thế trước khi đưa vào mô phỏng pha ngưng tụ hay không?
3. Câu hỏi nghiên cứu 3 (RQ3): Kỹ thuật mô phỏng ngẫu nhiên toàn cục Gibbs Ensemble Monte Carlo (GEMC-NVT) dựa trên hàm thế ab initio mới có thể tái lập giản đồ pha lỏng-hơi, áp suất hơi bão hòa, enthalpy bay hơi và các điểm tới hạn ($T_c, P_c, \rho_c$) với độ lệch tương đối đạt chuẩn quy chuẩn thực nghiệm quốc tế hay không?
4. Giả thuyết khoa học (H1): Việc phân bố các vị trí tương tác ngoài hạt nhân (gồm 2 tâm nguyên tử, 1 tâm trọng khối và 2 tâm giả phân đoạn) kết hợp hàm giảm chấn Damping sẽ triệt tiêu hiện tượng kỳ dị năng lượng phân tán tầm ngắn, cho phép hàm thế phản ánh trung thực bản chất tương quan electron cấp cao $\text{CCSD(T)}$.
5. Giả thuyết khoa học (H2): Mô hình mạng thần kinh nhân tạo truyền thẳng kết hợp phân tích thành phần chính (PCA-ANN) có thể đóng vai trò như một phương pháp phi tuyến xấp xỉ phổ quát để dự báo chính xác các hệ số nhiệt động của phương trình trạng thái virial từ các thông số tới hạn cơ bản.

Khung lý thuyết của luận án tích hợp chặt chẽ: lý thuyết tương quan electron nhóm cặp kích thích đơn, đôi và nhiễu loạn ba $\text{CCSD(T)}$; lý thuyết mô phỏng đại tập chính tắc Gibbs (GEMC); lý thuyết nhiệt động solvat hóa mô hình sàng lọc vật dẫn (COSMO/COSMO-SAC); và lý thuyết phương trình trạng thái virial, Peng–Robinson (PR-EOS) và Deiters (D-EOS). Phạm vi nghiên cứu khảo sát toàn diện 4 hệ khí thuần khiết ($\text{Ar}$, $\text{N}_2$, $\text{Cl}_2$, $\text{CO}$) và hệ hỗn hợp nhị phân ($\text{CO}-\text{Cl}_2$) trải dài trên dải nhiệt độ rộng từ trạng thái dưới điểm sôi đến vùng lân cận điểm tới hạn. Ý nghĩa đột phá thể hiện qua việc giảm thiểu sai số dự báo mật độ pha lỏng/hơi và áp suất hơi xuống mức sai số trung bình tương đối (MARE) dưới $5%$, chứng minh tính khả thi của việc thay thế một phần thực nghiệm tốn kém bằng tính toán mô phỏng lượng tử độ tin cậy cao.

Literature Review và Positioning

Lịch sử tính toán nhiệt động học chất lưu chứng kiến sự chuyển dịch mạnh mẽ từ các phương trình thực nghiệm giải tích sang mô phỏng máy tính dựa trên cấu trúc vi mô. Nghiên cứu tiên phong của Michels (1960) về hệ số virial của hiđro và đơteri, cùng công trình kinh điển tổng hợp dữ liệu hệ số virial của Dymond và Smith (1980) đã đặt nền móng cho việc sử dụng hệ số virial bậc hai như một thước đo trực tiếp đánh giá lực liên phân tử. Bước ngoặt phương pháp luận xuất hiện khi Panagiotopoulos (1987) phát minh ra kỹ thuật mô phỏng Gibbs Ensemble Monte Carlo (GEMC), cho phép tính toán trực tiếp cân bằng pha lỏng-hơi mà không cần mô phỏng tường minh bề mặt phân cách pha giữa hai trạng thái. Tiếp theo đó, Spyriouni (1998) đã mở rộng GEMC cho các hydrocacbon mạch dài ($n$-hexadecan), trong khi Martin và Siepmann (1998, 1999) hoàn thiện mô hình thế chuyển tiếp liên kết (TraPPE) để tính toán năng lượng tự do Gibbs cho họ ankan.

                  TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VÀ VỊ TRÍ LUẬN ÁN
                  
   Dymond & Smith (1980) ────► Panagiotopoulos (1987) ────► Martin & Siepmann (1998)
   (Dữ liệu Virial thực nghiệm)   (Kỹ thuật GEMC cân bằng pha)    (Hàm thế thực nghiệm TraPPE)
                                                                       │
   Leonhard & Deiters (2000-2002) ◄────────────────────────────────────┘
   (Thế Ab Initio cho Ne, Ar, N2)
             │
             ▼
   ══════════════════════════════════════════════════════════════════════════════
   VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN (Nguyễn Thành Được, 2020):
   • Mở rộng mô hình thế Ab Initio 5-vị trí cho phân tử dị hạch (CO) và Halogen (Cl2).
   • Tích hợp: CCSD(T)/aug-cc-pVmZ ──► Cầu phương Gauss 4D ──► GEMC-NVT ──► COSMO-SAC/ANN.
   ══════════════════════════════════════════════════════════════════════════════

Trong dòng nghiên cứu ứng dụng thế ab initio, Leonhard và Deiters (2000, 2002) đã đạt được những kết quả đáng kể khi sử dụng năng lượng tương tác lượng tử để mô phỏng GEMC cho neon, argon và nitơ. Song song đó, Garrison và Sandler (2002) chứng minh rằng năng lượng tương tác lượng tử có thể tính toán chính xác các tính chất nhiệt động học mà không cần dựa vào dữ liệu khớp thực nghiệm. Tuy nhiên, các tranh luận học thuật sâu sắc vẫn tồn tại xung quanh hai trường phái:

• Trường phái thế tham số hóa thực nghiệm (Empirical Force Fields): Cho rằng việc tinh chỉnh tham số thế năng (như Lennard-Jones 12-6) trực tiếp từ dữ liệu thực nghiệm mang lại tốc độ tính toán nhanh và độ khớp cao trong phạm vi hẹp của điều kiện khảo sát, nhưng hoàn toàn mất tính khái quát khi ngoại suy ra ngoài miền dữ liệu hoặc khi áp dụng cho các điều kiện tới hạn cực đoan.
• Trường phái thế nguyên lý ban đầu (Ab Initio Potentials): Khẳng định chỉ có tính toán cấu trúc điện tử lượng tử cấp cao mới phản ánh đúng bản chất vật lý của tương tác tán xạ, phân cực và đẩy Pauli, dù đòi hỏi chi phí tính toán rất lớn và độ phức tạp cao trong giải thuật khớp mặt thế năng (PES).

Luận án của NCS. Nguyễn Thành Được định vị chính xác tại giao điểm này, giải quyết triệt để hạn chế của các nghiên cứu quốc tế tiền nhiệm. Trong khi các nghiên cứu của Leonhard và Deiters (2000, 2002) chủ yếu giới hạn ở các phân tử đồng hạch đối xứng đơn giản ($\text{N}_2$, $\text{Ar}$), công trình này đã mở rộng biên giới nghiên cứu sang phân tử dị hạch có momen lưỡng cực vĩnh cửu ($\text{CO}$) và phân tử halogen có lớp vỏ electron cồng kềnh với độ phân cực cao ($\text{Cl}_2$). Bằng việc đối chiếu với hai nghiên cứu quốc tế tiêu biểu—nghiên cứu của Wormer (2005) về hệ số virial của các phân tử bất đối xứng và nghiên cứu của William Ferreira (2005) về khớp mặt thế bằng thuật toán di truyền—công trình này đã tạo ra một bước tiến vượt bậc: kết hợp thuật toán tối ưu hóa toàn cục di truyền (Genetic Algorithm) với thuật toán cục bộ Levenberg-Marquardt để xác định bộ tham số thế 5-vị trí đạt hệ số tương quan $R^2 > 0.999$, đồng thời tích hợp mô hình nhiệt động solvat hóa COSMO-SAC (Lin & Sandler, 2002) để khảo sát cân bằng lỏng-hơi của hỗn hợp nhị phân $\text{CO}-\text{Cl}_2$.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp trực tiếp vào việc kiểm chứng, bổ khuyết và mở rộng lý thuyết tương tác liên phân tử và cơ học thống kê chất lưu thông qua các khía cạnh:

1. Mở rộng lý thuyết nhóm cặp Coupled Cluster $\text{CCSD(T)}$: Chứng minh năng lực của lý thuyết $\text{CCSD(T)}$ kết hợp bộ hàm cơ sở tương quan mở rộng Dunning $\text{aug-cc-pVmZ}$ ($m = \text{D, T, Q, 5}$) trong việc nắm bắt toàn diện năng lượng tương quan electron của các cấu hình đime bất đối xứng.
2. Xây dựng mô hình hình học 5-vị trí (5-site potential representation): Đối với các phân tử hai nguyên tử ($\text{N}_2$, $\text{Cl}_2$, $\text{CO}$), mô hình không chỉ đặt các tâm thế tại 2 hạt nhân nguyên tử mà còn bổ sung 1 tâm tại trọng tâm phân tử (M) và 2 tâm phụ nằm tại trung điểm giữa hạt nhân và trọng tâm. Cấu trúc này cho phép mô tả chính xác sự phân bố dị hướng của mật độ điện tích và đám mây electron.
3. Phát triển dạng hàm thế biến tính Morse-Damping: Khắc phục nhược điểm sụp đổ năng lượng ở cự ly cực ngắn của thế phân tán London cổ điển bằng cách tích hợp hai hàm giảm chấn $D_{ij}^{M6}(r)$ và $D^{CC}(r)$: $$D_{ij}^{M6}(r) = \exp\left[-f_d \left(1 - \frac{r}{r_{m,ij}}\right)^2\right] \quad \text{khi } r < r_d$$
4. Chứng minh lý thuyết hiệu chỉnh lượng tử bậc nhất: Chỉ ra rằng đối với các phân tử nhẹ hoặc ở vùng nhiệt độ thấp, hệ số virial cổ điển $B_{\text{cl}}^0(T)$ phải được bổ sung các số hạng hiệu chỉnh lượng tử từ chuyển động tịnh tiến, quay và mômen quán tính: $$B_2(T) = B_{\text{cl}}^0(T) + B_r^1(T) + B_a^{1I}(T) + B_a^{1\mu}(T)$$

┌─────────────────────────────────────────────────────────────────────────────┐
│                    KHUNG PHÂN TÍCH TỔNG THỂ CỦA ĐỀ TÀI                      │
└─────────────────────────────────────────────────────────────────────────────┘
                                       │
      ┌────────────────────────────────┴────────────────────────────────┐
      ▼                                                                 ▼
[TẦNG VI MÔ: HÓA LƯỢNG TỬ]                                    [TẦNG TRUNG MÔ: MẶT THẾ & VIRIAL]
• CCSD(T)/aug-cc-pVmZ                                         • Mô hình 5-vị trí (Morse + Damping)
• Khử sai số BSSE                                             • Thuật toán Di truyền + Levenberg-Marquardt
• Tối ưu độ dài liên kết:                                     • Tích phân cầu phương Gauss-Legendre 4D
  r(C-O)=1.128 Å; r(N-N)=1.097 Å; r(Cl-Cl)=2.009 Å            • Đánh giá hệ số Virial B2(T)
      │                                                                 │
      └────────────────────────────────┬────────────────────────────────┘
                                       │
                                       ▼
                         [TẦNG VĨ MÔ: MÔ PHỎNG & NHIỆT ĐỘNG]
                         • GEMC-NVT (2 hộp, 3 loại dịch chuyển)
                         • Phương trình trạng thái (PR-EOS, D-EOS)
                         • Mô hình dung môi COSMO-SAC
                         • Mạng nơron PCA-ANN: I(5)-HL(6)-O(3)
                                       │
                                       ▼
                         [KẾT QUẢ ĐẦU RA NHIỆT ĐỘNG]
                         • Giản đồ pha VLE lỏng - hơi
                         • Thông số tới hạn: Tc, Pc, Vc, ρc
                         • Enthalpy bay hơi ΔHvap, Áp suất hơi Pvap

Khung phân tích độc đáo

Khung phân tích của đề tài là sự dung hợp đa ngành giữa ba trụ cột lý thuyết: Hóa học lượng tử cấu trúc điện tử, Cơ học thống kê Gibbs, và Học máy thống kê (Machine Learning).

Khung phân tích vận hành tuần tự qua 4 bước nghiêm ngặt:

• Bước 1 (Mức vi mô): Quét bề mặt năng lượng thế đa chiều thông qua việc giải phương trình Schrödinger phi tuyến cho các cấu hình định hướng không gian đặc thù: Tuyến tính (L - Linear), Song song (H - Parallel), Chữ T (T - T-shaped), và Chéo (X - Crossed).
• Bước 2 (Mức trung gian): Khớp năng lượng ab initio vào hàm thế giải tích bằng giải thuật tối ưu kép (Genetic Algorithm lọc cực tiểu toàn cục $\rightarrow$ Levenberg-Marquardt tinh chỉnh tham số).
• Bước 3 (Kiểm định thống kê): Kiểm định độ tin cậy của hàm thế qua tích phân cầu phương Gauss-Legendre 4 chiều để thu được hệ số virial $B_2(T)$ so sánh với thực nghiệm.
• Bước 4 (Mức vĩ mô): Thực thi mô phỏng GEMC-NVT trên hệ hai hộp pha lỏng và hơi để trích xuất các đường cong mật độ bão hòa, enthalpy bay hơi và thông số tới hạn thông qua quan hệ tỷ lệ tới hạn với số mũ $\beta = 0.325$.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được thiết lập theo thế giới quan thực chứng lượng hóa (positivism paradigm) kết hợp chủ nghĩa hiện thực phản biện (critical realism), lấy quy luật cơ học lượng tử và cơ học thống kê làm nền tảng khách quan giải thích các hiện tượng nhiệt động vĩ mô. Thiết kế nghiên cứu đa cấp độ (multi-level design) liên kết trực tiếp các thang đo thời gian và không gian: từ bán kính hạt nhân và mật độ điện tích nguyên tử ($10^{-10} \text{ m}$ / Angstrom) đến các đại lượng nhiệt động vĩ mô của hệ đa hạt ($10^{23}$ hạt).

Các thông số hình học phân tử được cố định chính xác dựa trên cấu trúc hình học tối ưu cân bằng:

• Phân tử $\text{CO}$: Độ dài liên kết thực nghiệm $r(\text{C}-\text{O}) = 1.128 \text{ \AA}$.
• Phân tử $\text{N}_2$: Độ dài liên kết thực nghiệm $r(\text{N}-\text{N}) = 1.097 \text{ \AA}$.
• Phân tử $\text{Cl}_2$: Độ dài liên kết thực nghiệm $r(\text{Cl}-\text{Cl}) = 2.009 \text{ \AA}$.
• Nguyên tử $\text{Ar}$: Mô hình hóa dưới dạng khối cầu đối xứng hoàn hảo.

                  CÁC CẤU HÌNH ĐỊNH HƯỚNG KHÔNG GIAN CỦA ĐIME
                  
       (L) Tuyến tính                  (H) Song song
   ───●──────●───────●──────●───       ───●──────────●───
                                          │          │   r
                                          │          │
                                       ───●──────────●───
   
       (T) Chữ T                       (X) Chéo góc
          ●                            ───●──────────●───
          │                                  ╲    ╱
        r │                                    ╳    r
          │                                  ╱    ╲
   ───●───┴───●───                     ───●──────────●───

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các giao thức chuẩn hóa cao độ với sự kiểm soát sai số chặt chẽ:

┌─────────────────────────────────────────────────────────────────────────────┐
│             TIẾN TRÌNH THỰC THI PHƯƠNG PHÁP NGHIÊN CỨU CHI TIẾT             │
└─────────────────────────────────────────────────────────────────────────────┘
  [1] GAUSSIAN 03W
      └─ Tính toán CCSD(T) / aug-cc-pVmZ
      └─ Hiệu chỉnh sai số chồng chập bộ hàm (BSSE) theo Boys-Bernardi
  [2] TỐI ƯU HÓA THAM SỐ THẾ (Auto2Fit 3.0 & C/C++)
      └─ Thuật toán Di truyền (Genetic Algorithm) ──► Tìm vùng đáy toàn cục
      └─ Thuật toán Levenberg-Marquardt ──► Tinh chỉnh chính xác cực tiểu
  [3] TÍCH PHÂN CẦU PHƯƠNG GAUSS-LEGENDRE 4 CHIỀU
      └─ Giải hệ phương trình tích phân hệ số Virial B2(T)
      └─ Hiệu chỉnh lượng tử bậc nhất: tịnh tiến, quay, mômen quán tính
  [4] MÔ PHỎNG TOÀN CỤC GEMC-NVT
      └─ Thiết lập 2 hộp (Hộp I: Lỏng, Hộp II: Hơi)
      └─ 3 bước dịch chuyển Markov: Dịch chuyển hạt - Đổi thể tích - Trao đổi hạt
  [5] MÔ HÌNH DUNG MÔI COSMO-SAC & HỌC MÁY PCA-ANN
      └─ Trích xuất mật độ điện tích bề mặt sigma (σ-profile)
      └─ Huấn luyện mạng truyền thẳng I(5)-HL(6)-O(3) với thuật toán lan truyền ngược

Độ tin cậy và giá trị nội tại (internal validity) được đảm bảo thông qua kỹ thuật tam giác giác hóa phương pháp (methodological triangulation): kết quả mật độ pha và áp suất cân bằng từ GEMC-NVT được đối chiếu chéo đồng thời với dữ liệu thực nghiệm tiêu chuẩn của Viện Tiêu chuẩn và Công nghệ Quốc gia Hoa Kỳ (NIST), phương trình trạng thái cải tiến Peng-Robinson (PR-EOS), phương trình Deiters (D-EOS) và mô hình solvat hóa COSMO-SAC.

Data và phân tích

Hệ thống công cụ phần mềm chuyên dụng được lập trình và tích hợp đồng bộ:

• Gaussian 03W™: Thực hiện các phép tính phiếm hàm mật độ và tương quan electron cấp cao $\text{CCSD(T)}$.
• Auto2Fit 3.0 & Trình biên dịch C/C++: Lập trình giải thuật di truyền và Levenberg-Marquardt tối ưu hóa phi tuyến tính các tham số thế tương tác.
• GEMC-NVT Engine (C/C++): Chương trình mô phỏng đại tập Gibbs chuyên dụng, thực thi hàng triệu bước lặp Monte Carlo để đạt trạng thái cân bằng nhiệt động và cân bằng hóa thế giữa hai pha.
• ThermoC Online Tool: Hỗ trợ tính toán đối sánh phương trình trạng thái nhiệt động học chất lưu.

Trong mô hình học máy, mạng thần kinh nhân tạo ANN được tích hợp với phân tích thành phần chính (PCA) với cấu hình mạng $\text{I(5)-HL(6)-O(3)}$:

• Lớp đầu vào (5 nút): Gồm 5 đặc trưng nhiệt động học tới hạn $[P_c, V_c, T_c, T_L, T_U]$ được trực giao hóa qua biến đổi PCA nhằm loại bỏ hoàn toàn hiện tượng đa cộng tuyến.
• Lớp ẩn (6 nơron): Sử dụng hàm kích hoạt phi tuyến tính dạng sigmoid, huấn luyện bằng giải thuật lan truyền ngược có giám sát (backpropagation).
• Lớp đầu ra (3 nút): Dự báo chính xác các hệ số trạng thái $a, b, c$ của phương trình virial dạng logarit: $$\ln B_2(T) = a + b \exp\left(\frac{c}{T}\right)$$ Các chỉ số đánh giá thống kê được tính toán nghiêm ngặt: căn bậc hai sai số bình phương trung bình (RMSE), sai số tương đối tuyệt đối (ARE, %) và sai số trung bình tương đối (MARE, %): $$\text{MARE} = \frac{1}{n}\sum_{i=1}^n \left|\frac{x_{\text{origin}} - x_{\text{pred}}}{x_{\text{origin}}}\right| \times 100%$$

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Xây dựng thành công bề mặt thế năng lượng tử 5-vị trí có độ chính xác cao: Các bề mặt thế năng ab initio của các đime $\text{Ar}-\text{Ar}$, $\text{N}_2-\text{N}_2$, $\text{Cl}_2-\text{Cl}_2$ và $\text{CO}-\text{CO}$ đã tái hiện sắc nét lực hút phân tán và lực đẩy xen phủ tại mọi góc định hướng. Trích dẫn nguyên văn khẳng định phương pháp luận của luận án: "Phương pháp mô phỏng toàn cục Gibbs Ensemble Monte Carlo (GEMC) không giống như các phương pháp mô phỏng chính xác khác ở chỗ nó là một phương pháp mô phỏng ngẫu nhiên, phương pháp này không đi giải phương trình tường minh mà đúng hơn là nhận các kết quả bằng mô phỏng các hạt riêng lẻ và ghi lại các tham số của trạng thái trung bình của chúng." Độ khớp giữa năng lượng ab initio và năng lượng tính từ hàm thế đạt hệ số xác định $R^2 > 0.999$.

2. Xác định chính xác hệ số virial bậc hai $B_2(T)$ qua hiệu chỉnh lượng tử: Luận án nhấn mạnh vai trò sống còn của việc đánh giá virial: "Hệ số virial bậc hai B2(T) chỉ phụ thuộc vào sự tương tác cặp giữa các hạt. Do đó các hệ số virial thể hiện tầm quan trọng trong việc đánh giá độ tin cậy của các hàm thế cặp tương tác phân tử." Các giá trị $B_2(T)$ tính toán cho argon, nitơ, clo và cacbon monoxit trên dải nhiệt độ từ $100 \text{ K}$ đến $1000 \text{ K}$ hoàn toàn trùng khớp với dữ liệu thực nghiệm chuẩn của Dymond và Smith, khẳng định tính đúng đắn của bộ tham số thế trước khi bước vào mô phỏng pha ngưng tụ.

3. Tái lập hoàn hảo giản đồ cân bằng lỏng-hơi bằng GEMC-NVT: Mô phỏng GEMC-NVT sử dụng hàm thế 5-vị trí mới đã dựng nên toàn bộ đường cong bao pha lỏng-hơi cho $\text{N}_2$ và $\text{CO}$. Mật độ pha lỏng bão hòa $\rho_L$ và mật độ pha hơi $\rho_V$ ở các mức nhiệt độ khác nhau cho thấy độ sai lệch MARE cực thấp:

   • Đối với nitơ ($\text{N}_2$): $\text{MARE}(\rho_L) = 1.24%$, $\text{MARE}(\rho_V) = 3.15%$.
   • Đối với cacbon monoxit ($\text{CO}$): $\text{MARE}(\rho_L) = 1.48%$, $\text{MARE}(\rho_V) = 3.62%$.
   • Tái hiện chính xác các hàm phân bố xuyên tâm cặp nguyên tử $g(r_{\text{N}-\text{N}})$, $g(r_{\text{C}-\text{C}})$, $g(r_{\text{O}-\text{O}})$, $g(r_{\text{C}-\text{O}})$ và $g(r_{\text{M}-\text{M}})$, phản ánh cấu trúc trật tự tầm ngắn của chất lỏng bão hòa.

4. Khám phá và mô hình hóa cân bằng pha hỗn hợp nhị phân $\text{CO}-\text{Cl}_2$ qua COSMO-SAC: Lần đầu tiên tại Việt Nam, sự phân bố mật độ điện tích bề mặt sàng lọc ($\sigma$-profile) của phân tử đơn $\text{CO}$ và $\text{Cl}_2$ được trích xuất hoàn chỉnh từ hóa lượng tử để tính toán hệ số hoạt độ pha lỏng $\gamma_i$ theo mô hình nhiệt động NRTL cải tiến. Giản đồ cân bằng pha đẳng nhiệt ($P-x-y$) và đẳng áp ($T-x-y$) của hệ hỗn hợp $\text{CO}-\text{Cl}_2$ được thiết lập với độ chính xác tương đương thực nghiệm, cung cấp bức tranh nhiệt động toàn diện cho hệ phản ứng tổng hợp phosgene ($\text{COCl}_2$).

5. Hiệu năng vượt trội của mô hình học máy PCA-ANN $\text{I(5)-HL(6)-O(3)}$: Mạng nơron nhân tạo đã phục hồi xuất sắc các hệ số nhiệt động $a, b, c$ từ dữ liệu nén PCA với hệ số sai số $\text{RMSE} < 0.008$ và hệ số tương quan tuyến tính giữa giá trị tính toán và thực nghiệm đạt $R^2 = 0.9998$.

    GIẢN ĐỒ CÂN BẰNG PHA LỎNG - HƠI VÀ ĐIỂM TỚI HẠN DỰ BÁO TỪ GEMC-NVT
    
    Áp suất P / Nhiệt độ T
         ▲
         │                       ● Điểm tới hạn (Tc, Pc)
         │                     ╱   ╲
         │                   ╱       ╲
         │       Pha lỏng  ╱           ╲  Pha hơi
         │       (ρL)    ╱               ╲ (ρV)
         │              │                 │
         │              │                 │
         │              │                 │
         └──────────────┴─────────────────┴──────────────► Mật độ ρ (g/cm³)
                        ◄── Vùng hai pha ──►

Implications đa chiều

• Về mặt lý thuyết: Chứng minh rằng mô hình thế cặp 5-vị trí kết hợp hàm giảm chấn phi tuyến đủ năng lực đại diện cho tương tác liên phân tử trong mô phỏng đại tập Gibbs, mở ra hướng tiếp cận mô hình hóa pha ngưng tụ mà không cần dựa vào các tham số điều chỉnh thực nghiệm.
• Về mặt phương pháp luận: Thiết lập một quy trình chuẩn hóa (pipeline) có khả năng chuyển giao trực tiếp: Gaussian (Quantum) -> Auto2Fit (Global Fitting) -> Gauss Quadrature (Virial) -> GEMC Engine (Phase Equilibrium) -> COSMO/ANN (Thermodynamic Screening).
• Về mặt ứng dụng thực tiễn: Cung cấp bộ dữ liệu nhiệt động học chuẩn xác cao về tỷ trọng bão hòa, độ bay hơi tương đối và áp suất hơi của các khí công nghiệp độc hại ($\text{Cl}_2$, $\text{CO}$), phục vụ thiết kế tháp chưng cất, tháp hấp thụ và bồn chứa áp lực cao trong công nghệ hóa chất và lọc hóa dầu.
• Về mặt chính sách và an toàn công nghiệp: Hỗ trợ xây dựng các tiêu chuẩn định mức an toàn lưu trữ và vận chuyển khí hóa lỏng ở điều kiện áp suất cao dựa trên các thông số tới hạn chính xác.

Limitations và Future Research

Luận án ghi nhận một cách khách quan các giới hạn khoa học nội tại:

1. Bỏ qua đóng góp của tương tác ba vật thể (Three-body interactions): Mô hình mô phỏng hiện tại xây dựng trên giả định cộng tính đôi (pairwise additivity). Ở vùng mật độ chất lỏng rất cao gần điểm ba thể, đóng góp của thế đa thân Axilrod-Teller-Muto có thể chiếm từ $5-10%$ tổng năng lượng liên kết, gây ra độ lệch nhỏ ở vùng nhiệt độ cực thấp.
2. Giả định phân tử cứng (Rigid rotor approximation): Việc cố định độ dài liên kết phân tử ($r(\text{C}-\text{O}) = 1.128 \text{ \AA}$, $r(\text{N}-\text{N}) = 1.097 \text{ \AA}$, $r(\text{Cl}-\text{Cl}) = 2.009 \text{ \AA}$) chưa xét đến dao động dãn liên kết nội phân tử ở dải nhiệt độ rất cao.
3. Chi phí tính toán cấu trúc điện tử: Phương pháp $\text{CCSD(T)}$ với bộ hàm $\text{aug-cc-pVQZ}$ đòi hỏi tài nguyên tính toán theo cấp số nhân $O(N^7)$, tạo ra rào cản lớn khi mở rộng sang các hệ phân tử hữu cơ cồng kềnh hoặc chất lỏng ion.

Chương trình nghiên cứu tiếp nối trong tương lai định hướng vào 4 mũi nhọn:

• Tích hợp số hạng thế 3-thân phi cộng tính (non-additive three-body potentials) vào thuật toán GEMC để hoàn thiện độ chính xác ở vùng lân cận điểm ba thể.
• Mở rộng mô hình thế linh động (flexible potential models) cho phép liên kết hóa học co dãn theo hàm thế Morse nội phân tử.
• Ứng dụng các thuật toán học sâu hiện đại (Deep Message Passing Neural Networks) để biểu diễn trực tiếp mặt thế năng lượng tử nhiều chiều thay thế cho các hàm thế giải tích truyền thống.
• Mở rộng quy trình tính toán cho các hệ chất lỏng thu giữ carbon dioxide ($\text{CO}_2$), dung môi eutecti sâu (DES) và chất lỏng ion phục vụ công nghệ chuyển đổi năng lượng xanh.

Tác động và ảnh hưởng

Nghiên cứu tạo nên những tác động đa diện đối với cộng đồng khoa học và ứng dụng công nghệ:

• Tác động học thuật: Đặt cột mốc quan trọng cho phân ngành hóa lý thuyết và hóa tính toán tại Việt Nam; cung cấp một hệ phương pháp luận hoàn chỉnh kết hợp hóa lượng tử và cơ học thống kê đã được công bố trên các tạp chí chuyên ngành uy tín như Tạp chí Khoa học và Công nghệ (VAST), Tạp chí Hóa học (VAST), SmartScience.
• Chuyển đổi công nghệ và công nghiệp: Cung cấp thông số thiết kế chính xác cho ngành kỹ thuật hóa học, đặc biệt trong việc tính toán số đĩa lý thuyết, chiều cao tháp chưng cất tách khí không khí ($\text{N}_2/\text{Ar}/\text{O}_2$) và công nghệ tổng hợp hữu cơ sử dụng khí than quyển ($\text{CO}/\text{Cl}_2$).
• Lợi ích an toàn môi trường: Giảm thiểu tối đa việc phải tiến hành các thử nghiệm thực nghiệm nguy hiểm với khí độc $\text{CO}$ và $\text{Cl}_2$ ở nhiệt độ và áp suất cao, đóng góp thiết thực vào phát triển hóa học xanh và an toàn lao động.

Đối tượng hưởng lợi

• Nghiên cứu sinh và học viên cao học ngành Hóa lý, Vật lý tính toán: Thừa hưởng một khung quy trình mẫu mực về việc lập trình và tích hợp các công cụ mô phỏng từ ab initio đến Monte Carlo.
• Các nhà khoa học cao cấp và giảng viên đại học: Sử dụng các kết quả nghiên cứu và bộ tham số thế 5-vị trí như dữ liệu chuẩn để phát triển các lý thuyết phương trình trạng thái mới.
• Kỹ sư R&D trong công nghiệp Hóa chất, Luyện kim và Dầu khí: Ứng dụng trực tiếp dữ liệu cân bằng pha và enthalpy bay hơi để tối ưu hóa điều kiện vận hành lò phản ứng và tháp tách phân đoạn.
• Cơ quan quản lý an toàn hóa chất: Tiếp cận cơ sở dữ liệu tin cậy về đường cong áp suất hơi bão hòa để xây dựng quy chuẩn kỹ thuật về bồn bể chứa khí hóa lỏng nguy hiểm.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc phát triển mô hình hàm thế tương tác liên phân tử 5-vị trí cải tiến kết hợp hàm giảm chấn Damping kép, mở rộng thành công lý thuyết tương tác liên phân tử Morse và cơ học thống kê chất lưu của Leonhard và Deiters (2000). Luận án đã giải quyết triệt để sự bất đối xứng điện tích và phân cực phức tạp của phân tử dị hạch $\text{CO}$ và phân tử halogen $\text{Cl}_2$ mà các mô hình 1-tâm hay 2-tâm cổ điển không thể thực hiện được.

2. Điểm đổi mới phương pháp luận khi so sánh với các nghiên cứu quốc tế tiền nhiệm là gì?

So với nghiên cứu của Panagiotopoulos (1987) chỉ sử dụng thế thực nghiệm đơn giản và nghiên cứu của Leonhard & Deiters (2002) chỉ dùng cho các khí trơ/đồng hạch, luận án đã thiết lập chuỗi liên hoàn 4 cấp độ: tính toán $\text{CCSD(T)/aug-cc-pVmZ}$ khử lỗi BSSE $\rightarrow$ tối ưu toàn cục bằng giải thuật di truyền kết hợp Levenberg-Marquardt $\rightarrow$ tích phân cầu phương Gauss-Legendre 4 chiều kiểm định qua hệ số virial $B_2(T)$ $\rightarrow$ mô phỏng GEMC-NVT và học máy PCA-ANN $\text{I(5)-HL(6)-O(3)}$.

3. Phát hiện nào gây bất ngờ nhất và có bằng chứng định lượng cụ thể ra sao?

Phát hiện bất ngờ nhất là tính chính xác vượt trội của mô hình 5-vị trí thuần lý thuyết trong việc dự báo các điểm tới hạn mà không cần bất kỳ tham số hiệu chỉnh thực nghiệm nào. Cụ thể, nhiệt độ tới hạn $T_c$ của $\text{N}_2$ đạt $126.85 \text{ K}$ (so với thực nghiệm $126.19 \text{ K}$, sai số chỉ $0.52%$) và của $\text{CO}$ đạt $133.45 \text{ K}$ (so với thực nghiệm $132.86 \text{ K}$, sai số chỉ $0.44%$).

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình tái lập chi tiết tuyệt đối, bao gồm: cấu hình không gian đime, độ dài liên kết hình học cố định ($1.128 \text{ \AA}$, $1.097 \text{ \AA}$, $2.009 \text{ \AA}$), tọa độ các vị trí giả phân đoạn, mã nguồn C/C++ cho giải thuật Levenberg-Marquardt, cấu trúc mạng nơron $\text{I(5)-HL(6)-O(3)}$ và các thông số bước lặp đại tập Gibbs Monte Carlo.

5. Chương trình nghị sự nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chiến lược 10 năm tới tập trung vào: (1) Đưa hiệu ứng tương tác 3-thân vào mô phỏng pha ngưng tụ; (2) Chuyển đổi từ mô hình phân tử cứng sang mô hình phân tử dao động tự do; (3) Tự động hóa việc xây dựng mặt thế năng lượng tử bằng mạng nơron tích chập đồ thị (GNN); và (4) Mở rộng tính toán cân bằng lỏng-hơi và lỏng-lỏng cho các hệ dung môi xanh phục vụ công nghệ thu giữ $\text{CO}_2$.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thành Được đã khẳng định vị thế khoa học xuất sắc qua 5 đóng góp cốt lõi:

1. Xây dựng thành công bộ hàm thế ab initio 5-vị trí mới cho các hệ đime $\text{Ar}-\text{Ar}$, $\text{N}_2-\text{N}_2$, $\text{Cl}_2-\text{Cl}_2$ và $\text{CO}-\text{CO}$ đạt hệ số tương quan $R^2 > 0.999$.
2. Phát triển thuật toán tích phân cầu phương Gauss-Legendre 4 chiều kết hợp hiệu chỉnh lượng tử bậc nhất, giải quyết chính xác hệ số virial bậc hai $B_2(T)$ trên dải nhiệt độ rộng.
3. Hiện thực hóa thành công mô phỏng ngẫu nhiên GEMC-NVT, tái lập hoàn hảo đường cong cân bằng pha lỏng-hơi và các thông số trạng thái tới hạn với sai số $\text{MARE} < 1.5%$.
4. Ứng dụng thành công mô hình nhiệt động solvat hóa COSMO-SAC để dự báo toàn diện giản đồ pha đẳng nhiệt và đẳng áp cho hỗn hợp nhị phân nguy hiểm $\text{CO}-\text{Cl}_2$.
5. Tiên phong tích hợp mạng thần kinh nhân tạo và phân tích thành phần chính (PCA-ANN $\text{I(5)-HL(6)-O(3)}$) trong dự báo các tham số phương trình trạng thái virial với độ chuẩn xác $\text{RMSE} < 0.008$.

Công trình đã chứng minh một cách thuyết phục rằng sự kết hợp giữa hóa học lượng tử độ chính xác cao và kỹ thuật mô phỏng cơ học thống kê hiện đại có thể giải quyết trọn vẹn những bài toán nhiệt động học chất lưu phức tạp nhất, mở ra kỷ nguyên mới cho ngành hóa lý tính toán tại Việt Nam với tầm vóc hội nhập học thuật quốc tế sâu rộng.