<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Dear Mark,<BR>
Thank you very much for the quick reply. According to your suggesstion, I had a look at the structure or trajectory using ngmx. Only to find that there are some big hole in the center, which seems to indicate the system is very inhomogeneous. And the system volume decreased slowly. when I using cutoff instead of PME for columbic force, and change the cutoff from 1.4 to 1.0, the mdrun finished 300 ps successfully. Another 3ns is started to try whether it can succeed too. more help from these descriptions?<BR>
Best regards,<BR>
xiaowu<BR>
<BR>
Message: 2<BR>Date: Sun, 26 Apr 2009 21:41:30 +1000<BR>From: Mark Abraham <Mark.Abraham@anu.edu.au><BR>Subject: Re: [gmx-users] Re:The X-size of the box times triclinic skew<BR>factor<BR>To: Discussion list for GROMACS users <gmx-users@gromacs.org><BR>Message-ID: <49F4486A.3000305@anu.edu.au><BR>Content-Type: text/plain; charset=GB2312<BR> <BR>wuxiao wrote:<BR>> Dear Mark,<BR>> Thank you for your reply. The system I simulated was one Polymer <BR>> chain solvated by 200 water that had totally about 3000 atoms. I had <BR>> defined the system density as 1000 using editconf.<BR> <BR>That just does a scaling. It might not fix a gross problem, so do look<BR>at the box dimensions before and after.<BR> <BR>> This density was very <BR>> the experimental one. After enough minimization (convergent to 0.1 <BR>> kj/mol) and 3 ns of NVT MD, the NPT MD was performed as describe in the <BR>> original question. It happened to terminate with the fatal error after <BR>> about 300 ps. It should be noted that I ran the mpimdrun using the <BR>> following command:<BR>> mpirun -np 4 mpimdrun -s 03md.tpr -c poly_md.gro -v<BR>> I have a four-core CPU. Can you give more help according to these <BR>> descriptions?<BR> <BR>Your preparation seems good enough not to be the source of the problem.<BR> <BR>Have a look at the trajectory. What does the structure look like? What<BR>is happening to the box dimensions and/or density (look in the .log file<BR>and/or use g_energy)?.<BR> <BR>I believe that the error message you're getting is consistent with a<BR>drastic decrease in the box size. Since EM and NVT don't break and NPT<BR>does, the most likely hypothesis is that it started far too large.<BR> <BR>Mark<BR> <BR>> Best regards,<BR>> xiaowu<BR>> <BR>> wuxiao wrote:<BR>> > >Dear gmx users,<BR>> > >While doing a NPT molecular dynamics(MD), A fatal error was given:<BR>> > >The X-size of the box times triclinic skew factor is smaller than the<BR>> > >number of DD cells times the smallest allowed cell size.<BR>> > >However, either NVT or NVE MD can work well. Can you give any help to<BR>> > >deal with this problem?<BR>> <BR>> >It sounds like you might have started with a density that is too low for<BR>> >your box size.<BR>> <BR>> >A more complete description of your system and what you've done with it<BR>> >might help you get useful advice.<BR>> <BR>> Mark<BR>> <BR>> <BR>> ------------------------------------------------------------------------<BR><br /><hr />把MSN装进手机，更多聊天乐趣等你挖掘！ <a href='http://mobile.msn.com.cn/' target='_new'>立刻下载！</a></body>
</html>