24小时热门版块排行榜    

查看: 597  |  回复: 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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员AcTPy7 的主题更新
信息提示
请填处理意见