| 查看: 669 | 回复: 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 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
» 猜你喜欢
关于filecode,很负责任的告诉大家
已经有5人回复
关于豆爷回答的JTJC与%2F数量
已经有7人回复
2026国自然放榜时间
已经有3人回复
国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率
已经有7人回复
好奇怪的filecode
已经有4人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有5人回复













回复此楼