<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
<BR>Hi, All,<BR>
I am simulating an interface system composed of methane gas and liquid water in the NPT ensemble. The temperature is well controlled by using either the Berendsen or the Nose-Hoover T-couplings. As for pressure, the Berendsen P-coupling also works well. However, the Parrinello-Rahman P-coupling always produces a higher pressure than the expected. For example, I set ref_p=300 bar in the mdp file but obtain an average pressure as ~330 bar (over 0.5 ns). I test the tau_p in the range 0.1 ~ 10 ps, the case does not change. When I test the single-phase system, i.e., pure water, the Parrinello-Rahman P-coupling can indeed produce the expected pressure with the deviation less than 2 bar. What is the matter in the two-phase system?<BR>
<BR>
Thanks a lot.<BR>
Guang-Jun<BR><br /><hr />立刻下载 MSN 保护盾，保障 MSN 安全稳定！ <a href='http://im.live.cn/safe/' target='_new'>现在就下载！</a></body>
</html>