| 查看: 2831 | 回复: 3 | ||
xiaomeichen金虫 (小有名气)
|
[求助]
分享gromac的基本步骤，有木有好心人告知MC模拟 用towhee基本步骤
|
|
偶然看到这个，以下是进行分子动力学模拟的最基本步骤，来自Gromacs维基。具体的步骤可能会因为模拟目标的不同而存在差异，这里只做基本的概述。 1. 十分清楚要什么样的模拟系统性质和现象，明白模拟的目标。 2. 选择适当的工具以实现模拟目标。这需要熟知其他研究人员对相似模拟系统的模拟方法等资料，要好好读别人的论文，主要包括： -模拟使用的软件，当然还要考虑使用的力场，有时因为力场关系，需要更改使用的软件。 -力场描述了系统中原子的性质和相互作用，选择一个适合自己研究目标的力场至关重要，还是需要多看论文。根据sensenbobo血泪史和别人的经验，模拟核酸的话，可以采用amber力场，水分子可以使用TIP3P；模拟蛋白质的话，可以采用OPLS-AA，水分子使用TIP4P（重申一下，我本人没有跑过测试）。 3. 获取模拟分子的坐标文件（PDB），或者自己生成一个。自己生成很好啊，强烈鼓励，但是要十分明白有一定的危险性，要明白每一步模拟系统的现象是为什么。 4. 根据模拟软件，生成模拟系统初步坐标文件。 5. 获取或者生成系统的拓扑文件，比如Gromacs的pdb2gmx命令，或者使用网络服务器PRODRG（google一下这几个字母），amber的leap命令和antechamber，NAMD的psfgen或者VMD。请参考相关手册。 6. 给模拟系统添加一个盒子（盒子的翻译方法我十分不喜欢，箱子显得太大，边界显得太难理解，但是盒子听起来想到蛋糕，我们是严肃的科学家。娃哈哈。）在盒子中添加水分子，添加离子如NA和CL之类（添油加醋，生活就是这样进行的）。添加任何东西到你的系统中，就要相应的修改一下拓扑文件，以保持与系统对应，不然会出错哦。 7. 进行能量最优话。这样做是消去系统中原子之间距离太近造成的高能量。必要的时候，分子可以现在真空中能量最优话，然后再做溶剂环境中能量最优话。 8. 选择恰当的模拟参数模拟平衡模拟系统。需要慎重考虑力场的作用，必要的时候保持系统体积温度不变（NVT系综）进行位置限制性平衡系统，然后在保持系统温度压强不变（NPT系综）修正系统的密度。然后再选择正确的系综（NPT，NVT，NET，BT？？）进行数据收集。（真烦！） 9. 跑足够长时间的模拟，使系统达到理想的平衡。使用模拟软件的工具或者其他软件观察系统的性质，如温度，能量瞪瞪。 10. 使用适当的参数进行数据采集模拟，注意模拟衔接过程中各个系统性质保持不变，如果原子速度之类。这里还要考虑力场的作用，同时要考虑怎么测量描述模拟目标。 11. 模拟足够长的时间使模拟目标明显出现（可能也出现跟你期望相反的结果，出现率应该百分之五十，加上个人感情因素，可能会达到百分之八十。保重！）。 12. 分析模拟结果，解释得到结果。要好好利用视图工具，同时要十分注意漂亮的三维图说明不了什么问题。图表，比较，统计，一个都不能少。 就这样吧。  分割线有没有好心人高人 高尚的厉害的好心人帮我说一说蒙特卡洛模拟 MC模拟 最好是用towhee软件模拟的基本步骤，参照上面的那样，我觉得他说的很清楚 谢谢谢谢~~~ |
回复此楼» 本帖已获得的红花（最新10朵）
aken5555» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助懂医学的好心人
已经有11人回复
- 大家好，有好心人帮我计算四个分子的homo和lumo轨道能级图吗？急需，可以支付适当费用
已经有4人回复
- MATERIAL STUDIO画高分子晶体（聚乙烯PE）
已经有10人回复
- 如何用拟薄水铝石制备铝溶胶？？求具体步骤~~
已经有4人回复
- udf 重力方向变化
已经有17人回复
- 使用towhee进行GEMC模拟相平衡的MC运动定义问题..
已经有15人回复
- 【求助】英语太难了，哪位好心人能帮我翻译下下面关于旗袍的这段中文，翻译成英文，谢
已经有5人回复
- 【求助】哪位好心人帮我查查该物质的单晶结构有发表没？十分感谢
已经有7人回复
cuiyoutian88
金虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 890.1
- 散金: 27
- 帖子: 262
- 在线: 198.5小时
- 虫号: 1193813
- 注册: 2011-01-21
- 性别: GG
- 专业: 生物化学
2楼2012-02-05 11:10:55
五绝神鹰
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 240.9
- 散金: 5
- 帖子: 54
- 在线: 22.2小时
- 虫号: 3913553
- 注册: 2015-06-07
- 性别: GG
- 专业: 生物大分子结构与功能
3楼2016-12-07 17:26:55
aken5555
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3034.1
- 散金: 200
- 红花: 1
- 帖子: 387
- 在线: 57.4小时
- 虫号: 345700
- 注册: 2007-04-14
- 性别: GG
- 专业: 资源化工
4楼2020-02-07 16:32:20
