24小时热门版块排行榜    

查看: 2734  |  回复: 2

化工小巨人

新虫 (小有名气)

[交流] MS建模-分子动力学模拟参数与流程

各位MD大神,小弟刚开始接触MD模拟,最近看了很多资料,学习了MS的操作以及MD模拟的相关知识,也进行了前期的一些尝试计算,目前遇到几个问题希望大家能给予帮助。研究体系为三种有机物质A、B和C,均为液态,考察彼此之间的相互作用。操作在MS8.0上进行。
1.力场选择。所研究体系包括C、H、O、N、S元素,以及苯环、羰基、醛基、S=O,C-N、C=C。对于Compass和CVFF不知道该选哪一个?
2.流程:对分子A、B、C在Forcite中分别进行淬火,然后在Amorphous中建立盒子(建立盒子的时候有个Refinement options,是不是会在建好盒子时自动进行优化?),最后在Forcite Calculation中进行分子动力学计算。 该流程是否合适?
3.结构优化和动力学计算,怎么知道计算是否收敛或者结果准确?
4.一些计算细节:①淬火时,动力学优化系综选择是NVT,那么温度应该设为实验体系的温度吗?总模拟时间大于等于1000ps就会提示参数超上限,如果我要淬火1000ps,是不是分两次500ps即可?(第二次淬火起始速度选择current?)Thermostat选哪一个合适,个人觉得Berendsen和Andersen比较适合,选哪个?;Geometry Optimization中最大迭代次数一般设为多少合适?(文献没有提到)。 还有Optimize cell需要打钩吗? ②建盒子的时候密度怎么确定?根据实际体系估算?(混合物不太好估计)还是选择Ramp density from an initial value?   这里的温度是实验体系的温度吗?需要和单个分子淬火时的温度保持一致吗? Refinement options里面的动力学和最小化优化是对建好的盒子随即进行了优化吗?步数设为多少合适?5000合理吗? ③最后动力学计算静电力用Ewald,范德华用cutoff,这样可以吧?
5.关于静电力和范德华,在单个分子淬火优化和盒子建立时只有Atombased和Groupbased,那么选Atombased可以吧?(选项中没有Ewald) 那么静电力和范德华的cutoff需要分别设置吗?如果设置的不一样,在计算盒子边长的时候以范德华的cutoff为准可以吗?
问题有点多,希望大神可以帮帮忙,最近看了很多资料,这些细节问题还是不太懂!谢谢!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化工小巨人

新虫 (小有名气)

MD高手呢!大哥们帮忙解答问题呀
2楼2017-03-22 15:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

送你的太阳

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我也有相似的问题,请问楼主的问题解决了吗?

发自小木虫Android客户端
3楼2017-10-20 15:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 化工小巨人 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +15 低垂的野花 2026-07-31 20/1000 2026-08-01 14:10 by gltch
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 18/900 2026-08-01 13:56 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见