论文部分内容阅读
这篇论文介绍了一个一般的基于电多极展开势和Gay-Berne (GB)势的粗粒化分子力学模型.在这个粗粒化模型中,把有机溶液分子处理成为单轴的椭球体,用各向异性的GB函数表示单轴椭球体之间的相互作用,从而实现了范德瓦尔斯势的粗粒化.分子体系的电荷分布用电多极展开来表示,包括位于质心处的点电荷、偶极矩和四极矩.在基于Boltzmann分布对四种基本的分子构象进行了Monte Carlo取样之后,通过与全原子模型的范德瓦尔斯势比较得到了GB参数.又在对用量子化学方法计算得到的分子体系的静电势进行了电荷、偶极矩和四极矩的拟合之后,得到了电多极展开势参数.利用得到的粗粒化参数,基于此粗粒化模型,本文对CHCl3、CH2Cl2、C5H5N、DMF、DMSO、THF六种有机溶剂分子进行了分子动力学模拟,并将其结果同全原子模型的模拟结果进行了比较.计算得到的结果表明了用此粗粒化模型从整体上能够重复全原子模型的模拟结果,但是在某些细节的计算上与全原子模型的计算结果有些偏差,其原因可能是在我们目前的工作中仅仅考虑了单个位点的情况.为此,今后在对具有复杂结构的分子体系进行粗粒化模拟时还应该考虑合理地放置相互作用的位点,以及考虑合理地增加相应的相互作用位点的数目.通过粗粒化模型进行分子动力学模拟的主要优点是可以在保证能获得一定精度的条件下较大幅度地提高计算的效率,这主要是由于用粗粒化的GB-EMP相互作用代替了多对原子之间的相互作用,从而使自由度的数目得到了较大地减少.此外另外一个重要的因素就是由于忽略了键长、键角和二面角的高频运动,而使模拟的时间尺度增大了.还需要指出的是,由于在处理GB效应中引入的d w参数较好地控制了GB势的柔性,所以用在本文中介绍的粗粒化模型还可以被用来模拟更大的更复杂得分子体系.此外在考虑了能够更精确地描述静电相互作用势的电多极势效应之后,使基于此粗粒化模型模拟那些带高电量的生物大分子体系成为可能.