| 查看: 1456 | 回复: 1 | ||
[求助]
利用Gromacs在做纯溶剂的md.mdp时，有一个关于温度耦合的设定没有明白，求指点！ 已有1人参与
|
|
我利用Gromacs做离子液体的动力学模拟，在写md.mdp时，完全搬照网上的例子，这样写： title = try cpp = /usr/bin/cpp constraints = none integrator = md dt = 0.002 ; ps ! nsteps = 50000 ; total 100ps nstcomm = 1 nstxout = 500 ; collect data every 1ps nstvout = 0 nstfout = 0 nstlist = 10 ns_type = grid rlist = 1.0 coulombtype = PME rcoulomb = 1.0 vdwtype = cut-off rvdw = 1.4 fourierspacing = 0.12 fourier_nx = 0 fourier_ny = 0 fourier_nz = 0 pme_order = 4 ewald_rtol = 1e-5 optimize_fft = yes ;Berendsen temperature coupling is on Tcoupl = v-rescale tau_t = 0.1 0.1 tc-grps = protein non-protein ref_t = 300 300 ; Pressure coupling is on Pcoupl = parrinello-rahman Pcoupltype = isotropic 我现在有两个问题不明白： （1） 网上的教程中有这样一句话：“用于explicit solvation 的md.mdp文件内容（特殊注释：做真空模拟，去掉温度耦合中的 “sol”部分。在有补偿离子的模拟中，为离子加入相应的温度耦合参数）” 在我这里应该去掉哪部分？ （2） tc-grps = protein non-protein 表示的具体意思是什么，我看manul没有很懂，另外我模拟的体系中没有protein，所以是不是不用写这句？ 谢谢高人指点！！ |
回复此楼» 猜你喜欢
8月时间戳变的，举个手。玩一下，释放压力
已经有4人回复
-
UV压敏胶开发
已经有6人回复
-
有没有H口的？有收到消息的吗？
已经有3人回复
-
十年后又回来了，论文投稿求助
已经有3人回复
-
纯娱乐，不喜欢勿喷
已经有10人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有26人回复
-
面上提前没消息，有中的吗
已经有18人回复
-
娱乐
已经有4人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有16人回复
-
一个有机合成实验室都需要哪些设备？
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gromacs模拟温度耦合和压力耦合选择算法的问题，感谢！！！
已经有4人回复
- 怎么用gromacs模拟混合溶剂？跪求大牛指导。。。
已经有15人回复
- 在gromacs中运行mdrun命令时出错
已经有10人回复
fangsteel
木虫 (正式写手)
- 应助: 96 (初中生)
- 金币: 1215.1
- 散金: 32
- 红花: 22
- 帖子: 456
- 在线: 192.5小时
- 虫号: 2248690
- 注册: 2013-01-18
- 性别: GG
- 专业: 原子和分子物理
2楼2014-08-22 10:14:11