[论文解读] Pressure Correction in Classical Density Functional Theory: Hyper Netted Chain and Hard Sphere Bridge Functionals

Volodymyr P. Sergiievskyi, Guillaume Jeanmairet+2 more|arXiv (Cornell University)|Sep 4, 2015

Phase Equilibria and Thermodynamics参考文献 27被引用 3

一句话总结

本文提出了一种经典密度泛函理论（DFT）的后处理压力校正方法，以改进溶剂化自由能（SFE）的预测。结果表明，SFE计算中的误差源于HNC和硬球桥接泛函中系统压力不正确。该校正方法基于可压缩性路线压力与额外的ρkTΔV项，显著提升了预测精度，尤其在结合HNC/PC或HS-Bridge泛函时，对122种有机分子的SFE误差已达到与分子动力学模拟相当的水平。

ABSTRACT

Low accuracy of the Solvation Free Energy (SFE) calculation is a known problem of the numerical methods of the Integral Equation Theory of Liquids and the Classical Density Functional Theory (Classical DFT). Although functionals with empirical corrections can essentially improve the predictability of the methods, their universality is still a question. In our recent paper we connected the SFE calculation errors with the incorrect pressure in the Classical DFT and proposed the a posteriory correction to improve the results (J. Phys. Chem. Lett., 5, 1925-1942 ). This paper raised a discussion in the community. In particular, recently appeared a critical reply where pointed some thermodynamical inconsistencies of the derivations in our paper (J. Chem. Theory Comput., 11, 378-380). In the present work we re-derive the pressure correction in a more simple way and show that despite the inaccuracies during the derivation, the final form of the previously derived correction is correct. We also test the applicability of the proposed correction to the functionals which include a three- and many- body terms from the fundamental measure theory (FMT) for hard sphere fluid. We test all the functionals on a set of model systems and discuss the obtained results.

研究动机与目标

• 解决经典DFT中溶剂化自由能（SFE）预测长期存在的不准确问题，尤其针对HNC和硬球桥接泛函。
• 探究SFE误差的根本原因，识别出系统压力不正确为主要驱动因素，尤其在HNC中压力高达11.5 kBar。
• 开发并验证一种热力学一致的后处理压力校正方法，以提升SFE预测精度，且不改变底层泛函形式。
• 评估在参数化硬球直径时，热力学一致性（正确压力）与小溶质SFE预测精度之间的权衡。
• 提出一种通用的ρkTΔV校正项，以进一步降低SFE误差，尤其适用于具有显著偏摩尔体积的较大溶质。

提出的方法

• 采用更简单、更清晰的方法重新推导HNC近似下的压力校正，以验证先前压力校正公式的正确性。
• 将压力校正应用于HNC和硬球桥接（HS-Bridge）泛函，利用可压缩性路线计算系统压力。
• 对HS-Bridge泛函进行参数化，使其既能重现正确压力，也能重现甲烷的SFE，以探索热力学一致性与预测精度之间的权衡。
• 基于SFE误差与偏摩尔体积依赖关系的系统分析，引入额外的ρkTΔV校正项。
• 在包含122种有机分子的基准测试集中，测试校正后的泛函（HNC/PC+，HS/PRES+），评估其相对于分子动力学模拟的预测精度。
• 通过径向分布函数与溶剂结构分析，评估溶剂化壳层预测质量，超越SFE本身。

实验结果

研究问题

• RQ1 为何基于HNC的泛函尽管能正确预测溶剂结构的定性特征，却仍产生高度不准确的溶剂化自由能？
• RQ2 HNC和HS-Bridge泛函中不正确的系统压力在多大程度上导致SFE误差，尤其对具有大偏摩尔体积的溶质？
• RQ3 能否推导出一种通用的、后处理的压力校正方法，在不损害热力学一致性的情况下提升SFE预测精度？
• RQ4 硬球直径的选择如何影响系统压力与SFE预测？能否通过单一参数化同时实现正确压力与小溶质SFE的准确预测？
• RQ5 SFE预测中的系统性误差是否可通过添加ρkTΔV项得到纠正？该修正是否适用于多种类型溶质？

主要发现

• HNC近似导致系统压力高达11.5 kBar，这是溶剂化自由能被高估的主要原因，因为−PΔV项主导了误差。
• 后处理压力校正（HNC/PC）成功恢复了SFE对溶质尺寸的正确定性依赖关系，尤其在大硬球溶质中，HNC错误地预测了与体积相关的趋势。
• 将HS-Bridge泛函参数化为重现正确压力，可实现对大溶质的准确SFE与溶剂结构预测，但导致对小溶质（如甲烷）的SFE显著低估。
• 结合HS-Bridge与压力校正的混合泛函（HS/PRES+）未能改善小溶质的SFE预测，表明当前形式下压力一致性与小溶质精度之间存在根本性不兼容。
• 额外的ρkTΔV校正项能有效降低各类溶质的SFE误差，尤其对具有大偏摩尔体积的溶质效果显著；当应用于HNC/PC+时，误差降低至与分子动力学模拟相当的水平。
• 对于122种有机分子，采用ρkTΔV校正的HNC/PC+泛函实现了最佳的整体SFE预测精度，表明该修正是一种稳健且可推广的改进方法，适用于经典DFT的溶剂化自由能计算。