| 查看: 2494 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
gromacs 模拟NPT系宗出问题
|
|||
|
请教，各位，我在gromacs模拟NPT系宗的过程中，又需要限制一些原子的位置，使用freezegrps 或者 position restraint的时候，都会报错，请教各位遇到如此问题如何解决呢？ 多谢 追加问题： 我的模拟体系是中间有个膜（膜使用了position restraint），在Z方向上把溶液隔开，膜中间有空，可以使溶液通过孔来实现流动。现在要使用NPT系宗来调节密度，想要保持XY平面大小不变，通过调节z方向的长度来实现，设置如下： ; Pressure coupling = Pcoupl = berendsen Pcoupltype = semiisotropic ; Time constant (ps), compressibility (1/bar) and reference P (bar) = tau_p = 40 40 compressibility = 0 5e-5 ref_p = 0 10.0 refcoord_scaling = all NPT系宗模拟的结果是Z方向上的长度不断增大，最后模拟的盒子里水中出现了真空。这很明显不对的 按照我的理解，正确的结果应该是xy方向的盒子长度不变，而Z方向的盒子逐渐减小最后达到一个稳定的值。 请教各位，出现这样的错误，如何解决呢？ 非常感谢 ************************************************************************************************************************************ 我的问题解决了，问题出在使用 berendsen pressure coupling，并且在模拟体系里使用了 position restraint或者 freeze group，gromacs的手册里第32页有段话： “ If you allow full anisotropic deformations and use constraints you might have to scale more slowly or decrease your timestep to avoid errors from the constraint algorithms. It is important to note that although the Berendsen pressure control algorithm yields a simulation with the correct average pressure, it does not yield the exact NPT ensemble, and it is not yet clear exactly what errors this approximation may yield. ” 我把压力增大到1000bar的时候，体系正常压缩，密度达到合适终止此过程即可。 我再测试使用别的pressure coupling试试。 盒子的长度Z方向的变化与设想一致的： box length.wmf[ Last edited by jiaoyixiong on 2012-10-7 at 10:39 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : boxlength.wmf
2012-10-07 10:39:13, 47.88 K
» 收录本帖的淘帖专辑推荐
 gromacs |
» 本帖已获得的红花（最新10朵）
469187892» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
gromacs NPT 系综如何模拟？
已经有8人回复- gromacs，用pdb文件生成top文件时遇到的问题
已经有4人回复
- gromacs模拟分子在二氧化硅表面的吸附
已经有9人回复
- 如何使用gromacs进行多分子模拟？
已经有7人回复
- Gromacs处理二氧化硅为基底的体系
已经有3人回复
- 请教GROMACS安装
已经有8人回复
- Gromacs中加入离子
已经有5人回复
- 【转帖】GROMACS分析程序讨论，在gromacs中添加Amber力场
已经有8人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】gromacs出现了一个问题
已经有10人回复
- 【求助】有关gromacs得两个问题
已经有10人回复
5楼2012-10-06 19:13:26
2楼2012-10-06 13:14:14
4楼2012-10-06 18:56:24
469187892
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 909.7
- 散金: 15
- 帖子: 54
- 在线: 21.1小时
- 虫号: 1441709
- 注册: 2011-10-14
- 专业: 理论和计算化学
6楼2012-11-29 08:54:45