<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Dear Berk,<BR>
In fact, I calculated the potential and force from the analytical function of angles in degrees, which is multicentred Gaussian-based potentials. So I think it must be some other reasons that cause the issues. Anyway, thanks again for your reply.<BR>
Best wishes,<BR>Chaofu Wu, Dr.<BR>
<HR id=stopSpelling>
From: gmx3@hotmail.com<BR>To: gmx-users@gromacs.org<BR>Subject: RE: [gmx-users] a bit strange errors<BR>Date: Fri, 16 Oct 2009 09:39:07 +0200<BR><BR>
<STYLE>
.ExternalClass .ecxhmmessage P
{padding:0px;}
.ExternalClass body.ecxhmmessage
{font-size:10pt;font-family:Verdana;}
</STYLE>
Hi,<BR><BR>But there could be a very similar issue there for the conversion from radians to degrees.<BR>The force has to be in units of kJ/mol/degree, not kJ/mol/radian.<BR><BR>Berk<BR><BR>
<HR id=ecxstopSpelling>
From: xiaowu759@hotmail.com<BR>To: gmx-users@gromacs.org<BR>Subject: RE: [gmx-users] a bit strange errors<BR>Date: Fri, 16 Oct 2009 15:35:19 +0800<BR><BR>
<STYLE>
.ExternalClass .ecxhmmessage P
{padding:0px;}
.ExternalClass body.ecxhmmessage
{font-size:10pt;font-family:Verdana;}
</STYLE>
Dear Berk,<BR> Unfortunately, this time it is angle (table_a0.xvg) that cause the problem. The angles are taken in degrees. I think it is not due to the wrong unit. Anyway, thanks very much for your reply.<BR> Best wishes,<BR>Chaofu Wu, Dr.<BR>
<HR id=ecxecxstopSpelling>
From: gmx3@hotmail.com<BR>To: gmx-users@gromacs.org<BR>Subject: RE: [gmx-users] a bit strange errors<BR>Date: Fri, 16 Oct 2009 09:21:18 +0200<BR><BR>
<STYLE>
.ExternalClass .ecxhmmessage P
{padding:0px;}
.ExternalClass body.ecxhmmessage
{font-size:10pt;font-family:Verdana;}
</STYLE>
Hi,<BR><BR>Last week another user reported very similar problems (also 164% deviations, I think).<BR>In that case it seemed that the user did an Angstrom to nm conversion in the table file<BR>on x without scaling f(x). Could this also be your issue?<BR><BR>Berk<BR><BR>
<HR id=ecxecxecxstopSpelling>
From: xiaowu759@hotmail.com<BR>To: gmx-users@gromacs.org<BR>Subject: RE: [gmx-users] a bit strange errors<BR>Date: Fri, 16 Oct 2009 14:55:40 +0800<BR><BR>
<STYLE>
.ExternalClass .ecxhmmessage P
{padding:0px;}
.ExternalClass body.ecxhmmessage
{font-size:10pt;font-family:Verdana;}
</STYLE>
Dear Mark,<BR> Thank you very much for the quick reply to my post. According to your suggestion, I check the table and reformat it with equal spacing between entries. I also check the entries calculated from the analytical functions and no wrong values are found. However, the same warnings are given, which lead to terminate the simulation. Are there any other reasons responsible for it?<BR>Best wishes,<BR>Chaofu Wu, Dr. <BR>> Date: Fri, 16 Oct 2009 14:18:17 +1100<BR>> From: Mark.Abraham@anu.edu.au<BR>> To: gmx-users@gromacs.org<BR>> Subject: Re: [gmx-users] a bit strange errors<BR>> <BR>> wuxiao wrote:<BR>> > Dear GMXers,<BR>> > Recently, I has been performing coarse-grained MD simulations based <BR>> > on the tabulated potentials using GROMACS. In the initialization, some <BR>> > warnings are given as follows <BR>> > /WARNING: For the 999 non-zero entries for table 0 in table_b0.xvg the <BR>> >
forces deviate on average 164% from minus the numerical derivative of <BR>> > the potential/<BR>> <BR>> This looks like your table has either wrong values or wrong formatting.<BR>> It does have to have equal spacing between entries.<BR>> <BR>> > It proceeds untill at about 17 ns before it is interupted with some errors<BR>> > /WARNING: 1-4 interaction between 148 and 151 at distance 13.125 which <B>> > is larger than the 1-4 table size 2.000 nm/<BR>> > I think the duartion time is long enough to equilibrated the system, <BR>> > which is confirmed by the energy plot with time. Why it happens so? What <BR>> > should I do to cope with this problem? Thanks a lot for any reply to <BR>> <BR>> Perhaps much of the table is fine, and the problem only arises when you<BR>> get a force being looked up in an erroneous part.<BR>> <BR>> Mark<BR>> _______________________________________________<BR>> gmx-users mailing
list gmx-users@gromacs.org<BR>> http://lists.gromacs.org/mailman/listinfo/gmx-users<BR>> Please search the archive at http://www.gromacs.org/search before posting!<BR>> Please don't post (un)subscribe requests to the list. Use the <BR>> www interface or send it to gmx-users-request@gromacs.org.<BR>> Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR><BR></B>
<HR>
<B>使用新一代 Windows L iv e Messenger 轻松交流和共享！ <A href="http://www.windowslive.cn/Messenger/">立刻下载！</A></B> <BR>
<HR>
What can you do with the new Windows Live? <A href="http://www.microsoft.com/windows/windowslive/default.aspx">Find out</A> <BR>
<HR>
更多热辣资讯尽在新版MSN首页！ <A href="http://cn.msn.com/">立刻访问！</A> <BR>
<HR>
Express yourself instantly with MSN Messenger! <A href="http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/">MSN Messenger</A> <br /><hr />Messenger保护盾2.0，更安全可靠的Messenger聊天！ <a href='http://im.live.cn/safe/' target='_new'>现在就下载！</a></body>
</html>