Fragment Quantum Mechanical Method for Large-Sized Ion-Water Clusters

Jinfeng Liu; Lian-Wen Qi; John Z H Zhang; Xiao He

doi:10.1021/acs.jctc.7b00149

Fragment Quantum Mechanical Method for Large-Sized Ion-Water Clusters

J Chem Theory Comput. 2017 May 9;13(5):2021-2034. doi: 10.1021/acs.jctc.7b00149. Epub 2017 Apr 13.

Authors

Jinfeng Liu¹, Lian-Wen Qi¹, John Z H Zhang^{2

3

4}, Xiao He^{2

3}

Affiliations

¹ State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University , Nanjing, 210009, China.
² School of Chemistry and Molecular Engineering, East China Normal University , Shanghai, 200062, China.
³ NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai, 200062, China.
⁴ Department of Chemistry, New York University , New York, New York 10003, United States.

PMID: 28379695
DOI: 10.1021/acs.jctc.7b00149

Abstract

Fragmentation methods have been widely studied for computing quantum mechanical (QM) energy of medium-sized water clusters, but less attention has been paid to large-sized ion-water clusters, in which many-body QM interaction is more significant, because of the charge-transfer effect between ions and water molecules. In this study, we utilized electrostatically embedded generalized molecular fractionation (EE-GMF) method for full QM calculation of the large-sized ion-water clusters (up to 15 Na⁺ and 15 Cl^- ions solvated with 119 water molecules). Through systematic validation using different fragment sizes, we show that, by using distance thresholds of 6 Å for both the two-body and three-body QM interactions, the EE-GMF method is capable of providing accurate ground-state energies of large-sized ion-water clusters at different ab initio levels (including HF, B3LYP, M06-2X, and MP2) with significantly reduced computational cost. The deviations of EE-GMF from full system calculations are within a few kcal/mol. The result clearly shows that the calculated energies of the ion-water clusters using EE-GMF are close to converge after the distance thresholds are larger than 6 Å for both the two-body and three-body QM interactions. This study underscores the importance of the three-body interactions in ion-water clusters. The EE-GMF method can also accurately reproduce the relative energy profiles of the ion-water clusters.