24小时热门版块排行榜    

查看: 408  |  回复: 2

jiaoyixiong

荣誉版主 (职业作家)

[求助] gromacs 模拟的一个条件判断

请教大家一个技术性问题:
我现在用gromacs想模拟一个体系,采用NPT系宗,semiisotropic 加压的方法给模拟体系施加一个侧向压力,当模拟的盒子在xy平面的边长大于20nm的时候,施加侧正向压力,使盒子压缩,当压缩到20nm的时候,改为侧向负压力,使盒子扩张当盒子扩张到50nm的时候,再压缩。如此循环模拟,我要做5个这样的循环,如何在gromacs中设置呢?
非常感谢各位指教
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jiaoyixiong(金币+10): 有帮助 非常感谢,我是需要在体系中施加压力的作用而产生变化。 2012-03-02 10:38:12
回答不了在gromacs里面怎么实现你的想法, 但是看你描述, 你想让体系在20nm和50nm之间振荡. 如果不是很关注在指定的压力下体系的变化(猜想这是你用NPT的原因?), 是不是可以用NVT, 将几步模拟串起来跑? 比如

DO I  = 1, 5
(1) NVT, x=y=20nm, run for 1ns for equilibrium, followed by  another 10 ns for statistics;

(2) NVT, x=y=50nm, run for 1ns for equilibrium, followed by another 10 ns for statistics;

END DO

我在lammps就常采用类似的做法.
Life, Love, Laugh.
2楼2012-03-02 09:11:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

qphll

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jiaoyixiong(金币+10): 有帮助 thanks a lot 2012-03-02 13:58:43
嗯, 如果你一定要知道在不同压力下盒子的响应(体积变化), 那么可以设定一系列的压力值, 将NPT串起来跑.

更进一步, 如果你只是希望看到盒子在x, y方向的变化, 那么可以在模拟过程中, 将z方向fix.

只是想法而已, 供你参考一下, 并不一定是最优的方案. 但是原则上我觉得你的想法是可以实现的.
Life, Love, Laugh.
3楼2012-03-02 10:46:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiaoyixiong 的主题更新
信息提示
请填处理意见