| 查看: 1116 | 回复: 5 | ||
jackwangee金虫 (小有名气)
|
[求助]
gromacs NPT restrain 两个分子已有1人参与
|
|
做石墨烯水体系NPT 想restrain 两层石墨烯 用如下方法restrain: 1. genrestr -f tcn.gro -o strong_posre.itp -fc 418 418 418 2. ; Strong position restraints for MOL #ifdef STRONG_POSRES #include "strong_posre.itp" #endif in itp 3.define = -DSTRONG_POSRES in mdp 但是strong_posre.itp 内的原子个数必须和力场参数itp内原子个数相同。 做出来之后 。。只能限制其中一层. 同时我也改变力场。。将两层石墨烯的原子的全部变成一个itp文件。但是结果是 石墨烯跑散了。 请问还有其他办法可以限制两层石墨烯吗？ |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
gromacs 模拟NPT系宗出问题
已经有8人回复- gromacs NPT 系综如何模拟？
已经有8人回复
jackwangee
金虫 (小有名气)
| up |
2楼2014-03-20 09:55:45
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+5, 感谢指导！ 2014-03-21 16:37:42
jackwangee: 金币+10, ★★★★★最佳答案 2014-04-07 12:12:49
感谢参与，应助指数 +1
月只蓝: 金币+5, 感谢指导！ 2014-03-21 16:37:42
jackwangee: 金币+10, ★★★★★最佳答案 2014-04-07 12:12:49
|
帅哥，你的思路没错，就是你需要再仔细读一下手册，关于 position_restraints 的设置问题。 不管你有多少层石墨烯，你把石墨烯的参数都弄在一个itp文件中，比如取名就叫 graphene.itp 打个比方说你有100个原子，那么在 graphene.itp 文件里应该是这样的： [moleculetype] ; molname nrexcl graphene 1 [atoms] ; id type resnr residu atom cgnr charge mass 1 CA 1 GRA C4 1 0.0 12.011 2 CA 1 GRA C2 2 0.0 12.011 3 CA 1 GRA C4 3 0.0 12.011 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 96 CA 1 GRA C2 96 0.0 12.011 97 CA 1 GRA C4 97 0.0 12.011 98 CA 1 GRA C2 98 0.0 12.011 99 CA 1 GRA C4 99 0.0 12.011 100 CA 1 GRA C2 100 0.0 12.011 [ bonds ] ; i j funct bo kb [angles] ; ai aj ak funct [dihedrals] ; ai aj ak al funct ; Include Position restraint file #ifdef POSRES #include "graphene_pr.itp" #endif ；OK graphene.itp 到此结束，最后引用的graphene_pr.itp 文件里面就是你想限制位置的那些原子序列，这参数设置就放在这里，如果你的模拟过程中不需要限制，就在MDP文件中设置即可，如果需要限制，就在MDP文件中加上一条 define = -DPOSRES 就可以了。 打个比方我想限制第1个和第100个原子的位置，那么就在graphene_pr.itp文件中写： [ position_restraints ] ; atom type fx fy fz 1 1 1000000 1000000 1000000 100 1 1000000 1000000 1000000 ok，如果你还想限制其他的原子，就再加上即可。 祝科研顺利 |
3楼2014-03-20 10:37:02
4楼2014-03-20 10:41:10
jackwangee
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1168
- 散金: 3
- 红花: 2
- 帖子: 124
- 在线: 185.6小时
- 虫号: 517335
- 注册: 2008-03-03
- 性别: GG
- 专业: 理论和计算化学
5楼2014-03-21 09:21:50
6楼2014-03-21 10:09:40