School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
School of Marine Science and Environment, Dalian Ocean University, Dalian 116023, People's Republic of China.
J Phys Chem A. 2020 Jul 16;124(28):5963-5978. doi: 10.1021/acs.jpca.0c03941. Epub 2020 Jul 2.
Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH-water system were proposed. The difference between these two models is whether there is charge transfer between OH and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH(HO) ( = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH. Molecular dynamics simulations on the aqueous OH solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.
基于原子键电负性均衡方法与分子力学(ABEEM/MM)融合,提出了两种 OH-水体系的变电荷模型。这两种模型的区别在于 OH 和其第一层水分子之间是否存在电荷转移。利用这两种 ABEEM/MM 模型、OPLS/AA 力场、OPLS-SMOOTH/AA 力场和 QM 方法研究了 OH(HO)(=1-8、10、15、23)团簇的结构、电荷分布、电荷转移和结合能。结果表明,两种 ABEEM/MM 模型都可以像 QM 方法一样搜索到所有稳定结构,并且结构和电荷分布与 QM 计算结果吻合良好。第一水合壳层中的结构、电荷转移和氢键强度与 OH 的配位数密切相关。利用 ABEEM/MM-I 模型在 298 和 278 K 下对水溶液中的 OH 进行分子动力学模拟。MD 结果表明,在 298 K 时,三、四和五配位 OH 的比例分别为 29.6%、67.1%和 3.4%,而在 278 K 时,二、三、四和五配位 OH 的比例分别为 10.8%、44.9%、39.2%和 4.9%;第一层的平均氢键长度和氢键角随温度降低而增加。