基于粗粒化模型对有机溶剂的分子动力学模拟

在基于boltzmann分布对四种基本构象进行montecarlo取样后,通过与全原子模型的范德华势比较得到了gay-berne(gb)参数.又在对用量化计算得到的分子体系的电势进行电荷、偶极矩和四极矩的拟合后,得到了电多极展开势(emp)参数.利用得到的粗粒化参数,基于粗粒化模型,对chcl3及四氢呋喃(thf)两种有机溶剂进行了分子动力学模拟(mds),并将结果同全原子模拟进行了比较.计算结果表明用粗粒化模型从整体上能重复全原子模型的模拟结果,但在某些细节的计算与全原子模型有偏差,其原因可能是目前工作仅考虑了单位点情况,为此今后在对具有复杂结构的分子进行粗粒化模拟时还应考虑合理放置及增加相互作用位点.