Li Jicun, Wang Feng
Department of Chemistry and Biochemistry, University of Arkansas , Fayetteville, Arkansas 72701, United States.
J Phys Chem B. 2016 Sep 1;120(34):9088-96. doi: 10.1021/acs.jpcb.6b06102. Epub 2016 Aug 11.
Simple nonpolarizable molecular mechanics potentials were developed for Li(+), Na(+), and K(+) with the adaptive force matching (AFM) method using the second order Møller-Plesset perturbation theory (MP2) with the frozen core approximation as reference. The effects of different choices of core orbitals and basis sets in the MP2 calculations were investigated for Na(+) and Li(+). For Na(+), correlating the 2s2p electrons in MP2 changes its hydration free energy by 18 kJ/mol, which is surprisingly large, constituting to about 5% of the intrinsic hydration free energy of the ion. Whereas correlating the 2s2p electrons with the aug-cc-pCVTZ basis set leads to the best agreement with experiments, with the aug-cc-pVTZ basis set, a better hydration free energy will be obtained if the 2s2p are kept as frozen core orbitals. Even with nonpolarizable energy expressions, the AFM derived ion potentials predict the experimental hydration free energies of the various salts within 2% of experimental values, suggesting the robustness of the fitting procedure. However, the 2% agreement can only be achieved if the core correlation is modeled appropriately in the MP2 reference calculations.
使用二阶Møller-Plesset微扰理论(MP2)并采用冻结核心近似作为参考,通过自适应力匹配(AFM)方法为Li(+)、Na(+)和K(+)开发了简单的非极化分子力学势。研究了MP2计算中不同核心轨道和基组选择对Na(+)和Li(+)的影响。对于Na(+),在MP2中关联2s2p电子会使其水合自由能变化18 kJ/mol,这一变化惊人地大,约占离子本征水合自由能的5%。虽然使用aug-cc-pCVTZ基组关联2s2p电子与实验结果的一致性最佳,但使用aug-cc-pVTZ基组时,如果将2s2p保持为冻结核心轨道,则会得到更好的水合自由能。即使使用非极化能量表达式,AFM导出的离子势预测的各种盐的实验水合自由能与实验值的偏差在2%以内,这表明拟合过程的稳健性。然而,只有在MP2参考计算中对核心关联进行适当建模,才能实现2%的一致性。