[论文解读] To Pair or Not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water

Niamh O’Neill, Benjamin X. Shi+3 more|arXiv (Cornell University)|Nov 2, 2023

Spectroscopy and Quantum Chemical Studies被引用 2

一句话总结

本研究利用RPA和MP2方法构建基于波函数的机器学习势能，用于模拟水溶液中的NaCl，精确捕捉了溶剂化离子间的平均势能（PMF）。RPA和MP2对离子对态与溶剂共存态的能量预测几乎简并（能量差小于0.2 kcal/mol），解决了电解质模拟中长期存在的分歧。

ABSTRACT

The extent of ion pairing in solution is an important phenomenon to rationalise transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wavefunction based machine learning potentials with the Random Phase Approximation (RPA) and second order Moller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.

研究动机与目标

• 为解决计算模型在水溶液NaCl离子对行为方面缺乏共识的问题。
• 开发能够准确捕捉电子相关效应（对离子-溶剂相互作用至关重要）的波函数基机器学习势。
• 使用高阶电子结构方法基准测试水溶液中NaCl的平均势能（PMF）。
• 比较RPA、MP2与多种DFT泛函在预测溶剂化结构与离子对结合能方面的性能。
• 通过分析离子-水与离子-离子相互作用能，提供可靠的离子对结合能参考数据。

提出的方法

• 基于RPA和MP2电子结构计算训练NaCl在水中的机器学习势。
• 采用包含95个水分子的模型箱，通过周期性边界条件模拟溶剂化Na⁺与Cl⁻离子。
• 通过伞形采样沿离子间距坐标计算平均势能（PMF）。
• 使用参考电子结构方法在真空中计算离子-离子结合能曲线，以评估内在相互作用强度。
• 通过能量分解量化离子-水相互作用能：EIW = Esys − EW − EG。
• 通过分析水-水、离子-水与水-氢的径向分布函数（RDFs），探究溶剂化结构。

实验结果

研究问题

• RQ1 RPA与MP2方法在预测水溶液中NaCl离子对与溶剂共存态相对稳定性方面表现如何？
• RQ2 使用高阶波函数方法与DFT泛函计算的Na⁺与Cl⁻离子间PMF在定量上存在多大差异？
• RQ3 与RPA和MP2相比，不同DFT泛函（如revPBE-D3、r2SCAN、optB88-vdW）在描述离子-水与离子-离子相互作用方面表现如何？
• RQ4 核量子效应在多大程度上影响NaCl/水体系中的PMF与溶剂化结构？
• RQ5 通过离子-水与离子-离子相互作用能的分解，能获得关于离子对本质的哪些新见解？

主要发现

• RPA与MP2势能预测离子对态（CIP）与溶剂共存态（SSIP）能量几乎简并，能量差仅为0.2 kcal/mol。
• CIP势阱深度为RPA的27.9 kcal/mol与MP2的27.8 kcal/mol，表明存在浅但稳定的离子对势阱。
• Na⁺的离子-水相互作用能约为-173.8 kcal/mol（RPA）与-172.3 kcal/mol（MP2），而Cl⁻的相互作用为弱稳定化，约+4.8 kcal/mol。
• DFT泛函如revPBE-D3因离域化误差导致对Na⁺与Cl⁻的过度结合，Cl⁻-水相互作用能高达+12.68 kcal/mol。
• RPA与MP2模型中水的自扩散系数与实验值高度一致，优于经典力场与DeepMD势。
• RDF分析表明，RPA与MP2均能准确再现真实的溶剂化壳层，其第一壳层配位数与实验预期一致。