| 查看: 646 | 回复: 0 | ||
qq364085265金虫 (著名写手)
|
[求助]
高密度聚乙烯动力学模拟遇到问题，求大神！
|
|
小弟正在用MS做高密度聚乙烯在150度高温下降解的模拟，需要观察它什么时候断键，先断那根键等等。可是算到现在动力学模拟都跑了快100步了还没出现断键，后来又做一次300度的，结构跑到20多步发现一个碳氢键断了，但是继续跑两步后又连上了，很无语，。。，。，观测了一下键长，感觉在高温环境下键长变化不大，20%的幅度都很难达到。以下是我建模的过程，望高手指点。 1、建立乙烯重复单元并构建聚合度为10 的均聚物，扭转角设为random； 2、优化分子链（在Discover模块中：非键参数选择Vdw&Coulomb，计算方式选择Group based，并选用compass力场计算，最后用Discover的minimizer优化能量） 3、用AC模块建立无定形晶包（放三条链进去） 4、按步骤2的方法再优化一次晶包 5、用Dmol3跑分子动力学模拟，交换关联函数选择GGA的BLYP，选择NVT系宗，温度调为423K，跑500步，步长为1fs 6、键计算选择键长公差为0.6-1.12    |
回复此楼
