| 查看: 689 | 回复: 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    |
回复此楼» 猜你喜欢
面上意见出来了
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有16人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
有没有仍没收到信息的
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
10