SPC、SPC／E和TIP3P模型水的结构和性质.pdfVIP

· 92 · SPC、SPC／E和TIP3P模型水的结构和性质冰 谈荣 日。邵建群 (1．贵州大学理学院 贵州省光电子技术与应用重点实验室，贵州 贵阳55(X)25；2．南昌市第十中学，江西 南昌330000) 摘要：本文采用分子动力学方法分别计算了SPC、SPC／E和TIP3P~-gj密度、自扩散系数等热力学性质。通过对径向分布 函数的分析发现，三种模型水的g(r)没有明显的差异，且温度越低，径向分布函数的峰值越大，说明水分子的结构越明显， 分子间的氢键越强，分子极性也越大。 关键词：水模型 热力学性质 分子动力学 中图分类号：0561．1，0561．2 文献标识码：A 文章编号：1002—6886(2011)04-0092—03 StructuresandPropertiesofSPC，SPC／E，TIP3P Water TAN Rongri，SHAO Jianqun Abstract：ThethermodynamicpropertiesofSPC，SPC／EandTIP3Pwater，includingdensityandserf—diffusioncoefficient，are studiedusingmoleculardynamicssimulations．TheradialdistributionfunctionscalculatedindicatethattheRDFsofthethreewatersrae almostidentica1．Inaddition，thefirstpeakoftheRDFsbecomeshigherandhigherwiththetemperaturedecreasing．Itsuggeststhatthe structuresofwatermoleculesraemorepronouncedandtheactionofthehydrogenbondgetsstrongerandstrongerwhichcausesthepo— larityofwatermoleculesincreasing． Keywords：watermodel；thermodynamicproperties；moleculradynamicssimulations 了两个只带电量的粒子。不同的势模型 下计算得到 0 引言 的热力学性质经常不相同，适用于不同的体系。在模拟 x射线散射…、中子衍射 以及核磁共振 等实验方 中，可以根据实际需要选择合适的水模型。 法曾被用来研究过液态水的结构。近些年来，随着计算机 本文主要阐述用分子动力学方法计算SPC、SPC／E和 的飞速发展，分子动力学等理论模拟方法已经成为一种重 TIP3P模型水的密度、自扩散系数等热力学性质随温度的 要而又有效的研究微观结构和热力学性质的手段。在液 变化关系，并比较它们与实验值的差异。同时还通过径向 态水的分子动力学模拟中的一个关键问题是如何选择一 分布函数和配位数分析其微观结构。 个合适的势模型。势模型既可以从量子力学从头计算 ( 1 模型和MD计算细节 Abinitio)获得，如MCY模型，也可以由半经验方法获得， 如 SPC(simplepointchrage)H、SPC／E(extendedsimple 在MD模拟中，电荷决定分子的极性，短程相互作用 pointcharge) 、TIP3P(transferableintermoleculrapotential 由Lennrad—Jones(简写为 LI)势决定。其表达式可写 3points) 、TIP4P 以及TIPSP 等水模型。最简单的水 为 分子模型就是键长、键角等参数都固定的刚性模型，如 [()一()。卜 (1) SPC和T1P3P模型。这两种模型的原子质量和电量都在 同一个质点上，