24小时热门版块排行榜    

查看: 1430  |  回复: 1

zxy964777869

铜虫 (正式写手)

[求助] NVE系综的初始构型

我的一个体系,跑完NVT达到温度平均值300K,跑完NPT,压力平均值为1bar,RMSD值为800多。我用NPT模拟产生的gro文件,作为NVE模拟的初始文件,跑完之后发现在NVE的整个模拟过程中,压力呈一条直线,至为-240bar.
这是NVE的输出信息
Energy                      Average   Err.Est.       RMSD  Tot-Drift
-------------------------------------------------------------------------------
Temperature                 300.006    0.00092    5.29188  0.0049759  (K)

Energy                      Average   Err.Est.       RMSD  Tot-Drift
-------------------------------------------------------------------------------
Pres. DC                   -214.126          0          0 4.45431e-11  (bar)
我怀疑是NPT模拟中压力是上下震动的(从RMSD看出),虽然压力平均值为1bar,但它产生的npt.gro文件实际上不是压力为1bar时的构型,体积不是压力为1bar时的应有体积,所以导致NVE模拟中压力不是1bar。
不知各位高手是否遇到过这种情况,可否指点一二,金币不多,聊表心意,谢谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sheva7love

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
chaizhm: 金币+2, 谢谢~ 2012-07-13 18:04:37
NPT压力在平衡值附近上下波动是正常的,其平衡后的构型可以反映该压力下体系的密度和体积了。
我觉得不是初始构型的问题,应该是NVE模拟中的问题。
修身齐家治国平天下
2楼2012-07-13 10:18:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxy964777869 的主题更新
信息提示
请填处理意见